| 1 | 1/1 | 返回列表 |
| 查看: 580 | 回復: 0 | ||
[求助]
使用matlab進行Newton-Raphson迭代是出現虛數的問題
|
|
首先說明,我是利用牛頓迭代法求解非線性方程組。所涉及的關鍵代碼如下: sta_fir=Fa+F0-sum(sta_Kn.*deltan.^(1.5).*sina); sta_sec=Fr-sum(sta_Kn.*deltan.^(1.5).*(cos(phai)).*cosa); sta_thir=M-dm/2*sum(sta_Kn.*deltan.^(1.5).*cos(phai).*sina); static=[sta_fir;sta_sec;sta_thir]; 。。。 findsym(static); findsym(sina); [gen,n]=Newtonraphson(static,x0); 但是在實際運行的時候因為方程組中1.5次方的緣故會導致迭代出現負數,產生一系列后果。 我在調試中嘗試過將方程在牛頓迭代函數中重新定義并添加判斷來消除負數,但是這樣又會導致原方程組Jacobi矩陣改變,進而出現不收斂的情況。 因此特來此求助,希望有高手能告知妥善的解決方法,不勝感激。 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求調劑 +6 | chuanzhu川燭 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 求調劑一志愿海大,0703化學學碩304分,有大創(chuàng)項目,四級已過 +5 | 幸運哩哩 2026-03-22 | 9/450 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 考研調劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
[考研] 工科0856求調劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 278求調劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 306求0703調劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境總分308求調劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 310求調劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 085700資源與環(huán)境308求調劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 329求調劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
[考研] 材料學碩297已過四六級求調劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 一志愿西安交通大學 學碩 354求調劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 261求B區(qū)調劑,科研經歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 工科材料085601 279求調劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 0854可跨調劑,一作一項核心論文五項專利,省、國級證書40+數一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調劑 +4 | Jy? 2026-03-16 | 4/200 |
|