24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者ymgy将赠送您 200 个金币

ymgy

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 81.5
帖子: 7
在线: 9.9小时
虫号: 2022759
注册: 2012-09-23
专业: 可再生与替代能源利用中的

[求助] 多参数拟合，求助1stopt代算

1stopt程序，麻烦大神帮我跑一边程序吧，目前我没有正版软件，有点急用，非常感谢！！！

Constant t=20 ;
Parameter e=[1,10], a(t)=[0.01,100], p(t)=[0.01,100], c(t)=[0.01,100];
Variable x,y[OutPut];
StartProgram [Pascal];
Procedure MainModel;
var i, j: integer;
temD1, temD2, er, ei, k, n: double;
Begin
for j := 0 to DataLength - 1 do begin
      temD1 := 0;
      temD2 := 0;
      for i := 1 to t do begin
         temD1 := temD1 + (sqr(p)*(sqr(a)-sqr(x[j])))/(sqr(sqr(a)-sqr(x[j]))+sqr(c)*sqr(x[j]));
         temD2 := temD2 + (sqr(p)*c*x[j])/(sqr(sqr(a)-sqr(x[j]))+sqr(c)*sqr(x[j]));
      end;
      er := e + temD1;
      ei := temD2;
      k := ei/sqrt(2*(er+sqrt(er*er+ei*ei)));
      n := sqrt((er+sqrt(er*er+ei*ei))/2);
      y[j] :=(sqr(n-1)+k*k)/(sqr(n+1)+k*k);
end;
End;
EndProgram;
Data;
12.000       0.150
11.538       0.161
11.111       0.171
10.714       0.183
10.345       0.192
10.000       0.201
9.677       0.214
9.375       0.225
9.091       0.236
8.824       0.246
8.571       0.256
8.333       0.266
8.108       0.278
7.895       0.290
7.692       0.300
7.500       0.308
7.317       0.314
7.143       0.325
6.977       0.336
6.818       0.346
6.667       0.353
6.522       0.356
6.383       0.362
6.250       0.372
6.122       0.380
6.000       0.386
5.882       0.392
5.769       0.397
5.660       0.402
5.556       0.407
5.455       0.410
5.357       0.414
5.263       0.416
5.172       0.417
5.085       0.422
5.000       0.427
4.918       0.431
4.839       0.434
4.762       0.437
4.688       0.440
4.615       0.443
4.545       0.440
4.478       0.445
4.412       0.453
4.348       0.456
4.286       0.457
4.225       0.462
4.167       0.472
4.110       0.483
4.054       0.486
4.000       0.489
3.947       0.492
3.896       0.494
3.846       0.496
3.797       0.498
3.750       0.500
3.704       0.502
3.659       0.503
3.614       0.505
3.571       0.506
3.529       0.508
3.488       0.515
3.448       0.544
3.409       0.573
3.371       0.580
3.333       0.581
3.297       0.583
3.261       0.584
3.226       0.586
3.191       0.587
3.158       0.589
3.125       0.590
3.093       0.592
3.061       0.594
3.030       0.595
3.000       0.596
2.970       0.597
2.941       0.599
2.913       0.600
2.885       0.602
2.857       0.604
2.830       0.605
2.804       0.607
2.778       0.608
2.752       0.609
2.727       0.611
2.703       0.612
2.679       0.613
2.655       0.614
2.632       0.616
2.609       0.617
2.586       0.618
2.564       0.620
2.542       0.621
2.521       0.623
2.500       0.624
2.479       0.626
2.459       0.627
2.439       0.628
2.419       0.630
2.400       0.631
2.381       0.633
2.362       0.634
2.344       0.635
2.326       0.637
2.308       0.638
2.290       0.639
2.273       0.640
2.256       0.642
2.239       0.643
2.222       0.644
2.206       0.644
2.190       0.644
2.174       0.645
2.158       0.645
2.143       0.646
2.128       0.646
2.113       0.647
2.098       0.648
2.083       0.650
2.069       0.651
2.055       0.652
2.041       0.654
2.027       0.655
2.013       0.657
2.000       0.659
1.987       0.660
1.974       0.662
1.961       0.663
1.948       0.665
1.935       0.666
1.923       0.667
1.911       0.668
1.899       0.670
1.887       0.671
1.875       0.672
1.863       0.673
1.852       0.674
1.840       0.674
1.829       0.667
1.818       0.627
1.807       0.588
1.796       0.580
1.786       0.581
1.775       0.581
1.765       0.581
1.754       0.581
1.744       0.582
1.734       0.583
1.724       0.583
1.714       0.583
1.705       0.583
1.695       0.583
1.685       0.583
1.676       0.584
1.667       0.586
1.657       0.586
1.648       0.585
1.639       0.585
1.630       0.585
1.622       0.584
1.613       0.582
1.604       0.577
1.596       0.567
1.587       0.553
1.579       0.536
1.571       0.520
1.563       0.508
1.554       0.503
1.546       0.502
1.538       0.504
1.531       0.509
1.523       0.514
1.515       0.521
1.508       0.527
1.500       0.534
1.493       0.540
1.485       0.546
1.478       0.550
1.471       0.553
1.463       0.556
1.456       0.557
1.449       0.559
1.442       0.561
1.435       0.563
1.429       0.566
1.422       0.570
1.415       0.575
1.408       0.580
1.402       0.583
1.395       0.585
1.389       0.586
1.382       0.586
1.376       0.586
1.370       0.587
1.364       0.588
1.357       0.589
1.351       0.590
1.345       0.591
1.339       0.593
1.333       0.592
1.327       0.590
1.322       0.588
1.316       0.585
1.310       0.579
1.304       0.573
1.299       0.570
1.293       0.570
1.288       0.572
1.282       0.571
1.277       0.569
1.271       0.569
1.266       0.567
1.261       0.569
1.255       0.567
1.250       0.565
1.245       0.563
1.240       0.562
1.235       0.558
1.230       0.556
1.224       0.555
1.220       0.547
1.215       0.543
1.210       0.535
1.205       0.531
1.200       0.523

回复此楼

» 猜你喜欢

求助中医院临床中医中标的青年标书学习，谢谢已经有5人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
物理学I论文润色/翻译怎么收费? 已经有244人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有1人回复
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有19人回复
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计已经有1人回复

1楼 2020-08-27 02:05:47

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

相关版块跳转我要订阅楼主 ymgy 的主题更新

返回列表

不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)

普通表情龙兔虎猫