24小时热门版块排行榜    

查看: 840  |  回复: 1

plf_jason

铜虫 (初入文坛)

[求助] 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

求各位帮忙看一下,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +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
2楼2020-01-12 20:45:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 plf_jason 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +41 医学老男孩 2026-08-13 94/4700 2026-08-18 10:20 by 医学老男孩
[论文投稿] 投稿咨询 +3 wwm09 2026-08-17 4/200 2026-08-18 09:52 by xs74101122
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 5/250 2026-08-18 08:03 by 基诺咪客
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见