xcxd.m

十几个数字信号处理实验课的程序

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