版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (2284)
>虫友互识 (148)
>导师招生 (131)
>考研 (114)
>休闲灌水 (91)
>考博 (28)
>论文投稿 (21)
>基金申请 (16)
>文献求助 (16)
>教师之家 (14)
>硕博家园 (14)
>找工作 (14)
>博后之家 (12)
>公派出国 (7)
>论文道贺祈福 (6)
>计算模拟 (4)

返回列表

无树可用

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 409.6
帖子: 27
在线: 45.9小时
虫号: 2816186
注册: 2013-11-20
专业: 系统科学与系统工程

[求助] 求高版本1stopt求解非线性方程组已有1人参与

代码如下：
Parameters x(1:89);
Constants a=80,b=1,qw(1:2)=[9,10],rho(1:2)=[0.5,0.5],cqmax(1:2)=[50,50],wqsmax=50,c(1:2)=[20,20];
Function c1 - x57 + x59 - a + x55*b + b*(x55 + x56 + x69) - x31*((x1 - x55)/(x5^2 + (x1 - x55)^2)^(1/2) + 1) - x32*((x2 - x55)/(x6^2 + (x2 - x55)^2)^(1/2) + 1) - x39*((x1 - x55 + cqmax1)/((x1 - x55 + cqmax1)^2 + x9^2)^(1/2) - 1) - x40*((x2 - x55 + cqmax1)/((x2 - x55 + cqmax1)^2 + x10^2)^(1/2) - 1)=0;
c2 - x58 + x60 - a + x56*b + b*(x55 + x56 + x69) - x37*((x3 - x56)/(x7^2 + (x3 - x56)^2)^(1/2) + 1) - x38*((x4 - x56)/(x8^2 + (x4 - x56)^2)^(1/2) + 1) - x45*((x3 - x56 + cqmax2)/((x3 - x56 + cqmax2)^2 + x11^2)^(1/2) - 1) - x46*((x4 - x56 + cqmax2)/((x4 - x56 + cqmax2)^2 + x12^2)^(1/2) - 1)=0;
x23 - x61=0;
x24 - x62=0;
x29 - x63=0;
x30 - x64=0;
x17 + x31*((x1 - x55)/(x5^2 + (x1 - x55)^2)^(1/2) + 1) - rho1*(c1 - x21) + x47*(x69 - qw1 + (x1*(x69 - qw1)^2)/(x13^2 + x1^2*(x69 - qw1)^2)^(1/2)) + x39*((x1 - x55 + cqmax1)/((x1 - x55 + cqmax1)^2 + x9^2)^(1/2) - 1)=0;
x18 + x32*((x2 - x55)/(x6^2 + (x2 - x55)^2)^(1/2) + 1) - rho2*(c1 - x22) + x48*(x69 - qw2 + (x2*(x69 - qw2)^2)/(x14^2 + x2^2*(x69 - qw2)^2)^(1/2)) + x40*((x2 - x55 + cqmax1)/((x2 - x55 + cqmax1)^2 + x10^2)^(1/2) - 1)=0;
x19 + x37*((x3 - x56)/(x7^2 + (x3 - x56)^2)^(1/2) + 1) - rho1*(c2 - x21) + x53*(x69 - qw1 + (x3*(x69 - qw1)^2)/(x15^2 + x3^2*(x69 - qw1)^2)^(1/2)) + x45*((x3 - x56 + cqmax2)/((x3 - x56 + cqmax2)^2 + x11^2)^(1/2) - 1)=0;
x20 + x38*((x4 - x56)/(x8^2 + (x4 - x56)^2)^(1/2) + 1) - rho2*(c2 - x22) + x54*(x69 - qw2 + (x4*(x69 - qw2)^2)/(x16^2 + x4^2*(x69 - qw2)^2)^(1/2)) + x46*((x4 - x56 + cqmax2)/((x4 - x56 + cqmax2)^2 + x12^2)^(1/2) - 1)=0;
x17 + x33*((x3 - x56)/(x7^2 + (x3 - x56)^2)^(1/2) + 1) + x49*(x69 - qw1 + (x3*(x69 - qw1)^2)/(x15^2 + x3^2*(x69 - qw1)^2)^(1/2)) + x41*((x3 - x56 + cqmax2)/((x3 - x56 + cqmax2)^2 + x11^2)^(1/2) - 1)=0;
x18 + x34*((x4 - x56)/(x8^2 + (x4 - x56)^2)^(1/2) + 1) + x50*(x69 - qw2 + (x4*(x69 - qw2)^2)/(x16^2 + x4^2*(x69 - qw2)^2)^(1/2)) + x42*((x4 - x56 + cqmax2)/((x4 - x56 + cqmax2)^2 + x12^2)^(1/2) - 1)=0;
x19 + x35*((x1 - x55)/(x5^2 + (x1 - x55)^2)^(1/2) + 1) + x51*(x69 - qw1 + (x1*(x69 - qw1)^2)/(x13^2 + x1^2*(x69 - qw1)^2)^(1/2)) + x43*((x1 - x55 + cqmax1)/((x1 - x55 + cqmax1)^2 + x9^2)^(1/2) - 1)=0;
x20 + x36*((x2 - x55)/(x6^2 + (x2 - x55)^2)^(1/2) + 1) + x52*(x69 - qw2 + (x2*(x69 - qw2)^2)/(x14^2 + x2^2*(x69 - qw2)^2)^(1/2)) + x44*((x2 - x55 + cqmax1)/((x2 - x55 + cqmax1)^2 + x10^2)^(1/2) - 1)=0;
x1*rho1 - x25 - x23=0;
x2*rho2 - x26 - x24=0;
x3*rho1 - x29 - x27=0;
x4*rho2 - x30 - x28=0;
x31*(x5/(x5^2 + (x1 - x55)^2)^(1/2) - 1) - x23=0;
x32*(x6/(x6^2 + (x2 - x55)^2)^(1/2) - 1) - x24=0;
x33*(x7/(x7^2 + (x3 - x56)^2)^(1/2) - 1) - x25=0;
x34*(x8/(x8^2 + (x4 - x56)^2)^(1/2) - 1) - x26=0;
x35*(x5/(x5^2 + (x1 - x55)^2)^(1/2) - 1) - x27=0;
x36*(x6/(x6^2 + (x2 - x55)^2)^(1/2) - 1) - x28=0;
x37*(x7/(x7^2 + (x3 - x56)^2)^(1/2) - 1) - x29=0;
x38*(x8/(x8^2 + (x4 - x56)^2)^(1/2) - 1) - x30=0;
x23 + x39*(x9/((x1 - x55 + cqmax1)^2 + x9^2)^(1/2) - 1)=0;
x24 + x40*(x10/((x2 - x55 + cqmax1)^2 + x10^2)^(1/2) - 1)=0;
x25 + x41*(x11/((x3 - x56 + cqmax2)^2 + x11^2)^(1/2) - 1)=0;
x26 + x42*(x12/((x4 - x56 + cqmax2)^2 + x12^2)^(1/2) - 1)=0;
x27 + x43*(x9/((x1 - x55 + cqmax1)^2 + x9^2)^(1/2) - 1)=0;
x28 + x44*(x10/((x2 - x55 + cqmax1)^2 + x10^2)^(1/2) - 1)=0;
x29 + x45*(x11/((x3 - x56 + cqmax2)^2 + x11^2)^(1/2) - 1)=0;
x30 + x46*(x12/((x4 - x56 + cqmax2)^2 + x12^2)^(1/2) - 1)=0;
x47*(x13/(x13^2 + x1^2*(x69 - qw1)^2)^(1/2) - 1) - x23*(x69 - qw1)=0;
x48*(x14/(x14^2 + x2^2*(x69 - qw2)^2)^(1/2) - 1) - x24*(x69 - qw2)=0;
x49*(x15/(x15^2 + x3^2*(x69 - qw1)^2)^(1/2) - 1) - x25*(x69 - qw1)=0;
x50*(x16/(x16^2 + x4^2*(x69 - qw2)^2)^(1/2) - 1) - x26*(x69 - qw2)=0;
x51*(x13/(x13^2 + x1^2*(x69 - qw1)^2)^(1/2) - 1) - x27*(x69 - qw1)=0;
x52*(x14/(x14^2 + x2^2*(x69 - qw2)^2)^(1/2) - 1) - x28*(x69 - qw2)=0;
x53*(x15/(x15^2 + x3^2*(x69 - qw1)^2)^(1/2) - 1) - x29*(x69 - qw1)=0;
x54*(x16/(x16^2 + x4^2*(x69 - qw2)^2)^(1/2) - 1) - x30*(x69 - qw2)=0;
x5 - x9 + x21 - x65 + x13*(x69 - qw1)=0;
x6 - x10 + x22 - x66 + x14*(x69 - qw2)=0;
x7 - x11 + x21 - x67 + x15*(x69 - qw1)=0;
x8 - x12 + x22 - x68 + x16*(x69 - qw2)=0;
qw1 - x3 - x69 - x1=0;
qw2 - x4 - x69 - x2=0;
x5 - x1 + x55 - (x5^2 + (x1 - x55)^2)^(1/2)=0;
x6 - x2 + x55 - (x6^2 + (x2 - x55)^2)^(1/2)=0;
x7 - x3 + x56 - (x7^2 + (x3 - x56)^2)^(1/2)=0;
x8 - x4 + x56 - (x8^2 + (x4 - x56)^2)^(1/2)=0;
x1 + x9 - x55 + cqmax1 - ((x1 - x55 + cqmax1)^2 + x9^2)^(1/2)=0;
x2 + x10 - x55 + cqmax1 - ((x2 - x55 + cqmax1)^2 + x10^2)^(1/2)=0;
x3 + x11 - x56 + cqmax2 - ((x3 - x56 + cqmax2)^2 + x11^2)^(1/2)=0;
x4 + x12 - x56 + cqmax2 - ((x4 - x56 + cqmax2)^2 + x12^2)^(1/2)=0;
x13 - x1*(x69 - qw1) - (x13^2 + x1^2*(x69 - qw1)^2)^(1/2)=0;
x14 - x2*(x69 - qw2) - (x14^2 + x2^2*(x69 - qw2)^2)^(1/2)=0;
x15 - x3*(x69 - qw1) - (x15^2 + x3^2*(x69 - qw1)^2)^(1/2)=0;
x16 - x4*(x69 - qw2) - (x16^2 + x4^2*(x69 - qw2)^2)^(1/2)=0;
x55 + x57 - (x55^2 + x57^2)^(1/2)=0;
x56 + x58 - (x56^2 + x58^2)^(1/2)=0;
x59 - x55 + cqmax1 - (x59^2 + (x55 - cqmax1)^2)^(1/2)=0;
x60 - x56 + cqmax2 - (x60^2 + (x56 - cqmax2)^2)^(1/2)=0;
x61 + x65 - (x61^2 + x65^2)^(1/2)=0;
x62 + x66 - (x62^2 + x66^2)^(1/2)=0;
x63 + x67 - (x63^2 + x67^2)^(1/2)=0;
x64 + x68 - (x64^2 + x68^2)^(1/2)=0;
x71 - x70 + x72 + x73 - a - x13*x74 - x14*x75 - x15*x76 - x16*x77 + x69*b + b*(x55 + x56 + x69) + x86*(x1 + (x1^2*(x69 - qw1))/(x13^2 + x1^2*(x69 - qw1)^2)^(1/2)) + x87*(x2 + (x2^2*(x69 - qw2))/(x14^2 + x2^2*(x69 - qw2)^2)^(1/2)) + x88*(x3 + (x3^2*(x69 - qw1))/(x15^2 + x3^2*(x69 - qw1)^2)^(1/2)) + x89*(x4 + (x4^2*(x69 - qw2))/(x16^2 + x4^2*(x69 - qw2)^2)^(1/2))=0;
x72 + x78*((x1 - x55)/(x5^2 + (x1 - x55)^2)^(1/2) + 1) - x21*rho1 + x86*(x69 - qw1 + (x1*(x69 - qw1)^2)/(x13^2 + x1^2*(x69 - qw1)^2)^(1/2)) + x82*((x1 - x55 + cqmax1)/((x1 - x55 + cqmax1)^2 + x9^2)^(1/2) - 1)=0;
x73 + x79*((x2 - x55)/(x6^2 + (x2 - x55)^2)^(1/2) + 1) - x22*rho2 + x87*(x69 - qw2 + (x2*(x69 - qw2)^2)/(x14^2 + x2^2*(x69 - qw2)^2)^(1/2)) + x83*((x2 - x55 + cqmax1)/((x2 - x55 + cqmax1)^2 + x10^2)^(1/2) - 1)=0;
x72 + x80*((x3 - x56)/(x7^2 + (x3 - x56)^2)^(1/2) + 1) - x21*rho1 + x88*(x69 - qw1 + (x3*(x69 - qw1)^2)/(x15^2 + x3^2*(x69 - qw1)^2)^(1/2)) + x84*((x3 - x56 + cqmax2)/((x3 - x56 + cqmax2)^2 + x11^2)^(1/2) - 1)=0;
x73 + x81*((x4 - x56)/(x8^2 + (x4 - x56)^2)^(1/2) + 1) - x22*rho2 + x89*(x69 - qw2 + (x4*(x69 - qw2)^2)/(x16^2 + x4^2*(x69 - qw2)^2)^(1/2)) + x85*((x4 - x56 + cqmax2)/((x4 - x56 + cqmax2)^2 + x12^2)^(1/2) - 1)=0;
- x74 - x76 - rho1*(x1 + x3)=0;
- x75 - x77 - rho2*(x2 + x4)=0;
x78*(x5/(x5^2 + (x1 - x55)^2)^(1/2) - 1) - x74=0;
x79*(x6/(x6^2 + (x2 - x55)^2)^(1/2) - 1) - x75=0;
x80*(x7/(x7^2 + (x3 - x56)^2)^(1/2) - 1) - x76=0;
x81*(x8/(x8^2 + (x4 - x56)^2)^(1/2) - 1) - x77=0;
x74 + x82*(x9/((x1 - x55 + cqmax1)^2 + x9^2)^(1/2) - 1)=0;
x75 + x83*(x10/((x2 - x55 + cqmax1)^2 + x10^2)^(1/2) - 1)=0;
x76 + x84*(x11/((x3 - x56 + cqmax2)^2 + x11^2)^(1/2) - 1)=0;
x77 + x85*(x12/((x4 - x56 + cqmax2)^2 + x12^2)^(1/2) - 1)=0;
x86*(x13/(x13^2 + x1^2*(x69 - qw1)^2)^(1/2) - 1) - x74*(x69 - qw1)=0;
x87*(x14/(x14^2 + x2^2*(x69 - qw2)^2)^(1/2) - 1) - x75*(x69 - qw2)=0;
x88*(x15/(x15^2 + x3^2*(x69 - qw1)^2)^(1/2) - 1) - x76*(x69 - qw1)=0;
x89*(x16/(x16^2 + x4^2*(x69 - qw2)^2)^(1/2) - 1) - x77*(x69 - qw2)=0;
x69 + x70 - (x69^2 + x70^2)^(1/2)=0;
x71 - x69 + wqsmax - (x71^2 + (x69 - wqsmax)^2)^(1/2)=0;

