|
@ -1,4 +1,4 @@ |
|
|
filedata='maggio11_1731'; |
|
|
|
|
|
|
|
|
filedata='maggio11_1731'; |
|
|
str=['load ' filedata '.dat'];eval(str); |
|
|
str=['load ' filedata '.dat'];eval(str); |
|
|
str=['x=' filedata '(:,1);'];eval(str); |
|
|
str=['x=' filedata '(:,1);'];eval(str); |
|
|
str=['y=' filedata '(:,2);'];eval(str); |
|
|
str=['y=' filedata '(:,2);'];eval(str); |
|
@ -24,4 +24,9 @@ f=(0:1:length(x)-1)*1/(length(x)-1)*1/((x(2)-x(1))*3e-5); |
|
|
|
|
|
|
|
|
Rw=-imag(fft(y)./fft(ysI)); |
|
|
Rw=-imag(fft(y)./fft(ysI)); |
|
|
figure(2) |
|
|
figure(2) |
|
|
plot(f,Rw,f,Rw*0);xlim([0 2000]);ylim([-0.05 3]) |
|
|
|
|
|
|
|
|
plot(f,Rw,f,Rw*0);xlim([0 2000]);ylim([-0.05 3]) |
|
|
|
|
|
|
|
|
|
|
|
if 0, |
|
|
|
|
|
AA=[f' Rw]; str=['save fft.dat AA -ascii;'];eval(str); |
|
|
|
|
|
end; |
|
|
|
|
|
|