基于頻率抽樣法二維非遞歸濾波器設(shè)計及其簡單應(yīng)-資料下載頁
2025-06-27 21:17本頁面
【正文】 1),39。o39。,wdl,Hdr)。axis([0,1,])。title(39。Frequency Samples:M=4039。)xlabel(39。frequency in pi units39。)。ylabel(39。Hr(k)39。)subplot(2,2,2)。stem(l,h)。axis([1,M,])title(39。Impulse Response39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。subplot(2,2,3)。plot(ww/pi,Hr,wl(1:21)/pi,Hrs(1:21),39。o39。)。axis([0,1,])。title(39。Amplitude Response39。)xlabel(39。frequency in pi units39。)。ylabel(39。Hr(w)39。)subplot(2,2,4)。plot(w/pi,db)。axis([0,1,120,10])。gridtitle(39。MagnitudeResponse39。)。xlabel(39。frequencyinpiunits39。)。ylabel(39。Decibels39。)。%最優(yōu)設(shè)計法T1=clear all。M=40。alpha=(M1)/2。 l=0:M1。wl=(2*pi/M)*l。Hrs=[ones(1,5),zeros(1,29),ones(1,4)]。Hdr=[1,1,0,0]。wdl=[0,1]。k1=0:floor((M1)/2)。k2=floor((M1)/2)+1:M1。angH=[alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(Mk2)]。H=Hrs.*exp(j*angH)。h=real(ifft(H,M))。[db,mag,pha,grd,w]=freqz_m(h,1)。[Hr,ww,a,L]=Hr_Type2(h)。subplot(2,2,1)。plot(wl(1:21)/pi,Hrs(1:21),39。o39。,wdl,Hdr)。axis([0,1,])。title(39。Frequency Samples:M=4039。)xlabel(39。frequency in pi units39。)。ylabel(39。Hr(k)39。)subplot(2,2,2)。stem(l,h)。axis([1,M,])title(39。Impulse Response39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。subplot(2,2,3)。plot(ww/pi,Hr,wl(1:21)/pi,Hrs(1:21),39。o39。)。axis([0,1,])。title(39。Amplitude Response39。)xlabel(39。frequency in pi units39。)。ylabel(39。Hr(w)39。)subplot(2,2,4)。plot(w/pi,db)。axis([0,1,120,10])。gridtitle(39。Magnitude Response39。)。xlabel(39。frequency in pi units39。)。ylabel(39。Decibels39。)。％增加兩個過渡帶clear all。M=60。alpha=(M1)/2。 l=0:M1。wl=(2*pi/M)*l。Hrs=[ones(1,7),zeros(1,43),ones(1,6)]。Hdr=[1,1,0,0]。wdl=[0,1]。k1=0:floor((M1)/2)。k2=floor((M1)/2)+1:M1。angH=[alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(Mk2)]。H=Hrs.*exp(j*angH)。h=real(ifft(H,M))。[db,mag,pha,grd,w]=freqz_m(h,1)。[Hr,ww,a,L]=Hr_Type2(h)。subplot(2,2,1)。plot(wl(1:31)/pi,Hrs(1:31),39。o39。,wdl,Hdr)。axis([0,1,])。title(39。Frequency Samples:M=6039。)xlabel(39。frequency in pi units39。)。ylabel(39。Hr(k)39。)subplot(2,2,2)。stem(l,h)。axis([1,M,])title(39。Impulse Response39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。subplot(2,2,3)。plot(ww/pi,Hr,wl(1:31)/pi,Hrs(1:31),39。o39。)。axis([0,1,])。title(39。Amplitude Response39。)xlabel(39。frequency in pi units39。)。ylabel(39。Hr(w)39。)subplot(2,2,4)。plot(w/pi,db)。axis([0,1,120,10])。gridtitle(39。MagnitudeResponse39。)。xlabel(39。frequencyinpiunits39。)。ylabel(39。Decibels39。)。％產(chǎn)生20階具有最小二乘線性相位FIR濾波器f=[0 1]。 m=[1 1 1 1 0 0 1 1 1 1]。b=firls(20, f , m, 39。Hilbert39。)。[h0,w]=freqz(b,1,512,2)。plot(f,m,39。39。,w, abs(h0))。h1=fsamp2(f,f,m,[5 5])。figure(2)freqz2(h1,[32 32])。 ff=1/512:1/512:1。h2=fsamp2(ff,ff,h039。,[5 5])。figure(3)freqz2(h2,[32 32])。 f=imread(39。39。)。％利用濾波器h1濾波fh1=filter2(h1,f)。figure(4)imshow(fh1)。 ％利用濾波器h2濾波fh2=filter2(h2,f)。figure(5)imshow(fh2)。22
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1