Solving Integral in matlab

2023 4 17
2 个回答
更新时间：2023 4 17
6 次查看（30 天）

0 个投票

I need to solve this integral in matlab
Where a=1, b=0.71, c=0.53.
I have tried this in matlab
syms s
a=1;
b=0.71;
c=0.53;
f = 1/((s+a^2)*sqrt(((s+a)^2)*((s+b)^2)*((s+c)^2)));
F = int(f, 0, inf)
but I am unable to solve the problem

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

回答（2 个）

Dyuman Joshi
Dyuman Joshi 2023-4-17
编辑：Dyuman Joshi 2023-4-17
Ran in:

0 个投票

Ran in:
There are cases when int is unable to compute the value of definite integral. In such cases, use vpa to obtain the approximate numerical value
syms s
a = 1;
b = 0.71;
c = 0.53;
f = 1/((s+a^2)*sqrt(((s+a)^2)*((s+b)^2)*((s+c)^2)));
F = int(f, 0, inf);
val1 = vpa(F)
val1 = 
0.67917589608295134074024322188586
You can use vpaintegral directly as well -
val2 = vpaintegral(f,0,inf)
val2 = 
0.679176
Torsten
Torsten 2023-4-17
编辑：Torsten 2023-4-17
Ran in:

0 个投票

Ran in:
syms s
a=sym('1');
b=sym('71/100');
c=sym('53/100');
f = 1/((s+a^2)*(s+a)*(s+b)*(s+c));
F = a*b*c/2*int(f)
F = 
limit(F-subs(F,s,0),s,Inf)
ans = 

类别

在 帮助中心 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

标签

编辑：

Torsten
Torsten
2023-4-17

Community Treasure Hunt

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

Start Hunting!

Translated by