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

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dzheng93

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 829.1
红花: 2
帖子: 282
在线: 224.3小时
虫号: 4368115
注册: 2016-01-19
专业: 结构工程

[求助] 求助！为什么1stopt一直显示函数表达式有误已有2人参与

function y=p2*log(p3*0.00157)-p2*log(p3*0.00157*10^(y/p1)-x)-((2.303*8.3114*298.15)/(4*96494))*log(1-p4)+((2.303*8.3114*298.15)/(4*96494))*log(1-(p3*0.00157*10^(y/p1)-x)/(p3*0.00157/p4))
Parameters p1,p2,p3,p4;
Variable x,y;
Data；
//x,y
-0.0000152 0.0016
-0.0000089 0.0014
-0.0000045 0.0009
-0.0000012 0.0003
-0.0000016 -0.0004
-0.0000041 -0.0011
-0.0000064 -0.0019
-0.0000086 -0.0026
-0.0000105 -0.0034
-0.0000123 -0.0042
-0.0000141 -0.0050
-0.0000157 -0.0058
-0.0000173 -0.0066
-0.0000188 -0.0074
-0.0000203 -0.0083
-0.0000217 -0.0091
-0.0000231 -0.0099
-0.0000244 -0.0108
-0.0000258 -0.0116
-0.0000271 -0.0125
-0.0000283 -0.0133
-0.0000296 -0.0142
-0.0000305 -0.0151
-0.0000318 -0.0159
-0.0000329 -0.0168
-0.0000341 -0.0177
-0.0000352 -0.0185
-0.0000364 -0.0194
-0.0000375 -0.0203
-0.0000386 -0.0211
-0.0000397 -0.0220
-0.0000408 -0.0229
-0.0000418 -0.0238
-0.0000429 -0.0246
-0.0000440 -0.0255
-0.0000450 -0.0264
-0.0000460 -0.0273
-0.0000471 -0.0282
-0.0000481 -0.0290
-0.0000491 -0.0299
-0.0000501 -0.0308
-0.0000512 -0.0317
-0.0000522 -0.0326
-0.0000532 -0.0335
-0.0000542 -0.0343
-0.0000552 -0.0352
-0.0000562 -0.0361
-0.0000572 -0.0370
-0.0000581 -0.0379
-0.0000592 -0.0388
-0.0000602 -0.0397
-0.0000611 -0.0405
-0.0000621 -0.0414
-0.0000631 -0.0423
-0.0000641 -0.0432
-0.0000650 -0.0441
-0.0000661 -0.0450
-0.0000671 -0.0459
-0.0000681 -0.0467
-0.0000690 -0.0476
-0.0000700 -0.0485
-0.0000710 -0.0494
-0.0000720 -0.0503
-0.0000730 -0.0512
-0.0000739 -0.0521
-0.0000750 -0.0529
-0.0000760 -0.0538
-0.0000769 -0.0547
-0.0000778 -0.0556
-0.0000788 -0.0565
-0.0000798 -0.0574
-0.0000807 -0.0583
-0.0000817 -0.0592
-0.0000828 -0.0600
-0.0000837 -0.0609
-0.0000847 -0.0618
-0.0000857 -0.0627
-0.0000866 -0.0636
-0.0000875 -0.0645
-0.0000885 -0.0654
-0.0000895 -0.0663
-0.0000905 -0.0672
-0.0000915 -0.0680
-0.0000924 -0.0689
-0.0000935 -0.0698
-0.0000944 -0.0707
-0.0000955 -0.0716
-0.0000964 -0.0725
-0.0000974 -0.0734
-0.0000983 -0.0743
-0.0000993 -0.0751
-0.0001003 -0.0760
-0.0001013 -0.0769
-0.0001023 -0.0778
-0.0001032 -0.0787
-0.0001042 -0.0796
-0.0001053 -0.0804
-0.0001063 -0.0813
-0.0001073 -0.0822
-0.0001082 -0.0831
-0.0001092 -0.0840
-0.0001102 -0.0849
-0.0001112 -0.0858
-0.0001123 -0.0866
-0.0001133 -0.0875
-0.0001143 -0.0884
-0.0001153 -0.0893
-0.0001163 -0.0902
-0.0001174 -0.0911
-0.0001185 -0.0919
-0.0001195 -0.0928
-0.0001206 -0.0937
-0.0001216 -0.0946
-0.0001227 -0.0954
-0.0001238 -0.0963
-0.0001247 -0.0972
-0.0001259 -0.0981
-0.0001270 -0.0989
-0.0001280 -0.0998
-0.0001290 -0.1007
-0.0001301 -0.1016
-0.0001312 -0.1025
-0.0001323 -0.1033
-0.0001335 -0.1042
-0.0001345 -0.1051
-0.0001356 -0.1060
-0.0001367 -0.1068
-0.0001378 -0.1077
-0.0001389 -0.1086
-0.0001400 -0.1094
-0.0001411 -0.1103
-0.0001423 -0.1112
-0.0001434 -0.1121
-0.0001445 -0.1129
-0.0001457 -0.1138
-0.0001468 -0.1147
-0.0001480 -0.1155
-0.0001490 -0.1164
-0.0001503 -0.1173
-0.0001514 -0.1181
-0.0001526 -0.1190
-0.0001538 -0.1199
-0.0001550 -0.1207
-0.0001562 -0.1216
-0.0001574 -0.1224
-0.0001586 -0.1233
-0.0001598 -0.1242
-0.0001609 -0.1250
-0.0001621 -0.1259
-0.0001634 -0.1267
-0.0001646 -0.1276
-0.0001659 -0.1285
-0.0001671 -0.1293
-0.0001683 -0.1302
-0.0001696 -0.1310
-0.0001708 -0.1319
-0.0001721 -0.1327
-0.0001734 -0.1336
-0.0001747 -0.1344
-0.0001759 -0.1353
-0.0001772 -0.1362
-0.0001784 -0.1370
-0.0001798 -0.1379
-0.0001811 -0.1387
-0.0001824 -0.1396
-0.0001837 -0.1404
-0.0001851 -0.1412
-0.0001864 -0.1421
-0.0001877 -0.1429
-0.0001891 -0.1438
-0.0001904 -0.1446
-0.0001918 -0.1455
-0.0001932 -0.1463
-0.0001945 -0.1472
-0.0001959 -0.1480
-0.0001973 -0.1488
-0.0001986 -0.1497
-0.0002000 -0.1505
-0.0002014 -0.1514
-0.0002028 -0.1522
-0.0002042 -0.1530
-0.0002056 -0.1539
-0.0002070 -0.1547
-0.0002085 -0.1555
-0.0002098 -0.1564
-0.0002113 -0.1572
-0.0002127 -0.1581
-0.0002141 -0.1589
-0.0002154 -0.1597
-0.0002168 -0.1606
-0.0002182 -0.1614
-0.0002197 -0.1622
-0.0002212 -0.1631
-0.0002227 -0.1639
-0.0002242 -0.1647
-0.0002256 -0.1656
-0.0002271 -0.1664
-0.0002286 -0.1672
-0.0002301 -0.1680
-0.0002316 -0.1689

