24小时热门版块排行榜

返回列表

狂奔大蜗牛

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 10
帖子: 2
在线: 42分钟
虫号: 1909714
注册: 2012-07-25
专业: 传动机械学

[求助] 使用matlab进行Newton-Raphson迭代是出现虚数的问题

首先说明，我是利用牛顿迭代法求解非线性方程组。所涉及的关键代码如下：
sta_fir=Fa+F0-sum(sta_Kn.*deltan.^(1.5).*sina);
sta_sec=Fr-sum(sta_Kn.*deltan.^(1.5).*(cos(phai)).*cosa);
sta_thir=M-dm/2*sum(sta_Kn.*deltan.^(1.5).*cos(phai).*sina);
static=[sta_fir;sta_sec;sta_thir];
。。。
findsym(static);
findsym(sina);
[gen,n]=Newtonraphson(static,x0);
但是在实际运行的时候因为方程组中1.5次方的缘故会导致迭代出现负数，产生一系列后果。
我在调试中尝试过将方程在牛顿迭代函数中重新定义并添加判断来消除负数，但是这样又会导致原方程组Jacobi矩阵改变，进而出现不收敛的情况。
因此特来此求助，希望有高手能告知妥善的解决方法，不胜感激。

回复此楼

» 猜你喜欢

1楼 2012-07-25 12:44:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主狂奔大蜗牛的主题更新

返回列表