| 5 | 1/1 | 返回列表 |
| 查看: 722 | 回復: 5 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
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 ] |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經驗: +1059 |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經驗: +1059 |
|
MATLAB程序,復制進一個m文件,F(xiàn)5運行即可,根據(jù)自己的需要輸入已知參數(shù),即返回相應的解: 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))); |

鐵桿木蟲 (著名寫手)
工程和科學數(shù)值計算咨詢
|
你這是非線性方程組,最好用Newton法求解。設 三個函數(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 當然,你要給迭代初值。你這里涉及到sqrt函數(shù),所以給初值要小心。應該有經驗。 上面是一般方法。實現(xiàn)需要編程。主要的工作是求解3X3的逆矩陣。 |
鐵桿木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 北科281學碩材料求調劑 +16 | tcxiaoxx 2026-03-20 | 18/900 |
|
|---|---|---|---|---|
|
[考研]
材料調劑
5+3
|
想要一壺桃花水 2026-03-25 | 8/400 |
|
|
[考研] 297求調劑 +3 | 田洪有 2026-03-26 | 3/150 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 調劑310 +3 | 溫柔的晚安 2026-03-25 | 4/200 |
|
|
[考研] 考研調劑 +6 | 來好運來來來 2026-03-21 | 7/350 |
|
|
[考研] 290分調劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 材料調劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 340求調劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 生物學學碩求調劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 材料專碩找調劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 一志愿河北工業(yè)大學0817化工278分求調劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 一志愿山東大學藥學學碩求調劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 求調劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
[考研] 環(huán)境學碩288求調劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
[考研] 328求調劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 293求調劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 260求調劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 一志愿北京化工大學070300 學碩336求調劑 +5 | vv迷 2026-03-21 | 8/400 |
|