| 1 | 1/1 | 返回列表 |
| 查看: 770 | 回復(fù): 0 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者ymgy將贈(zèng)送您 100 個(gè)金幣 | ||
ymgy銅蟲(chóng) (初入文壇)
|
[求助]
麻煩幫忙跑一下1stopt的程序,萬(wàn)分感謝!
|
|
|
1stopt的程序,代碼已經(jīng)編好,可以直接運(yùn)行,只要復(fù)制到程序里即可。請(qǐng)大家?guī)兔ε芤幌?stopt的程序,麻煩把結(jié)果截圖以及擬合參數(shù)的數(shù)值發(fā)給我,謝謝大神們,萬(wàn)飛感謝!。 程序1; //************************************* 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; 1.2 0.9521 1.3 0.9318 1.4 0.8852 1.5 0.8678 1.6 0.8794 1.7 0.8972 1.8 0.9069 1.9 0.9116 2 0.9148 2.2 0.92 2.4 0.9228 2.6 0.9238 2.8 0.9242 3 0.9241 3.2 0.9243 3.4 0.9245 3.6 0.9246 3.8 0.9247 4 0.9248 4.2 0.9248 4.4 0.9249 4.6 0.9249 4.8 0.9249 4 0.9248 4.2 0.9248 4.4 0.9249 4.6 0.9249 4.8 0.9249 8 0.9269 8.5 0.9272 9 0.9277 9.5 0.9282 10 0.9286 10.5 0.9293 11 0.9298 11.5 0.9283 12 0.9224 12.5 0.9118 程序2; //************************************* 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(c)+sqr(x[j])); temD2 := temD2 +(sqr(p)*c)/(sqr(x[j])*(x[j])+sqr(c)*(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; 1.2 0.9521 1.3 0.9318 1.4 0.8852 1.5 0.8678 1.6 0.8794 1.7 0.8972 1.8 0.9069 1.9 0.9116 2 0.9148 2.2 0.92 2.4 0.9228 2.6 0.9238 2.8 0.9242 3 0.9241 3.2 0.9243 3.4 0.9245 3.6 0.9246 3.8 0.9247 4 0.9248 4.2 0.9248 4.4 0.9249 4.6 0.9249 4.8 0.9249 4 0.9248 4.2 0.9248 4.4 0.9249 4.6 0.9249 4.8 0.9249 8 0.9269 8.5 0.9272 9 0.9277 9.5 0.9282 10 0.9286 10.5 0.9293 11 0.9298 11.5 0.9283 12 0.9224 12.5 0.9118 程序3; //************************************* 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 程序4; //************************************* 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(c)+sqr(x[j])); temD2 := temD2 +(sqr(p)*c)/(sqr(x[j])*(x[j])+sqr(c)*(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.52@月只藍(lán) |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工 298 調(diào)劑 +9 | 小西笑嘻嘻 2026-03-03 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +3 | 是樂(lè)渝哇 2026-03-04 | 3/150 |
|
|
[考研] 成績(jī)276,專(zhuān)業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
[考研] 0855機(jī)械工程324求調(diào)劑 +3 | 鄭SK 2026-03-02 | 7/350 |
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
[考研] 301求調(diào)劑 +3 | 李LJR 2026-03-04 | 3/150 |
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專(zhuān)業(yè)的同學(xué)報(bào)考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | 杰尼龜aaa 2026-03-03 | 4/200 |
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
[基金申請(qǐng)] 請(qǐng)問(wèn)大家,研究風(fēng)險(xiǎn)與應(yīng)對(duì)措施那里, 大家都怎么寫(xiě)呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
[考研] 281求調(diào)劑 +5 | 2026計(jì)算機(jī)_誠(chéng)心 2026-03-01 | 8/400 |
|
|
[考研] 材料類(lèi)求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
[考研] 299求調(diào)劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|
|
[考研]
|
LYidhsjabdj 2026-02-28 | 4/200 |
|