【答案】应助回帖

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有108人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

有限差分法化之后化简为线性方程组已经有50人回复
求教如何用matlab解带约束的非线性方程组已经有4人回复
matlab 非线性微分方程求解已经有3人回复
用matlab求解一个非线性方程组的解已经有4人回复
大家能给看看这个方程属于什么类型，应该怎么求解？已经有10人回复
用fortran程序遗传算法解非线性方程组已经有7人回复
MATLAB求解非线性方程组已经有5人回复
求Matlab解非线性方程代码已经有12人回复
拟牛顿法求解复杂（带根号的）非线性方程组已经有17人回复
对流扩散方程迭代求解发散已经有11人回复
第5章_解线性方程组的直接方法.ppt 已经有13人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
求助matlab解非线性规划。求代码。已经有13人回复
求助matlab---fsolve解非线性方程组已经有6人回复
请问牛顿法和高斯牛顿法有什么不同？非常感谢已经有6人回复
matlab求解非线性方程组已经有16人回复
求一个mathmetica的程序已经有10人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】求教matlab解非线性方程组已经有9人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序已经有10人回复

1楼2014-05-08 20:31:28

reko34

木虫 (正式写手)

应助: 117 (高中生)
金币: 4794.2
散金: 22
红花: 12
帖子: 780
在线: 417.4小时
虫号: 3055206
注册: 2014-03-15
性别: GG
专业: 机械测试理论与技术

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
bsn_sky: 金币+50, ★★★★★最佳答案 2014-05-10 08:28:09

后面那个function文件要单独存成另一个m文件，名字也用fx，内容改成
function y=fx(x,x0,x1,a,b)
fx=x^5+x^3+x+a-b*(x+x0+x1);

主程序我不知道常数值，迭代初始值要根据对具体问题的了解定，所以你自己试吧

5楼2014-05-09 20:24:51

查看全部 5 个回答

bsn_sky

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 88.5
帖子: 76
在线: 21.7小时
虫号: 2194764
注册: 2012-12-19
专业: 金属结构材料

第二行应该是2到Nz-1

2楼2014-05-08 21:11:35

895701777

木虫 (著名写手)

进城农民工

应助: 36 (小学生)
金币: 1384.8
散金: 706
红花: 12
帖子: 1547
在线: 287.9小时
虫号: 2052402
注册: 2012-10-10
性别: GG
专业: 新闻理论

【答案】应助回帖

感谢参与，应助指数 +1

牛顿迭代很简单的，复杂的是高维的牛顿迭代

天道酬勤

3楼2014-05-09 09:23:57

reko34

木虫 (正式写手)

应助: 117 (高中生)
金币: 4794.2
散金: 22
红花: 12
帖子: 780
在线: 417.4小时
虫号: 3055206
注册: 2014-03-15
性别: GG
专业: 机械测试理论与技术

【答案】应助回帖

感谢参与，应助指数 +1

白天没机会碰MATLAB，急的话先试，有问题的话等我晚上试，但注释的2个东西要LZ根据实际情况自己定。
MATLAB解N维牛顿法有优势的，我之前用EXCEL编的，可能会有语法搞起来。

a=?;
b=?;
x=?; %迭代初始向量

Jacob=zeros(n,n);
y=zeros(n,1);

while (stopcriterion) %迭代停止条件自己定
pdx=5*x.^4+3*x.^2+1-b;
for i=1:n
Jacob(i,i)=pdx(i);
end
y(1)=x(1)^5+x(1)^3+x(1)+a-b*(x(1)+x(2));
for i=2:n-1
y(i)=fx(x(i),x(i-1),x(i+1),a,b);
Jacob(i-1,i)=-b;
Jacob(i,i-1)=-b;
end
y(n)=x(n)^5+x(n)^3+x(n)+a-b*(x(n)+x(n-1));
Jacob(n-1,n)=-b;
Jacob(n,n-1)=-b;
dx=y/Jacob;
x=x+dx;
end
%向量x是结果

function fx(x,x0,x1,a,b)
fx=x^5+x^3+x+a-b*(x+x0+x1);
end

4楼2014-05-09 12:49:16