资源描述:
《北京交通大学数字信号处理matlab习题问题详解(钱my).doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、数字信号处理课后Matlab习题解答学院:电子信息工程学院专业:通信工程指导老师:钱满义2015年6月26日目录课后Matlab习题解答11、第一章2M1-12M1-23M1-35M1-47M1-58M1-692、第二章10M2-110M2-213M2-314M2-416M2-5:18M2-6193、第四章21M4-121M4-222M4-323M4-425M4-526M4-627M4-729M4-830M4-9324、第五章33M5-133M5-233M5-335M5-437M5-537M5-637M5-738M5-838M5-938M5-10385、第六章38M6-1
2、38M6-2386、第七章42M7-142M7-243M7-345M7-446M7-5461、第一章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(k,h1,'r.',t,g1,'r');xlabel('t');ylabel('g1(t)');subplot(4,1,2);plot(k
3、,h2,'g.',t,g2,'g');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)');图像:M1-2(1)代码:b=[1];a=[1-1.8450.850586]
4、;w=linspace(0,2*pi,200);y=filter(b,a,w);plot(w,y);ylabel('单位脉冲响应');xlabel('Normalizedfrequency');%b=[1];%a=[1-1.8450.850586];h=impz(b,a,21);%figure(1);stem(0:20,h);ylabel('单位脉冲响应');xlabel('Normalizedfrequency');图像:(2)代码:b=[1];a=[1-1.850.85];w=linspace(0,2*pi,200);y=filter(b,a,w);plot(w,y)
5、;ylabel('单位脉冲响应');xlabel('Normalizedfrequency');b=[1];a=[1-1.850.85];h=impz(b,a,21);figure(1);stem(0:20,h);ylabel('单位脉冲响应');图像:M1-3(1)代码:x=firls(511,[00.40.4041],[1100]);k=linspace(0,1,512);stem(k,x);图像:(2)代码:x=firls(511,[00.40.4041],[1100]);k=linspace(0,1,512);stem(k,x);k1=0.4*pi;y1=x.*c
6、os(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);k3=0.9*pi;y3=x.*cos(k3*k);C=abs(y3);subplot(4,1,3);plot(C);k4=pi;y4=x.*cos(k4*k);D=abs(y4);subplot(4,1,4);plot(D);图像:M1-4代码:b=[11];a=2;w=linspace(0,pi,512);H1=freqz(b,a,w);plot(w/pi,unwrap
7、(abs(H1)),'r');holdonb1=[1.8];a1=[1-0.8];w=linspace(0,pi,512);H2=freqz(b1,a1,w);plot(w/pi,unwrap(abs(H2)),'g');holdonb2=conv(b,b1);a2=conv(a,a1);w=linspace(0,pi,512);H3=freqz(b2,a2,w);plot(w/pi,unwrap(abs(H3),'b'));ylabel('幅度');xlabel('Normalizedfrequency');图像:M1-5(