| 查看: 569 | 回复: 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 |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有54人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有3人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有5人回复