DSP实验报告B13011025分析

Q2.28运行程序P2.7，对序列h[n]和x[n]求卷积，生成y[n]，并用FIR滤波器h[n]对输入x[n]滤波，求得y1[n]。y[n]和y1[n]有差别吗？为什么要使用对x[n]补零后得到的x1[n]作为输入来产生y1[n]？

答：y[n] 和 y1[n] 的差别为 - 没有差别

将x[n]补零后得到 x1[n]作为输入,产生y1[n]的原因是 – 对于长度N1和N2的两个序列，转化率返回得到的序列长度N1+N2-1。与此相反，滤波器接收一个输入信号和系统规范，返回的结果是相同的长度作为输入信号。因此，为了从转化率和滤波器得到直接比较的结果，有必要供应滤波器的输入已经填充为长度L(x)+L(h)-1。

Q3.1在程序P3.1中，计算离散时间傅里叶变换的原始序列是什么？MATLAB命令pause的作用是什么？

2?z?1答：离散时间傅里叶变换的原始序列：； ?11?0.6zpause 命令的作用：暂停程序，直至用户按任意键，程序才可以开始。

Q3.2运行程序 P3.1，求离散时间傅里叶变换的实部、虚部以及幅度和相位谱。离散时间傅里叶变换是?的周期函数吗？若是，周期是多少？描述这四个图表示的对称性。

答：DTFT 是关于 ?的周期函数么？是周期函数 周期是 - 2π

四个图形的对称性为：实部是偶对称，虚部是奇对称，幅度谱是偶对称相位谱是奇对称。 Q3.4修改程序 P3_1 重做Q3.2的程序如下: clf;

w = -4*pi:8*pi/511:4*pi;

num = [1 3 5 7 9 11 13 15 17];den = [1]; h = freqz(num, den, w); subplot(2,1,1)

plot(w/pi,real(h));grid

title('Real part of H(e^{j\\omega})') xlabel('\\omega /\\pi'); ylabel('Amplitude'); subplot(2,1,2)

plot(w/pi,imag(h));grid

title('Imaginary part of H(e^{j\\omega})') xlabel('\\omega /\\pi'); ylabel('Amplitude'); pause

subplot(2,1,1)

plot(w/pi,abs(h));grid

title('Magnitude Spectrum |H(e^{j\\omega})|') xlabel('\\omega /\\pi'); ylabel('Amplitude'); subplot(2,1,2)

plot(w/pi,angle(h));grid

title('Phase Spectrum arg[H(e^{j\\omega})]') xlabel('\\omega /\\pi');

ylabel('Phase in radians'); 修改程序后的运行结果为: