Matlab - ¿Cómo puedo guardar resultados de una operación en forma de vector en .dat ?

clc

clear all

% Programa que sirve para construir la matriz dinámica y diagonalizarla det

% det(A - wI)=0,

fprintf('Operación matricial \n\n')

disp('Cálculo de la inversa de una matriz')

        fprintf('Ingresar una Matriz A \n\n')

natom=3;

beta=1.0d0;

masa=1.0d0;

md=zeros(natom,natom);

numdiv=50

VDOS=0;

Xfinal=4

pasos=Xfinal/numdiv

for z=1:numdiv

  zreal=z*pasos

    w2=(i*(1/natom))+zreal

  freq(z)=zreal

    md(1,1)=beta; %Matriz dinámica

for j=2:natom

    md(j,j)=2.0d0*beta;

    md(j-1,j)=-beta;

    md(j,j-1)=-beta;

end

   md(natom,natom)= beta;

mI=eye(natom); %Matriz identidad

md2=masa*w2*mI;

mm=md-md2;

y=imag(inv(mm))

temporal=0;

for t=1:1:natom

    temporal= temporal+y(t,t);

end

suma=temporal

    VDOS(z)= suma * (1/pi);

end

for z1=1:numdiv

    w1(z1)=sqrt(freq(z1))

    VDOS1(z1)=w1(z1)*VDOS(z1)

end

z1=1:numdiv

plot(w1(z1),VDOS1(z1))

export(mat2dataset([w1',VDOS'],'Varnames',{'w1';'VDOS'}),'Delimiter','\t','file','datos.dat')

export(mat2dataset([w1',VDOS']),'file','datos.dat')

clear all

% Programa que sirve para construir la matriz dinámica y diagonalizarla det

% det(A - wI)=0,

fprintf('Operación matricial \n')

disp('Cálculo de la inversa de una matriz')

fprintf('Ingresar una Matriz A \n\n')

natom=3;

beta=1.0d0;

masa=1.0d0;

md=zeros(natom,natom);

numdiv=50;

VDOS=0;

Xfinal=4;

pasos=Xfinal/numdiv;

for z=1:numdiv

  zreal=z*pasos;

    w2=(i*(1/natom))+zreal;

  freq(z)=zreal;

    md(1,1)=beta; %Matriz dinámica

for j=2:natom

    md(j,j)=2.0d0*beta;

    md(j-1,j)=-beta;

    md(j,j-1)=-beta;

end

md(natom,natom)= beta;

mI=eye(natom); %Matriz identidad

md2=masa*w2*mI;

mm=md-md2;

y=imag(inv(mm));

temporal=0;

for t=1:1:natom

    temporal= temporal+y(t,t);

end

suma=temporal;

VDOS(z)= suma * (1/pi);

end

for z1=1:numdiv

    w1(z1)=sqrt(freq(z1));

    VDOS1(z1)=w1(z1)*VDOS(z1);

end

z1=1:numdiv;

plot(w1(z1),VDOS1(z1))

M=[w1' VDOS'];

nombre='archivos.dat';

fid=fopen(nombre,'w');

fprintf(fid,'  w1   VDOS   \n');

for i=1:size(M,1)

    for j=1:size(M,2)

     fprintf(fid,'%7.4f\t',M(i,j));

    end

    fprintf(fid,'\n');

end

fclose(fid);

open(nombre)

>> archivos_dat

Operación matricial

Cálculo de la inversa de una matriz

Ingresar una Matriz A