1、数字信号处理课后 Matlab 习题解答学院：电子信息工程学院专业：通信工程指导老师：钱满义2015 年 6 月 26 日第 2 页目录课后 Matlab 习题解答 .11、第一章 .2M1-1.2M1-2.3M1-3.5M1-4.7M1-5.8M1-6.92、第二章 .10M2-1.10M2-2.13M2-3.14M2-4.16M2-5:.18M2-6.193、第四章 .21M4-1.21M4-2.22M4-3.23M4-4.25M4-5.26M4-6.27M4-7.29第 3 页M4-8.30M4-9.324、第五章 .33M5-1.33M5-2.33M5-3.35M5-4.37M5-5.
2、37M5-6.37M5-7.38M5-8.38M5-9.38M5-10.385、第六章 .38M6-1.38M6-2.386、第七章 .42M7-1.42M7-2.43M7-3.45M7-4.46M7-5.46第 4 页第 5 页1、第一章M1-1代码：f=10; t=-0.2:0.001:0.2; g1=cos(6.*pi.*t); g2=cos(14.*pi.*t); g3=cos(26.*pi.*t); k=-0.2:1/f:0.2; h1=cos(6.*pi.*k); h2=cos(14.*pi.*k); h3=cos(26.*pi.*k); subplot(4,1,1); plot(
3、k,h1,r.,t,g1,r); xlabel(t); ylabel(g1(t); subplot(4,1,2); plot(k,h2,g.,t,g2,g); 第 6 页xlabel(t);ylabel(g2(t); subplot(4,1,3); plot(k,h3,b.,t,g3,b); xlabel(t); ylabel(g3(t); subplot(4,1,4);plot(t,g1,r,t,g2,g,t,g3,b,k,h1,r.,k,h2,g.,k,h3,b.) xlabel(t); ylabel(g(t); legend(g1(t),g2(t),g3(t);图像：第 7 页M1-2(
4、1)代码：b=1;a=1 -1.845 0.850586;w=linspace(0,2*pi,200); y=filter(b,a,w);plot(w,y);ylabel(单位脉冲响应);xlabel(Normalized frequency);%b=1;%a=1 -1.845 0.850586;h=impz(b,a,21);%figure(1);stem(0:20,h);ylabel(单位脉冲响应);xlabel(Normalized frequency);图像：第 8 页(2)代码：b=1;a=1 -1.85 0.85; w=linspace(0,2*pi,200); y=filter(b
5、,a,w); plot(w,y); ylabel(单位脉冲响应); xlabel(Normalized frequency); b=1; a=1 -1.85 0.85;h=impz(b,a,21); figure(1);stem(0:20,h); ylabel(单位脉冲响应);第 9 页图像：M1-3(1)代码：x=firls(511,0 0.4 0.404 1,1 1 0 0); k=linspace(0,1,512); stem(k,x);图像：第 10 页(2)代码：x=firls(511,0 0.4 0.404 1,1 1 0 0); k=linspace(0,1,512); stem(k,x);k1=0.4*pi;y1=x.*cos(k1*k);A=abs(y1);subplot(4,1,1);plot(A);k2=0.8*pi;y2=x.*cos(k2*k);B=abs(y2);subplot(4,1,2);plot(B);
