数字信号处理matlab实验第三章报告分解.docxVIP

第三章离散时间信号的频域分析一、实验目的通过matlab仿真进一步研究离散时间系统的时域特性。二、实验要求学习并调试本章所给的例子。回答书后给出的问题。三、实验过程Q.2运行程序P3.1,求离散时间傅里叶变换的实部、虚部以及幅度和相位谱。离散时间傅里叶变换是w的周期函数吗?若是周期是多少？描述这四个图形表示的对称性。clf;w = -4*pi:8*pi/511:4*pi;num = [2 1];den = [1 -0.6];h = freqz(num, den, w);% Plot the DTFTsubplot(2,2,1)plot(w/pi,real(h));gridtitle(H(e^{j\omega})的实部)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,2)plot(w/pi,imag(h));gridtitle(H(e^{j\omega})的虚部)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,3)plot(w/pi,abs(h));gridtitle(|H(e^{j\omega})|幅度谱)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,4)plot(w/pi,angle(h));gridtitle([H(e^{j\omega})]相位谱)xlabel(\omega /\pi);ylabel(以弧度为单位的相位);由图得出是周期函数，周期为2，实部是和幅度谱是实偶对称，虚部和相位谱是实奇对称Q.4修改程序P3.1，计算如下有限长序列的离散时间傅里叶变换：g[n]=[1 3 5 7 9 11 13 15 17]并重做习题Q.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);% Plot the DTFTsubplot(2,2,1)plot(w/pi,real(h));gridtitle(H(e^{j\omega})的实部)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,2)plot(w/pi,imag(h));gridtitle(H(e^{j\omega})的虚部)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,3)plot(w/pi,abs(h));gridtitle(|H(e^{j\omega})|幅度谱)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,4)plot(w/pi,angle(h));gridtitle([H(e^{j\omega})]相位谱)xlabel(\omega /\pi);ylabel(以弧度为单位的相位);1.由图得出是周期函数，周期为2，实部是和幅度谱是实偶对称，虚部和相位谱是实奇对称。2.要计算一个系统相频特性，就要用到反正切函数，计算机中反正切函数规定，在一、二象限中的角度为0～pi，三四象限的角度为0～-pi。若一个角度从0变到2pi，但实际得到的结果是0～pi，再由-pi～0，在w=pi处发生跳变。Q.6通过加入合适的注释语句和程序语句，修改程序P3.2，对程序生成的图形的两个轴加标记。哪个参数控制时移量？clf;w = -pi:2*pi/255:pi; wo = 0.4*pi; D = 10;num = [1 2 3 4 5 6 7 8 9];h1 = freqz(num, 1, w);h2 = freqz([zeros(1,D) num], 1, w);subplot(2,2,1)plot(w/pi,abs(h1));gridtitle(原序列的幅度谱)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,2)plot(w/pi,abs(h2));gridtitle(时移后序列的幅度谱)xlabel(\omega /\pi);ylabel(振幅);subplot(2,2,3)plot(w/pi,angle(h1));gridtitle(原序列的相位谱)xlabel(\omega /\pi);ylabel(以弧度为单位的相位);subplot(2,2,4)plot(w/pi,angle(h2));gridtitle(时移后序列的相位谱)xlabel(\omega /\pi);ylabel(以弧度为单位的相位);参数D控制位移量Q.8选取不同的时移值重做习题Q.6令Q.6中的位移量D为20，程序结果如图：Q.10通过加入合适的注释语句和程序语句，修