weilinhy 发表于 2010-3-21 06:08

请求积分问题

有一个积分关于x的积分 主要是积分函数含有积分下限
H=0.68; s=1; t=30;
syms x;
f=int(power(x,H-0.5)*power(x-s,H-1.5),x,s,t);

运行结果
Warning: Explicit integral could not be found.
> In sym.int at 58

但是个积分在Mathematica里面可以得到具体数值
很奇怪

如果用
double(int(power(x,H-1.5)*power(x-s,H-0.5),x,s,t));
也不行 有waring

请位高手帮忙看看了
谢谢

rocwoods 发表于 2010-3-22 09:10


H = 0.68;
t = 30;
IntResult = @(s) quadl(@(x) power(x,H-0.5).*power(x-s,H-1.5),s,t);
%当s=1时
s = 1;
IntResult(s)
手头没有MATLAB，楼主运行看看

weilinhy 发表于 2010-3-22 22:31

谢谢 recwoods
运行有结果 但是
Warning: Minimum step size reached; singularity possible.
> In quadl at 95
In @(s) quadl(@(x) power(x,H-0.5).*power(x-s,H-1.5),s,t)
请问可以改进么
多谢了

ChaChing 发表于 2010-3-22 23:34

ff=@(x) power(x,H-0.5).*power(x-s,H-1.5)
ezplot(ff)
的确在s=1处singularity

weilinhy 发表于 2010-3-23 00:37

谢谢chachinq论坛高手好多 拜谢
页: [1]
查看完整版本: 请求积分问题