| 24小時(shí)熱門版塊排行榜 |
| 6 | 1/1 | 返回列表 |
| 查看: 2165 | 回復(fù): 5 | ||
[求助]
非線性方程組的解法
|
|
如果用fortran解下面兩個(gè)方程組: 1. 5+3.2ln(x)+ln(1-x)-2x+5x^2+6x^3-3y=0; 2. 0.3x+0.7y=0.9 方程中既有l(wèi)n,又有3次冪。 應(yīng)該怎么編程呢?我是用別人的程序,現(xiàn)在需要加入這個(gè)。請(qǐng)高手指點(diǎn),多謝! |
銀蟲(chóng) (小有名氣)
|
將上述兩個(gè)方程化為: y=(5+3.2ln(x)+ln(1-x)-2x+5x^2+6x^3)/3; x=(0.9-0.7y)/0.3 然后將(x0,y0)賦初值采用上述式子迭代計(jì)算,得到的本次結(jié)果(xn,yn)作為下一次結(jié)果的初值。循環(huán)迭代至前后兩次的結(jié)果滿足精度要求,可以取誤差的二范數(shù)或者無(wú)窮范數(shù)等作為判斷。 該方法稱為迭代法,具體的fortran算法,可以參考gauss-賽德?tīng)柕惴。網(wǎng)上有很多這樣的程序,也可以參見(jiàn)徐士良的《fortran常用算法程序集》這本書(shū)。 以上僅介紹計(jì)算的主要過(guò)程,愿對(duì)你有幫助。 不足之處,請(qǐng)見(jiàn)諒…… |

|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
至尊木蟲(chóng) (正式寫(xiě)手)
|
解題思路:先對(duì)上述方程組第2式進(jìn)行變形,把1式中y用x代替 最終得:x^22.4=exp(5x-35x^2-42x^3-8)/((1-x))^7 用迭代法VB程序如下: Dim x, x1 As Single x = 0.1 For i = 1 To 1000 x1 = ((Exp(5 * x - 35 * x ^ 2 - 42 * x ^ 3 - 8)) / ((1 - x) ^ 7)) ^ (1 / 22.4) x = x1 Next i Print x 解得:x=0.3031055 則 y=1.155811929(x帶入方程組2式手算得) |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
.6lL 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +6 | d如愿上岸 2026-03-12 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 妮妮ninicgb 2026-03-15 | 11/550 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[論文投稿]
投稿問(wèn)題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|