| 查看: 741 | 回复: 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 |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有204人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解带积分的非线性方程组，牛顿迭代过程中的奇异矩阵如何修正
已经有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人回复
材料廖
木虫 (正式写手)
- 应助: 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
夜的宁静
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 604.9
- 散金: 90
- 帖子: 93
- 在线: 74.4小时
- 虫号: 2307122
- 注册: 2013-02-28
- 性别: GG
- 专业: 机械动力学
很完美，太感谢啦 5楼2015-05-20 08:37:09
20