南邮2013DSP实验报告（英文mitra版） - 图文

Q2.28 程序P2_7产生的序列 y[n] and y1[n] 如下所示：

① y[n] 和 y1[n] 的差别为--- 它们无差别。

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

程序P3_1计算离散时间傅里叶变换的原始序列为---H(e

j?Q3.1

2?z?1)= 1?0.6z?1

pause 命令的作用为- 不加参数,直接用pause的话,就是程序暂停,直至用户按任意一个按键。如果加参数，比如pause(1.5)就是程序暂停1.5秒。

Q3.2 程序 P3_1 运行结果为:

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

四个图形的对称性为：实部是2?周期和偶对称； 虚部是2?周期和奇对称；幅度

是2?周期和偶对称；相位是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('H(e^{j\\omega})的实部') xlabel('\\omega /\\pi');

ylabel('振幅'); subplot(2,1,2)

plot(w/pi,imag(h));grid title('H(e^{j\\omega})的虚部') xlabel('\\omega /\\pi'); ylabel('振幅'); pause

subplot(2,1,1)

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

title('|H(e^{j\\omega})|幅度谱') xlabel('\\omega /\\pi'); ylabel('振幅'); subplot(2,1,2)

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

title('相位谱arg[H(e^{j\\omega})]') xlabel('\\omega /\\pi'); ylabel('以弧度为单位的相位');

修改程序后的运行结果为: