fsolve循环使用中的问题

hebut

还是昨天的八元二次方程组循环求解（http://forum.vibunion.com/thread-106169-1-1.html），每次循环方程组里的参数会发生改变

现在的问题是：每次开始运行，第一次循环得到的结果是正确的，但是第二次第三次和以后的循环要不就是提示

1. Norm of First-order Trust-region
   
2. Iteration Func-count f(x) step optimality radius
   
3. 0 9 1.82858e-014 0 1
   
4. Optimization terminated: first-order optimality is less than options.TolFun.

复制代码

或者提示，初值不对，需要重新设置初值。


下面是主程序

1. 
   
2. x0=[0;0;0;0;0;0;0;0];
   
3. options=optimset('LevenbergMarquardt','on','FunValCheck','on','TolFun',1e-12,'TolX',1e-12,'Display','iter');   
   
4. [x,fval] = fsolve(@myfun,x0,options);

复制代码

下面是方程子函数

1. 
   
2. function F = myfun(x)
   
3. load spectral.mat;
   
4. F = [Gapap_1-Gapap_2-2.....;
   
5. Gapbp_1_re-Gapbp_2_re-Ganbp_1.....;
   
6. Gapbp_1_im-Gapbp_......;
   
7. Gbpbp_1-Gbpbp_2-2*Ganbp_1_re*x(5.....;
   
8. Gapap_2-Gapap_3-2*Ganap_2_re*x(1)+2*Ganap_3_....;
   
9. Gapbp_2_re-Gapbp_3_re-Ganbp_2_re*x(....;
   
10. Gapbp_2_im-Gapbp_3_im-Ganbp_2_im*x(..........;
    
11. Gbpbp_.....];

复制代码

函数太长，没写全
spectral.mat;是主程序中保存的变量值