sin00.m

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

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

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)');