Plotting IIR filter error

Question

Nathan Jaqua 2019-11-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/488870-plotting-iir-filter-error

回答： Star Strider 2019-11-2

采纳的回答： Star Strider

I have an IIR filter I am trying to plot, but I am not sure what I am doing wrong here. Can someone please help?

My code:

%% This is a first part of an exercise.

% Here you are implementing the transfer function H(z).

% Pass different values of theta and r to have differnt outcome.

clc; clear all;

%% Part -1 simulation

% initialize r

r = 0.99 ;

%initialize theta

theta = pi/3 ;

% initialize omega vector. Let, ω to start from 0 to pi with % 200 intermediate values.

w = 200 ;

% Now, define z = e^(j*w).

z = exp(1i*w);

% implement H

H = 1./(1 - (2*r*cos(theta)*(z.^-1))+((r^2)*(z.^-2)));

% Now, plot the filter magnitude response in dB scale w.r.t w.

figure;plot(w, );grid on;

xlabel('normalized frequency'); ylabel('amplitude in dB')

title('Frequency and phase response at r = and theta = ');

2 个评论
显示无隐藏无

Star Strider 2019-11-2

What result are you supposed to get, and what result are you getting?

Nathan Jaqua 2019-11-2

I have not even been able to get anything to show up in the plot. I need it to plot properly.

请先登录，再进行评论。

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

Answer 1

Star Strider 2019-11-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/488870-plotting-iir-filter-error#answer_399505

在 MATLAB Online 中打开

You are plotting one point, so you need to plot a marker instead of a line.

I am not certain what you are doing, so consider this example:

% Now, define z = e^(j*w).
z = @(w) exp(1i*w);
 
% implement H
H = @(w) 1./(1 - (2*r*cos(theta)*(z(w).^-1))+((r^2)*(z(w).^-2)));
 
% Now, plot the filter magnitude response in dB scale w.r.t w.  
w = linspace(-100*pi, 100*pi, 250);
figure;plot(w, abs(H(w)));grid on;
xlabel('normalized frequency'); ylabel('amplitude in dB')
title('Frequency and phase response at r = and theta = ');