24小时热门版块排行榜    

查看: 595  |  回复: 0
【悬赏金币】回答本帖问题,作者ymgy将赠送您 200 个金币

ymgy

铜虫 (初入文坛)

[求助] 多参数拟合,求助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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ymgy 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +3 ahbd 2026-03-14 3/150 2026-03-16 13:05 by Iveryant
[考研] 304求调剂 +4 素年祭语 2026-03-15 4/200 2026-03-16 09:42 by 闲人终南山
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 求调剂! +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 328,0703考生求调剂,一志愿为东北师范大学 +4 观素律 2026-03-09 5/250 2026-03-14 01:24 by JourneyLucky
[考研] 一志愿郑大070303,338分,求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
信息提示
请填处理意见