matlab迭代法代码
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
matlab 迭代法代码
1、％用不动点迭代法求方程x-e A x+4=0的正根与负根，误差限是
10A-6% disp(' 不动点迭代法 ');
n0=100;
p0=-5;
for i=1:n0
p=exp(p0)-4;
if abs(p-p0)<=10(6)
if p<0
disp('|p-p0|=')
disp(abs(p-p0))
disp(' 不动点迭代法求得方程的负根为 :')
disp(p);
break;
else
disp(' 不动点迭代法无法求出方程的负根 .')
end
else
p0=p;
end
end
if i==n0 disp(n0) disp(' 次不动点迭代后无法求出方程的负根')
end
p1=1.7;
for i=1:n0
pp=exp(p1)-4;
if abs(pp-p1)<=10(6)
if pp>0
disp('|p-p1|=')
disp(abs(pp-p1))
disp(' 用不动点迭代法求得方程的正根为 ')
disp(pp);
else
disp(' 用不动点迭代法无法求出方程的正根 ');
end
break;
else
p1=pp;
end
end
if i==n0
disp(n0)
disp(' 次不动点迭代后无法求出方程的正根 ') end
2、％用牛顿法求方程x-e A x+4=0的正根与负根，误差限是disp(' 牛顿法')
n0=80;
p0=1;
for i=1:n0 p=p0-(p0-exp(p0)+4)/(1-exp(p0));
if abs(p-p0)<=10(6)
disp('|p-p0|=')
disp(abs(p-p0))
disp(' 用牛顿法求得方程的正根为 ')
disp(p);
break;
else
p0=p;
end
end
if i==n0
disp(n0)
disp(' 次牛顿迭代后无法求出方程的解
p1=-3;
for i=1:n0 p=p1-(p1-exp(p1)+4)/(1-exp(p1)); 10A-6
') end
if abs(p-p1)<=10A(-6) disp('|p-p1|=') disp(abs(p-p1)) disp(' 用牛顿法求得方程的负根为 ') disp(p); break; else
p1=p;
end
end
if i==n0
disp(n0)
disp(' 次牛顿迭代后无法求出方程的解') end