rawdata.m

来自「基于Matlab编写的用于对数据进行插值分析和傅立叶变换的程序」· M 代码 · 共 30 行

clear
clc
load('C:\D_001.dat')
R=D_001(:,3);
B=D_001(:,1);

Bf=1./B;
Bf2=Bf(1:9501);
Rf2=R(1:9501);

figure
subplot(211)
plot(Bf2,Rf2)


n=11
dim=2^n
f=fft(Rf2,dim)/dim;
f2=abs(f(1:dim/2));


nf=0:dim-1;
incf=1/1.9;
freq=incf*nf(1:dim/2);
subplot(212)
plot(freq,f2,'r-');
axis([0 100 0 0.3])
xlabel('Freq  n=11');
ylabel('Intensity (arb. units)');