matlab用fsolve得不到解啊，低手问题

Question

方程组m文件如下：
function q=myfunn(p)
a=p(1);
b=p(2);
c=p(3);
d=p(4);
e=p(5);
syms x
q(1)=int(x^1*exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)/int(exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)+0.0472;
q(2)=int(x^2*exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)/int(exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)-0.9442;
q(3)=int(x^3*exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)/int(exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)+0.1419;
q(4)=int(x^4*exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)/int(exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)-2.5563;
q(5)=int(x^5*exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)/int(exp(a*x+b*x^2+c*x^3+d*x^4+e*x^5),x,-3,3)+0.5687;
两个定积分之商等于一个常数，求解
fsolve('myfunn',[10,10,10,10,10]',optimset('Display','off'))

结果显示10个Warning: Explicit integral could not be found.
外加??? Undefined function or method 'full' for input arguments of type 'sym'.

Error in ==> trustnleqn at 28
Fvec = full(Fvec);

Error in ==> fsolve at 378
[x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=...
请问这是我写的m文件的问题还是初值的问题还是这道题不能用fsolve？
高手如果觉得讲不清能不能留个q？小辈谢谢了

Answer 1

楼主是余海洋狼么？我觉得楼主不是低手呢，在百度混的都半斤八两哈哈。
我觉得你的方法挺牛的，只不过你的函数的返回值给q的是sym类型的，
在你的函数最后加个q=double(q);就能运行了
但是积分那步int函数会说算不出精确值，然后一直蹦Explicit integral could not be found.
感觉没办法，除非修改matlab源函数把它屏蔽了。。。
fsolve得运行停长时间做循环运算的，
我觉得它的option不如改成optimset('Display','on')
这样你每次循环能看到优化的cost function f(x)的收敛程度
我在主函数用的是fsolve('myfunn',[1,1,1,1,1]',optimset('Display','on'))
它会显示这种的:
Norm of First-order Trust-region
Iteration Func-count f(x) step optimality radius
0 6 65895.2 146 1
.....................
.....................
56 247 2.2224e-18 4.56686e-06 4.89e-07 0.202

Equation solved.
fsolve completed because the vector of function values is near zero
as measured by the default value of the function tolerance, and
the problem appears regular as measured by the gradient.
<stopping criteria details>
ans =
-0.0123
-0.5313
-0.0275
0.0025
0.0032

期间会有一堆因为符号积分运算不能给出具体表达式的警告Warning: Explicit integral could not be found.
请楼主无视它。。。