sin00.asv

基于Matlab编写的用于对数据进行插值分析和傅立叶变换的程序

clear
x=0.1:0.01:2;
y=(12*x/sinh(3*x))*exp(-10*x);
figure
plot(x,y,'r-');

n=16
dim=2^n;
inc=1.9/(2^n-1);

xf=0.1:inc:2;

t=interp1(x,y,xf);
offset=sum(t)/dim;
f=fft(t-offset,dim)/dim;
f2=2*abs(f(1:dim/2));

figure
plot(xf,t,'r-');
xlabel('X:time');
ylabel('Y(arb. units)');

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