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

登录注册

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

24小时热门版块排行榜

返回列表

W57786

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 58.8
帖子: 21
在线: 2.7小时
虫号: 14723246
注册: 2019-03-25
专业: 控制理论与方法

[求助] 1stOpt共享模式运行以下代码已有1人参与

因为1stOpt1.5不支持共享模式，烦请哪位大神帮我运行以下代码，新人金币不多
因为每个代码有个4结果所以还请标明对应情况的运行结果，麻烦了。
以下格式供参考
（代码1:1-3-a
<运行结果>）
（代码1:1-4-b
<运行结果>

....
（代码2:2-4-a
<运行结果>）

（附件中代码与以下代码相同）
代码1：

CODE:

//请运行1-3-a,1-4-b,2-3-c,2-4-d Parameters a01=[0,0.001],a10=[0,0.01],a02,a03,a20,a30; ConstStr //v<119 //--------------------------1-------------------// k1=a02/2/a03, fx1=(1/a01)*ln(abs((2*a03*x/3.6+a02-a01)/(2*a03*x/3.6+a02+a01))), //b^2 >4ac fy1=(1/a01)*ln(abs((2*a03*y/3.6+a02-a01)/(2*a03*y/3.6+a02+a01))), sx1=(1/2/a03)*ln(abs((a03*(x/3.6)^2)+a02*(x/3.6)+(a02^2-a01^2)/4/a03)), sy1=(1/2/a03)*ln(abs((a03*(y/3.6)^2)+a02*(y/3.6)+(a02^2-a01^2)/4/a03)), f1_119=(1/a01)*ln(abs((2*a03*119/3.6+a02-a01)/(2*a03*117/3.6+a02+a01))), s1_119=(1/2/a03)*ln(abs((a03*(119/3.6)^2)+a02*(119/3.6)+(a02^2-a01^2)/4/a03)), //--------------------------2-------------------// //k2=a02/2/a03, //fx2=(2/a01)*atan((2*a03*x/3.6+a02)/a01), //b^2 <4ac //fy2=(2/a01)*atan((2*a03*y/3.6+a02)/a01), //sx2=(1/2/a03)*ln(abs((a03*(x/3.6)^2)+a02*(x/3.6)+(a02^2+a01^2)/4/a03)), //sy2=(1/2/a03)*ln(abs((a03*(y/3.6)^2)+a02*(y/3.6)+(a02^2+a01^2)/4/a03)), //f2_119=(2/a01)*atan((2*a03*119/3.6+a02)/a01), //s2_119=(1/2/a03)*ln(abs((a03*(119/3.6)^2)+a02*(119/3.6)+(a02^2+a01^2)/4/a03)), //v>119 //--------------------------3-------------------// //k3=a20/2/a30, //fx3=(1/a10)*ln(abs((2*a30*x/3.6+a20-a10)/(2*a30*x/3.6+a20+a10))), //b^2 >4ac //fy3=(1/a10)*ln(abs((2*a30*y/3.6+a20-a10)/(2*a30*y/3.6+a20+a10))), //sx3=(1/2/a30)*ln(abs((a30*(x/3.6)^2)+a20*(x/3.6)+(a20^2-a10^2)/4/a30)), //sy3=(1/2/a30)*ln(abs((a30*(y/3.6)^2)+a20*(y/3.6)+(a20^2-a10^2)/4/a30)), //f3_119=(1/a10)*ln(abs((2*a30*119/3.6+a20-a01)/(2*a30*117/3.6+a20+a01))), //s3_119=(1/2/a30)*ln(abs((a30*(119/3.6)^2)+a20*(119/3.6)+(a20^2-a10^2)/4/a30)); //--------------------------4-------------------// k4=a20/2/a30, fx4=(2/a10)*atan((2*a30*x/3.6+a20)/a10), //b^2 <4ac fy4=(2/a10)*atan((2*a30*y/3.6+a20)/a10), sx4=(1/2/a30)*ln(abs((a30*(x/3.6)^2)+a20*(x/3.6)+(a20^2+a10^2)/4/a30)), sy4=(1/2/a30)*ln(abs((a30*(y/3.6)^2)+a20*(y/3.6)+(a20^2+a10^2)/4/a30)), f4_119=(2/a10)*atan((2*a30*119/3.6+a20)/a10), s4_119=(1/2/a30)*ln(abs((a30*(119/3.6)^2)+a20*(119/3.6)+(a20^2+a10^2)/4/a30)); SharedModel; Variable x,y,z,z1; //--------------------------a-------------------// //Function z=if(y<119,(fy1-fx1),if((x<119)and(y>119),(f1_119-fx1+fy3-f3_119),(fy3-fx3))); //z1=if(y<119,(sy1-sx1-k1*(fy1-fx1)),if((x<119)and(y>119),(s1_119-sx1-k1*(f1_119-fx1)+sy3-s3_119-k3*(fy3-f3_119)),(sy3-sx3-k3*(fy3-fx3)))); //--------------------------b-------------------// Function z=if(y<119,(fy1-fx1),if((x<119)and(y>119),(f1_119-fx1+fy4-f4_119),(fy4-fx4))); z1=if(y<119,(sy1-sx1-k1*(fy1-fx1)),if((x<119)and(y>119),(s1_119-sx1-k1*(f1_119-fx1)+sy4-s4_119-k4*(fy4-f4_119)),(sy4-sx4-k4*(fy4-fx4)))); //--------------------------c-------------------// //Function z=if(y<119,(fy2-fx2),if((x<119)and(y>119),(f2_119-fx2+fy3-f3_119),(fy3-fx3))); //z1=if(y<119,(sy2-sx2-k2*(fy2-fx2)),if((x<119)and(y>119),(s2_119-sx2-k2*(f2_119-fx2)+sy3-s3_119-k3*(fy3-f3_119)),(sy3-sx3-k3*(fy3-fx3)))); //--------------------------d-------------------// //Function z=if(y<119,(fy2-fx2),if((x<119)and(y>119),(f2_119-fx2+fy4-f4_119),(fy4-fx4))); //z1=if(y<119,(sy2-sx2-k2*(fy2-fx2)),if((x<119)and(y>119),(s2_119-sx2-k2*(f2_119-fx2)+sy4-s4_119-k4*(fy4-f4_119)),(sy4-sx4-k4*(fy4-fx4)))); Data; x y z z1 0 45 26 164 0 80 47 535 0 160 98 2260 0 200 132 3974 0 250 193 7819 0 300 289 15213 45 80 21 368 45 160 72 2096 45 200 106 3810 45 250 167 7655 45 300 263 15049 80 160 51 1728 80 200 85 3442 80 250 146 7287 80 300 242 14681 160 200 34 1714 160 250 95 5590 200 250 61 3845 200 300 157 11239 250 300 96 7394

