【正文】 0。Fp=4000。Fs=3500。wp1=tan(pi*Fp/Ft)。 %高通到低通濾波器參數轉換ws1=tan(pi*Fs/Ft)。wp=1。ws=wp1*wp/ws1。[n13,wn13]=cheb1ord(wp,ws,1,50,39。s39。)。 %求模擬的低通濾波器階數和截止頻率[b13,a13]=cheby1(n13,1,wn13,39。s39。)。 %求S域的頻率響應的參數[num,den]=lp2hp(b13,a13,wn13)。 %將S域低通參數轉為高通的[num13,den13]=bilinear(num,den,)。 %利用雙線性變換實現S域到Z域轉換[h,w]=freqz(num13,den13)。[y,fs,nbits]=wavread (39。C:\program files\matlab71\work\speech\39。)。n = length (y) 。 %求出語音信號的長度s=y。 S=fft(s)。 %傅里葉變換z13=filter(num13,den13,s)。sound(z13)。m13=fft(z13)。 %求濾波后的信號axes()。plot(abs(S),39。g39。)。title(39。濾波前信號的頻譜39。)。axes()。plot(abs(m13),39。r39。)。title(39。濾波后信號的頻譜39。)。axes()。plot(s)。 title(39。濾波前信號的波形39。)。axes()。plot(z13)。title(39。濾波后的信號波形39。)。合成信號低通濾波f1=20。f2=200。f3=500。t=(1:100)/2000。x1=sin(2*pi*t*f1)。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。x3=sin(2*pi*t*f3)。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。n=[1:100]。t=n/2000X=fft(x,512)。w=(0:255)/256*1000。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。Ft=8000。Fp=1000。Fs=1200。As=100 。Ap=1。wp=2*pi*Fp/Ft。ws=2*pi*Fs/Ft。fp=2*Fp*tan(wp/2)。fs=2*Fs*tan(ws/2)。[n11,wn11]=buttord(wp,ws,1,50,39。s39。)。 [b11,a11]=butter(n11,wn11,39。s39。)。 [num11,den11]=bilinear(b11,a11,)。 [h,w]=freqz(num11,den11)。 z11=filter(num11,den11,x)。m11=fft(z11)。 axes()。plot(abs(x),39。g39。)。title(39。濾波前信號的頻譜39。)。axes()。plot(abs(m11),39。r39。)。title(39。濾波后信號的頻譜39。)。axes()。plot(x)。 title(39。濾波前信號的波形39。)。axes()。plot(z11)。title(39。濾波后的信號波形39。)。合成信號帶通濾波f1=20。f2=200。f3=500。t=(1:100)/2000。x1=sin(2*pi*t*f1)。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。x3=sin(2*pi*t*f3)。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。n=[1:100]。t=n/2000X=fft(x,512)。w=(0:255)/256*1000。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。Ft=8000。Fp1=1200。Fp2=3000。Fs1=1000。Fs2=3200。As=100。Ap=1。wp1=tan(pi*Fp1/Ft)。 wp2=tan(pi*Fp2/Ft)。ws1=tan(pi*Fs1/Ft)。ws2=tan(pi*Fs2/Ft)。w=wp1*wp2/ws2。bw=wp2wp1。wp=1。ws=(wp1*wp2w.^2)/(bw*w)。[n12,wn12]=buttord(wp,ws,1,50,39。s39。)。 [b12,a12]=butter(n12,wn12,39。s39。)。 [num2,den2]=lp2bp(b12,a12,sqrt(wp1*wp2),bw)。[num12,den12]=bilinear(num2,den2,)。[h,w]=freqz(num12,den12)。 z12=filter(num12,den12,x)。m12=fft(z12)。 axes()。plot(abs(x),39。g39。)。title(39。濾波前信號的頻譜39。)。axes()。plot(abs(m12),39。r39。)。title(39。濾波后信號的頻譜39。)。axes()。plot(x)。 title(39。濾波前信號的波形39。)。axes()。plot(z12)。title(39。濾波后的信號波形39。)。合成信號高通濾波f1=20。f2=200。f3=500。t=(1:100)/2000。x1=sin(2*pi*t*f1)。 %繪制x(t)的圖形x2=sin(2*pi*t*f2)。x3=sin(2*pi*t*f3)。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。n=[1:100]。t=n/2000X=fft(x,512)。w=(0:255)/256*1000。x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*t*f3)。Ft=8000。Fp=4000。Fs=3500。wp1=tan(pi*Fp/Ft)。 %高通到低通濾波器參數轉換ws1=tan(pi*Fs/Ft)。wp=1。ws=wp1*wp/ws1。[n13,wn13]=cheb1ord(wp,ws,1,50,39。s39。)。 %求模擬的低通濾波器階數和截止頻率[b13,a13]=cheby1(n13,1,wn13,39。s39。)。 %求S域的頻率響應的參數[num,den]=lp2hp(b13,a13,wn13)。 %將S域低通參數轉為高通的[num13,den13]=bilinear(num,den,)。 %利用雙線性變換實現S域到Z域轉換[h,w]=freqz(num13,den13)。z13=filter(num13,den13,x)。sound(z13)。m13=fft(z13)。 %求濾波后的信號axes()。plot(abs(x),39。g39。)。title(39。濾波前信號的頻譜39。)。axes()。plot(abs(m13),39。r39。)。title(39。濾波后信號的頻譜39。)。axes()。plot(x)。 title(39。濾波前信號的波形39。)。axes()。plot(z13)。title(39。濾波后的信號波形39。)。28