Thread Subject: bilinear transformation

Subject: bilinear transformation

From: Ana

Date: 4 Nov, 2009 22:33:03

Message: 1 of 4

hi ,
If i understood correctly bilinear transformation is used in order to transform from s domain to the z domain. Now why then when using;
example;

 fp = 2e03;
 fs = 20e03;
 s = tf('s');
 num = 1;
 den = (s^2 + 1.414*s + 1); % butterworth 2nd order
 [numd,dend] = bilinear(num,den,fs,fp)
it is telling me that function rdivide cannot be classified by tf, i guess numd and dend will give basically the transfer function in the zdomain.

Any ideas highly appreciated.Thanks beforhand

Ana

Subject: bilinear transformation

From: James Pau

Date: 4 Nov, 2009 22:42:02

Message: 2 of 4

"Ana " <filter_world@yahoo.com> wrote in message <hcsvev$gbb$1@fred.mathworks.com>...
> hi ,
> If i understood correctly bilinear transformation is used in order to transform from s domain to the z domain. Now why then when using;
> example;
>
> fp = 2e03;
> fs = 20e03;
> s = tf('s');
> num = 1;
> den = (s^2 + 1.414*s + 1); % butterworth 2nd order
> [numd,dend] = bilinear(num,den,fs,fp)
> it is telling me that function rdivide cannot be classified by tf, i guess numd and dend will give basically the transfer function in the zdomain.
>
> Any ideas highly appreciated.Thanks beforhand
>
> Ana

Hi Ana,

I think it's just the syntax you're using is incorrect.
Try:

fp = 2e03;
 fs = 20e03;
 s = tf('s');
 num = 1;
 den = [2 1.414 1]; % butterworth 2nd order
 [numd,dend] = bilinear(num,den,fs,fp)

Hope that helps,
James

Subject: bilinear transformation

From: Ana

Date: 6 Nov, 2009 19:16:03

Message: 3 of 4

Thanks alot james. Ill give it a good try.

Regards,
Ana

Subject: bilinear transformation

From: Ana

Date: 7 Nov, 2009 18:15:05

Message: 4 of 4

hi there,

 %fp = 2*pi*2e03; For prewarping
 fs = 20e03;
 s = tf('s');
 num = 1;
 den = [1 1.414 1 ];
 
 [numd,dend] = bilinear(num,den,fs)
 
 
z= [1:0.1:20]';
m = filter(numd,dend,z)

freqz(m)

It seems that it is working now, if i understood well the command freqz is used to give the response of discrete time systems, like digital systems does this therefore mean that when obtianing the plot from the code above, I am getting the frequency response of a digital filter in this case the butterworth? Thanks alot for any help.

Regards,
Ana

Tags for this Thread

Add a New Tag:

Separated by commas
Ex.: root locus, bode

What are tags?

A tag is like a keyword or category label associated with each thread. Tags make it easier for you to find threads of interest.

Anyone can tag a thread. Tags are public and visible to everyone.

rssFeed for this Thread
 

MATLAB Central Terms of Use

NOTICE: Any content you submit to MATLAB Central, including personal information, is not subject to the protections which may be afforded information collected under other sections of The MathWorks, Inc. Web site. You are entirely responsible for all content that you upload, post, e-mail, transmit or otherwise make available via MATLAB Central. The MathWorks does not control the content posted by visitors to MATLAB Central and, does not guarantee the accuracy, integrity, or quality of such content. Under no circumstances will The MathWorks be liable in any way for any content not authored by The MathWorks, or any loss or damage of any kind incurred as a result of the use of any content posted, e-mailed, transmitted or otherwise made available via MATLAB Central. Read the complete Terms prior to use.

Contact us at files@mathworks.com