| 查看: 563 | 回复: 0 | |||
| 【悬赏金币】回答本帖问题，作者ymgy将赠送您 200 个金币 | |||
[求助]
多参数拟合，求助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 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有118人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复