| 查看: 908 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
夜的宁静金虫 (小有名气)
|
[求助]
matlab牛顿迭代的问题，出现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('迭代次数过多，防止死循环终止'); else fprintf('迭代到%d次时得到结果%f\n',k,x) end  360反馈意见截图16570207506364.png |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有66人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解带积分的非线性方程组，牛顿迭代过程中的奇异矩阵如何修正
已经有0人回复
- matlab实现牛顿迭代法求非线性方程组程序运行出现问题
已经有8人回复
- 求助用matlab解多元指数型方程组
已经有10人回复
- matlab 牛顿迭代法
已经有3人回复
- matlab对含贝塞尔函数的超越方程求根
已经有5人回复
- 牛顿迭代法接非线性方程组求解
已经有1人回复
- 求助如何用matlab求非线性方程组的解
已经有8人回复
- 自己的写的matlab跑不动，哪位大神帮我用牛顿迭代编一下
已经有4人回复
- 求用牛顿迭代法解非线性方程组的matlab程序
已经有4人回复
- 牛顿迭代法解方程。
已经有11人回复
- Matlab用牛顿法求解非线性方程组问题
已经有8人回复
- MATLAB牛顿-SOR迭代法解n维方程组
已经有0人回复
- 求：牛顿迭代法matlab程序及输出结果.求指点
已经有1人回复
- matlab 求解非线性方程组 牛顿解法
已经有1人回复
- 在matlab中用牛顿迭代法
已经有1人回复
- 使用matlab进行Newton-Raphson迭代是出现虚数的问题
已经有0人回复
- [请教]用牛顿迭代整这个问题
已经有7人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
夜的宁静
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 604.9
- 散金: 90
- 帖子: 93
- 在线: 74.4小时
- 虫号: 2307122
- 注册: 2013-02-28
- 性别: GG
- 专业: 机械动力学
送红花一朵 
很完美，太感谢啦 5楼2015-05-20 08:37:09
材料廖
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 4889.5
- 散金: 200
- 红花: 11
- 帖子: 891
- 在线: 295.1小时
- 虫号: 1794189
- 注册: 2012-05-04
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2015-05-18 11:45:01
夜的宁静
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 604.9
- 散金: 90
- 帖子: 93
- 在线: 74.4小时
- 虫号: 2307122
- 注册: 2013-02-28
- 性别: GG
- 专业: 机械动力学
|
还是不知道怎么改，我是初学，因为用到现学的，所以基础不怎么样。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('迭代次数过多，防止死循环终止'); else fprintf('迭代到%d次时得到结果%f\n',k,x) end |
3楼2015-05-18 22:23:32
材料廖
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 4889.5
- 散金: 200
- 红花: 11
- 帖子: 891
- 在线: 295.1小时
- 虫号: 1794189
- 注册: 2012-05-04
- 性别: GG
- 专业: 金属材料的合金相、相变及
|
试试这个看对不对 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('迭代次数过多，防止死循环终止'); else fprintf('迭代到%d次时得到结果%f\n',k,x) end |
» 本帖已获得的红花（最新10朵）
4楼2015-05-19 15:56:37
50