頻譜校正與濾波器設(shè)計matlab程序

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)頻譜校正與濾波器設(shè)計matlab程序1、數(shù)字濾波器Matlab程序代碼如下：wc=(400/1000)*pi; w1=boxcar(81); w2=hamming(81);w3=hanning(81);n=1:1:81;hd=sin(wc*(n-41)./(pi*(n-41); hd(41)=wc/pi;h1=hd.*w1; h2=hd.*w2;h3=hd.*w3;mag1,rad=freqz(h1); mag2,rad=freqz(h2);mag3,rad=freqz

2、(h3);figure(1); plot(rad,20*log10(abs(mag1);xlabel(Normalized Frequency(rad);ylabel(Normaliaed Magnitude(dB);axis(0,3,-80,0);title(矩形窗濾波器);grid on;figure(2);plot(rad,20*log10(abs(mag2);xlabel(Normalized Frequency(rad);ylabel(Normaliaed Magnitude(dB);axis(0,3,-80,0);title(hamming窗濾波器);grid on;figure(

3、3);plot(rad,20*log10(abs(mag3);xlabel(Normalized Frequency(rad);ylabel(Normaliaed Magnitude(dB);axis(0,3,-80,0);title(hanning窗濾波器);grid on;輸出結(jié)果如下：2、頻譜校正（比值校正法）Matlab程序如下：functionZ=bizhi(x,fs,N,nx1,nx2,method)nx,mx=size(x);if mx=1, x=x;endnx,mx=size(x);if mxA(index+1); df=indsecL.*A(index-1)./(Amax+A

4、(index-1)-(1-indsecL).*A(index+1)./(Amax+A(index+1); Z(1)=(index-1-df)*ddf; Z(2)=Amax/sinc(df); Z(3)=(phmax+pi*df); end if (WindowType=2) indsecL=A(index-1)A(index+1); df=indsecL.*(2*A(index-1)-Amax)./(Amax+A(index-1)-(1-indsecL).*( 2*A(index+1)-Amax)./(Amax+A(index+1); Z(1)=(index-1-df)*ddf; Z(2)=(

5、1-df2)*Amax/sinc(df); Z(3)=(phmax+pi*df); end Z(3)=mod(Z(3),2*pi); Z(3)=Z(3)-(Z(3)pi)*pi;設(shè)計例子如下：Fs=1000;n=0:1/Fs:1;xn=10+15*sin(10*2*pi*n)+randn(size(n);nfft=1024;Xn = fftshift(fft(xn,nfft);Xn(513)=Xn(513)/2;A=2*abs(Xn)/Fs;subplot(2,1,1),plot(xn)subplot(2,1,2)stem(-nfft/2:nfft/2 - 1*Fs/nfft,A)xlabel(頻率),ylabel(幅值)P=Xn.*conj(Xn)/nfft;xlim(0 500)Z=bizhi(xn,Fs,nfft,8,12,2)輸出結(jié)果如下這樣得到Z為10