xdxc.m

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

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
figure
subplot(2,1,1);stem(abs(fft(xc,16)));title('16点的fft[xc]');
subplot(2,1,2);stem(abs(fft(xd,16)));title('16点的fft[xd]');
figure
subplot(2,1,1);stem(angle(fft(xc,16)));title('16点的fft[xc]的相位谱');
subplot(2,1,2);stem(angle(fft(xd,16)));title('16点的fft[xd]的相位谱');