How can i plot zero and its frequency responsE

3 ビュー (過去 30 日間)
PAUL
PAUL 2014 年 3 月 21 日
編集済み: Azzi Abdelmalek 2014 年 3 月 23 日
How can i plot zero and its frequency response
b = conv([1 -1.08 0.3132],[1 -1.315 0.6182])
c =b;
a = [0 0 0 0.0686]
[z,p,K] = tf2zp(c,a);
figure(1)
zplane(a,c)
xlabel('Re(z)');
ylabel('Im(z)');
title('POLE DIAGRAM OF THE TRANSFER FUNCTION H(z)')
w = -pi:(pi/500):pi;
H = freqz(a,c,w);
figure(2)
subplot(2,1,1);
plot(w, abs(H),'r','LineWidth',2);
grid on;
xlabel('frequency in pi units');
ylabel('Magntiude');
title('MAGNITUDE FREQUENCY RESPONSE OF THE TRANSFER FUNCTION H(z)')
subplot(2,1,2);
plot(w, angle(H),'r', 'LineWidth',2);
grid on;
xlabel('frequency in pi units');
ylabel('phase');
  1 件のコメント
Azzi Abdelmalek
Azzi Abdelmalek 2014 年 3 月 23 日
Paul, why are you changing your question?

サインインしてコメントする。

採用された回答

Azzi Abdelmalek
Azzi Abdelmalek 2014 年 3 月 21 日
編集済み: Azzi Abdelmalek 2014 年 3 月 21 日
The numerator should be
a = 0.0686
To find poles and zeros
[z,p,K] = tf2zp(a,c); % not tf2zp(c,a)
Add
b = conv([1 -1.08 0.3132],[1 -1.315 0.6182])
c =b;
a = 0.0686
[z,p,K] = tf2zp(a,c);
figure(1)
zplane(a,c)
xlabel('Re(z)');
ylabel('Im(z)');
title('POLE DIAGRAM OF THE TRANSFER FUNCTION H(z)')
w = -pi:(pi/500):pi;
H = freqz(a,c,w);
figure(2)
subplot(2,1,1);
plot(w, abs(H),'r','LineWidth',2);
grid on;
xlabel('frequency in pi units');
ylabel('Magntiude');
title('MAGNITUDE FREQUENCY RESPONSE OF THE TRANSFER FUNCTION H(z)')
subplot(2,1,2);
plot(w, angle(H),'r', 'LineWidth',2);
grid on;
xlabel('frequency in pi units');
ylabel('phase');
  3 件のコメント
Azzi Abdelmalek
Azzi Abdelmalek 2014 年 3 月 21 日
Paul commented
Hi Azzi Abdelmalek Can you please help me with this one.I'm new to MAT LAB.
Azzi Abdelmalek
Azzi Abdelmalek 2014 年 3 月 21 日
The frequency interval is not well chosen. It should be
w = -2*pi:(pi/500):2*pi

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeSystem Identification Toolbox についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by