24小时热门版块排行榜    

查看: 1424  |  回复: 8

kimileegdut

捐助贵宾 (小有名气)

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

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

» 猜你喜欢

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

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

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的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
kimileegdut: 金币+5 2015-06-11 10:17:43
当然是采用牛顿拉斐逊方法,先给出一组初值,对方程组进行线性化,以一阶线性化的步长为未知数进行求解,然后将计算出的步长加到原来的初值上,再次对方程组线性化,求解步长......,直至收敛。
5楼2015-06-11 06:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kimileegdut

捐助贵宾 (小有名气)

引用回帖:
5楼: Originally posted by peterflyer at 2015-06-11 06:41:17
当然是采用牛顿拉斐逊方法,先给出一组初值,对方程组进行线性化,以一阶线性化的步长为未知数进行求解,然后将计算出的步长加到原来的初值上,再次对方程组线性化,求解步长......,直至收敛。

你说的是牛顿法吗?如果是的话,牛顿法求解非线性方程组的时候要求雅可比矩阵,我的方程太复杂了,求雅可比矩阵不太可行
6楼2015-06-11 10:12:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

引用回帖:
6楼: Originally posted by kimileegdut at 2015-06-11 10:12:24
你说的是牛顿法吗?如果是的话,牛顿法求解非线性方程组的时候要求雅可比矩阵,我的方程太复杂了,求雅可比矩阵不太可行...

如果你的 {A1}和{A2}可以解析表达为{A3},雅可比就可以做,所以Newton法应该可以。

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

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

7楼2015-06-11 11:06:59
已阅   回复此楼   关注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的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

引用回帖:
8楼: Originally posted by kimileegdut at 2015-06-11 16:50:49
谢谢,又见到您了,前辈!矩阵和列向量相加那个是我打错了,矩阵应该改为列向量,即{K}. *是两个矩阵相乘,应该是我没表达清楚。雅可比矩阵的话我再尝试一下。另外,我再想问下,这种情况能否用Monte Carlo 法?...

孤陋寡闻,不懂Monte Carlo

[ 发自手机版 http://muchong.com/3g ]
9楼2015-06-11 18:40:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kimileegdut 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +4 超级伊昂大王 2026-03-24 4/200 2026-03-24 14:31 by JourneyLucky
[考研] 一志愿哈工大,085400,320,求调剂 +3 gdlf9999 2026-03-24 3/150 2026-03-24 10:08 by 搏击518
[考研] 资源与环境 调剂申请(333分) +6 holy J 2026-03-21 6/300 2026-03-24 09:51 by hengsmith
[考研] 298-一志愿中国农业大学-求调剂 +11 手机用户 2026-03-17 12/600 2026-03-23 23:51 by 热情沙漠
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +5 niko- 2026-03-22 5/250 2026-03-23 22:01 by fuyu_
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 289材料与化工(085600)B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
信息提示
请填处理意见