| 5 | 1/1 | 返回列表 |
| 查看: 723 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
ldoop
|
[求助]
高手們,請幫忙解一個三元方程組啊,謝謝了
|
|||
|
b*x -a -b-sqrt((b*x-a-b)^2-4*ab)=b_0*x_0 -a_0 -b_0-sqrt((b_0*x_0-a_0-b_0)^2-4*a_0*b_0 b*x -a -b+sqrt((b*x-a-b)^2-4ab)=b0*x0 -a0 -b0+sqrt((b0*x0-a0-b_0)^2-4*a_0*b_0 x*b/(a+b+2*sqrt(a*b))=b_0*x_0/(a_0+b_0+2*sqrt(a_0*b_0)) 其中a_0,b_0,x_0為已知常數(shù),求未知數(shù) a b x [ Last edited by ldoop on 2013-5-17 at 23:42 ] |
鐵桿木蟲 (著名寫手)
工程和科學(xué)數(shù)值計算咨詢
|
你這是非線性方程組,最好用Newton法求解。設(shè) 三個函數(shù) f1=b*x -a -b-sqrt((b*x-a-b)^2-4*ab)-b_0*x_0 +a_0 +b_0+sqrt((b_0*x_0-a_0-b_0)^2+4*a_0*b_0 f2=b*x -a -b+sqrt((b*x-a-b)^2-4ab)-b0*x0 +a0 +b0-sqrt((b0*x0-a0-b_0)^2+4*a_0*b_3=0 f3=x*b/(a+b+2*sqrt(a*b))-b_0*x_0/(a_0+b_0+2*sqrt(a_0*b_0)) 記3X3矩陣 A=∂(f1,f2,f3)/∂(x,a,b) let dX=(dx,da,db), F=(f1,f2,f3) 解出 dX = -inverse(A) F 然后迭代 x=x+dx a=a+da b=b+db 當(dāng)然,你要給迭代初值。你這里涉及到sqrt函數(shù),所以給初值要小心。應(yīng)該有經(jīng)驗。 上面是一般方法。實現(xiàn)需要編程。主要的工作是求解3X3的逆矩陣。 |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |
|
MATLAB程序,復(fù)制進(jìn)一個m文件,F(xiàn)5運行即可,根據(jù)自己的需要輸入已知參數(shù),即返回相應(yīng)的解: function solvenole314 clear all;clc format long global a_0 b_0 x_0 a_0=input('請輸入已知常數(shù)a_0='); b_0=input('請輸入已知常數(shù)b_0='); x_0=input('請輸入已知常數(shù)x_0='); aint=input('請輸入未知數(shù)a的猜想值='); bint=input('請輸入未知數(shù)b的猜想值='); xint=input('請輸入未知數(shù)x的猜想值='); options=optimset('TolFun',1e-9,'MaxIter',1000); x1=fsolve(@fun299792,[aint bint xint],options) %-------------------------------- function y=fun299792(z) global a_0 b_0 x_0 a=z(1); b=z(2); x=z(3); y(1)=b*x-a-b-sqrt((b*x-a-b)^2-4*a*b)-(b_0*x_0-a_0-b_0-sqrt((b_0*x_0-a_0-b_0)^2-4*a_0*b_0)); y(2)=b*x-a-b+sqrt((b*x-a-b)^2-4*a*b)-(b_0*x_0-a_0-b_0+sqrt((b_0*x_0-a_0-b_0)^2-4*a_0*b_0)); y(3)=x*b/(a+b+2*sqrt(a*b))-(b_0*x_0/(a_0+b_0+2*sqrt(a_0*b_0))); |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |

鐵桿木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我!! +3 | 祁祺祺 2026-03-25 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 329求調(diào)劑 +4 | 星野? 2026-03-26 | 4/200 |
|
|
[考研] 321求調(diào)劑 +4 | 璞玉~~ 2026-03-25 | 4/200 |
|
|
[考研] 277跪求調(diào)劑 +3 | 1915668 2026-03-27 | 4/200 |
|
|
[考研] 292求調(diào)劑 +12 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 13/650 |
|
|
[考研] 269專碩求調(diào)劑 +9 | 金恩貝 2026-03-21 | 9/450 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|