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@月只蓝

回复此楼

» 猜你喜欢

澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有256人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有120人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有71人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有5人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复

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

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

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

返回列表

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

普通表情龙兔虎猫

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

[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +7	IDs3scOF0tjC 2026-09-28	8/400	2026-10-01 19:08 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4	CfXuS1rDhLYN 2026-09-28	4/200	2026-10-01 09:23 by Equinoxhua
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	cqwQDCxMcL3I 2026-09-29	4/200	2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3	cqwQDCxMcL3I 2026-09-29	4/200	2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +10	zju2000 2026-09-26	13/650	2026-09-30 10:26 by yexuqing
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4	IDs3scOF0tjC 2026-09-28	4/200	2026-09-29 19:27 by cqwQDCxMcL3I
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3	ZYXdhzDAy9ZX 2026-09-28	4/200	2026-09-29 16:37 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3	tqUTOxClQUMF 2026-09-28	3/150	2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4	tqUTOxClQUMF 2026-09-28	4/200	2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	tqUTOxClQUMF 2026-09-28	3/150	2026-09-29 15:04 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6	GDBe8tDZqE8z 2026-09-28	6/300	2026-09-29 14:41 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6	IDs3scOF0tjC 2026-09-28	6/300	2026-09-29 14:13 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4	IDs3scOF0tjC 2026-09-28	5/250	2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4	GDBe8tDZqE8z 2026-09-28	4/200	2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3	ZYXdhzDAy9ZX 2026-09-28	3/150	2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3	ZYXdhzDAy9ZX 2026-09-28	3/150	2026-09-29 06:28 by ZvyPK8n6Nfki
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3	GDBe8tDZqE8z 2026-09-28	3/150	2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3	CfXuS1rDhLYN 2026-09-28	4/200	2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3	IDs3scOF0tjC 2026-09-28	6/300	2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	CfXuS1rDhLYN 2026-09-28	3/150	2026-09-28 19:08 by 9lS3ad5oOymn