應(yīng)用matlab對圖像信號進(jìn)行頻譜分析及濾波(編輯修改稿)
2025-07-20 17:21 本頁面
【文章內(nèi)容簡介】 階數(shù)以及3dB截止頻率[b,a]=butter(Nm,Wc)。 %計(jì)算濾波器的分子分母系數(shù)H=freqz(b,a,f*2*pi/fs)。%計(jì)算濾波器頻率響應(yīng)mag=abs(H)。pha=angle(H)。%幅度響應(yīng)和相位響應(yīng)mag1=20*log((mag+eps)/max(mag))。 %將幅頻響應(yīng)轉(zhuǎn)化為dB的形式%%用IIR濾波器對帶噪序列進(jìn)行濾波并求頻率響應(yīng)x6=filter(b,a,x4)。%用IIR濾波Xk6=fft(x6)。Xf6=fftshift(Xk6)。x7=zeros(M,N)。for i=1:M for j=1:N x7(i,j)=x6(N*(i1)+j)。 endend%%設(shè)計(jì)FIR濾波器并分析相關(guān)指標(biāo)wc=280*2/fs。 %6dB截止頻率280kHzfx=[0 wc wc 1]。m=[1 1 0 0]。%理想頻幅響應(yīng)b1=fir2(40,fx,m,hamming(41))。%計(jì)算FIR濾波器多項(xiàng)式系數(shù)（不聲明窗默認(rèn)為Hamming窗）H1=freqz(b1,1,f*2*pi/fs)。%計(jì)算濾波器頻率響應(yīng)mag2=abs(H1)。pha1=angle(H1)。%幅度響應(yīng)和相位響應(yīng)mag3=20*log((mag2+eps)/max(mag2))。 %將幅頻響應(yīng)轉(zhuǎn)化為dB的形式%%用FIR濾波器對帶噪序列進(jìn)行濾波并求頻率響應(yīng)x8=filter(b1,1,x4)。%用FIR進(jìn)行濾波grd=grpdelay(b1,1,f*2*pi/fs)。%計(jì)算群延時(shí)K=round(grd(1))。%修正因子(冒號操作做索引時(shí)，需要整型數(shù)操作)x8=[x8((K+1):L),x8(1:K)]。%對群延遲進(jìn)行修正Xk8=fft(x8)。Xf8=fftshift(Xk8)。x9=zeros(M,N)。for i=1:M for j=1:N x9(i,j)=x8(N*(i1
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1