代码2：

CODE:

//请运行1-3-a,1-4-b,2-3-c,2-4-d Parameters a01=[0,0.1],a10=[0,0.1],a02,a03,a20,a30; ConstStr //v<106.7 //--------------------------1-------------------// //k1=a02/2/a03, //fx1=(1/a01)*ln(abs((2*a03*x/3.6+a02-a01)/(2*a03*x/3.6+a02+a01))), //b^2 >4ac //fy1=(1/a01)*ln(abs((2*a03*y/3.6+a02-a01)/(2*a03*y/3.6+a02+a01))), //sx1=(1/2/a03)*ln(abs((a03*(x/3.6)^2)+a02*(x/3.6)+(a02^2-a01^2)/4/a03)), //sy1=(1/2/a03)*ln(abs((a03*(y/3.6)^2)+a02*(y/3.6)+(a02^2-a01^2)/4/a03)), //f1_106=(1/a01)*ln(abs((2*a03*106.7/3.6+a02-a01)/(2*a03*117/3.6+a02+a01))), //s1_106=(1/2/a03)*ln(abs((a03*(106.7/3.6)^2)+a02*(106.7/3.6)+(a02^2-a01^2)/4/a03)), //-------------------------2--------------------// k2=a02/2/a03, fx2=(2/a01)*atan((2*a03*x/3.6+a02)/a01), //b^2 <4ac fy2=(2/a01)*atan((2*a03*y/3.6+a02)/a01), sx2=(1/2/a03)*ln(abs((a03*(x/3.6)^2)+a02*(x/3.6)+(a02^2+a01^2)/4/a03)), sy2=(1/2/a03)*ln(abs((a03*(y/3.6)^2)+a02*(y/3.6)+(a02^2+a01^2)/4/a03)), f2_106=(2/a01)*atan((2*a03*106.7/3.6+a02)/a01), s2_106=(1/2/a03)*ln(abs((a03*(106.7/3.6)^2)+a02*(106.7/3.6)+(a02^2+a01^2)/4/a03)), //v>106.7 //------------------------3----------------------// //k3=a20/2/a30, //fx3=(1/a10)*ln(abs((2*a30*x/3.6+a20-a10)/(2*a30*x/3.6+a20+a10))), //b^2 >4ac //fy3=(1/a10)*ln(abs((2*a30*y/3.6+a20-a10)/(2*a30*y/3.6+a20+a10))), //sx3=(1/2/a30)*ln(abs((a30*(x/3.6)^2)+a20*(x/3.6)+(a20^2-a10^2)/4/a30)), //sy3=(1/2/a30)*ln(abs((a30*(y/3.6)^2)+a20*(y/3.6)+(a20^2-a10^2)/4/a30)), //f3_106=(1/a10)*ln(abs((2*a30*106.7/3.6+a20-a01)/(2*a30*117/3.6+a20+a01))), //s3_106=(1/2/a30)*ln(abs((a30*(106.7/3.6)^2)+a20*(106.7/3.6)+(a20^2-a10^2)/4/a30)); //-----------------------4-----------------------// k4=a20/2/a30, fx4=(2/a10)*atan((2*a30*x/3.6+a20)/a10), //b^2 <4ac fy4=(2/a10)*atan((2*a30*y/3.6+a20)/a10), sx4=(1/2/a30)*ln(abs((a30*(x/3.6)^2)+a20*(x/3.6)+(a20^2+a10^2)/4/a30)), sy4=(1/2/a30)*ln(abs((a30*(y/3.6)^2)+a20*(y/3.6)+(a20^2+a10^2)/4/a30)), f4_106=(2/a10)*atan((2*a30*106.7/3.6+a20)/a10), s4_106=(1/2/a30)*ln(abs((a30*(106.7/3.6)^2)+a20*(106.7/3.6)+(a20^2+a10^2)/4/a30)); SharedModel; Variable x,y,z,z1; //-----------------------a------------------------// //Function z=if(y<106.7,(fy1-fx1),if((x<106.7)and(y>106.7),(f1_106-fx1+fy3-f3_106),(fy3-fx3))); //z1=if(y<106.7,(sy1-sx1-k1*(fy1-fx1)),if((x<106.7)and(y>106.7),(s1_106-sx1-k1*(f1_106-fx1)+sy3-s3_106-k3*(fy3-f3_106)),(sy3-sx3-k3*(fy3-fx3)))); //-----------------------b------------------------// //Function z=if(y<106.7,(fy1-fx1),if((x<106.7)and(y>106.7),(f1_106-fx1+fy4-f4_106),(fy4-fx4))); //z1=if(y<106.7,(sy1-sx1-k1*(fy1-fx1)),if((x<106.7)and(y>106.7),(s1_106-sx1-k1*(f1_106-fx1)+sy4-s4_106-k4*(fy4-f4_106)),(sy4-sx4-k4*(fy4-fx4)))); //-----------------------c------------------------ // //Function z=if(y<106.7,(fy2-fx2),if((x<106.7)and(y>106.7),(f2_106-fx2+fy3-f3_106),(fy3-fx3))); //z1=if(y<106.7,(sy2-sx2-k2*(fy2-fx2)),if((x<106.7)and(y>106.7),(s2_106-sx2-k2*(f2_106-fx2)+sy3-s3_106-k3*(fy3-f3_106)),(sy3-sx3-k3*(fy3-fx3)))); //------------------------d-------------------------// Function z=if(y<106.7,(fy2-fx2),if((x<106.7)and(y>106.7),(f2_106-fx2+fy4-f4_106),(fy4-fx4))); z1=if(y<106.7,(sy2-sx2-k2*(fy2-fx2)),if((x<106.7)and(y>106.7),(s2_106-sx2-k2*(f2_106-fx2)+sy4-s4_106-k4*(fy4-f4_106)),(sy4-sx4-k4*(fy4-fx4)))); Data; x y z z1 300 250 16 1227 300 200 32 2230 300 160 42 2733 300 80 63 3440 300 45 73 3611 300 0 85 3682 250 200 16 1003 250 160 26 1506 250 80 47 2213 250 45 57 2384 250 0 69 2455 200 160 10 503 200 80 31 1210 200 45 41 1381 200 0 53 1452 160 80 21 707 160 45 31 878 160 0 43 949 80 45 10 171 80 0 22 242 45 0 12 71

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有233人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

