課程設(shè)計---基于雙線性變換法的iir數(shù)字低通濾波器設(shè)計-資料下載頁
2025-01-17 04:06本頁面
【正文】 和操作功能，其程序編寫也比其他語言要簡單。參考文獻[1] 蔣志凱. 數(shù)字濾波與卡爾曼濾波[M].北京：中國科學(xué)技術(shù)出版社，1993[2] 施陽．等．MATLAB語言工具箱[M]..西安：西北工業(yè)大學(xué)出版社,1999[3] 樓順天、李博菡. 基于MATLAB的系統(tǒng)分析與設(shè)計－信號處理[M].西安：西安電子科技大學(xué)出版社，1998.[4] 蒙以正. [M].北京：科學(xué)出版社，1999.[5] 胡廣書. 數(shù)字信號處理：理論、算法與實現(xiàn)[M].北京：清華大學(xué)出版社，西安：1997.[6] 倪養(yǎng)華．?dāng)?shù)字信號處理與實現(xiàn)．上海：上海交通大學(xué)出版社,1998[7] 施陽. MATLAB語言工具箱．西安：西北工業(yè)大學(xué)出版[8] 邊肇祺. 模式識別（第二版）.北京： 清華大學(xué)出版社,1988，25～35[9] 胡廣書. 數(shù)字信號處理—理論、算法與實現(xiàn)[M]. 北京:清華大學(xué)出版社,1997.附錄MATLAB程序clear。clc。close all。[S,Fs,Bit] = wavread(39。39。)。 %讀取音頻信號sound(S,Fs)。 fnoise1 = 10000。 % 頻率為fnoise的正弦噪聲信號fnoise2 = 50。 % 頻率為fnoise的正弦噪聲信號N = length(S)。T= N/Fs。t=T/N:T/N:T。for j=1:N Noise1(j)=sin(2*pi*fnoise1*t(j))。 Noise2(j)=sin(2*pi*fnoise2*t(j))。end% 對于低通濾波器：S1 = S+*Noise139。fp = 3000。fs = 4000。Rp = 2。Rs = 45。 % 求數(shù)字頻率wp = 2*pi*fp/Fs。ws = 2*pi*fs/Fs。% 預(yù)畸變處理Wp = 2*Fs*tan(wp/2)。Ws = 2*Fs*tan(ws/2)。% 開始設(shè)計模擬低通濾波器[N,Wn] = buttord(Wp,Ws,Rp,Rs,39。s39。)。[B,A] = butter(N,Wn,39。s39。)。[b,a] = bilinear(B,A,Fs)。 %impinvar 沖激響應(yīng)不變法 [h,w] = freqz(b,a,1000,39。whole39。)。h = h(1:501)。w = w(1:501)。mag = abs(h)。db = 10*log10((mag+eps)/max(mag))。subplot(1,2,1)。plot(w/pi,mag)。title(39。幅頻響應(yīng)39。)。axis([0,1,0,])。subplot(1,2,2)。plot(w/pi,db)。title(39。幅頻響應(yīng)(dB)39。)。axis([0,1,40,1])。set(gca,39。XTickMode39。,39。manual39。,39。XTick39。,[0,1])。set(gca,39。YTickMode39。,39。manual39。,39。YTick39。,[40,20,3,0])。grid。 X1 = abs(fft(S))。Z1 = (0:150000)/150000。subplot(3,1,1)。plot(Z1,X1(1:150001))。 X2 = abs(fft(S1))。Z2 = (0:150000)/150000。subplot(3,1,2)。plot(Z2,X2(1:150001))。 S2 = filter(b,a,S1)。sound(S2,Fs)。X3=abs(fft(S2))。Z3 = (0:150000)/150000。subplot(3,1,3)。plot(Z3,X3(1:150001))。14
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號-1