{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

Kaiser 设计FIR低通滤&a

Kaiser 设计FIR低通滤&a

Info icon This preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
DSP ? ? ? ? Kaiser ? ? FIR ? ? ? ? ? ( ? ? ? ? ) F S = 5000Hz F Pass = 600Hz F Stop = 800Hz R P = 2dB R S = 60dB Matlab ? ? : >> Rp = 2; >> Rs = 60; >> Fp = 600; >> Fstop = 800; >> Fs = 5000; >> devs = [(10^(Rp/20)-1)/(10^(Rp/20)+1) 10^(-Rs/20)]; >> alfa = -min(20*log10(devs)); >> if alfa>50, beta = .5842*((alfa-21).^(.4))+0.07886*(alfa-21);
Image of page 1

Info icon This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document Right Arrow Icon
Image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: else beta = 0; end >> n = ceil((alfa-8)/(2.285*(Fstop-Fp)*2*pi/Fs)); >> n = max(1,n); >> h = fir1(n,(Fp+(Fstop-Fp)/2)*2/Fs,kaiser(n+1,beta)); >> H = fft(h,2048); >> k = 0:2047; >> plot(k/1024,20*log10(abs(H)));grid on; ? ? ? ? ? ? ? ? : ? ? 91...
View Full Document

{[ snackBarMessage ]}