How to manually replicate the bode() gain plot from a transfer function

Question

Shane Palmer 2020-6-4

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/542288-how-to-manually-replicate-the-bode-gain-plot-from-a-transfer-function

评论： Shane Palmer 2020-6-9

采纳的回答： Star Strider

在 MATLAB Online 中打开

Hello,

I have been trying to replicate a bode gain plot from a given transfer function (see code below).

I am not seeing where I went wrong, but my manual plots looks like this:

But the actual bode plots should look like this:

I thought that I was setting up my amplitude equation exactly the same as how the bode() function operates, but apparently not. Can you see where I made the mistake?

Here is my code:

syms C R_1 R_2 R_3 L s V_o i omega
%Equations
R_1 = 10;
R_2 = 10;
R_3 = 10;
L = 0.001; 
C = 2*10^-6;
Z_1 = R_1+L*s;
Z_2 = 1/(C*s)+R_2;
Z_3 = R_3;
Z_23 = 1/(1/(Z_2)+1/(Z_3));
Z_tot = Z_1+Z_23;
V_i = Z_1*i+V_o;
V_o = (V_i/(Z_1))/((1/((1/(C*s))+R_2))+1/R_3+1/(Z_1));
Transfer_func(s) = vpa(simplify(V_o/V_i),5);
amp = (sqrt((5.903*10^24)^2+(1.1806*10^20*omega)^2))/(sqrt((1.1806*10^25+(-2.3613*10^16*omega^2))^2+(9.4447*10^20*omega)^2))
amp_dB = 20*log(amp)
fplot(amp_dB,[0.001 1*10^6])
figure
subplot(2,1,1)
fplot(20*log(abs(Transfer_func)), [0.001 1*10^6])
legend('|H( j\omega )|')
title('Gain (dB)')
grid
subplot(2,1,2)
fplot(180/pi()*angle(Transfer_func), [0.001 1*10^6])
grid
title('Phase (degrees)')
xlabel('Frequency (rad/s)')
legend('\phi( j\omega )')

H = tf([118059162071741125000 5902958103587056517120000],[23611832414348225 944473296573929026712 11805916207174113034240000]);
figure
bode(H,{0.001,1*10^6})

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2020-6-4

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/542288-how-to-manually-replicate-the-bode-gain-plot-from-a-transfer-function#answer_446312

在 MATLAB Online 中打开

‘But the actual bode plots should look like this:’

It does.

Harkening back to yesterday:

syms C R_1 R_2 R_3 L s V_o i omega
%Equations
R_1 = 10;
R_2 = 10;
R_3 = 10;
L = 0.001; 
C = 2*10^-6;
Z_1 = R_1+L*s;
Z_2 = 1/(C*s)+R_2;
Z_3 = R_3;
Z_23 = 1/(1/(Z_2)+1/(Z_3));
Z_tot = Z_1+Z_23;
V_i = Z_1*i+V_o;
V_o = (V_i/(Z_1))/((1/((1/(C*s))+R_2))+1/R_3+1/(Z_1));
Transfer_func(s) = vpa(simplify(V_o/V_i), 5);
Transfer_func(omega) = subs(Transfer_func, {s},{1j*omega});
figure
subplot(2,1,1)
fplot(20*log10(abs(Transfer_func)), [0.001 1E6])
set(gca, 'XScale','log')
legend('|H( j\omega )|')
title('Gain (dB)')
grid
subplot(2,1,2)
fplot(180/pi*angle(Transfer_func), [0.001 1E6])
set(gca, 'XScale','log')
grid
title('Phase (degrees)')
xlabel('Frequency (rad/s)')
legend('\phi( j\omega )')

produces:

Note the slight changes between your code and mine.

7 个评论
显示 5更早的评论隐藏 5更早的评论

Star Strider 2020-6-4

You likely don’t need that, because your code works. When you want to do something new, search the documentation and see what is closest, then adapt it to your problem. That’s what I do.

The difference between your code and mine initially had to do with substituting

for s, and I experimented with the subs function for that. That was just experimentation on my part, and it worked. (Reference: How to substitute jw for s in a transfer function in Matlab? for others reading this.)

With this code, I used the code you posted in your earlier Question (with my changes, along with the correct decibel transformation that was correct in part of your code, although not in the fplot call) and setting the frequency scale to 'log'. (Setting the 'XScale' to 'log' I remember from reading Axis Properties.) I doubt if there are any shortcuts, other than searching the documentation. Also, when you have the time, surf Answers for topics that interest you, and read them. I’ve learned (and am still learning) a lot from Answers over the years by doing just that. Knowledge about MATLAB (and likely everything else) comes with experience using it.

Star Strider 2020-6-9

I suspect Shane Palmer is encountering this for the first time and is experimenting in order to understand how it all works. (I did something similar when I first encountered them, so I understand.)

Shane Palmer 2020-6-9

Star Strider is correct, I am experimenting and learning how to work with matlab and transfer functions simultaneously. I originally started the code symbolically, and then went back in afterwards to add the input values. I don't have a solid answer for you on whether the H is 100% correct, but the manual plot of the bode charts match the matlab bode() command for the charts, so I would say yes as far as the process goes.

请先登录，再进行评论。

How to manually replicate the bode() gain plot from a transfer function

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

采纳的回答

7 个评论
显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

How to manually replicate the bode() gain plot from a transfer function

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

采纳的回答

7 个评论 显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

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

7 个评论
显示 5更早的评论隐藏 5更早的评论