一个简单的频率谱瀑布图

closeall;clearall;A1=3;A2=4;f=100;%选取频率为100Hzfori=1:12t(i,:)=0:0.001:2;%时间间隔为0.001，说明采样频率为1000Hzforj=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).*sin(2*pi*f*t);R=randn(12,length(t));Y=X1+X2+R;Y2=X1+X2;%频域分析l=length(Y);fori=1:12r(i,:)=fft(Y(i,:),1000)/l;%。。。。。采样点个数1000r1(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;holdon;figure;waterfall(abs(r));%。。。。。。。。对称频谱瀑布图figure;waterfall(abs(r1));%。。。。。。。。对称反转频谱瀑布图figure;waterfall(abs(rr));%。。。。。。。。消除对称瀑布图figure;waterfall(abs(rr1));