24小时热门版块排行榜

返回列表

无树可用

银虫 (初入文坛)

应助: 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;

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有102人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复

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的回帖