i want to find out area under the curve using quad command ,the curve is generated after running the programme

Question

sibananda mohanty 2012-11-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/53026-i-want-to-find-out-area-under-the-curve-using-quad-command-the-curve-is-generated-after-running-the

clc
clear all
s=tf('s')
alpha=0.05
mu=0.05;
ztap=0;
ztaa=sqrt(3*mu/(8*(1+mu)));
omegaa=(1/(1+mu));
a1=1-alpha;
a2=2*ztaa*omegaa+2*ztap+2*ztaa*omegaa*mu;
a3=omegaa^2+4*ztap*ztaa*omegaa+mu*omegaa^2+1;
a4=2*ztap*omegaa^2+2*ztaa*omegaa;
a5=omegaa^2;
A=s^2+2*ztaa*omegaa*s+omegaa^2;
D=a1*s^4+a2*s^3+a3*s^2+a4*s+a5;
G=A/D;
bodemag(G)

After running this programme a curve is generated (bode plot) i want find out area under this curve by using quad syntax i tried to solve this many times i m very new to matlab please help me to find out the answer

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

Grzegorz Knor 2012-11-7

Please edit your code.

请先登录，再进行评论。

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

Answer 1

Grzegorz Knor 2012-11-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/53026-i-want-to-find-out-area-under-the-curve-using-quad-command-the-curve-is-generated-after-running-the#answer_64616

在 MATLAB Online 中打开

If you draw yourself Bode plot (see code), you receive vectors 'x' and 'y', which can be numerically integrate.

[HMag, HPhase, HOmega] = bode(G);
HMag   = squeeze(HMag);
HPhase = squeeze(HPhase);
figure
semilogx(HOmega, 20*log10(HMag),'r')
xlabel('Frequency (rad/s)')
ylabel('Magnitude (dB)')
title('Bode Diagram')