循环使用fsolve的问题
本帖最后由 hebut 于 2011-10-8 11:06 编辑对一个八元二次方程组进行循环求解；方程组里有一些数是每次循环求的的
每次进行一次循环，主程序求出这些已知量，然后保存成mat文件，在子程序中load使用
现在的问题是：不论设置初值是多少，fsolve都几乎不做迭代求解会提示已经找到根
提示如下：
First-Order Norm of
IterationFunc-count Residual optimality Lambda step
0 9 1.10716e-011 1.09e-011 0.01
1 18 1.10716e-011 1.09e-011 0.001 1.54531e-009
Optimization terminated: the relative change in the sum-of-squares of the functions is less than options.TolFun
First-Order Norm of
IterationFunc-count Residual optimality Lambda step
0 9 1.91064e-011 0 0.01
Optimization terminated: the first-order optimality measure is less than 1e-4 times options.TolFun.
First-Order Norm of
IterationFunc-count Residual optimality Lambda step
0 9 1.80203e-009 0 0.01
Optimization terminated: the first-order optimality measure is less than 1e-4 times options.TolFun.
First-Order Norm of
IterationFunc-count Residual optimality Lambda step
0 9 3.46033e-007 0 0.01
Optimization terminated: the first-order optimality measure is less than 1e-4 times options.TolFun.
First-Order Norm of
IterationFunc-count Residual optimality Lambda step
0 9 8.02442e-007 0 0.01
Optimization terminated: the first-order optimality measure is less than 1e-4 times options.TolFun.
First-Order Norm of
IterationFunc-count Residual optimality Lambda step
0 9 1.01221e-006 0 0.01
Optimization terminated: the first-order optimality measure is less than 1e-4 times options.TolFun.
First-Order Norm of
IterationFunc-count Residual optimality Lambda step
0 9 4.86868e-006 0 0.01
Optimization terminated: the first-order optimality measure is less than 1e-4 times options.TolFun.
待求函数如下：
其中x（1）到x（8）为未知数，其余量为已知数，从主程序中获得的
F = [Gapap_1-Gapap_2-2*Ganap_1_re*x(1)+2*Ganap_2_re*x(1)+(Ganan_1-Ganan_2)*x(1)*x(1)+2*(Ganap_1_im-Ganap_2_im)*x(2)+(Ganan_1-Ganan_2)*x(2)*x(2)-2*Gbnap_1_re*x(3)+2*Gbnap_2_re*x(3)+2*Ganbn_1_re*x(1)*x(3)-2*Ganbn_2_re*x(1)*x(3)-2*Ganbn_1_im*x(2)*x(3)+2*Ganbn_2_im*x(2)*x(3)+Gbnbn_1*x(3)*x(3)-Gbnbn_2*x(3)*x(3)+2*(Gbnap_1_im-Gbnap_2_im+Ganbn_1_im*x(1)-Ganbn_2_im*x(1)+Ganbn_1_re*x(2)-Ganbn_2_re*x(2))*x(4)+(Gbnbn_1-Gbnbn_2)*x(4)*x(4);
Gapbp_1_re-Gapbp_2_re-Ganbp_1_re*x(1)+Ganbp_2_re*x(1)+Ganbp_1_im*x(2)-Ganbp_2_im*x(2)-Gbnbp_1_re*x(3)+Gbnbp_2_re*x(3)+Gbnbp_1_im*x(4)-Gbnbp_2_im*x(4)-Ganap_1_re*x(5)+Ganap_2_re*x(5)+Ganan_1*x(1)*x(5)-Ganan_2*x(1)*x(5)+Ganbn_1_re*x(3)*x(5)-Ganbn_2_re*x(3)*x(5)+Ganbn_1_im*x(4)*x(5)-Ganbn_2_im*x(4)*x(5)+Ganap_1_im*x(6)-Ganap_2_im*x(6)+Ganan_1*x(2)*x(6)-Ganan_2*x(2)*x(6)-Ganbn_1_im*x(3)*x(6)+Ganbn_2_im*x(3)*x(6)+Ganbn_1_re*x(4)*x(6)-Ganbn_2_re*x(4)*x(6)-Gbnap_1_re*x(7)+Gbnap_2_re*x(7)+Ganbn_1_re*x(1)*x(7)-Ganbn_2_re*x(1)*x(7)-Ganbn_1_im*x(2)*x(7)+Ganbn_2_im*x(2)*x(7)+Gbnbn_1*x(3)*x(7)-Gbnbn_2*x(3)*x(7)+Gbnap_1_im*x(8)-Gbnap_2_im*x(8)+Ganbn_1_im*x(1)*x(8)-Ganbn_2_im*x(1)*x(8)+Ganbn_1_re*x(2)*x(8)-Ganbn_2_re*x(2)*x(8)+Gbnbn_1*x(4)*x(8)-Gbnbn_2*x(4)*x(8);
Gapbp_1_im-Gapbp_2_im-Ganbp_1_im*x(1)+Ganbp_2_im*x(1)-Ganbp_1_re*x(2)+Ganbp_2_re*x(2)-Gbnbp_1_im*x(3)+Gbnbp_2_im*x(3)-Gbnbp_1_re*x(4)+Gbnbp_2_re*x(4)+Ganap_1_im*x(5)-Ganap_2_im*x(5)+Ganan_1*x(2)*x(5)-Ganan_2*x(2)*x(5)-Ganbn_1_im*x(3)*x(5)+Ganbn_2_im*x(3)*x(5)+Ganbn_1_re*x(4)*x(5)-Ganbn_2_re*x(4)*x(5)+Ganap_1_re*x(6)-Ganap_2_re*x(6)-Ganan_1*x(1)*x(6)+Ganan_2*x(1)*x(6)-Ganbn_1_re*x(3)*x(6)+Ganbn_2_re*x(3)*x(6)-Ganbn_1_im*x(4)*x(6)+Ganbn_2_im*x(4)*x(6)+Gbnap_1_im*x(7)-Gbnap_2_im*x(7)+Ganbn_1_im*x(1)*x(7)-Ganbn_2_im*x(1)*x(7)+Ganbn_1_re*x(2)*x(7)-Ganbn_2_re*x(2)*x(7)+Gbnbn_1*x(4)*x(7)-Gbnbn_2*x(4)*x(7)+(Gbnap_1_re-Gbnap_2_re-Ganbn_1_re*x(1)+Ganbn_2_re*x(1)+Ganbn_1_im*x(2)-Ganbn_2_im*x(2)-Gbnbn_1*x(3)+Gbnbn_2*x(3))*x(8);
Gbpbp_1-Gbpbp_2-2*Ganbp_1_re*x(5)+2*Ganbp_2_re*x(5)+Ganan_1*x(5)*x(5)-Ganan_2*x(5)*x(5)+2*Ganbp_1_im*x(6)-2*Ganbp_2_im*x(6)+Ganan_1*x(6)*x(6)-Ganan_2*x(6)*x(6)-2*Gbnbp_1_re*x(7)+2*Gbnbp_2_re*x(7)+2*Ganbn_1_re*x(5)*x(7)-2*Ganbn_2_re*x(5)*x(7)-2*Ganbn_1_im*x(6)*x(7)+2*Ganbn_2_im*x(6)*x(7)+Gbnbn_1*x(7)*x(7)-Gbnbn_2*x(7)*x(7)+2*(Gbnbp_1_im-Gbnbp_2_im+Ganbn_1_im*x(5)-Ganbn_2_im*x(5)+Ganbn_1_re*x(6)-Ganbn_2_re*x(6))*x(8)+(Gbnbn_1-Gbnbn_2)*x(8)*x(8);
Gapap_2-Gapap_3-2*Ganap_2_re*x(1)+2*Ganap_3_re*x(1)+(Ganan_2-Ganan_3)*x(1)*x(1)+2*(Ganap_2_im-Ganap_3_im)*x(2)+(Ganan_2-Ganan_3)*x(2)*x(2)-2*Gbnap_2_re*x(3)+2*Gbnap_3_re*x(3)+2*Ganbn_2_re*x(1)*x(3)-2*Ganbn_3_re*x(1)*x(3)-2*Ganbn_2_im*x(2)*x(3)+2*Ganbn_3_im*x(2)*x(3)+Gbnbn_2*x(3)*x(3)-Gbnbn_3*x(3)*x(3)+2*(Gbnap_2_im-Gbnap_3_im+Ganbn_2_im*x(1)-Ganbn_3_im*x(1)+Ganbn_2_re*x(2)-Ganbn_3_re*x(2))*x(4)+(Gbnbn_2-Gbnbn_3)*x(4)*x(4);
Gapbp_2_re-Gapbp_3_re-Ganbp_2_re*x(1)+Ganbp_3_re*x(1)+Ganbp_2_im*x(2)-Ganbp_3_im*x(2)-Gbnbp_2_re*x(3)+Gbnbp_3_re*x(3)+Gbnbp_2_im*x(4)-Gbnbp_3_im*x(4)-Ganap_2_re*x(5)+Ganap_3_re*x(5)+Ganan_2*x(1)*x(5)-Ganan_3*x(1)*x(5)+Ganbn_2_re*x(3)*x(5)-Ganbn_3_re*x(3)*x(5)+Ganbn_2_im*x(4)*x(5)-Ganbn_3_im*x(4)*x(5)+Ganap_2_im*x(6)-Ganap_3_im*x(6)+Ganan_2*x(2)*x(6)-Ganan_3*x(2)*x(6)-Ganbn_2_im*x(3)*x(6)+Ganbn_3_im*x(3)*x(6)+Ganbn_2_re*x(4)*x(6)-Ganbn_3_re*x(4)*x(6)-Gbnap_2_re*x(7)+Gbnap_3_re*x(7)+Ganbn_2_re*x(1)*x(7)-Ganbn_3_re*x(1)*x(7)-Ganbn_2_im*x(2)*x(7)+Ganbn_3_im*x(2)*x(7)+Gbnbn_2*x(3)*x(7)-Gbnbn_3*x(3)*x(7)+Gbnap_2_im*x(8)-Gbnap_3_im*x(8)+Ganbn_2_im*x(1)*x(8)-Ganbn_3_im*x(1)*x(8)+Ganbn_2_re*x(2)*x(8)-Ganbn_3_re*x(2)*x(8)+Gbnbn_2*x(4)*x(8)-Gbnbn_3*x(4)*x(8);
Gapbp_2_im-Gapbp_3_im-Ganbp_2_im*x(1)+Ganbp_3_im*x(1)-Ganbp_2_re*x(2)+Ganbp_3_re*x(2)-Gbnbp_2_im*x(3)+Gbnbp_3_im*x(3)-Gbnbp_2_re*x(4)+Gbnbp_3_re*x(4)+Ganap_2_im*x(5)-Ganap_3_im*x(5)+Ganan_2*x(2)*x(5)-Ganan_3*x(2)*x(5)-Ganbn_2_im*x(3)*x(5)+Ganbn_3_im*x(3)*x(5)+Ganbn_2_re*x(4)*x(5)-Ganbn_3_re*x(4)*x(5)+Ganap_2_re*x(6)-Ganap_3_re*x(6)-Ganan_2*x(1)*x(6)+Ganan_3*x(1)*x(6)-Ganbn_2_re*x(3)*x(6)+Ganbn_3_re*x(3)*x(6)-Ganbn_2_im*x(4)*x(6)+Ganbn_3_im*x(4)*x(6)+Gbnap_2_im*x(7)-Gbnap_3_im*x(7)+Ganbn_2_im*x(1)*x(7)-Ganbn_3_im*x(1)*x(7)+Ganbn_2_re*x(2)*x(7)-Ganbn_3_re*x(2)*x(7)+Gbnbn_2*x(4)*x(7)-Gbnbn_3*x(4)*x(7)+(Gbnap_2_re-Gbnap_3_re-Ganbn_2_re*x(1)+Ganbn_3_re*x(1)+Ganbn_2_im*x(2)-Ganbn_3_im*x(2)-Gbnbn_2*x(3)+Gbnbn_3*x(3))*x(8);
Gbpbp_2-Gbpbp_3-2*Ganbp_2_re*x(5)+2*Ganbp_3_re*x(5)+Ganan_2*x(5)*x(5)-Ganan_3*x(5)*x(5)+2*Ganbp_2_im*x(6)-2*Ganbp_3_im*x(6)+Ganan_2*x(6)*x(6)-Ganan_3*x(6)*x(6)-2*Gbnbp_2_re*x(7)+2*Gbnbp_3_re*x(7)+2*Ganbn_2_re*x(5)*x(7)-2*Ganbn_3_re*x(5)*x(7)-2*Ganbn_2_im*x(6)*x(7)+2*Ganbn_3_im*x(6)*x(7)+Gbnbn_2*x(7)*x(7)-Gbnbn_3*x(7)*x(7)+2*(Gbnbp_2_im-Gbnbp_3_im+Ganbn_2_im*x(5)-Ganbn_3_im*x(5)+Ganbn_2_re*x(6)-Ganbn_3_re*x(6))*x(8)+(Gbnbn_2-Gbnbn_3)*x(8)*x(8)]; 你好，我现在也碰到了相同问题，
http://forum.vibunion.com/forum.php?mod=viewthread&tid=129365
请问你以前是如何解决的？
页:
[1]