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