 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 724 | 回復(fù): 5 | |||||
ldoop
|
[求助]
高手們，請(qǐng)幫忙解一個(gè)三元方程組啊，謝謝了
|
|
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 ] |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有198人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求解一元三次方程
已經(jīng)有11人回復(fù)
- 用matlab解帶參數(shù)的一元三次方程--求解
已經(jīng)有3人回復(fù)
- 求Microsoft Math3.0高手，幫畫一張一元三次方程圖，最后再往郵箱發(fā)一份。
已經(jīng)有10人回復(fù)
- 求Microsoft Math3.0高手，幫畫一張一元三次方程圖
已經(jīng)有9人回復(fù)
- 求高手幫我編一個(gè)一元三次方程求解的小程序！重謝哦！
已經(jīng)有5人回復(fù)
- 【求助】關(guān)于MATLAB解三次方程
已經(jīng)有7人回復(fù)
- 【求助】求助：vb編程中用牛頓迭代解三次方程為什么只得到一個(gè)根？【已完成】
已經(jīng)有15人回復(fù)
- 【求助】如何用Fortran解一元三次方程【已解決】
已經(jīng)有6人回復(fù)
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 鼓勵(lì)交流 2013-05-18 14:43:15
ldoop: 金幣+10, ★有幫助, 謝謝回復(fù)，我是想求得解析解而非數(shù)值解，不過(guò)還是感謝了，麻煩幫忙研究下是否能求得解析解 2013-05-20 16:26:21
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 鼓勵(lì)交流 2013-05-18 14:43:15
ldoop: 金幣+10, ★有幫助, 謝謝回復(fù)，我是想求得解析解而非數(shù)值解，不過(guò)還是感謝了，麻煩幫忙研究下是否能求得解析解 2013-05-20 16:26:21
|
MATLAB程序，復(fù)制進(jìn)一個(gè)m文件，F(xiàn)5運(yùn)行即可，根據(jù)自己的需要輸入已知參數(shù)，即返回相應(yīng)的解： function solvenole314 clear all;clc format long global a_0 b_0 x_0 a_0=input('請(qǐng)輸入已知常數(shù)a_0='); b_0=input('請(qǐng)輸入已知常數(shù)b_0='); x_0=input('請(qǐng)輸入已知常數(shù)x_0='); aint=input('請(qǐng)輸入未知數(shù)a的猜想值='); bint=input('請(qǐng)輸入未知數(shù)b的猜想值='); xint=input('請(qǐng)輸入未知數(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))); |
2樓2013-05-18 09:43:50
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
★ ★ ★
ldoop(fegg7502代發(fā)): 金幣+3, 鼓勵(lì)交流 2013-05-18 14:46:48
ldoop(fegg7502代發(fā)): 金幣+3, 鼓勵(lì)交流 2013-05-18 14:46:48
| 一個(gè)計(jì)算的例子如附件。 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 附圖1.jpg
2013-05-18 09:47:31, 52.57 K
3樓2013-05-18 09:47:58
pippi6
鐵桿木蟲(chóng) (著名寫手)
工程和科學(xué)數(shù)值計(jì)算咨詢
- 應(yīng)助: 413 (碩士)
- 貴賓: 0.002
- 金幣: 7116.5
- 散金: 15
- 紅花: 63
- 帖子: 1639
- 在線: 798.9小時(shí)
- 蟲(chóng)號(hào): 2469437
- 注冊(cè): 2013-05-14
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
ben_ladeng: 金幣+2, 歡迎常來(lái) 2013-05-19 23:09:22
ldoop: 金幣+10, ★有幫助, 謝謝回復(fù)，我是想求得解析解而非數(shù)值解，不過(guò)還是感謝了，麻煩幫忙研究下是否能求得解析解 2013-05-20 17:02:46
感謝參與，應(yīng)助指數(shù) +1
ben_ladeng: 金幣+2, 歡迎常來(lái) 2013-05-19 23:09:22
ldoop: 金幣+10, ★有幫助, 謝謝回復(fù)，我是想求得解析解而非數(shù)值解，不過(guò)還是感謝了，麻煩幫忙研究下是否能求得解析解 2013-05-20 17:02:46
|
你這是非線性方程組，最好用Newton法求解。設(shè) 三個(gè)函數(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)驗(yàn)。 上面是一般方法。實(shí)現(xiàn)需要編程。主要的工作是求解3X3的逆矩陣。 |
4樓2013-05-19 07:30:01
ldoop 
鐵桿木蟲(chóng) (著名寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 3485.8
- 散金: 3151
- 紅花: 16
- 沙發(fā): 1
- 帖子: 2225
- 在線: 285.5小時(shí)
- 蟲(chóng)號(hào): 533053
- 注冊(cè): 2008-03-25
- 專業(yè): 群落生態(tài)學(xué)
5樓2013-05-20 17:05:40
pippi6
鐵桿木蟲(chóng) (著名寫手)
工程和科學(xué)數(shù)值計(jì)算咨詢
- 應(yīng)助: 413 (碩士)
- 貴賓: 0.002
- 金幣: 7116.5
- 散金: 15
- 紅花: 63
- 帖子: 1639
- 在線: 798.9小時(shí)
- 蟲(chóng)號(hào): 2469437
- 注冊(cè): 2013-05-14
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
【答案】應(yīng)助回帖
★
csgt0: 金幣+1, 歡迎交流 2013-05-21 15:36:49
csgt0: 金幣+1, 歡迎交流 2013-05-21 15:36:49
|
不要迷戀解析解。實(shí)際中遇到的絕大多數(shù)方程都沒(méi)有解析解。從實(shí)用觀點(diǎn)出發(fā)，只要能解就好，沒(méi)有解析解也沒(méi)有關(guān)系。解析解如果太復(fù)雜，也沒(méi)有什么意義。20多年前，我的一個(gè)同事給我一個(gè)解析解，用maple做的，解的表達(dá)式寫了7頁(yè)紙。從計(jì)算速度、簡(jiǎn)明性上和精度上說(shuō)，都不如數(shù)值解。另外，三元線性方程是有公式解，但也是很復(fù)雜。還要判斷根的情況。我看都不如newton法。一招制敵。會(huì)了這個(gè)，基本不需要其他。一隅之見(jiàn)，希望有幫助。 |
6樓2013-05-21 07:03:08
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ldoop 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +6 | 孅華 2026-03-22 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 295求調(diào)劑 +4 | 1428151015 2026-03-27 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +8 | PolarBear11 2026-03-26 | 8/400 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 08開(kāi)頭275求調(diào)劑 +4 | 拉誰(shuí)不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +4 | 李軼男003 2026-03-20 | 4/200 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|