第 PAGE 3頁數字信號處理課后習題解答學院：電子信息工程學院專業：通信工程指導老師：錢滿義2015年6月26日目錄 TOC \o "1-3" \h \z \u 課后習題解答 1、第一章 M1-1 M1-2 M1-3 M1-4 M1-5 M1-6 2、第二章 M2-1 M2-2 M2-3 M2-4 M2-5: M2-6 3、第四章 M4-1 M4-2 M4-3 M4-4 M4-5 M4-6 M4-7 M4-8 M4-9 4、第五章 M5-1 M5-2 M5-3 M5-4 M5-5 M5-6 M5-7 M5-8 M5-9 M5-10 5、第六章 M6-1 M6-2 6、第七章 M7-1 M7-2 M7-3 M7-4 M7-5 461、第一章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); (4,1,1); plot(k,h1,'r。',t,g1,'r'); ('t'); ('g1(t)'); (4,1,2); plot(k,h2,'g。',t,g2,'g'); ('t');('g2(t)'); (4,1,3); plot(k,h3,'b。',t,g3,'b'); ('t'); ('g3(t)'); (4,1,4);plot(t,g1,'r',t,g2,'g',t,g3,'b',k,h1,'r。',k,h2,'g。',k,h3,'b。') ('t'); ('g(t)'); ('g1(t)','g2(t)','g3(t)');圖像：M1-2(1)代碼：b=[1];a=[1 -1。
845 0。];w=(0,2*pi,200); y=(b,a數字信號處理matlab實驗答案,w);plot(w,y);('單位脈沖響應');(' ');%b=[1];%a=[1 -1。845 0。];h=impz(b,a,21);%(1);stem(0:20,h);('單位脈沖響應');(' ');圖像：(2)代碼：b=[1];a=[1 -1。85 0。85]; w=(0,2*pi,200); y=(b,a,w); plot(w,y); ('單位脈沖響應'); (' '); b=[1]; a=[1 -1。85 0。85];h=impz(b,a,21); (1);stem(0:20,h); ('單位脈沖響應');圖像：M1-3(1)代碼：x=firls(511,[0 0。4 0。404 1],[1 1 0 0]); k=(0,1,512); stem(k,x);圖像：(2)代碼：x=firls(511,[0 0。
4 0。404 1],[1 1 0 0]); k=(0,1,512); stem(k,x);k1=0。4*pi;y1=x。*cos(k1*k);A=abs(y1);(4,1,1);plot(A);k2=0。8*pi;y2=x。*cos(k2*k);B=abs(y2);(4,1,2);plot(B);k3=0。9*pi;y3=x。*cos(k3*k);C=abs(y3);(4,1,3);plot(C);k4=pi;y4=x。*cos(k4*k);D=abs(y4);(4,1,4);plot(D);圖像：M1-4代碼：b=[1 1];a=2;w=(0,pi,512); H1=freqz(b,a,w); plot(w/pi,(abs(H1)),'r'); hold on b1=[1。8]; a1=[1 -0。8]; w=(0,pi數字信號處理matlab實驗答案,512); H2=freqz(b1,a1,w); plot(w/pi,(abs(H2)),'g'); hold on b2=conv(b,b1); a2=conv(a,a1); w=(0,pi,512); H3=freqz(b2,a2,w); plot(w/pi,(abs(H3),'b'));('幅度'); (' ');圖像：M1-5(1)代碼：b1=[0。
0534 0。0534]; b2=[1 2。0332 2。0169 -1。0166 1]; a1=[1 -0。683]; a2=[1 -1。4661 0。7957]; b=conv(b1,b2); a=conv(a1,a2); w=(0,pi,200); H=freqz(b,a,w);(2,1,1); plot(w/pi,abs(H)); ('幅度'); (' '); (2,1,2); plot(w,angle(H)); ('相位'); (' ');圖像：(2)代碼：b1=[1 -2 1];b2=[1 -2 1]; a1=[1 -1。499 0。8482];a2=[1 -1。5548 0。6493]; b=conv(b1,b2); a=conv(a1,a2); w=(0,pi,200); H=freqz(b,a,w); (2,1,1); plot(w/pi,abs(H)); ('幅度'); (' '); (2,1,2); plot(w,angle(H)); ('相位'); (' ');圖像：M1-6代碼：b=[1 2 0。
99]; a=[1 1。55 0。6]; w=(0,pi,512); H=freqz(b,a,w); plot(w,(angle(H)),'r'); hold on; b1=[1。1 1。99 0。9]; a1=[1 1。55 0。6]; w=(0,pi,512); Hmin=freqz(b1,a1,w); plot(w,(angle(Hmin)),'g');hold on; b2=[0。9 1。99 1。1]; a2=[1 1。55 0。6]; w=(0,pi,512); Hmax=freqz(b2,a2,w); plot(w,(angle(Hmax)),'m'); ('相位'); ('omega');圖像：2、第二章M2-1（1）代碼：g=[1 -3 4 2 0 -2]; h=[3 0 1 -1 2 1]; l=(g); L=2*l-1; GE=fft(g,L); HE=fft(h,L); y1=ifft(GE。*HE); for n=1:lifn+l