li.m

来自「十几个数字信号处理实验课的程序」· M 代码 · 共 29 行

for i=1:4 
    xc(i)=i;
end
for i=5:8
    xc(i)=9-i;
end

for i=1:4 
    xd(i)=5-i;
end
for i=5:8
    xd(i)=i-4;
end
n=0:7;
w=0:pi/511:pi;
Xc=xc*(exp(-j).^(n'*w)); 
Xd=xd*(exp(-j).^(n'*w)); 

figure
subplot(2,2,1);plot(w/pi,abs(Xc));title('DTFT[xc]的幅度谱');
grid on
subplot(2,2,3);plot(w/pi,abs(Xd));title('DTFT[xd]的幅度谱');
grid on


subplot(2,2,2);plot(w/pi,angle(Xc));title('DTFT[xc]的相位谱');
grid 
subplot(2,2,4);plot(w/pi,angle(Xd));title('DTFT[xd]的相位谱');
grid