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

回复此楼

» 猜你喜欢

津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名已经有0人回复
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有222人回复
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）已经有6人回复
【修正版】单环用磁约束低速电子实验方案（简化版）已经有0人回复
桂林理工大学物理学专业招收调剂，还有三个名额！！！已经有22人回复
考博自荐已经有4人回复
山东大学第二批博士研究生招生已经有0人回复
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复

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

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

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

返回列表

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

普通表情龙兔虎猫