回复此楼

» 猜你喜欢

计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值已经有7人回复
2026年循环经济功能材料国际会议（ICFMCE 2026）已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有99人回复
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）已经有0人回复
哈尔滨理工大学物理系招收物理学考研调剂已经有5人回复
0702一志愿吉大B区求调剂已经有5人回复
求调剂已经有0人回复
0702一志愿吉大B区求调剂有论文已经有0人回复
请问还有没有用Latex写文章的小伙伴们？已经有0人回复
光学工程学硕调剂信息已经有9人回复

1楼 2015-11-10 09:57:28

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

lipenggg

专家顾问 (正式写手)

专家经验: +65
应助: 122 (高中生)
金币: 4841.3
红花: 43
帖子: 722
在线: 435.7小时
虫号: 3929201
注册: 2015-06-17
性别: GG
专业: 流体力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
无树可用: 金币+10, ★有帮助 2015-11-12 15:06:21

Constants  a=80,b=1,qw(1:2)=[9,10],rho(1:2)=[0.5,0.5],cqmax(1:2)=[50,50],wqsmax=50,c(1:2)=[20,20];
改为：
Constant a=80,
         b=1,
         qw1=9,
         qw2=10,
         rho1=0.5,
         rho2=0.5,
         cqmax1=50,
         cqmax2=50,
         wqsmax=50,
         c1=20,
         c2=20;
