MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   数学 (https://www.labfans.com/bbs/forumdisplay.php?f=15)
-   -   在高数中可积,但matlab不能积,大家帮忙看看 (https://www.labfans.com/bbs/showthread.php?t=10700)

seather 2010-02-01 21:14

在高数中可积,但matlab不能积,大家帮忙看看
 
syms t;
syms s;
f=s^0.9*(t-s)^(-0.1);
int(f,s,0,t)
上面是程序。

s^0.9*(t-s)^(-0.1)在[0,t]上对s积分。
高数的书上有结果,但是matlab积不出来。
提示:Warning: Explicit integral could not be found.

ans =

piecewise([t = 1, (5*gamma(9/10)^2)/(8*gamma(4/5))], [Otherwise,
int(s^(9/10)/(t - s)^(1/10), s = 0..t)])

adolphice 2010-02-26 22:01

回复: 在高数中可积,但matlab不能积,大家帮忙看看
 
试了下,你的可以计算出来啊,要是看做0到T的符号积分,可以算出以下
>> syms t;
>> syms s;
>> f=s^0.9*(t-s)^(-0.1);
>> int(f,s,0,t)

ans =

5/8/(1/t)^(9/10)*t^(9/10)*gamma(9/10)*2^(1/5)/pi^(1/2)*sin(2/5*pi)*gamma(3/5)
不知道你问的是这个意思吗?

adolph 2010-02-26 22:02

回复: 在高数中可积,但matlab不能积,大家帮忙看看
 
试了下,你的可以计算出来啊,要是看做0到T的符号积分,可以算出以下
>> syms t;
>> syms s;
>> f=s^0.9*(t-s)^(-0.1);
>> int(f,s,0,t)

ans =

5/8/(1/t)^(9/10)*t^(9/10)*gamma(9/10)*2^(1/5)/pi^(1/2)*sin(2/5*pi)*gamma(3/5)
不知道你问的是这个意思吗?


所有时间均为北京时间。现在的时间是 21:25

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.