hc327.m

来自「这个是西安电子科技大学的陈怀琛教授编写的《数字信号处理教程——MATLAB释义与」· M 代码 · 共 24 行

%《数字信号处理教程——MATLAB释义与实现》第三章例3.2.7程序hc327
% 信号序列的DTFT
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
n = 0:10; x = (0.8*exp(j*pi/3)).^n;         % 给定输入序列x
k = -200:200; w = (pi/100)*k;               % 很密的频率下标和频率向量
X = x * (exp(-j*pi/100)) .^ (n'*k);         % x的DTFT计算
x1 = (0.8).^n;                              % 给定输入序列x的实部
X1 = x1 * (exp(-j*pi/100)) .^ (n'*k);       % x1的DTFT计算
subplot(4,2,1); stem(n,real(x));grid		% 以下绘制时间序列
title('复序列'),ylabel('实部')
subplot(4,2,3); stem(n,imag(x));grid
ylabel('虚部')
subplot(4,2,2); stem(n,real(x1));grid
title('实序列')
subplot(4,2,4); stem(n,imag(x1));grid
subplot(4,2,5); plot(w,abs(X));grid		    % 以下绘制频谱图
ylabel('幅度')
subplot(4,2,7); plot(w,angle(X));grid
ylabel('相角')
subplot(4,2,6); plot(w,abs(X1));grid
subplot(4,2,8); plot(w,angle(X1));grid
set(gcf,'color','w')                        % 置图形背景色为白