feizhouqi.m

fft做谱分析各种连续信号 用MATLAB实现

%求xa(t)=exp(-0.01t)cos2t+2exp(-0.015t)sin2.1t  t>0  的频谱 
T0=[0.6 0.15 0.15 0.15];  %各次计算采用的 T 编成向量T0
N0=[256,256,256,2048];    %各次计算采用的 N 编成向量N0
for r=1:4
T=T0(r);N=N0(r);n=0:N-1;  %根据计算顺序选用T和N
    D=2*pi/(N*T);             %频率分辨率 D=2*pi/(N*T)
    x=exp(-0.01*n*T).*cos(2*n*T)+2*exp(-0.015*n*T);  %  序列 x
    Xa=T*fftshift(fft(x));                            %  求x的FFT变换  移动到绝对位置
    [r,Xa(1)]                                         %检查各次循环所得奈奎斯特频率处的副特性
    k=floor(-(N-1)/2:(N-1)/2);                      %奈奎斯特频率区间下标向量
   subplot(2,2,r),plot(k*D,abs(Xa))                %画幅频特性
 end