【正文】 其次，通過這次的課設，讓我初次接觸了 MATLAB，學會了用它來分析和處理信號，也讓我為以后在學習和生活中需要運用它解決問題打下了良好的基礎。這次的課程設計讓我對以后的設計有了進一步的了解，也讓我知道平時的學習范圍太過 于狹隘，課外的學習必不可少。如果在 S平面上去逼近，就得到模擬濾波器；如果在 z平面上去逼近，就得到數字濾波器。 窗函數法分析： 窗函數法：相位響應有嚴格 的線性 ,不存在穩(wěn)定性問題 , 設計簡單。)。 subplot(2,2,4)。plot(z1)。 title(39。)。plot(abs(Y1))。 x=filter(b,a,z1)。stop39。 wc=[2*fc1/Fs,2* fc2/Fs]。fc2=3200 。 Y1=fft(y1)。 圖形分析如下： ②選用 cheby1 程序設計如下： clear。 title(39。濾波前信號波形 39。)。plot(abs(X))。 title(39。 figure(2)。 figure(1)。 [N,ws]=buttord(wc,wp,Ap,As)。Ap=1。fp2=3000。d:\39。)。 subplot(2,2,4)。plot(z1)。 title(39。)。plot(abs(Y1))。 x=filter(b,a,z1)。high39。 wc=2*fc/Fs。fp=3000 。d:\39。)。 subplot(2,2,4)。plot(z1)。 title(39。)。plot(abs(Y1))。 x=filter(b,a,z1)。high39。 wc=2*fc/Fs。fp=5000 。d:\39。)。 subplot(2,2,4)。plot(z1)。 title(39。)。plot(abs(Y1))。 x=filter(b,a,z1)。 [n,wp]=cheb1ord(wc,wb,Ap,As)。 。 fp=1000。close all [z1,fs,bits]=wavread(39。濾波前信號波形 39。)。 subplot(2,2,3)。axis([0,1000,0,])。濾波前信號頻譜 39。 subplot(2,2,1)。 freqz(b,a,512,Fs)。 wp=2*fp/Fs。As=100。) y1=z1(1:8192)。 sound(x,fs,bits)。plot(x)。 title(39。 title(39。)。plot(abs(Y1))。 freqz(b,1) x=fftfilt(b,z1)。 ws =[(wp1+wc1)/2/pi,(wp2+wc2)/2/pi]。 wc2=2*pi*fc2/Fs。Fs=8000 。fc1=1000 。) y1=z1(1:8192)。 sound(x,fs,bits)。plot(x)。 title(39。濾波后信號頻譜 39。 subplot(2,2,2)。axis([0,1000,0,])。 X=fft(x,8192)。,wn)。 wn= kaiser(N,beta)。 wp=2*pi*fp/Fs。As=100。) y1=z1(1:8192)。 sound(x,fs,bits)。plot(x)。 title(39。濾波后信號頻譜 39。 subplot(2,2,2)。axis([0,1000,0,])。 X=fft(x,8192)。 b=fir1(N,ws,wn)。 beta=*()。Fs=8000。 fp=1000。close all [z1,fs,bits]=wavread(39。 Y1=fft(y1)。 圖像輸出如圖 1 圖 1 語 音信號時域波形 ②對語音信號進行頻譜分析，在 MATLAB中，可 z1=wavread(39。 2.．語音信號的頻譜分析 ①首先畫出語音信號的時域波形 z1=wavread(39。首先，一個線性相位的模擬濾波器經雙線性變換后得到非線性相位的數字濾波器，不再保持原有的線性相位了；其次，這種非線性關系要求模擬濾波器的幅頻響應必須是分段常數型的，即某一頻率段的幅頻 響應近似 等于某一常數（這正是一般典型的低通、高通、帶通、帶阻型濾波器的響應特性），不然變換所產生的數字濾波器幅頻響應相對于原模擬濾波器的幅頻響應會有畸變，如圖 3所示。這個關系如式（ 4）所示，重寫如下： 上式表明， S平面上Ω與 Z平面的ω成非線性的正切關系，如圖 2所示。因此，穩(wěn)定的模擬濾波器經雙線性變換后所得的 數字濾波器也一定是穩(wěn)定的。將式（ 1）寫成 將此關系解析延拓到整個 S平面和 S1平面，令 jΩ =s， jΩ 1=s1，則得 再將 S1平面通過以下標準變換關系映射到 Z平面 z=es1T 從而得到 S平面和 Z平面的單值映射關系為： (2) (3) 式（ 2）與式（ 3）是 S平面與 Z平面之間的單值映射關系，這種變換都是兩個線性函數之比，因此稱為雙線性變換 式（ 1）與式（ 2）的 雙線性變換符合映射變換應滿足的兩點要求。也就是說，第一步先將整個 S 平面壓縮映射到 S1 平面的 π /T～π /T 一條橫帶里；第二步再通過標準變換關系z=es1T將此橫帶變換到整個 Z 平面上去。 驗算技術指標是否滿足要求 ，為了計算數字濾波器在頻域中的特性，可調用 freqz子程序，如果不滿足要求，可根據具體情況，調整窗函數類型或長