u1_3.m

来自「matlab 实现的离散信号分析」· M 代码 · 共 18 行

%单位脉冲序列
n=1:50;
x=zeros(1,50);
x(1)=1;
close all;
subplot(3,1,1);
stem(x);
title('单位冲击信号序列');
k=-25:25;
X=x*(exp(-j*pi/12.5)).^(n'*k);
magX=abs(X);
subplot(3,1,2);
stem(magX);
title('单位冲击信号的幅度谱');
angX=angle(X);
subplot(3,1,3);
stem(angX);
title('单位冲击信号的相位谱');