一个简单的频率谱瀑布图(共1页).doc

《一个简单的频率谱瀑布图(共1页).doc》由会员分享，可在线阅读，更多相关《一个简单的频率谱瀑布图(共1页).doc（1页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上close all;clear all;A1=3;A2=4;f=100;%选取频率为100Hzfor i=1:12 t(i,:)=0:0.001:2;%时间间隔为0.001，说明采样频率为1000Hz for j=1:length(t(i,:) X1(i,j)=A1*exp(-0.9*t(i,j).*sin(2*pi*f*t(i,j); X2(i,j)=A2*exp(-1.5*t(i,j).*sin(2*pi*f*t(i,j); endend%N=length(t);%X1=A1*exp(-0.9*t).*sin(2*pi*f*t);%X2=A2*exp(-1.5*t

2、).*sin(2*pi*f*t);R=randn(12,length(t);Y=X1+X2+R;Y2=X1+X2;%频域分析l=length(Y);for i=1:12 r(i,:)=fft(Y(i,:),1000)/l;%。采样点个数1000 r1(i,:)=fftshift(r(i,:);%。反转折叠 rr(i,:)=r(i,1:500);%。消除对称 rr1(i,:)=r1(i,1:500);%。消除对称endfigure;plot(abs(r(1,:);%。单个频谱测试图%z=linspace(0,fs/2,l);%figure(1);plot(z,abs(r); grid; hold on;figure;waterfall(abs(r);%。对称频谱瀑布图figure;waterfall(abs(r1);%。对称反转频谱瀑布图figure;waterfall(abs(rr);%。消除对称瀑布图figure;waterfall(abs(rr1);专心-专注-专业