HOW TO DESIGN AN IIR Low PASS FILTER WITH MATLAB
21 次查看(过去 30 天)
显示 更早的评论
I have asked this type of questions many times but can you give me the matlab code for designing a lowpass IIR filter.
thankyou in advance
0 个评论
回答(2 个)
Wayne King
2012-5-10
Hi Raj, you can either use basic functions like butter(), cheby1()
For example:
Lowpass filter for data sampled at 10 kHz, passband below 1 kHz
Wc = (2*1e3)/1e4;
[B,A] = butter(10,Wc);
% view magnitude response
fvtool(B,A,'Fs',1e4)
Or something like:
[B,A] = cheby1(10,0.5,Wc);
% view magnitude response
fvtool(B,A,'Fs',1e4)
Or you can use the fdesign workflow. Using fdesign.lowpass
Fs = 1e4;
d = fdesign.lowpass('N,F3dB',10,1000,Fs);
Hd = design(d,'butter');
fvtool(Hd)
There are a number of specification strings for fdesign.lowpass that support IIR designs. After you specify a filter, you can use
designmethods(d)
to see which design methods are supported.
2 个评论
Wayne King
2012-5-10
You can easily design an FIR filter, or use the impulse response from the filters above, and then use linear prediction on those coefficients (or impulse response). The problem you are going to have is choosing the order. It won't take a very large AR order at all before you start getting a peaky response, which isn't going to match your filter's frequency response very well at all.
Wc = (2*1e3)/1e4;
[B,A] = butter(10,Wc);
h = impz(B,A);
A1 = lpc(h,2);
fvtool(1,A1,'Fs',1e4);
Or starting with an FIR filter
d = fdesign.lowpass('Fp,Fst,Ap,Ast',1e3,1.1e3,0.5,40,Fs);
Hd = design(d);
Afilt = lpc(Hd.Numerator,1);
fvtool(1,Afilt,'Fs',1e4)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Digital Filter Design 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!