24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者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@月只蓝

回复此楼

» 猜你喜欢

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

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

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

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

返回列表

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

普通表情龙兔虎猫

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3	1uaZt2kAZj2S 2026-10-30	14/700	2026-11-01 07:01 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3	X5bHxOULxFqb 2026-10-30	17/850	2026-11-01 06:48 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4	X5bHxOULxFqb 2026-10-30	21/1050	2026-11-01 06:46 by YXM1FDCSE8Bn
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	xvJS1BYDxvBS 2026-10-31	8/400	2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	xvJS1BYDxvBS 2026-10-31	8/400	2026-11-01 05:43 by YXM1FDCSE8Bn
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	1uaZt2kAZj2S 2026-10-31	6/300	2026-11-01 05:41 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3	1uaZt2kAZj2S 2026-10-31	6/300	2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	1uaZt2kAZj2S 2026-10-30	15/750	2026-11-01 05:22 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	1uaZt2kAZj2S 2026-10-31	7/350	2026-11-01 04:04 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3	1uaZt2kAZj2S 2026-10-31	3/150	2026-11-01 04:04 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3	1uaZt2kAZj2S 2026-10-31	8/400	2026-11-01 03:59 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	1uaZt2kAZj2S 2026-10-31	6/300	2026-11-01 03:56 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	1uaZt2kAZj2S 2026-10-31	5/250	2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	1uaZt2kAZj2S 2026-10-30	12/600	2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3	1uaZt2kAZj2S 2026-10-31	4/200	2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4	X5bHxOULxFqb 2026-10-30	12/600	2026-11-01 00:15 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	X5bHxOULxFqb 2026-10-30	13/650	2026-10-31 17:12 by hKe7U1kad5Sg
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3	yyh@ 2026-10-28	4/200	2026-10-31 12:47 by 可淡不可忘
[基金申请] 论文的挂名怎么评判贡献？基金申请中，排名第一的挂名文章 +3	biily 2026-10-25	3/150	2026-10-30 18:33 by 旁观2020
[考博] 申博求助 20+3	Katrinae 2026-10-26	4/200	2026-10-27 09:56 by 北京莱茵编辑