how to see the frequency response of difference equation

22 views (last 30 days)
Hi I am having the difference equation
y[n] = -0.9y[n-1] +x [n] -x [n-1]
I want to see its frequency response. What should i do
should i first convert it into Z-Form and then proceede Any ideas how to go abt it

Accepted Answer

Wayne King
Wayne King on 15 Dec 2011
You can do the following:
W = -pi:(2*pi)/8192:pi;
[H,W] = freqz([1 -1],[1 0.9],W);
plot(W,abs(H));
  4 Comments
Boom
Boom on 7 Aug 2016
How did you get [1 0.9]? Why 1? and why 0.9? I see the equation has a negative 0.9, but not positive 0.9. Thanks
Mike Shane
Mike Shane on 27 Jul 2023
I hope by now you've figured it out, but for those who have the same question here is a simple explanation
y[n] = -0.9y[n-1] +x [n] -x [n-1] is the origianl expression
you have to put y terms on one side and x terms on one side, as below
y[n] + 0.9y[n-1] = x[n] -x[n-1]
now you can take the coefficients of y and x correctly, otherwise you are making a mistake while taking the coefficients
after that you'll get the below,
y has [1 0.9]
x has [1 -1]

Sign in to comment.

More Answers (4)

Wayne King
Wayne King on 15 Dec 2011
Hi, here you go. If you want to visualize the magnitude:
fvtool([1 -1],[1 0.9]);
If you want the complex-valued response:
[H,W] = freqz([1 -1],[1 0.9]);
plot(W,abs(H));
  2 Comments
moonman
moonman on 15 Dec 2011
Hi king thanks
i have plotted it and changed the scale as well in fvtool but not getting the desired option
I have uploaded question here
can u kindly review it
http://www.4shared.com/photo/Y3A4z6qx/question.html
Thanks a lot

Sign in to comment.


Wayne King
Wayne King on 15 Dec 2011
Hi, that's because your graphs are showing the frequency response over [-pi,pi). You can use the 'whole' option in freqz()
[H,W] = freqz([1 -1],[1 0.9],'whole');
plot(W,abs(H));

moonman
moonman on 15 Dec 2011
When i plot by this
[H,W] = freqz([1 -1],[1 0.9],'whole');
W = W-pi;
plot(W,abs(H));
I sm getting the peak magnitude of 20 where as in all options peak is 20. why its so
Anyhow as per ur direction Option 'A' is correct am i right?
  4 Comments
Wayne King
Wayne King on 15 Dec 2011
the book is wrong, the frequency response in magnitude should be 20 at pi radians/sample and -pi radians/sample
moonman
moonman on 15 Dec 2011
Thanks king for explaining me
So the answer is 'A'. Am i right
but one of my frd has solved it and he is saying answer is D
and in email he wrote me
answer would be D as the value of H(0) = 0
Is he wrong.Plz tell me

Sign in to comment.


moonman
moonman on 15 Dec 2011
Thanks a lot King People on forum say that i only accept the answer of Wayne King U always explain in such a nice manner and move with me in question till the time i m not satisfied Thanks a lot and may God give u more knowledge so that u can enlighten more people

Community Treasure Hunt

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

Start Hunting!