24小时热门版块排行榜

返回列表

plf_jason

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 112.5
帖子: 47
在线: 30.7小时
虫号: 3497553
注册: 2014-10-25
专业: 结构工程

[求助] 1stopt求助已有1人参与

代码如下：
parameter b(1:10)[0,],c(1:10)[0,],d;
variable x,y1,y2,y3,y4,y5,y6,y7,y8;
sharedmodel;
function y1=4.966934+sum(i=1:10)((exp(d*0.007262)*b*(2*pi*x*c)^2/(1+(2*pi*x*c)^2)));
y2=sum(i=1:10)((exp(d*0.007262)*b*(2*pi*x*c)/(1+(2*pi*x*c)^2));
y3=5.818023+sum(i=1:10)((exp(d*0.007763)*b*(2*pi*x*c)^2/(1+(2*pi*x*c)^2)));
y4=sum(i=1:10)(exp(d*0.007763)*b*(2*pi*x*c)/(1+(2*pi*x*c)^2));
y5=6.398997+sum(i=1:10)((exp(d*0.03222)*b*(2*pi*x*c)^2/(1+(2*pi*x*c)^2)));
y6=sum(i=1:10)(exp(d*0.03222)*b*(2*pi*x*c)/(1+(2*pi*x*c)^2));
y7=8.018632+sum(i=1:10)((exp(d*0.14)*b*(2*pi*x*c)^2/(1+(2*pi*x*c)^2)));
y8=sum(i=1:10)(exp(d*0.14)*b*(2*pi*x*c)/(1+(2*pi*x*c)^2));
data;
1 13.6301 2.67121 46.169 10.707 50.517 11.059 65.853 14.171
1.25893 13.73917658 2.73016 50.108 11.359 54.29 11.265 75.379 16.476
1.58489 13.8190748 2.78912 50.088 11.295 55.407 11.83 72.225 15.115
1.99526 13.93969992 2.84808 50.383 11.504 55.145 12.259 71.24 16.01
2.51189 14.02105722 2.90704 50.947 11.78 55.561 12.807 72.351 15.368
3.16228 14.13315203 2.966 51.718 12.14 56.391 13.491 72.99 15.622
3.98107 14.24598972 3.02496 52.551 12.542 57.092 13.969 73.469 15.878
5.01187 14.6397 3.0839 53.158 13.217 57.387 14.518 73.648 16.01
6.30957 14.926244 3.3251336 54.25 13.445 58.669 15.107 74.385 16.407
7.94328 15.502236 3.5412557 56.122 14.017 59.681 15.682 75.12 16.686
10 15.9782 3.73542 57.183 14.544 61.208 16.31 75.905 17.091
12.5893 16.558436 3.9107003 59.207 15.237 62.038 16.958 77.91 17.709
15.8489 17.038644 4.0702904 60.965 15.695 63.687 17.633 79.321 18.571
19.9526 17.5723 4.21734 62.572 16.351 65.337 18.365 81.143 18.95
25.1189 17.903276 4.3548752 64.611 17.089 67.014 19.188 82.936 19.701
31.6228 18.2495 4.48618 66.282 17.922 68.728 20.134 83.78 20.35

点运行的时候出现错误提示：error: incompatible types: 'double' and 'boolean'
error: '.' expected but ';' found

求各位帮忙看一下，谢谢！

回复此楼

» 猜你喜欢

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

1楼 2020-01-11 13:08:15

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
独孤神宇: 金币+2, 鼓励交流 2020-01-13 08:12:10

第2个人公式左右括号都不匹配吧，修正后结果参考如下，总体效果不是太好，而且是多解：

Root of Mean Square Error (RMSE): 11.7182174871856
Sum of Squared Residual: 17576.5274978538
Correlation Coef. (R): 0.895839102004841
R-Square: 0.80252769668084
Adjusted R-Square: 0.520846619934372
Determination Coef. (DC): -26.0952519944772
F-Statistic: -0.253386716752826

Parameter Best Estimate
-------------------- -------------
b1 0.314991587740832
b2 8.52829896941523
b3 0.00145758610491548
b4 1.17582158304953
b5 0.48549825956673
b6 2.20762814740227
b7 1.72240185032738
b8 31.1284028597247
b9 3.57002654711787E-8
b10 964376.045463489
c1 0.57895425196262
c2 0.0158441065831744
c3 0.585511625329171
c4 0.0341885758959382
c5 6.94462999348029E-8
c6 1.54257713355745E-7
c7 1.94640242988118E-7
c8 0.579522487829381
c9 0.00381917073197707
c10 5.46858624798338E-8
d 4.85333043063028

赞一下