|
正好这么多老师在 我借个地方 问1问
我自己用了几个简单的函数帮助自己理解傅立叶变换 请你们帮我看看对不对
===分割线==
%输入数据
fn=5000; %取点的个数,同时也是测试的频率,因此要大于500的两倍
x=-0.05:1/fn:0.05; %激励条件x
f1=50; %1阶固有频率
f2=200; %2阶固有频率
f3=500; %3阶固有频率
y1=1*sin(2*pi*f1*x); %1阶响应方程
y2=2*sin(2*pi*f2*x); %2阶响应方程
y3=3*sin(2*pi*f3*x); %3阶响应方程
y=y1+y2+y3; %固有响应方程
%分别绘制3阶的波形和总的响应
subplot(411);
plot(x,y1);
title('y1');
subplot(412);
plot(x,y2);
title('y2');
subplot(413);
plot(x,y3);
title('y3');
subplot(414);
plot(x,y);
title('y');
%定义点个数,频率数,快速傅立叶变换,绘制频谱图,观察3阶频率
N=2^12;
f=[[0:(N/2)-1]*fn/N];
FFTy=fft(y,N);
figure;
plot(f,abs(FFTy(1:(N/2))));
xlim([0,1250])
===分割线===
然后通过拾点 可以找到几个响应最大的点 y轴的大小是不是也体现了贡献的大小呢? |
|