24小时热门版块排行榜    

查看: 1607  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见