24小时热门版块排行榜    

查看: 645  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见