回复此楼

» 猜你喜欢

数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
物理学I论文润色/翻译怎么收费? 已经有66人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有1人回复
投稿chemical physical letters不送审？已经有2人回复
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生已经有22人回复
杂志要求提供official language-editing certificate 已经有4人回复
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组已经有15人回复

1楼 2016-10-11 16:33:42

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

lu_yu_lan

新虫 (初入文坛)

应助: 3 (幼儿园)
金币: 101
散金: 5
红花: 1
帖子: 43
在线: 26.4小时
虫号: 4935903
注册: 2016-08-22

【答案】应助回帖

引用回帖:

7楼: Originally posted by lipenggg at 2016-10-12 09:45:47
迭代数: 39
计算用时(时:分:秒:微秒): 00:01:33:462
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 1.29237480386255E-15
残差平方和(SSE): 3.340 ...

似乎与这个结果稍有差距吧。

结果（p1,p2,p3,p4，目标函数值）：
-6.553501528181543e-003 6.553501528181543e-003 244620310527.4335 -2.258610642562602e-040 7.415825341858829e-032

赞一下

回复此楼

9楼2016-10-12 16:58:43

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

查看全部 22 个回答

dzheng93

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 829.1
红花: 2
帖子: 282
在线: 224.3小时
虫号: 4368115
注册: 2016-01-19
专业: 结构工程

