[求助] Warning: Explicit integral could not be found.怎么解决?
function z=my_fun(x)syms w
z=int(abs((-x(1)*w^2+x(3))/(-i*w^3-x(1)*w^2+i*x(2)*w+x(3))-(-i*(w-3.14)^3+i*x(2)*(w-3.14))/(-i*(w-3.14)^3-x(1)*(w-3.14)^2+i*x(2)*(w-3.14)+x(3))),w)
这个函数有绝对值计算,有复数,积分结果计算不出来,请问怎么办?谢谢! 积分不存在，改用数值积分
回复
function中去掉绝对值abs后，我算了一下，可以算出一个结果，但相当复杂：任取x=;
my_fun(x);
结果经化简后为：
z=
sum(1/125000*_R*log(w^2+(-7/5429687500*i*_R^2+329/1737500*_R-157/50-36/139*i)*w+(7/17375000000+1099/542968750000*i)*_R^2+(-51653/173750000+21/868750*i)*_R+1413/3475*i+324/139),_R = RootOf(7*_Z^3+875000*RootOf(_Z^2+1)*_Z^2+6250000000*_Z-234375000000000*RootOf(_Z^2+1)))-w
转化成手写形式:
pretty(z)
/ -----
|\ / 2
| ) |1/125000 _R log(w
|/ \
| -----
\_R = %1
/ 2 329 157 36 \
+ |- 7/5429687500 i _R+ ------- _R - --- - --- i| w
\ 1737500 50 139/
/ 1099 \ 2 / 51653 21 \
+ |7/17375000000 + ------------ i| _R+ |---------- + ------ i| _R
\ 542968750000/ \ 173750000 868750/
\
1413 324 \|
+ ---- i + ---)|| - w
3475 139 /|
|
/
3 2 2
%1 := RootOf(7 _Z+ 875000 RootOf(_Z+ 1) _Z+ 6250000000 _Z
2
- 234375000000000 RootOf(_Z+ 1))
这个结果几乎还是等于不可积，所以还是按照Happy教授的建议，改用数值积分较好． 不过有点奇怪，当我取x=; 时，竟然可以算出解析结果！
z =
-21/145*i*log(9+25*w^2)+42/145*atan(5/3*w)+11/232*i*log(-3*w^2+4*w^4+1)+11/116*atan(4*w^3-w)+11/116*atan(2*w)+7/58*7^(1/2)*atanh(1/7*(4*w+i)*7^(1/2))-w-21/145*i*log(2500*w^2-15700*w+25549)+42/145*atan(5/3*w-157/30)+11/232*i*log(365047500*w^2-744541100*w+562918826+6250000*w^4-78500000*w^3)+11/116*atan(-942/25*w^2+4*w^3+73322/625*w-1885884/15625)+11/116*atan(2*w-157/25)+7/58*7^(1/2)*atanh(1/175*(100*w-314+25*i)*7^(1/2))
也许参数x选得合适，还是可以解析计算的．．．
页:
[1]