应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab牛顿迭代的问题，出现Matrix dimensions must agre的错误
51/1返回列表
查看: 702  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

夜的宁静

金虫 (小有名气)

[求助] 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

matlab牛顿迭代的问题，出现Matrix dimensions must agre的错误
360反馈意见截图16570207506364.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
Tomakeeverydayccount.
1楼 2015-05-17 19:33:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

材料廖

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
1.你的x0没有更新
2.报的错，调试时候发现feval(df,x0)是空矩阵
赞一下(1人) 回复此楼
2楼2015-05-18 11:45:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜的宁静

金虫 (小有名气)
引用回帖:
2楼: Originally posted by 材料廖 at 2015-05-18 11:45:01
1.你的x0没有更新
2.报的错，调试时候发现feval(df,x0)是空矩阵

还是不知道怎么改，我是初学，因为用到现学的，所以基础不怎么样。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
赞一下 回复此楼
Tomakeeverydayccount.
3楼2015-05-18 22:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

材料廖

木虫 (正式写手)
引用回帖:
3楼: Originally posted by 夜的宁静 at 2015-05-18 22:23:32
还是不知道怎么改，我是初学，因为用到现学的，所以基础不怎么样。x0没有更新什么意思
之前的赋值可能有些错误，改成这样还是不行，麻烦再看一下，十分感谢function =newton()
clc;
syms x;
f=@(x)2*log10(5.4 ...

试试这个看对不对
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
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

夜的宁静
4楼2015-05-19 15:56:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜的宁静

金虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by 材料廖 at 2015-05-19 15:56:37
试试这个看对不对
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 ...

很完美，太感谢啦
回复此楼
Tomakeeverydayccount.
5楼2015-05-20 08:37:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夜的宁静 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