牛顿迭代法求解非线性方程组 matlab 200
用牛顿迭代法或二分法求下面非线性方程组的根{█(f_1(x,y)=x+2y-3=0@〖f〗_2(x,y)=2x^2+y^2-5=0)┤解：将f1中x的值代入f2中得，g(...
用牛顿迭代法或二分法求下面非线性方程组的根
{█(f_1 (x,y)=x+2y-3=0@〖 f〗_2 (x,y)=2x^2+y^2-5=0)┤
解：将f1中x的值代入f2中得，
g(y)=9y2-24y+13=0
g^'(y)=18y-24
则由牛顿迭代法：
设g'(y_(n-1))≠0
有牛顿迭代式yn=yn-1-(g(y_(n-1)))/(g'(y_(n-1))) , (n=1,2,3,…)
因为g(1)g(2)＜0
故令初值y0=2，在matlab中计算迭代结果，过程如下
g=inline('9*y^2-24*y+13');
gy=inline('18*y-24 ');
a=2; e=0.0001;n=1;
while abs(y-a)>e
y=a-(g(a)/gy(a));
a= y;n=n+1;
end
a,n
结果是：
a=1.9167
n=2
又令初值y0=1，在matlab中计算迭代结果，过程如下
g=inline('9*y^2-24*y+13');
gy=inline('18*y-24 ');
a=1; e=0.0001;n=1;
while abs(y-a)>e
y=a-(g(a)/gy(a));
a= y;n=n+1;
end
a,n
结果是：
a=0.6667
n=2
想不通哪里错了，答案是错误的，y应该是1.9107和0.7560 展开
{█(f_1 (x,y)=x+2y-3=0@〖 f〗_2 (x,y)=2x^2+y^2-5=0)┤
解：将f1中x的值代入f2中得，
g(y)=9y2-24y+13=0
g^'(y)=18y-24
则由牛顿迭代法：
设g'(y_(n-1))≠0
有牛顿迭代式yn=yn-1-(g(y_(n-1)))/(g'(y_(n-1))) , (n=1,2,3,…)
因为g(1)g(2)＜0
故令初值y0=2，在matlab中计算迭代结果，过程如下
g=inline('9*y^2-24*y+13');
gy=inline('18*y-24 ');
a=2; e=0.0001;n=1;
while abs(y-a)>e
y=a-(g(a)/gy(a));
a= y;n=n+1;
end
a,n
结果是：
a=1.9167
n=2
又令初值y0=1，在matlab中计算迭代结果，过程如下
g=inline('9*y^2-24*y+13');
gy=inline('18*y-24 ');
a=1; e=0.0001;n=1;
while abs(y-a)>e
y=a-(g(a)/gy(a));
a= y;n=n+1;
end
a,n
结果是：
a=0.6667
n=2
想不通哪里错了，答案是错误的，y应该是1.9107和0.7560 展开
展开全部
g = inline('9*y^2-24*y+13');
gy = inline('18*y-24 ');
a = 2;
e1=0.0001;n=1;
y = a-(g(a)/gy(a));
z = abs(y-a);
while z>e1
y = a-(g(a)/gy(a));
z = abs(y-a);
a = y;n=n+1;
end
a,n
本回答被网友采纳
已赞过
已踩过<
你对这个回答的评价是？
系科仪器
2024-08-02 广告
科仪器致力于为微纳薄膜领域提供精益级测量及控制仪器，包括各种光谱椭偏、激光椭偏、反射式光谱等，从性能参数、使用体验、价格、产品可靠性及工艺拓展性等多个维度综合考量，助客户提高研发和生产效率，以及带给客户更好的使用体验。...
点击进入详情页
本回答由系科仪器提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
