| 查看: 2324 | 回复: 13 | |||
yuanhongjuam木虫 (初入文坛)
|
[求助]
求解非线性方程组的代码 已有3人参与
|
|
求解非线性方程组的代码，什么软件可以解这样的13个方程、13个未知数的方程组？x1到x13是要求的未知数。 9.52*10^6*x13-5.48*10^12*x1-2.14*10^11*x1*x13+2.79*10^(-4)*x2*x9-2.44*10^7*x1*x12+1.96*10^5*x3*x13-1.81*10^5*x1*x12+6.35*10^6*x4*x3=0; 2.14*10^11*x1*x13-2.79*10^(-4)*x2*x9+6.23*10^7*x4*x13-17.6*x2*x10-2.07*10^(-3)*x2*x12+3.32*10^11*x5*x13-1.13*10^(-3)*x2+1.19*10^7*x13=0; 2.44*10^7*x1*x12-1.96*10^5*x3*x13-6.11*10^3*x3*x13+2.94*10^3*x5*x9=0; 1.81*10^5*x1*x12-6.35*10^6*x4*x13-6.23*10^7*x4*x13+17.6*x2*x10=0; 6.11*10^3*x3*x13-2.94*10^3*x5*x9+2.07*10^(-3)*x2*x12-3.32*10^11*x5*x13-3.76*10^10*x5*x13+1.16*10^(-2)*x6*x9=0; 3.76*10^10*x5*x13-1.16*10^(-2)*x6*x9-1.10*10^7*x6*x12+3.09*10^10*x7*x13=0; 1.10*10^7*x6*x12-3.09*10^10*x7*x13-1.72*10^7*x7*x12+4.48*10^5*x8*x13=0; 1.72*10^7*x7*x12-4.48*10^5*x8*x13-1.35*10^6*x8*x12=0; 2.14*10^11*x1*x13-2.79*10^(-4)*x2*x9+6.11*10^3*x3*x13-2.94*10^3*x5*x9+3.76*10^10*x5*x13-1.16*10^(-2)*x6*x9-1.45*x9*x12+8.54*10^5*x10*x13=0; 6.23*10^7*x4*x13-17.6*x2*x10+1.45*x9*x12-8.54*10^5*x10*x13-25.2*x10*x12+1.01*10^8*x11*x13=0; 25.2*x10*x12-1.01*10^8*x11*x13-1.7*10^10*x11+1.49*10^7*x13=0; 2*4.47*10^7*(x13)^2-3.03*10^4*(x12)^2-2.44*10^7*x1*x12+1.96*10^5*x3*x13-1.81*10^5*x1*x12+6.35*10^6*x4*x13-2.07*10^(-3)*x2*x12+3.32*10^11*x5*x13-1.45*x9*x12+8.54*10^5*x10*x13-25.2*x10*x12-1.01*10^8*x11*x13-1.10*10^7*x6*x12+3.09*10^10*x7*x13-1.72*10^7*x7*x12+4.48*10^5*x8*x13-1.35*10^6*x8*x12=0; x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12+x13=1; |
回复此楼» 收录本帖的淘帖专辑推荐
 常微分方程 |
