| 1 | 1/1 | 返回列表 |
| 查看: 774 | 回復(fù): 0 | |||
| 【懸賞金幣】回答本帖問題,作者ymgy將贈送您 100 個金幣 | |||
ymgy銅蟲 (初入文壇)
|
[求助]
麻煩幫忙跑一下1stopt的程序,萬分感謝。
|
||
|
1stopt的程序,代碼已經(jīng)編好,可以直接運(yùn)行,只要復(fù)制到程序里即可。請大家?guī)兔ε芤幌?stopt的程序,麻煩把結(jié)果截圖以及擬合參數(shù)的數(shù)值發(fā)給我,謝謝大神們,萬飛感謝!。 程序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ā)表 | |
|---|---|---|---|---|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +3 | mysdl 2026-03-11 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0857環(huán)境調(diào)劑 +5 | 熠熠_11 2026-03-10 | 5/250 |
|
|
[考研] 材料專碩288分求調(diào)劑 一志愿211 +3 | 在家想你 2026-03-11 | 3/150 |
|
|
[考研] 0856材料與化工353分求調(diào)劑 +11 | NIFFFfff 2026-03-09 | 11/550 |
|
|
[考研] 求調(diào)劑材料專碩293 +6 | 段_(:з」∠)_ 2026-03-10 | 6/300 |
|
|
[考研] 收調(diào)劑 +7 | 調(diào)劑的考研學(xué)生 2026-03-10 | 7/350 |
|
|
[考研] 311求調(diào)劑 +3 | 牛乳糖的卡卡 2026-03-10 | 3/150 |
|
|
[考研] 環(huán)境調(diào)劑 +5 | 曉看天暮看云 2026-03-09 | 5/250 |
|
|
[基金申請] 面上項目還需要AI說明嗎? +3 | liyundong 2026-03-08 | 3/150 |
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
[考研]
|
程晴之 2026-03-06 | 6/300 |
|
|
[考研] 268求調(diào)劑 +4 | 劉合華 2026-03-05 | 4/200 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財 2026-03-05 | 5/250 |
|
|
[考研] 271求調(diào)劑 +7 | 月色c 2026-03-05 | 8/400 |
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼?~+123456 2026-03-05 | 5/250 |
|
|
[考研] 紡織、生物、化學(xué)、材料等專業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|