1楼2019-03-29 19:37:23

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

你的代码真是看的头晕，你不能整理一下再传上来么。

赞一下

回复此楼

数值计算

2楼2019-03-30 10:53:44

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

W57786

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 58.8
帖子: 21
在线: 2.7小时
虫号: 14723246
注册: 2019-03-25
专业: 控制理论与方法

引用回帖:

2楼: Originally posted by 独孤神宇 at 2019-03-30 10:53:44
你的代码真是看的头晕，你不能整理一下再传上来么。

第一次发帖，以为代码是自动排版，附件是代码文件

赞一下

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : code1_jiasu.mff

2019-03-31 09:51:55, 1.02 K

附件 2 : code2_jiansu.mff

2019-03-31 09:52:00, 1 K

3楼2019-03-31 09:52:35

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

jiasu
1-3-1/
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 8.10239134498161
残差平方和(SSR): 2625.94982028932
相关系数(R): 0.999986099670877
相关系数之平方(R^2): 0.999972199534974
修正R平方(Adj. R^2): 0.999964786232726
确定系数(DC): 0.999970170677951
F统计(F-Statistic): 51888.5147101725

参数最佳估算
-------------------- -------------
a01 2.98983559425701E-5
a10 2.74282980436971E-6
a02 7.31640885242362E-9
a03 -4.83081641055329E-10
a20 -0.0128232639363781
a30 4.88668055460729E-5

