| 5 | 1/1 | 返回列表 |
| 查看: 725 | 回復(fù): 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è)作家)
![]() |
專家經(jīng)驗: +1059 |
|
MATLAB程序,復(fù)制進一個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 |

鐵桿木蟲 (著名寫手)
工程和科學(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 當然,你要給迭代初值。你這里涉及到sqrt函數(shù),所以給初值要小心。應(yīng)該有經(jīng)驗。 上面是一般方法。實現(xiàn)需要編程。主要的工作是求解3X3的逆矩陣。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +3 | 阿姨復(fù)古不過 2026-03-27 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0856材料化工調(diào)劑 總分330 +10 | zhubinhao 2026-03-27 | 10/500 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 機械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|