[求助] Warning: Explicit integral could not be found.怎么解决?

lovv

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)

这个函数有绝对值计算,有复数,积分结果计算不出来,请问怎么办?谢谢!

happy

积分不存在，改用数值积分

xjzuo

function中去掉绝对值abs后，我算了一下，可以算出一个结果，但相当复杂：
任取x=[1 2 3];
         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教授的建议，改用数值积分较好．

xjzuo

不过有点奇怪，当我取x=[0.1 0.2 0.3]; 时，竟然可以算出解析结果！
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选得合适，还是可以解析计算的．．．