想求出残差最小的四个参数值

赞一下(1人)

回复此楼

2楼2016-10-11 16:34:53

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

lu_yu_lan

新虫 (初入文坛)

应助: 3 (幼儿园)
金币: 101
散金: 5
红花: 1
帖子: 43
在线: 26.4小时
虫号: 4935903
注册: 2016-08-22

【答案】应助回帖

感谢参与，应助指数 +1

用OpenLu求解：

CODE:

!!!using["luopt","math"]; //使用命名空间

f(p1,p2,p3,p4 :: x,y)=

{

   sum[(p2*lg(p3*0.00157)-p2.*lg(p3*(0.00157).*(10.0).^(y./p1)-x)-((2.303*8.3114*298.15)/(4*96494))*lg(1-p4)+((2.303*8.3114*298.15)/(4*96494)).*lg(1.0-(p3*(0.00157).*(10.0).^(y./p1)-x)./(p3*0.00157/p4))-y).^2.0]

};

main(:xy: x,y)=

{

   xy=matrix{   //存放实验数据

"

 -0.0000152         0.0016

 -0.0000089         0.0014

 -0.0000045         0.0009

 -0.0000012         0.0003

 -0.0000016         -0.0004

 -0.0000041         -0.0011

 -0.0000064         -0.0019

 -0.0000086         -0.0026

 -0.0000105         -0.0034

 -0.0000123         -0.0042

 -0.0000141         -0.0050

 -0.0000157         -0.0058

 -0.0000173         -0.0066

 -0.0000188         -0.0074

 -0.0000203         -0.0083

 -0.0000217         -0.0091

 -0.0000231         -0.0099

 -0.0000244         -0.0108

 -0.0000258         -0.0116

 -0.0000271         -0.0125

 -0.0000283         -0.0133

 -0.0000296         -0.0142

 -0.0000305         -0.0151

 -0.0000318         -0.0159

 -0.0000329         -0.0168

 -0.0000341         -0.0177

 -0.0000352         -0.0185

 -0.0000364         -0.0194

 -0.0000375         -0.0203

 -0.0000386         -0.0211

 -0.0000397         -0.0220

 -0.0000408         -0.0229

 -0.0000418         -0.0238

 -0.0000429         -0.0246

 -0.0000440         -0.0255

 -0.0000450         -0.0264

 -0.0000460         -0.0273

 -0.0000471         -0.0282

 -0.0000481         -0.0290

 -0.0000491         -0.0299

 -0.0000501         -0.0308

 -0.0000512         -0.0317

 -0.0000522         -0.0326

 -0.0000532         -0.0335

 -0.0000542         -0.0343

 -0.0000552         -0.0352

 -0.0000562         -0.0361

 -0.0000572         -0.0370

 -0.0000581         -0.0379

 -0.0000592         -0.0388

 -0.0000602         -0.0397

 -0.0000611         -0.0405

 -0.0000621         -0.0414

 -0.0000631         -0.0423

 -0.0000641         -0.0432

 -0.0000650         -0.0441

 -0.0000661         -0.0450

 -0.0000671         -0.0459

 -0.0000681         -0.0467

 -0.0000690         -0.0476

 -0.0000700         -0.0485

 -0.0000710         -0.0494

 -0.0000720         -0.0503

 -0.0000730         -0.0512

 -0.0000739         -0.0521

 -0.0000750         -0.0529

 -0.0000760         -0.0538

 -0.0000769         -0.0547

 -0.0000778         -0.0556

 -0.0000788         -0.0565

 -0.0000798         -0.0574

 -0.0000807         -0.0583

 -0.0000817         -0.0592

 -0.0000828         -0.0600

 -0.0000837         -0.0609

 -0.0000847         -0.0618

 -0.0000857         -0.0627

 -0.0000866         -0.0636

 -0.0000875         -0.0645

 -0.0000885         -0.0654

 -0.0000895         -0.0663

 -0.0000905         -0.0672

 -0.0000915         -0.0680

 -0.0000924         -0.0689

 -0.0000935         -0.0698

 -0.0000944         -0.0707

 -0.0000955         -0.0716

 -0.0000964         -0.0725

 -0.0000974         -0.0734

 -0.0000983         -0.0743

 -0.0000993         -0.0751

 -0.0001003         -0.0760

 -0.0001013         -0.0769

 -0.0001023         -0.0778

 -0.0001032         -0.0787

 -0.0001042         -0.0796

 -0.0001053         -0.0804

 -0.0001063         -0.0813

 -0.0001073         -0.0822

 -0.0001082         -0.0831

 -0.0001092         -0.0840

 -0.0001102         -0.0849

 -0.0001112         -0.0858

 -0.0001123         -0.0866

 -0.0001133         -0.0875

 -0.0001143         -0.0884

 -0.0001153         -0.0893

 -0.0001163         -0.0902

 -0.0001174         -0.0911

 -0.0001185         -0.0919

 -0.0001195         -0.0928

 -0.0001206         -0.0937

 -0.0001216         -0.0946

 -0.0001227         -0.0954

 -0.0001238         -0.0963

 -0.0001247         -0.0972

 -0.0001259         -0.0981

 -0.0001270         -0.0989

 -0.0001280         -0.0998

 -0.0001290         -0.1007

 -0.0001301         -0.1016

 -0.0001312         -0.1025

 -0.0001323         -0.1033

 -0.0001335         -0.1042

 -0.0001345         -0.1051

 -0.0001356         -0.1060

 -0.0001367         -0.1068

 -0.0001378         -0.1077

 -0.0001389         -0.1086

 -0.0001400         -0.1094

 -0.0001411         -0.1103

 -0.0001423         -0.1112

 -0.0001434         -0.1121

 -0.0001445         -0.1129

 -0.0001457         -0.1138

 -0.0001468         -0.1147

 -0.0001480         -0.1155

 -0.0001490         -0.1164

 -0.0001503         -0.1173

 -0.0001514         -0.1181

 -0.0001526         -0.1190

 -0.0001538         -0.1199

 -0.0001550         -0.1207

 -0.0001562         -0.1216

 -0.0001574         -0.1224

 -0.0001586         -0.1233

 -0.0001598         -0.1242

 -0.0001609         -0.1250

 -0.0001621         -0.1259

 -0.0001634         -0.1267

 -0.0001646         -0.1276

 -0.0001659         -0.1285

 -0.0001671         -0.1293

 -0.0001683         -0.1302

 -0.0001696         -0.1310

 -0.0001708         -0.1319

 -0.0001721         -0.1327

 -0.0001734         -0.1336

 -0.0001747         -0.1344

 -0.0001759         -0.1353

 -0.0001772         -0.1362

 -0.0001784         -0.1370

 -0.0001798         -0.1379

 -0.0001811         -0.1387

 -0.0001824         -0.1396

 -0.0001837         -0.1404

 -0.0001851         -0.1412

 -0.0001864         -0.1421

 -0.0001877         -0.1429

 -0.0001891         -0.1438

 -0.0001904         -0.1446

 -0.0001918         -0.1455

 -0.0001932         -0.1463

 -0.0001945         -0.1472

 -0.0001959         -0.1480

 -0.0001973         -0.1488

 -0.0001986         -0.1497

 -0.0002000         -0.1505

 -0.0002014         -0.1514

 -0.0002028         -0.1522

 -0.0002042         -0.1530

 -0.0002056         -0.1539

 -0.0002070         -0.1547

 -0.0002085         -0.1555

 -0.0002098         -0.1564

 -0.0002113         -0.1572

 -0.0002127         -0.1581

 -0.0002141         -0.1589

 -0.0002154         -0.1597

 -0.0002168         -0.1606

 -0.0002182         -0.1614

 -0.0002197         -0.1622

 -0.0002212         -0.1631

 -0.0002227         -0.1639

 -0.0002242         -0.1647

 -0.0002256         -0.1656

 -0.0002271         -0.1664

 -0.0002286         -0.1672

 -0.0002301         -0.1680

 -0.0002316         -0.1689

"

    },

    x=xy(all:0), //x取矩阵的列

    y=xy(all:1), //y取矩阵的列

    Opt[@f,optmode,100]     //Opt函数全局优化

};