» 猜你喜欢
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有235人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有15人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有3人回复
-
经费如何使用，花不完怎么办？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助大神解一个非线性方程组
已经有8人回复
- MATLAB中采用同fsolve求解复杂非线性方程组的问题
已经有1人回复
- 用1stOpt解出非线性方程组的解作为初值用Matlab求解，解不出来
已经有4人回复
- 关于用Broyden法解非线性方程组的问题
已经有2人回复
- fsolve函数与inline结合求解方程组（不需要写m文件）
已经有4人回复
- matlab求解一非线性方程组
已经有0人回复
- 用matlab解非线性方程组编程
已经有2人回复
- 求助matlab解非线性方程组并作图的问题，需要具体程序
已经有1人回复
- Matlab：求高人指点用matlab求解非线性方程组，解出来的值不收敛，提前结束
已经有0人回复
- 求助matlab解非线性方程。求代码
已经有5人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 求助关于用matlab求解三元二次非线性方程组的问题
已经有1人回复
- matlab求解非线性方程组 代码
已经有2人回复
- matlab 求解 二元非线性方程组
已经有0人回复
- 【分享】C常用程序算法集
已经有154人回复
15212223817
铁杆木虫 (著名写手)
- 应助: 15 (小学生)
- 金币: 7181.2
- 散金: 15
- 红花: 34
- 沙发: 6
- 帖子: 2865
- 在线: 685.2小时
- 虫号: 3097080
- 注册: 2014-03-28
- 专业: 控制理论与方法
2楼2017-08-22 09:17:47
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
1stOpt就可以很容易解决的，不过楼主的方程组不加服务约束的话是多解： 1： x1: 0 x2: 0 x3: 0.816647152090798 x4: 0 x5: 0 x6: 0 x7: 0 x8: 0 x9: 0 x10: 0.183352847908732 x11: 0 x12: 0 x13: 0 2： x1: 0 x2: 0 x3: 0.241538466580999 x4: 0 x5: 0.154368837115183 x6: 0 x7: 0 x8: 0.100107642087444 x9: 0 x10: 0.503985054216372 x11: 0 x12: 6.94727786959747E-29 x13: 0 3： x1: 0 x2: 0 x3: 0.560704408540411 x4: 1.16276785327303E-22 x5: 0.022221922230131 x6: 0.184928574582382 x7: 0.0586524090149989 x8: 0.0546504785595244 x9: 0 x10: 0.118842207072553 x11: 0 x12: 0 x13: 0 |
» 本帖已获得的红花（最新10朵）
yuanhongjuam3楼2017-08-22 09:23:12
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2017-08-22 09:24:33
yuanhongjuam
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2018.7
- 帖子: 35
- 在线: 394.1小时
- 虫号: 1457854
- 注册: 2011-10-24
- 性别: MM
- 专业: 理论和计算化学
5楼2017-08-22 09:29:57
yuanhongjuam
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2018.7
- 帖子: 35
- 在线: 394.1小时
- 虫号: 1457854
- 注册: 2011-10-24
- 性别: MM
- 专业: 理论和计算化学
6楼2017-08-22 09:30:00
梦在潇湘
铁杆木虫 (著名写手)
- 应助: 252 (大学生)
- 金币: 5778.8
- 散金: 2717
- 红花: 44
- 帖子: 2792
- 在线: 862小时
- 虫号: 1031369
- 注册: 2010-05-29
- 专业: 光学
7楼2017-08-22 09:31:21
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
8楼2017-08-22 09:31:23
yuanhongjuam
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 2018.7
- 帖子: 35
- 在线: 394.1小时
- 虫号: 1457854
- 注册: 2011-10-24
- 性别: MM
- 专业: 理论和计算化学
9楼2017-08-22 09:36:16
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
|
多解，每次计算结果不一样正常，即使范围约束在[0，1]之间，下面1stOpt10次计算10次结果： No. x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 1 0 0 0.0889221226824879 0 0 0 0.619075143009344 0.109599740672744 0.109691347939043 0.0727116456963814 0 0 0 2 0 0 0 0.584989159602328 0 0 0.0452828646877163 0.211765122616297 0.0122695987795364 0.145693254314122 0 1.02803604918793E-33 0 3 0 0 0 0 6.02099577131991E-8 0.351041849119781 0.13308452253865 0 0 0.515873568131611 0 1.71621373573786E-31 0 4 0 0 0.176361481099356 0 7.35585878988614E-5 0.434691371459343 0 0.0469852799261609 7.76505353833031E-25 0.341888308927241 1.11900267287777E-34 0 0 5 0 0 0 7.89809158858893E-9 0 0 0.864870318595261 0 0 0.135129673506648 0 3.6870871699034E-26 0 6 0 8.38315393517137E-20 0.021687121144077 0 8.89329081307484E-5 0.347615637426877 0.160065986721588 0 0 0.470542321799327 0 1.06193168104163E-25 0 7 0 2.36642962188249E-16 0.00574335874420293 0 0.606750024187809 0.152642263123701 0.0763600690555254 0.157164968595071 0 0.00133931629369072 0 0 0 8 0 0 0.201064715408182 0 0 0.798935284591817 0 0 8.35220037254079E-16 0 0 1.60675976936197E-24 0 9 0 1.44360273101539E-18 5.72011305382658E-6 1.11189021665092E-18 0.266722866682437 0.0336612414045716 0 0.0130473007419822 0 0.686562871057955 0 0 0 10 0 5.28987274661098E-18 0 0.435261093020827 0.00501720330021524 0.0232884993186285 6.14698548529995E-9 0 0 0.536433198213344 2.30495144996429E-27 0 1.61720219804766E-26 |
10楼2017-08-22 09:58:51
