hc451a.m

某本matlab书籍中的代码

%《数字信号处理教程——MATLAB释义与实现》第四章例4.5.1程序hc451a
% 连续信号频谱计算中的T,N,L选择过程的进一步讨论
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
%T0=[2,1,0.5,0.1,0.1,0.1];   % 各次计算拟采用的T,编成向量T0
%L0=[10,10,10,10,20,40];     % 各次计算拟采用的N,编成向量N0
%for r=1:6                   % 循环计算六次
T=0.1;N=400;      % 根据计算顺序选用T及N
DOmega=2*pi/(N*T);              % 频率分辨率
n=0:N-1;x=exp(-0.1*n*T);    % 给出样本序列x
X=T*fftshift(fft(x));       % 求x的FFT，移到对称位置
X(1)                   % 显示乃奎斯特边界处的频谱X
k=floor(-(N-1)/2:(N-1)/2);  % 频率下标
subplot(2,2,3),plot(k*DOmega,abs(X))    % 在位置r处绘图
ylabel('X(\Omega)'),xlabel('\Omega')
axis([-10,10,0,10])           % 统一坐标，以便比较
set(gcf,'color','w')              % 置图形背景色为白
%end
T=0.1;N=400;      % 根据计算顺序选用T及N
DOmega=2*pi/(N*T);              % 频率分辨率
n=0:N-1;x1=[exp(-0.1*0.1*[1:100]),zeros(1,400-100)];    % 给出样本序列x
X1=T*fftshift(fft(x1));       % 求x的FFT，移到对称位置
[r,X1(1)]                    % 显示乃奎斯特边界处的频谱X
k=floor(-(N-1)/2:(N-1)/2);  % 频率下标
subplot(2,2,4),plot(k*DOmega,abs(X1))    % 在位置r处绘图
ylabel('X1(\Omega)'),xlabel('\Omega')
axis([-10,10,0,10])           % 统一坐标，以便比较
subplot(2,2,1),plot(n,x)    % 在位置r处绘图
ylabel('x(n)'),xlabel('n')
subplot(2,2,2),plot(n,x1)    % 在位置r处绘图
ylabel('x1(n)'),xlabel('n')