| 查看: 2327 | 回复: 0 | ||
[求助]
Matlab：求高人指点用matlab求解非线性方程组，解出来的值不收敛，提前结束
|
|
fsolve解非线性方程组，只迭代了8次就退出了，得到的结果进行复算不等于0，存在较大的误差。（是不是需要修改算法？） 运行程序得到的结果提示如下： x = 1.0e+003 * Columns 1 through 4 0.008229968721187 0.000000000008083 0.000006858621379 1.206234368839563 Column 5 0.000029820218392 fval = -0.003404951183545 0.109922506812515 0.055200484333056 0.009207767273016 -0.000813641543979 exitflag = -2 output = iterations: 8 funcCount: 54 algorithm: 'trust-region dogleg' firstorderopt: 3.163833174781097e+003 message: [1x775 char] ans = No solution found. fsolve stopped because the relative size of the current step is less than the selected value of the step size tolerance squared, but the vector of function values is not near zero as measured by the selected value of the function tolerance. Stopping criteria details: fsolve stopped because the relative norm of the current step, 2.142692e-011, is less than max(options.TolX^2,eps) = 2.500000e-011. However, the sum of squared function values, r = 1.522709e-002, exceeds sqrt(options.TolFun) = 1.000000e-025. Optimization Metric Options relative norm(step) = 2.14e-011 max(TolX^2,eps) = 3e-011 (selected) r = 1.52e-002 sqrt(TolFun) = 1.0e-025 (selected) 程序代码如下：（其中调用了mafun函数，mafun函数中包括五个方程组，Iph，Io，Rs，Rp，Vt是要求的参数，已经赋初值，可以保证初值选择正确，因为所有解的范围是可以估计到的） options = optimset('Display','iter','TolFun',1e-50,'TolX',0.5e-5); [x,fval,exitflag,output] = fsolve(@mafun,x0,options) output.message format long, x, 如果我将程序中的optimset括号中的“0.5e-5”改为“0.5e-20”，则得到的结果为（好像说是变成奇异方程了）： No solution found. fsolve stopped because the problem appears to be locally singular. Stopping criteria details: fsolve stopped because the trust-region radius, 1.248209e-016, is less than 2*eps. Optimization Metric Options trust-region radius = 1.25e-016 2*eps = 4.4e-016 上面的问题已经说的很详细了，希望高手能帮忙解决。谢谢了。 只有那么多Q币，全给了。 [ Last edited by mataohk on 2012-9-30 at 09:32 ] |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有7人回复
-
这样的filecode谁见过
已经有18人回复
-
关于代码变化问题，想知道的进来
已经有17人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
求助 MATLAB解方程组-fslove
已经有7人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab求解方程组“”
已经有4人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
