基于頻域抽樣法的fir數(shù)字帶阻濾波器-資料下載頁
2025-06-22 21:15本頁面
【正文】 譜圖39。)。xlabel(39。頻率/HZ39。)。ylabel(39。幅度39。)程序2wsl=*pi。%低阻帶邊緣wsh=*pi。%高阻帶邊緣wpl=*pi。%低通帶邊緣wph=*pi。%高通帶邊緣delta=(wplwsl)。%過度帶N=ceil(2*pi*3/delta)1。%抽樣點數(shù)al=(N1)/2。wl=(2*pi/N)。 %抽樣間隔k=0:N1。T1=。 T2=。%過渡帶樣本點%建立符幅特性樣本序列Hrs=[ones(1,ceil(*pi/wl)+1),T2,T1,zeros(1,ceil(*pi/wl)),T1,T2,ones(1,ceil(*pi/wl)),T2,T1,zeros(1,ceil(*pi/wl)),T1,T2,ones(1,ceil(*pi/wl)+1)]。wdl=[0 1]。k1=0:floor((N1)/2)。k2=floor((N1)/2)+1:N1。 %k分為兩部分k1=0~(N1)/2,k2=(N1)/2~N1angH=[al*(2*pi)/N*k1,al*(2*pi)/N*(Nk2)]。%建立相位特性樣本序列 求H(k)的相角θ(k)H=Hrs.*exp(j*angH)。%建立頻率特性樣本序列 求H(k)h=real(ifft(H))。%傅立葉反變換求脈沖序列由H(k)傅立葉反變換得到h%沖激響應圖figure(1)subplot(211)。stem(k,h)。title(39。沖激響應39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。grid。%幅頻曲線圖Hf=abs(H)。w=k*wl/pi。subplot(212)。plot(w,Hf,39。*b39。)axis([0 1 ])。title(39。幅頻響應39。)。xlabel(39。頻率/π39。)。ylabel(39。Hr(w)39。)。set(gca,39。xtickmode39。,39。manual39。,39。xtick39。,wdl)。 %設置所要顯示數(shù)據(jù)的位置set(gca,39。ytickmode39。,39。manual39。,39。ytick39。,[0 1])。 %給這些數(shù)據(jù)加標簽，用來在指定范圍內標注網(wǎng)格線grid。%相頻響應曲線pha=unwrap(angle(H))。figure(2)subplot(211)。plot(w,pha)。 %相頻曲線title(39。相頻響應39。)。xlabel(39。39。)。ylabel(39。39。)。grid。fs=15000。[c,f3]=freqz(h,1)。f3=f3/pi*fs/2。subplot(212)。plot(f3,20*log10(abs(c)))。title(39。頻譜特性39。)。xlabel(39。頻率/HZ39。)。ylabel(39。衰減 39。)。grid。t=(0:100)/fs。x=sin(2*pi*t*1000)+sin(2*pi*t*3000)+sin(2*pi*t*6000)。q=filter(h,1,x)。[a,f1]=freqz(x)。f1=f1/pi*fs/2。[b,f2]=freqz(q)。f2=f2/pi*fs/2。figure(3)subplot(211)。plot(f1,abs(a))。title(39。輸入信號頻譜圖39。)。xlabel(39。頻率/HZ39。)。ylabel(39。dB39。)subplot(212)。plot(f2,abs(b))。title(39。輸出信號頻譜圖39。)。xlabel(39。頻率/HZ39。)。ylabel(39。dB39。)15 / 16
環(huán)評公示相關推薦
鄂ICP備17016276號-1