24小时热门版块排行榜    

查看: 1366  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kimileegdut

捐助贵宾 (小有名气)

[求助] 求解非线性方程组 已有2人参与

如题,当未知数比较多,比如有49个的时候,用什么方法可以比较有效地求解?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

kimileegdut

捐助贵宾 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by pippi6 at 2015-06-11 11:06:59
如果你的 {A1}和{A2}可以解析表达为{A3},雅可比就可以做,所以Newton法应该可以。

另外,你的方程{ui}=+{Ai} 和{uix}=+{Ai}有误,矩阵不能和列向量相加,而且在方程里 你又写成*之类,很是混淆。...

谢谢,又见到您了,前辈!矩阵和列向量相加那个是我打错了,矩阵应该改为列向量,即{K}. [U3X]*[U3X]是两个矩阵相乘,应该是我没表达清楚。雅可比矩阵的话我再尝试一下。另外,我再想问下,这种情况能否用Monte Carlo 法?
8楼2015-06-11 16:50:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
要看具体问题!
一般的要用最优化方法!

» 本帖已获得的红花(最新10朵)

善恶到头终有报,人间正道是沧桑.
2楼2015-06-10 20:27:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kimileegdut

捐助贵宾 (小有名气)

引用回帖:
2楼: Originally posted by wurongjun at 2015-06-10 20:27:36
要看具体问题!
一般的要用最优化方法!

我要求解的方程组是这样的:
{([u1x]+0.5*[u3x]*[u3x])+0.3*([u2y]+0.5*[u3y]*[u3y])}*{u3xx}+0.7*([u1y]+[u2x]+[u3x][u3y])*{u3xy}
+{([u2y]+0.5*[u3y][u3y])+0.3*([u1x]+0.5*[u3x]*[u3x])}*
{u3yy}+0.004136*({p}-2*ks*{u3})

其中
[] 为49*49的对角矩阵
{}为49*1的列阵
{ui}=[K]+[D]{Ai}  ,i=1,2,3
{uix}=[KX]+[DX]{Ai},i=1,2,3
如此类推,[K],[KX],...,[D],[DX],...为系数矩阵,{p}也为已知列阵。
对角矩阵其对角元素和对应的列阵相同,其余元素均为0,比如,[u1x]的对角元素和{u1x}相同。
然后{A1}和{A2}又可以通过{A3}求得,现在要求解{A3}
3楼2015-06-10 21:01:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kimileegdut

捐助贵宾 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by wurongjun at 2015-06-10 20:27:36
要看具体问题!
一般的要用最优化方法!

打漏了一点,方程组最后是={0}
4楼2015-06-10 21:07:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见