24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

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

ymgy

铜虫 (初入文坛)

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

[求助] 麻烦帮忙跑一下1stopt的程序，万分感谢！！

1stopt的程序，代码已经编好，可以直接运行，只要复制到程序里即可。请大家帮忙跑一下1stopt的程序，麻烦把结果截图以及拟合参数的数值发给我，谢谢大神们，万飞感谢！。

程序1；
//*************************************
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;
1.2 0.9521
1.3 0.9318
1.4 0.8852
1.5 0.8678
1.6 0.8794
1.7 0.8972
1.8 0.9069
1.9 0.9116
2 0.9148
2.2 0.92
2.4 0.9228
2.6 0.9238
2.8 0.9242
3 0.9241
3.2 0.9243
3.4 0.9245
3.6 0.9246
3.8 0.9247
4 0.9248
4.2 0.9248
4.4 0.9249
4.6 0.9249
4.8 0.9249
4 0.9248
4.2 0.9248
4.4 0.9249
4.6 0.9249
4.8 0.9249
8 0.9269
8.5 0.9272
9 0.9277
9.5 0.9282
10 0.9286
10.5 0.9293
11 0.9298
11.5 0.9283
12 0.9224
12.5 0.9118

程序2；
//*************************************
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(c)+sqr(x[j]));
         temD2 := temD2 +(sqr(p)*c)/(sqr(x[j])*(x[j])+sqr(c)*(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;
1.2 0.9521
1.3 0.9318
1.4 0.8852
1.5 0.8678
1.6 0.8794
1.7 0.8972
1.8 0.9069
1.9 0.9116
2 0.9148
2.2 0.92
2.4 0.9228
2.6 0.9238
2.8 0.9242
3 0.9241
3.2 0.9243
3.4 0.9245
3.6 0.9246
3.8 0.9247
4 0.9248
4.2 0.9248
4.4 0.9249
4.6 0.9249
4.8 0.9249
4 0.9248
4.2 0.9248
4.4 0.9249
4.6 0.9249
4.8 0.9249
8 0.9269
8.5 0.9272
9 0.9277
9.5 0.9282
10 0.9286
10.5 0.9293
11 0.9298
11.5 0.9283
12 0.9224
12.5 0.9118

程序3；

//*************************************
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

程序4；
//*************************************
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(c)+sqr(x[j]));
         temD2 := temD2 +(sqr(p)*c)/(sqr(x[j])*(x[j])+sqr(c)*(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.52@月只蓝

回复此楼

» 猜你喜欢

🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队已经有0人回复
第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有180人回复
求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有13人回复

1楼2020-08-18 15:20:15

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

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

返回列表