my plot shows nothing! any suggestion would be appreciated. I know i can use bode() and tf, but i wanna do it in this way. thanks

Question

ocsse 2017-11-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/369843-my-plot-shows-nothing-any-suggestion-would-be-appreciated-i-know-i-can-use-bode-and-tf-but-i-wa

回答： Star Strider 2017-11-28

syms w
% Resistor = 1000 Ohm
ZR = 1e3;
% Inductor = 100 mH
ZL = 100e-3;
% Capacitor = 500 microF
ZC = 500e-6;
for w = 200*pi:2000*pi
      %transfer function
      H = 100 * ( 1j*w / ( 1j*w + ( 1 / (ZR*ZC) ) ) );
      %gain
      gain_dB = mag2db(real(H))
      %phase
      phase = angle(H);
end
%frequency
fr_Hz = linspace(100,1000);
% plot the gain
figure(1)
plot(fr_Hz,gain_dB)
title('Bode Plot Gain')
xlabel('Frequency (Hz)')
ylabel('Gain (dB)')
% plot the phasor
figure(2)
plot(fr_Hz, phase)
title('Bode Plot (Phasor)')
xlabel('Frequency (Hz)')
ylabel('Phase (\circ)')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2017-11-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/369843-my-plot-shows-nothing-any-suggestion-would-be-appreciated-i-know-i-can-use-bode-and-tf-but-i-wa#answer_293648

在 MATLAB Online 中打开

Save your variables to vectors by subscripting them:

w = 200*pi:2000*pi;
for k1 = 1:length(w)
        %transfer function
        H(k1) = 100 * ( 1j*w(k1) / ( 1j*w(k1) + ( 1 / (ZR*ZC) ) ) );
        %gain
        gain_dB(k1) = mag2db(real(H(k1)));
        %phase
        phase(k1) = angle(H(k1));
end