24小时热门版块排行榜    

查看: 1415  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298-一志愿中国农业大学-求调剂 +11 手机用户 2026-03-17 12/600 2026-03-23 23:51 by 热情沙漠
[考研] 求材料,环境专业调剂 +3 18567500178 2026-03-18 3/150 2026-03-23 23:50 by 热情沙漠
[考研] 一志愿华东理工大学081700,初试分数271 +4 kotoko_ik 2026-03-23 5/250 2026-03-23 23:23 by 呆呆师姐
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 招08考数学 +6 laoshidan 2026-03-20 14/700 2026-03-23 14:37 by 15614359529
[考研] 北科281学硕材料求调剂 +8 tcxiaoxx 2026-03-20 9/450 2026-03-23 12:16 by tcxiaoxx
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 307求调剂 +3 余意卿 2026-03-21 3/150 2026-03-23 10:32 by Iveryant
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
信息提示
请填处理意见