Cannot deduce a transfer function! Help!

1 次查看（过去 30 天）
Alex
Alex 2011-11-19
Hey everyone, I've tried creating a bandpass filter and finding it's transfer function but have had no luck. As you can see below, I tried making the filter and the line for 'd' runs just fine, but when I try to create the IIR filter 'HdIIR' from object d, it doesn't work. It has worked for other objects like 'd' I have created so I do not understand why this one doesn't work (albeit that my other objects used different parameters, but this should still be valid).
Once I get the filter to work, I would like to find the numerator and the denominator of the transfer function (S-domain preferable) of the filter but cannot find a method that will do so given either 'd' or 'HdIIR'.
Please help if you can!
-------------------------------------------------------------
clear all
close all
f_p1=40000; %Pass Frequency in Hz
f_p2=50000; %Stop Frequency in Hz
r_s=25; %Stopband Attenuation in dB
r_p=3; %Passband Ripple
n=10;
d = fdesign.bandpass('N,F3dB1,F3dB2,Ast1,Ap,Ast2',n,f_p1,f_p2,r_s,r_p,r_s,150e3);
HdIIR = design(d,'butter');
-----------------------------------------------------------------
??? Error using ==> fdesign.abstracttype.superdesign at 60
butter is an invalid design method.
Error in ==> fdesign.abstracttype.design at 11
varargout{1} = superdesign(this, varargin{:});
Error in ==> Demo at 13
HdIIR = design(d,'butter');
  1 个评论
Alex
Alex 2011-11-21
Does anyone have any insight as to why this filter will not compile?

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

Kaustubha Govind
Kaustubha Govind 2011-11-21
It appears that 'butter' is not a valid design method for the band-pass filter object.
Using designmethods on your filter design object, I got:
>> designmethods(d)
Design Methods for class fdesign.bandpass (N,F3dB1,F3dB2,Ast1,Ap,Ast2):
ellip
So 'ellip' seems to be the only available design method for your band-pass filter.
  6 个评论
显示 4更早的评论
Wayne King
Wayne King 2011-11-22
B are the numerator coefficients and A are the denominator coefficients. They are just the coefficients for the difference equation.
H is just a MATLAB convention for a handle to an object Hd is a filter object. You don't have to use H. You could use hd, or no h at all, I just use h when I want to remind myself that the output is a handle object.
Alex
Alex 2011-11-22
Hey Wayne, so what domain do I assume these coefficients correlate to? The S or Z domain?
Also, is there any convention on what the denominator coefficients correspond to? Example, is the numerator (A): s^5 + 2.9945s^4 + 8.0928s^3 + ... ?

请先登录，再进行评论。

Naz
Naz 2011-11-21
I always use
fdatool
it's pretty straight forward and you can see the result right away.

类别

在 Help Center 和 File Exchange 中查找有关 Filter Design 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by