
¿Cómo puedo guardar resultados de una operación en forma de vector en .dat ?
Publicado por Sarita (5 intervenciones) el 18/08/2016 07:05:06
Este es mi programa y lo que quiero es guardar los valores de W1 y VDOS, en forma de vectores en archivos .dat no sé como hacerlo, por favor agradecerÃa mucho su ayuda. saludos. Les dejo el programa.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
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))
Valora esta pregunta


0