====== 结果输出 ======

文件: 数据文件-1
No 目标 z 计算 z 目标 z1 计算 z1
1 26 27.0206646852816 164 168.879153186254
2 47 48.0367391697461 535 533.741560895029
3 98 98.2568372514507 2260 2247.94652284579
4 132 132.89874132105 3974 3988.82756079974
5 193 193.729089762663 7819 7819.24752066881
6 289 289.855614843433 15213 15218.8060481488
7 21 21.0160744844645 368 364.862407708775
8 72 71.2361725661691 2096 2079.06736965953
9 106 105.878076635768 3810 3819.94840761348
10 167 166.708425077381 7655 7650.36836748255
11 263 262.834950158151 15049 15049.9268949625
12 51 50.2200980817046 1728 1714.20496195076
13 85 84.8620021513035 3442 3455.08599990471
14 146 145.692350592917 7287 7285.50595977378
15 242 241.818875673687 14681 14685.0644872538
16 34 34.6419040695989 1714 1740.88103795395
17 95 95.4722525112121 5590 5571.30099782302
18 61 60.8303484416132 3845 3830.41995986907
19 157 156.956873522383 11239 11229.978487349
20 96 96.12652508077 7394 7399.55852747998

******************************************

1-4-b/

优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 4.3602971227511
残差平方和(SSR): 760.487639946861
相关系数(R): 0.999971859014608
相关系数之平方(R^2): 0.999943718821132
修正R平方(Adj. R^2): 0.99992871145491
确定系数(DC): 0.999914763688369
F统计(F-Statistic): 16459.4643127049

参数最佳估算
-------------------- -------------
a01 2.28964248672128E-5
a10 0.0041527624253373
a02 7.38721388614751E-5
a03 -5.49440426444815E-7
a20 -0.0155201603770716
a30 6.96683382123042E-5

====== 结果输出 ======

文件: 数据文件-1
No 目标 z 计算 z 目标 z1 计算 z1
1 26 27.0187018934527 164 168.889979608791
2 47 48.0497490599691 535 534.024930327789
3 98 97.3620816628271 2260 2256.89089327256
4 132 131.566961192452 3974 3976.3982674351
5 193 192.682450771443 7819 7825.72066559193
6 289 288.734142049596 15213 15216.2166030194
7 21 21.0310471665163 368 365.134950718998
8 72 70.3433797693743 2096 2088.00091366377
9 106 104.548259298999 3810 3807.50828782631
10 167 165.663748877991 7655 7656.83068598314
11 263 261.715440156144 15049 15047.3266234106
12 51 49.312332602858 1728 1722.86596294478
13 85 83.5172121324828 3442 3442.37333710732
14 146 144.632701711474 7287 7291.69573526415
15 242 240.684392989627 14681 14682.1916726916
16 34 34.2048795296248 1714 1719.50737416254
17 95 95.3203691086164 5590 5568.82977231937
18 61 61.1154895789916 3845 3849.32239815683
19 157 157.167180857144 11239 11239.8183355843
20 96 96.0516912781528 7394 7390.49593742745

赞一下

回复此楼

数值计算

4楼2019-03-31 10:37:40

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
W57786: 金币+15, ★★★★★最佳答案 2019-04-02 22:24:25

内容已删除

赞一下(1人)

回复此楼

数值计算

5楼2019-03-31 10:46:51

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

W57786

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 58.8
帖子: 21
在线: 2.7小时
虫号: 14723246
注册: 2019-03-25
专业: 控制理论与方法

引用回帖:

5楼: Originally posted by 独孤神宇 at 2019-03-31 10:46:51
jiasu
2-3-c/

均方差(RMSE): 8.00205895360656
残差平方和(SSR): 2561.3178998798
相关系数(R): 0.999991422836794
相关系数之平方(R^2): 0.999982845747155
修正R平方(Adj. R^2): 0.999978271322797
确定 ...

大神，上次代码有误，能不能帮我运行下新代码，还是15金币。

赞一下

回复此楼