參與評論

評論

登錄后參與評論

全部評論(2條)

• 

  如夢*粱 2010-12-10 00:00:00

  clear all close all fedge=[0 .175 .2 .4 .425 .575 .6 .8 .825 1]; mag=[0 0 1 1 0 0 1 1 0 0]; b=remez(45，fedge，mag); [h1，w1]=freqz(b，1，512); plot(w1/pi，20*log10(abs(h1))); grid on; axis([0，1，-80，20]); xlabel('歸一化頻率/π'); ylabel('幅度/dB');

  贊(2)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論

• 

  qxxheart 2018-04-09 23:58:19

  說明：（1）為了使濾波器階數(shù)盡可能低，每個濾波器的邊界頻率選擇原則是盡量使濾波器過渡帶寬盡可能寬。 （2）與信號產(chǎn)生函數(shù)mstg相同，采樣頻率Fs=10kHz。 （3）為了濾波器階數(shù)Z低，選用橢圓濾波器。（之后，再依次實現(xiàn)巴特沃斯、切比雪夫1、切比雪夫2數(shù)字濾波器） 按照圖2 所示的程序框圖編寫的實驗程序為exp1.m。 2、實驗程序清單 %實驗1程序exp1.m % IIR數(shù)字濾波器設(shè)計及軟件實現(xiàn) clear all;close all Fs=10000;T=1/Fs; %采樣頻率 %調(diào)用信號產(chǎn)生函數(shù)mstg產(chǎn)生由三路YZ載波調(diào)幅信號相加構(gòu)成的復(fù)合信號st st=mstg; %低通濾波器設(shè)計與實現(xiàn)========================================= fp=280;fs=450; wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指標(biāo)（低通濾波器的通、阻帶邊界頻） [N，wp]=ellipord(wp，ws，rp，rs); %調(diào)用ellipord計算橢圓DF階數(shù)N和通帶截止頻率wp [B，A]=ellip(N，rp，rs，wp); %調(diào)用ellip計算橢圓帶通DF系統(tǒng)函數(shù)系數(shù)向量B和A y1t=filter(B，A，st); %濾波器軟件實現(xiàn) % 低通濾波器設(shè)計與實現(xiàn)繪圖部分 figure(5); subplot(2，1，1); myplot(B，A); %調(diào)用繪圖函數(shù)myplot繪制損耗函數(shù)曲線 yt='y_1(t)'; subplot(2，1，2); tplot(y1t，T，yt); %調(diào)用繪圖函數(shù)tplot繪制濾波器輸出波形 %帶通濾波器設(shè)計與實現(xiàn)==================================================== fpl=440;fpu=560;fsl=275;fsu=900; wp=[2*fpl/Fs，2*fpu/Fs];ws=[2*fsl/Fs，2*fsu/Fs];rp=0.1;rs=60; [N，wp]=ellipord(wp，ws，rp，rs); %調(diào)用ellipord計算橢圓DF階數(shù)N和通帶截止頻率wp [B，A]=ellip(N，rp，rs，wp); %調(diào)用ellip計算橢圓帶通DF系統(tǒng)函數(shù)系數(shù)向量B和A y2t=filter(B，A，st); %濾波器軟件實現(xiàn) % 帶通濾波器設(shè)計與實現(xiàn)繪圖部分 figure(3); subplot(2，1，1); myplot(B，A); %調(diào)用繪圖函數(shù)myplot繪制損耗函數(shù)曲線 yt='y_2(t)'; subplot(2，1，2); tplot(y2t，T，yt); %調(diào)用繪圖函數(shù)tplot繪制濾波器輸出波形 %高通濾波器設(shè)計與實現(xiàn)================================================ fp=890;fs=600; wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指標(biāo)（低通濾波器的通、阻帶邊界頻） [N，wp]=ellipord(wp，ws，rp，rs); %調(diào)用ellipord計算橢圓DF階數(shù)N和通帶截止頻率wp [B，A]=ellip(N，rp，rs，wp，'high'); %調(diào)用ellip計算橢圓帶通DF系統(tǒng)函數(shù)系數(shù)向量B和A y3t=filter(B，A，st); %濾波器軟件實現(xiàn) % 高低通濾波器設(shè)計與實現(xiàn)繪圖部分 figure(4); subplot(2，1，1); myplot(B，A); %調(diào)用繪圖函數(shù)myplot繪制損耗函數(shù)曲線 yt='y_3(t)'; subplot(2，1，2); tplot(y3t，T，yt); %調(diào)用繪圖函數(shù)tplot繪制濾波器輸出波形 function myplot(B，A) %時域離散系統(tǒng)損耗函數(shù)繪圖 %B為系統(tǒng)函數(shù)分子多項式系數(shù)向量 %A為系統(tǒng)函數(shù)分母多項式系數(shù)向量 [H，W]=freqz(B，A，1000); m=abs(H); plot(W/pi，20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅度(dB)') axis([0，1，-80，5]);title('損耗函數(shù)曲線'); function tplot(xn，T，yn) %時域序列連續(xù)曲線繪圖函數(shù) % xn:信號數(shù)據(jù)序列，yn:繪圖信號的縱坐標(biāo)名稱（字符串） % T為采樣間隔 n=0:length(xn)-1;t=n*T; plot(t，xn); xlabel('t/s');ylabel(yn); axis([0，t(end)，min(xn)，1.2*max(xn)]) function st=mstg N=2000 Fs=10000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; fm1=fc1/10; fc2=Fs/20; fm2=fc2/10; fc3=Fs/40; fm3=fc3/10; xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); st=xt1+xt2+xt3; fxt=fft(st，N); subplot(3，1，1) plot(t，st);grid;xlabel('t/s');ylabel('s(t)'); axis([0，Tp/8，min(st)，max(st)]);title('(a) s(t)的波形') subplot(3，1，2) stem(f，abs(fxt)/max(abs(fxt))，'.');grid;title('(b) s(t)的頻譜') axis([0，F(xiàn)s/5，0，1.2]); xlabel('f/Hz');ylabel('幅度');

  贊(14)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論