24小时热门版块排行榜    

查看: 1323  |  回复: 6

kelly_ling

新虫 (初入文坛)

[求助] matlab解方程组 已有3人参与

前边发了一个,多谢朋友回复帮忙。不过发现方程整理错误,于是还请大家帮忙。小女子在此谢过了。求解方程组:
联立两个方程式:

matlab解方程组
1.jpg
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2014-06-17 12:35:37
还不知道如何发文本格式代码?
2楼2014-06-11 11:51:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2014-06-17 12:35:47
程序:
Y=-0.0113*X^3+0.6513*X^2-12.378*X+89.932;
Z=X^2+(Y-25)^2;
@SQRT(Z-29.13^2)/29.13-@ACOS(@SQRT(29.13^2/Z))=@ATAN2(Y-25,X);
@FREE(X);@FREE(Y);

结果:
  Feasible solution found.
  Infeasibilities:                        0.543468559044E-08
  Extended solver steps:                              28
  Total solver iterations:                          7131
  Elapsed runtime seconds:                         12.50
  Model Class:                                       NLP
  Total variables:                      3
  Nonlinear variables:                  3
  Integer variables:                    0
  Total constraints:                    3
  Nonlinear constraints:                3
  Total nonzeros:                       8
  Nonlinear nonzeros:                   6

                                               Variable                Value
                                                      Y        59.9840986965
                                                      X        2.81642703906
                                                      Z        1231.81942287

或者:
  Feasible solution found.
  Infeasibilities:                        0.247140451393E-05
  Extended solver steps:                               0
  Total solver iterations:                           122
  Elapsed runtime seconds:                          0.28
  Model Class:                                       NLP
  Total variables:                      3
  Nonlinear variables:                  3
  Integer variables:                    0
  Total constraints:                    3
  Nonlinear constraints:                3
  Total nonzeros:                       8
  Nonlinear nonzeros:                   6

                                               Variable                Value
                                                      Y       -92.4240557017
                                                      X        40.4130515088
                                                      Z        15421.6235922

                                                    Row         Slack or Surplus
                                                      1      -0.181707505931E-06
                                                      2      -0.247140451393E-05
                                                      3      -0.250989451445E-09

1stOpt对于这种非线性问题应该有更好的答案,你@dingd版主吧!我这是用LINGO 14做的。
Be water, my friend.
3楼2014-06-11 12:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王言哲

至尊木虫 (正式写手)


fegg7502: 金币+1, 鼓励交流 2014-06-17 12:35:54
不知道也没有虫友考虑不用matlab自带算法,而是自己编程计算的
哈哈
4楼2014-06-11 21:35:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kelly_ling

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by 王言哲 at 2014-06-11 21:35:55
不知道也没有虫友考虑不用matlab自带算法,而是自己编程计算的

谢谢,不过我看不太懂你的解法,不过还是谢谢了
5楼2014-06-17 11:42:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shikang999

新虫 (小有名气)

【答案】应助回帖

1、对于这种非线性求解问题,可考虑采用智能算法求解,下面的结果采用自己写的遗传算法得到结果
x=25.3397050727989
y=10.6314145851593
残差平方和=0.26623403859127=f(1)^2+f(2)^2,其中f(1)、f(2)分别表示你给的两个方程移项后等式一边只有0的情形,如 f(1)=-0.0113*x^3+0.6513*x^2-12.378*x+89.932-y

2、因为代码直接用程序导出的,可读性不是很强,所以就不上传代码了。遗传算法很简单,自己查一下书,编写应该没有什么问题!
我还是那个我,过去如此,未来亦如此!
6楼2014-06-17 12:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

1stOpt求解:
CODE:
ConstStr Z=X^2+(Y-25)^2;
Function
Y=-0.0113*X^3+0.6513*X^2-12.378*X+89.932;
SQRT(Z-29.13^2)/29.13-ArcCOS(SQRT(29.13^2/Z))=ArcTAN(Y-25/X);

y: 111.307236362773
x: -1.59015428059843
7楼2014-06-17 15:35:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kelly_ling 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:12 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见