 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 732 | 回復(fù): 4 | ||
夜的寧靜金蟲 (小有名氣)
|
[求助]
matlab牛頓迭代的問題，出現(xiàn)Matrix dimensions must agre的錯誤 已有1人參與
|
|
在用用牛頓迭代求2*log10(5.4*10^-4+2.51*10^-5/sqrt(x)) +1/sqrt(x)=0中的x 初值是0.03 function [x,fval]=newton() clc; syms x; f=@(x)2*log10(5.4*10^-4+2.51*10^-5/sqrt(x))+1/sqrt(x); df=@(x)diff(2*log10(5.4*10^-4+2.51*10^-5/sqrt(x))+1/sqrt(x)); t=1e-5; x0=0.03; x=x0+2*t; k=0; while abs(x-x0)>t k=k+1; x=x0; x=x0-feval(f,x0)/feval(df,x0);%這句話有錯，提示Matrix dimensions must agree. end if k==500 disp('迭代次數(shù)過多，防止死循環(huán)終止'); else fprintf('迭代到%d次時得到結(jié)果%f\n',k,x) end 360反饋意見截圖16570207506364.png |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有289人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求解帶積分的非線性方程組，牛頓迭代過程中的奇異矩陣如何修正
已經(jīng)有0人回復(fù)
- matlab實現(xiàn)牛頓迭代法求非線性方程組程序運行出現(xiàn)問題
已經(jīng)有8人回復(fù)
- 求助用matlab解多元指數(shù)型方程組
已經(jīng)有10人回復(fù)
- matlab 牛頓迭代法
已經(jīng)有3人回復(fù)
- matlab對含貝塞爾函數(shù)的超越方程求根
已經(jīng)有5人回復(fù)
- 牛頓迭代法接非線性方程組求解
已經(jīng)有1人回復(fù)
- 求助如何用matlab求非線性方程組的解
已經(jīng)有8人回復(fù)
- 自己的寫的matlab跑不動，哪位大神幫我用牛頓迭代編一下
已經(jīng)有4人回復(fù)
- 求用牛頓迭代法解非線性方程組的matlab程序
已經(jīng)有4人回復(fù)
- 牛頓迭代法解方程。
已經(jīng)有11人回復(fù)
- Matlab用牛頓法求解非線性方程組問題
已經(jīng)有8人回復(fù)
- MATLAB牛頓-SOR迭代法解n維方程組
已經(jīng)有0人回復(fù)
- 求：牛頓迭代法matlab程序及輸出結(jié)果.求指點
已經(jīng)有1人回復(fù)
- matlab 求解非線性方程組 牛頓解法
已經(jīng)有1人回復(fù)
- 在matlab中用牛頓迭代法
已經(jīng)有1人回復(fù)
- 使用matlab進行Newton-Raphson迭代是出現(xiàn)虛數(shù)的問題
已經(jīng)有0人回復(fù)
- [請教]用牛頓迭代整這個問題
已經(jīng)有7人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
材料廖
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 4889.5
- 散金: 200
- 紅花: 11
- 帖子: 891
- 在線: 295.1小時
- 蟲號: 1794189
- 注冊: 2012-05-04
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
2樓2015-05-18 11:45:01
夜的寧靜
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 604.9
- 散金: 90
- 帖子: 93
- 在線: 74.4小時
- 蟲號: 2307122
- 注冊: 2013-02-28
- 性別: GG
- 專業(yè): 機械動力學(xué)
|
還是不知道怎么改，我是初學(xué)，因為用到現(xiàn)學(xué)的，所以基礎(chǔ)不怎么樣。x0沒有更新什么意思 之前的賦值可能有些錯誤，改成這樣還是不行，麻煩再看一下，十分感謝function [x,fval]=newton() clc; syms x; f=@(x)2*log10(5.4*10^-4+2.51*10^-5/sqrt(x)) +1/sqrt(x); df=@(x)diff(2*log10(5.4*10^-4+2.51*10^-5/sqrt(x)) +1/sqrt(x)); t=1e-5; x=0.03; x0=x+2*t; k=0; while abs(x-x0)>t k=k+1; x0=x; x=x0-feval(f,x0)/feval(df,x0);%這句話有錯，但不知道錯在哪里 end if k==500 disp('迭代次數(shù)過多，防止死循環(huán)終止'); else fprintf('迭代到%d次時得到結(jié)果%f\n',k,x) end |
3樓2015-05-18 22:23:32
材料廖
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 4889.5
- 散金: 200
- 紅花: 11
- 帖子: 891
- 在線: 295.1小時
- 蟲號: 1794189
- 注冊: 2012-05-04
- 性別: GG
- 專業(yè): 金屬材料的合金相、相變及
|
試試這個看對不對 function x = newton() clc clear syms x f df f = 2*log10(5.4*10^-4+2.51*10^-5/sqrt(x)) +1/sqrt(x); df = diff(2*log10(5.4*10^-4+2.51*10^-5/sqrt(x)) +1/sqrt(x)); t=1e-5; x=0.03; x0=x+2*t; k=0; while abs(x-x0)>t k=k+1; x0=x; x=x0-double(subs(f,x0)/subs(df,x0)); % feval(f,x0)./feval(df,x0);%這句話有錯，但不知道錯在哪里 end if k==500 disp('迭代次數(shù)過多，防止死循環(huán)終止'); else fprintf('迭代到%d次時得到結(jié)果%f\n',k,x) end |
» 本帖已獲得的紅花（最新10朵）
夜的寧靜4樓2015-05-19 15:56:37
夜的寧靜
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 604.9
- 散金: 90
- 帖子: 93
- 在線: 74.4小時
- 蟲號: 2307122
- 注冊: 2013-02-28
- 性別: GG
- 專業(yè): 機械動力學(xué)
很完美，太感謝啦 5樓2015-05-20 08:37:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 夜的寧靜 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