破解版也可以运行。

赞一下

回复此楼

一团和气

2楼2015-11-11 13:45:18

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

无树可用

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 409.6
帖子: 27
在线: 45.9小时
虫号: 2816186
注册: 2013-11-20
专业: 系统科学与系统工程

引用回帖:

2楼: Originally posted by lipenggg at 2015-11-11 13:45:18
Constants  a=80,b=1,qw(1:2)=,rho(1:2)=,cqmax(1:2)=,wqsmax=50,c(1:2)=;
改为：
Constant a=80,
         b=1,
         qw1=9,
         qw2=10,
         rho1=0.5,
         rho2=0.5,
...

破解版跑出来的结果不合理

赞一下

回复此楼

3楼2015-11-12 14:34:29

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

河流123

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 298.6
红花: 1
帖子: 31
在线: 18.9小时
虫号: 4000668
注册: 2015-08-01
专业: 运筹与管理

引用回帖:

这样的方程组用什么算法解呢望指教

赞一下

回复此楼

4楼2018-12-30 13:19:13

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

相关版块跳转我要订阅楼主无树可用的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 本人考085602 化学工程专硕 +13	不知道叫什么！ 2026-03-15	15/750	2026-03-17 07:43 by laoshidan
[考研] 286求调剂 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[文学芳草园] 伙伴们，祝我生日快乐吧 +17	myrtle 2026-03-10	26/1300	2026-03-16 18:32 by 青橙Ln
[考研] 333求调剂 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 289求调剂 +5	步川酷紫123 2026-03-11	5/250	2026-03-15 00:45 by kruisytel
[考研] 288求调剂 +4	奇点0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 297一志愿上交085600求调剂 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 材料与化工（0856）304求B区调剂 +7	邱gl 2026-03-10	11/550	2026-03-14 12:18 by 邱gl
[考研] 297求调剂 +4	学海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 热情沙漠
[考研] 0856材料与化工309分求调剂 +6	ZyZy…… 2026-03-10	6/300	2026-03-14 00:38 by JourneyLucky
[考研] 336求调剂 +6	Iuruoh 2026-03-11	6/300	2026-03-13 22:06 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 304求调剂 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 四川大学085601材料工程专硕初试294求调剂 +4	祝我们好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 332求调剂 +3	Zz版 2026-03-13	3/150	2026-03-13 20:36 by 18595523086
[考研] 302求调剂 +6	负心者当诛 2026-03-11	6/300	2026-03-13 16:11 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8	努力奋斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[论文投稿] 投稿问题 5+4	星光灿烂xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考研] 08食品或轻工求调剂，本科发表3篇sci一区top论文，一志愿南师大食品科学与工程 +3	我是一个兵， 2026-03-10	3/150	2026-03-13 10:21 by Yuyi.

24小时热门版块排行榜

无树可用

[求助] 求高版本1stopt求解非线性方程组 已有1人参与

» 猜你喜欢

lipenggg

【答案】应助回帖

无树可用

河流123

[求助] 求高版本1stopt求解非线性方程组已有1人参与