Bilinear transformation Please Help

Discussion in 'The Projects Forum' started by oshaye3, Mar 28, 2011.

  1. oshaye3

    Thread Starter New Member

    Jan 28, 2011
    2
    0
    I am reading a text book and I wanted to use Matlab to get the same answer as they deerived it in the text. The answer that the Matlab displayed is not the same as the text, please this is the question, my code with result and the expected result from the text. The problem is getting the same Transfer function as it is. Note: Also den2=poles of the TF

    Ouestion: Using Bilinear Transformation obtain the transfer function H(z) of the TF below assuming cutt off frequency of 150Hz and sampling frequency of 1.28khz
    Transfer function:
    1
    -----------------
    s^2 + 1.414 s + 1

    Expected Results:

    0.0878 z^2 + 0.1756 z + 0.0878
    ----------------
    z^2 - 1.0048 z + 0.3591

    My Code:

    fsamp = 1280
    %%%%%%Normalize and Prewarp%%%%%%
    % Determine the normalized angular
    wd=linspace(0,pi,100);
    wp1 = 2* pi * ffc / fsamp
    % Prewarp the digital edge frequencies
    wpa1 = tan (wp1 / 2.0)
    [numd, dend] = bilinear (1,den2,fsamp,wpa1);
    H3=tf(numd,dend,0.5)
    H3=tf(numd,dend,0.5,'variable','z^-1')

    my Results:
    wpa1 =

    0.3857


    Transfer function:
    0.288 z^2 + 0.576 z + 0.288
    ---------------------------
    z^2 - 2.828 z + 2.98

    Sampling time: 0.5

    Transfer function:
    0.288 + 0.576 z^-1 + 0.288 z^-2
    -------------------------------
    1 - 2.828 z^-1 + 2.98 z^-2

    Any help! Many thanks! ​
     
Loading...