Inserting vertical and horizontal line in bodeplot

Question

Carl Overgaard 2019-11-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/494012-inserting-vertical-and-horizontal-line-in-bodeplot

评论： Carl Overgaard 2019-11-30

采纳的回答： Stephan

在 MATLAB Online 中打开

Hello

With below code I have made a bode plot. Then I have calculatted an angular frequency.

I want to insert a vertical line from my found angular frequency(found to be 0.041889), and also a horizontal line so I can determine the phase angel on the y-axis:

Can this be done?

IThe plan is that it should look like this:

My code:

sys1=tf([1805.1],[4637.61 136.2 1]) %transfer function
w=logspace(-5,2); %Define frequency range
[A,phi]=bode(sys1,w); %Amplitude ratio (A) and phase shift (phi)
figure
subplot(2,1,1)
loglog(w,A(:)) %Amplitude ratio vs frequency
grid on
ylabel('AR')
subplot(2,1,2)
semilogx(w,phi(:))
ylabel('Phase angle(Degrees)')
xlabel('Frequency(rad/s)')
legend({'Lower tank'},'Fontsize',14)  %add "title" to plot...

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

请先登录，再进行评论。

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

Answer 1

Stephan 2019-11-30

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/494012-inserting-vertical-and-horizontal-line-in-bodeplot#answer_404002

编辑：Stephan 2019-11-30

xline and yline should work for this purpose. Both functions were introduced with R2018b. If you use an earlier release, just define two points for every line you need and use the plot command combined with hold on and hold off including the line format as you wish it to look like.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Stephan 2019-11-30

编辑：Stephan 2019-11-30

在 MATLAB Online 中打开

you could interpolate:

phi1 = squeeze(phi);
ypoint = interp1(w,phi1,0.041889)
ypoint =
 -140.6024

the whole code would look like:

sys1=tf([1805.1],[4637.61 136.2 1]) %transfer function
w=logspace(-5,2); %Define frequency range
[A,phi]=bode(sys1,w); %Amplitude ratio (A) and phase shift (phi)
figure
subplot(2,1,1)
loglog(w,A(:)) %Amplitude ratio vs frequency
grid on
ylabel('AR')
subplot(2,1,2)
semilogx(w,phi(:))
ylabel('Phase angle(Degrees)')
xlabel('Frequency(rad/s)')
legend({'Lower tank'},'Fontsize',14)  %add "title" to plot...
hold on
phi1 = squeeze(phi);
ypoint = interp1(w,phi1,0.041889)
scatter(0.041889,ypoint,'or')
xline(0.041889,'--r')
yline(ypoint,'--r')