test.m

来自「可以进行傅立叶变换的matlab程序 保证可以使用」· M 代码 · 共 30 行

clear
close all
load D001.dat
R=D001(:,2);
B=D001(:,1);

%figure
%plot(B,R,'r.');

n=16
dim=2^n;
inc=(1-0.25)/(2^n-1)
Bf=0.25:inc:1;

t=interp1(1./B,R,Bf);
offset=sum(t)/dim;
f=fft(t-offset,dim)/dim;
f2=2*abs(f(1:dim/2));

figure
plot(Bf,t);
xlabel('1/B(1/T):time');
ylabel('R(arb. units)');

figure
freq=0:1/inc/(dim-1):1/inc/2;
plot(freq,f2,'r-');
axis([-100 100 0 1])
xlabel('Freq');
ylabel('Intensity (arb. units)');