北京交通大学数字信号处理matlab习题答案钱my.docx

数字信号处理课后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-138M6-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,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.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);图像：(2)代码：b=1;a=1 -1.85 0.85; w=linspace(0,2*pi,200); y=filter(b,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(单位脉冲响应);图像：M1-3(1)代码：x=firls(511,0 0.4 0.404 1,1 1 0 0); k=linspace(0,1,512); stem(k,x);图像：(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);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=1 1;a=2;w=linspace(0,pi,512); H1=freqz(b,a,w); plot(w/pi,unwrap(abs(H1),r); hold on b1=1.8; a1=1 -0.8; w=linspace(0,pi,512); H2=freqz(b1,a1,w); plot(w/pi,unwrap(abs(H2),g); hold on b2=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(Normalized frequency);图像：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=linspace(0,pi,200); H=freqz(b,a,w);subplot(2,1,1); plot(w/pi,abs(H); ylabel(幅度); xlabel(Normalized frequency); subplot(2,1,2); plot(w,angle(H); ylabel(相位); xlabel(Normalized frequency);图像：(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=linspace(0,pi,200); H=freqz(b,a,w); subplot(2,1,1); plot(w/pi,abs(H); ylabel(幅度); xlabel(Normalized frequency); subplot(2,1,2); plot(w,angle(H); ylabel(相位); xlabel(Normalized frequency);图像：M1-6代码：b=1 2 0.99; a=1 1.55 0.6; w=linspace(0,pi,512); H=freqz(b,a,w); plot(w,unwrap(angle(H),r); hold on; b1=1.1 1.99 0.9; a1=1 1.55 0.6; w=linspace(0,pi,512); Hmin=freqz(b1,a1,w); plot(w,unwrap(angle(Hmin),g);hold on; b2=0.9 1.99 1.1; a2=1 1.55 0.6; w=linspace(0,pi,512); Hmax=freqz(b2,a2,w); plot(w,unwrap(angle(Hmax),m); ylabel(相位); xlabel(omega);图像：2、第二章M2-1（1）代码：g=1 -3 4 2 0 -2; h=3 0 1 -1 2 1; l=length(g); L=2*l-1; GE=fft(g,L); HE=fft(h,L); y1=ifft(GE.*HE); for n=1:lifn+l=L y2(n)=y1(n)+y1(n+l); elsey2(n)=y1(n); endendy2 stem(0:l-1,y2) xlabel(k) ylabel(y(k) title(循环卷积)结果：y2 = 6.0000 -3.0000 17.0000 -2.0000 7.0000 -13.0000图像：（2）代码：k=0:5; x=cos(pi.*k./2); y=3.k; l=length(x); L=2*l-1;GE=fft(x,L); HE=fft(y,L); y1=ifft(GE.*HE); for n=1:l ifn+l=Ly2(n)=y1(n)+y1(n+l); elsey2(n)=y1(n); endendy2 stem(0:l-1,y2) xlabel(k) ylabel(y(k) title(循环卷积)结果：y2 = -71.0000 -213.0000 89.0000 267.0000 73.0000 219.0000图像：M2-2（1）代码：N=10; k=-N:N; x=cos(k.*pi./(2*N); W=linspace(-pi,pi,512); X=zeros(1,length(W); for k=-N:N X1=x(k+N+1).*exp(-j.*W.*k); X=X+X1; endplot(W,abs(X) xlabel(W); ylabel(abs(X);图像：（2）代码：N=10; k=-N:N; x=cos(k.*pi./(2*N); X_21=fft(x,21); L=-10:10; W=linspace(-pi,pi,1024); X=zeros(1,length(W); for k=-N:N X1=x(k+N+1).*exp(-j.*W.*k); X=X+X1; endplot(W,abs(X);hold on; plot(2*pi*L/21,fftshift(abs(X_21),o); xlabel(W); ylabel(abs(X);图像：M2-3代码：N=64; L=1024; f1=100;f2=120;fs=800; A=1;B1=1;B2=0.5;B3=0.25;B4=0.05; T=1/fs; ws=2*pi*fs; k=0:N-1; x1=A*cos(2*pi*f1*T*k)+B1*cos(2*pi*f2*T*k); x2=A*cos(2*pi*f1*T*k)+B2*cos(2*pi*f2*T*k); x3=A*cos(2*pi*f1*T*k)+B3*cos(2*pi*f2*T*k); x4=A*cos(2*pi*f1*T*k)+B4*cos(2*pi*f2*T*k); hf=(hamming(N); x1=x1.*hf; x2=x2.*hf; x3=x3.*hf; x4=x4.*hf; X1=fftshift(fft(x1,L); X2=fftshift(fft(x2,L); X3=fftshift(fft(x3,L); X4=fftshift(fft(x4,L); W=T*(-ws/2+(0:L-1)*ws/L)/(2*pi); subplot(2,2,1);plot(W,abs(X1); title(A=1,B=1); xlabel(W); ylabel(X1); subplot(2,2,2);图像：F2=120HZ时：F2=140HZ时：F2=160HZ时：M2-4（1）代码：W0=2*pi/15; W1=2.3*pi/15; N=64; k=0:N-1; x=cos(W0*k)+0.75*cos(W1*k); X=fft(x); plot(k/N,abs(X);grid on; title(64点FFT);图像：（2）代码：W0=2*pi/15; W1=2.3*pi/15; N=64; L=1024; k=0:N-1; x=cos(W0*k)+0.75*cos(W1*k); X=fft(x,L); plot(0:L-1)/N,abs(X); grid on; title(1024点FFT);图像：M2-5:代码：fs=100; ws=2*pi*fs; Ts=1/fs; N=fs; x=exp(-3*Ts*(0:N-1); y=fft(x,N); l=length(y); k=linspace(-ws/2,ws/2,l); plot(k,Ts*fftshift(abs(y),b:); hold on; w=linspace(-ws/2,ws/2,1024); y1=sqrt(1./(9+w.2); plot(w,y1,r) title(fs=100Hz时的频谱) legend(近似值,实际值);图像：fs=10HZ时：fs=50Hz时：fs=100HZ时:M2-6代码：Ts=0.5; N=4;N0=64; k=(-N/2:(N/2)*Ts;x=exp(-pi*(k).2); X=Ts*fftshift(fft(x,N0); w=-pi/Ts:2*pi/N0/Ts:(pi-2*pi/N0)/Ts; XT=(pi/pi)0.5*exp(-w.2/4/pi);subplot(2,1,1) plot(w/pi,abs(X),-o,w/pi,XT); xlabel(omega/pi); ylabel(X(jomega); legend(试验值,理论值); title(Ts=,num2str(Ts) N=,num2str(N); subplot(2,1,2) plot(w/pi,abs(X)-XT) ylabel(实验误差) xlabel(omega/pi);图像：Ts=1,N=2:Ts=0.5 ,N=2:Ts=0.5 ,N=4:3、第四章M4-1代码：WP=10;WS=2;Ap=1;As=40; wp=1/WP;ws=1/WS; w0=1;B=2; N,Wc=cheb1ord(wp,ws,Ap,As,s); num,den = cheby1(N,Ap,Wc,s);numt,dent = lp2hp(num,den,1); w=linspace(1,12,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h) ; grid ; xlabel(Frequency in rad/s); ylabel(Gain in dB)图像：M4-2代码：wp=1;ws=3.3182;Ap=1;As=32; w0=sqrt(48);B=2; N,Wc=buttord(wp,ws,Ap,As,s); num,den = butter(N,Wc,s); numt,dent = lp2bp(num,den,w0,B); w=linspace(2,12,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h) ; grid ; xlabel(Frequency in rad/s); ylabel(Gain in dB)图像：M4-3（1）代码：wp=1;ws=3.3182;Ap=1;As=32; w0=sqrt(48);B=2; N,Wc=cheb1ord(wp,ws,Ap,As,s); num,den = cheby1(N,Ap,Wc,s); numt,dent = lp2bp(num,den,w0,B); w=linspace(2,12,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h) ; grid ; xlabel(Frequency in rad/s); ylabel(Gain in dB)图像：（2）代码 ：wp=1;ws=3.3182;Ap=1;As=32; w0=sqrt(48);B=2; N,Wc=cheb2ord(wp,ws,Ap,As,s); num,den = cheby2(N,Ap,Wc,s); numt,dent = lp2bp(num,den,w0,B); w=linspace(2,12,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h) ; grid ; xlabel(Frequency in rad/s); ylabel(Gain in dB)图像：（3）代码：wp=1;ws=3.3182;Ap=1;As=32; w0=sqrt(48);B=2; N,Wc=ellipord(wp,ws,Ap,As,s); num,den = ellip(N,Ap,As,Wc,s); numt,dent = lp2bp(num,den,w0,B); w=linspace(2,12,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h) ; grid ; xlabel(Frequency in rad/s); ylabel(Gain in dB)图像：M4-4代码：Ap=1;As=10;wp1=6;wp2=13;ws1=9;ws2=11; B=ws2-ws1;w0=sqrt(ws1*ws2); wLp1=B*wp1/(w0*w0-wp1*wp1); wLp2=B*wp2/(w0*w0-wp2*wp2);wLp=max(abs(wLp1),abs(wLp2); N,Wc=buttord(wLp,1,Ap,As,s); num,den = butter(N,Wc,s ); numt,dent=lp2bs(num,den,w0,B); w=linspace(5,14,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h); w=wp1 ws1 ws2 wp2; set(gca,xtick,w);grid; h=freqs(numt,dent,w);A=-20*log10(abs(h);图像：M4-5代码：Ap=1;As=10;wp1=6;wp2=13;ws1=9;ws2=11; B=ws2-ws1;w0=sqrt(ws1*ws2); wLp1=B*wp1/(w0*w0-wp1*wp1); wLp2=B*wp2/(w0*w0-wp2*wp2); wLp=max(abs(wLp1),abs(wLp2); N,Wc=ellipord(wLp,1,Ap,As,s); num,den = ellip(N,Ap,As,Wc,s); numt,dent = lp2bp(num,den,w0,B); w=linspace(5,14,1000); h=freqs(numt,dent,w); plot(w,20*log10(abs(h); w=wp1 ws1 ws2 wp2; set(gca,xtick,w);grid; h=freqs(numt,dent,w);A=-20*log10(abs(h);图像：M4-6代码：Wp=0.1*pi;Ws=0.4*pi;Ap=1;As=25; Fs=1 wp=Wp*Fs;ws=Ws*Fs; N=buttord(wp,ws,Ap,As,s); wc=wp/(10(0.1*Ap)-1)(1/2/N); numa,dena=butter(N,wc,s); numd,dend=impinvar(numa,dena,Fs);w=linspace(0,pi,512); h=freqz(numd,dend,w); norm=max(abs(h); numd=numd/norm; plot(w/pi,20*log10(abs(h)/norm);w=WpWs; h=freqz(numd,dend,w); fprintf(Ap= %.4fn,-20*log10( abs(h(1);fprintf(As= %.4fn,-20*log10( abs(h(2);图像：代码：WP=0.1*pi;WS=0.4*pi;Ap=1;As=25; Fs=1; wp=WP*Fs;ws=WS*Fs; N=cheb2ord(wp,ws,Ap,As,s);wc=wp/(10(0.1*Ap)-1)(1/2/N); numa,dena=cheby2(N,As,wc,s);numd,dend=impinvar(numa,dena,Fs); w=linspace(0,pi,512); h=freqz(numd,dend,w); norm=max(abs(h); numd=numd/norm; plot(w/pi,20*log10(abs(h)/norm) w=wpws; h=freqz(numd,dend,w); fprintf(Ap= %.4fn,-20*log10( abs(h(1);fprintf(As= %.4fn,-20*log10( abs(h(2);图像：M4-7代码:Wp=0.2*pi; Ws=0.4*pi; Ap=1; As=15;Fs=1; wp=Wp*Fs; ws=Ws*Fs; N,wc=buttord(wp,ws,Ap,As,s); numa,dena=butter(N,wc,s); num,dend=impinvar(numa,dena,Fs); w=linspace(0,pi,1024); h=freqz(num,dend,w); norm=max(abs(h); plot(w/pi,20*log10(abs(h)/norm); xlabel(Normalized frequency); ylabel(Gain,dB);grid on; w=wpws; numd=num/norm; h=freqz(numd,dend,w);图像：M4-8(1)脉冲响应不变法代码：fp=2000;fs=10000; Ap=0.5; As=50; Fs=44100; wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; N,wc=buttord(wp,ws,Ap,As,s);numa,dena=butter(N,wc,s); num,dend=impinvar(numa,dena,Fs); w=linspace(0,2*pi,1024); h=freqz(num,dend,w); norm=max(abs(h); plot(w/pi,20*log10(abs(h)/norm);xlabel(Normalized frequency); ylabel(Gain,dB);grid on;图像：（2）双线性变换法代码：fp=2000; fs=10000;Ap=0.5; As=50;Fs=20000; wp=2*pi*fp/Fs; ws=2*pi*fs/Fs;N,wc=buttord(wp,ws,Ap,As,s); num,den=butter(N,wc,high); numc,dend=bilinear(num,den,Fs); w=linspace(0,2*pi,1024); h=freqz(numc,dend,w); norm=max(abs(h); plot(w/pi,20*log10(abs(h)/norm); xlabel(Normalized frequency); ylabel(Gain,dB);grid on;图像：M4-9代码：wp=0.8*pi;ws=0.6*pi;Ap=0.5;As=30; N,wc=buttord(wp,ws,Ap,As,s);num,den=butter(N,wc,high);numc,dend=bilinear(num,den,Fs);w=linspace(0,pi,2048); h=freqz(numc,dend,w); norm=max(abs(h); plot(w/pi,20*log10(abs(h)/norm); xlabel(Normalized frequency); ylabel(Gain,dB);grid on;图像：4、第五章M5-1M5-2(1)求增益响应代码：Wp=0.6*pi;Ws=0.4*pi;As=45;Ap=1; N=ceil(7*pi/(Wp-Ws); N=mod(N+1,2)+N; M=N-1; w=hamming(N); Wc=(Wp+Ws)/2; k=0:M; hd=-(Wc/pi)*sinc(Wc*(k-0.5*M)/pi); hd(0.5*M+1)=hd(0.5*M+1)+1;h=hd.*w; omega=linspace(0,pi,512); mag=freqz(h,1,omega); magdb=20*log10(abs(mag); plot(omega/pi,magdb); xlabel(Normalized frequency); ylabel(Gain, db); grid;图像：(2)代码：Wp=0.6*pi;Ws=0.4*pi;As=45;Ap=1; N=ceil(7*pi/(Wp-Ws); N=mod(N+1,2)+N; M=N-1; w=hamming(N); Wc=(Wp+Ws)/2; k=0:M; hd=-(Wc/pi)*sinc(Wc*(k-0.5*M)/pi); hd(0.5*M+1)=hd(0.5*M+1)+1; h=hd.*w; omega=linspace(0,pi,512); mag=freqz(h,1,omega); magdb=20*log10(abs(mag); %plot(omega/pi,magdb); xlabel(Normalized frequency); ylabel(Gain, db);grid; w=WpWs; H=freqz(h,1,w); fprintf(Ap=%.4fn,-20*log10(abs(H(1);fprintf(As=%.4fn,-20*log10(abs(H(2);结果：Ap=0.0182As=53.5634M5-3(1)blackman窗：代码：Wp=0.4*pi;Ws=0.6*pi;Ap=0.5;As=45; N=ceil(11.4*pi/(Ws-Wp); N=mod(N+1,2)+N; M=N-1; w=blackman(N); wc=(Ws+Wp)/2; k=0:M; hd=-(wc/pi)*sinc(wc*(k-0.5*M)/pi); h=hd.*w; n=linspace(0,pi,512); mag=freqz(h,1,n); magdb=20*log10(abs(mag); plot(n/pi,magdb); xlabel(normalized frequency);ylabel(Gain(dB); title(blackman); grid图像：(2)kaizer窗（a为通，阻带。Ftype代表返回滤波器的类型，f代表频带，Rs代表阻带荡）代码：wp=0.4*pi; ws=0.6*pi; As=45;Ap=0.5; Rs=0.01; f=0.4,0.6; a=1,0; dev=Rs*ones(1,length(a);M,Wc,beta,ftype = kaiserord(f,a,dev); M=mod(M,2)+M; h = fir1(M,Wc,ftype,kaiser(M+1,beta); omega=linspace(0,pi,512); mag=freqz(h,1,omega); plot(omega/pi,20*log10(abs(mag); xlabel(Normalized frequency); ylabel(Gain, db);legend(Kaiser);grid; fprintf(ftype=%sn,ftype); fprintf(Order of the filter=%.0fn,M); fprintf(Wc=%.4fn,Wc);fprintf(beta=%.4fn,beta);图像：M5-4M5-5M5-6代码：M=44; wp1=0.3*pi;wp2=0.5*pi; m=0:M/2; Wm=2*pi*m./(M+1);mtr1=ceil(wp1*(M+1)/(2*pi) ;mtr2=floor(wp2*(M+1)/(2*pi)+2; Ad=Wm=wp1; Ad(mtr1)=0.38; Ad(mtr2)=0.38; Hd=Ad.*exp(-j*0.5*M*Wm); Hd=Hdconj(fliplr(Hd(2:M/2+1); h=real(ifft(Hd); w=linspace(0,pi,1000);H=freqz(h,1,w); plot(w/pi,20*log10(abs(H);axis(0 1 -100 10); grid; xlabel(Normalized frequency);ylabel(Gain,dB);图像：M5-7M5-8M5-9M5-105、第六章M6-1M6-2代码：figure(1) a=0.1; k=50; n=100:200; x(n)=cos(0.3*pi*n)+0.5*cos(0.6*pi*n); y(n)=x(n)+a*x(n-k);stem(n,y(100:200);figure(2) t=-199:199; R=xcorr(y);stem(t,R);figure(3) clear y(n) hold on k=40:5:60; d=r g b c m; for e=1:length(k) y(n)=x(n)+a*x(n-k(e); R=xcorr(y); stem(t,R,d(e); endfigure(4) hold on clear y(n) k=0:1:60; for e=1:length(k) y(n)=x(n)+a*x(n-k(e); R=xcorr(y);stem(k(e),max(R); endfigure(5) clear y(n) hold on k=50; a=0:0.05:0.2; for e=1:length(a) y(n)=x(n)+a(e)*x(n-k); R=xcorr(y); stem(t,R,d(e); endfigure(6) hold on a=0:0.01:0.2; for e=1:length(a) y(n)=x(n)+a(e)*x(n-k); R=xcorr(y);stem(a(e),max(R); end图像：6、第七章M7-1代码：x = firls(511,0 0.1 0.1 0.5 0.5 1,1 1 1 0 0 0); f=0.1 0.4;a=1 0;dev=0.01 0.001; N,fo,ao,w = remezord(f,a,dev); h = remez(N,fo,ao,w);xd=filter(h,1,x);y=xd(1:4:end); w=linspace(0,pi,512); mag=freqz(h,1,w); subplot(3,1,1);plot(w/pi,20*log10(abs(mag); axis(0,1,-100,20);ylabel(Gain, dB);xlabel(Normalized frequency);legend(FIR filter, );grid; X=freqz(x,1,w);Y=freqz(y,1,w); subplot(3,1,2);plot(w/pi,(abs(X); xlabel(Normalized frequency);ylabel(magnitude);legend(Spectrum of x); grid; subplot(3,1,3);plot(w/pi,abs(Y); xlabel(Normalized frequency);ylabel(magnitude);legend(Spectrum of y);grid;图像：M7-2代码：x = firls(511,0 0.4 0.4 0.8 0.8 1,1 1 1 0 0 0); f=0.2 0.3;a=1 0;dev=0.01 0.001; xi=zeros(1,4*length(x); xi(1:4:end)=x; y=interp(x,4); w=linspace(0,pi,512); X=freqz(x,1,w);XI=freqz(xi,1,w);Y=freqz(y,1,w); subplot(3,1,1);plot(w/pi,20*log10(abs(X); axis(0,1,-100,20); ylabel(Gain, dB);xlabel(Normalized frequency); grid; subplot(3,1,2);plot(w/pi,(abs(XI); xlabel(Normalized frequency);ylabel(magnitude); grid; subplot(3,1,3);plot(w/pi,abs(Y);xlabel(Normalized frequency);ylabel(magnitude); grid;图像：M7-3代码：Fp=0.45;Fs=0.55;Wp=Fp*pi;Ws=Fs*pi; dp=0.001; N,fo,ao,w= remezord(Fp Fs,1 0,dpdp); n=mod(N+1,4);N=N+3-n; h= remez(N,fo,ao,w); h_half=zeros(1,N+1);h_half(1:2:N+1)=h(1:2:N+1); h_half(N/2+1)=1/2; w=linspace(0,pi,512);H_half=20*log10(abs(freqz(h_half,1,w); subplot(2,1,1);plot(h_half); subplot(2,1,2);plot(w/pi,H_half); grid;图像：M7-4代码：x = firls(511,0 0.5 0.5 1,1 1 1 0);w=linspace(-pi,pi,512); X=freqz(x,1,w); plot(w/pi,abs(X);图像：M7-5(1)代码： h0,h1,g0,g1 = firpr2chfb(13,.4);fvtool(h0,1,h1,1,g0,1,g1,1);图像：(2)代码： h0,h1,g0,g1 = firpr2chfb(13,.4); h2=conv(h0,h0); h3=conv(h1,h1); H0=freqz(h2,1,w);H1=freqz(h3,1,w); plot(w/pi,10*log10(abs(H0+H1)图像：第72页