[x1_k,x2_k]=meshgrid(-10:0.1:10,-10:0.1:10);
f1=(x1_k+3).*(x2_k.^3-7)+18;
f2=sin(x2_k.*exp(x1_k)-1);
plot3(x1_k,x2_k,f1);
plot3(x1_k,x2_k,f2)
clear;
clc;
eps=1e-8;
iters=1;
maxiters=50000;
x1=10*rand(1,1);
x2=10*rand(1,1);
x1_k=x1;
x2_k=x2;
f1=(x1_k+3)*(x2_k^3-7)+18;
f2=sin(x2_k*exp(x1_k)-1);
jacobi=[x2_k^3-7 3*(x1_k+3)*x2_k^2;cos(x2_k*exp(x1_k)-1)*x2_k*exp(x1_k) cos(x2_k*exp(x1_k)-1)*exp(x1_k)];
x_k=[x1_k;x2_k];
x_kk=x_k-jacobi\[f1;f2];
while (norm([f1,f2])>eps)
x_k=x_kk;
x1_k=x_k(1);
x2_k=x_k(2);
f1=(x1_k+3)*(x2_k^3-7)+18;
f2=sin(x2_k*exp(x1_k)-1);
jacobi=[x2_k^3-7 3*(x1_k+3)*x2_k^2;cos(x2_k*exp(x1_k)-1)*x2_k*exp(x1_k) cos(x2_k*exp(x1_k)-1)*exp(x1_k)];
x_kk=x_k-jacobi\[f1;f2];
iters=iters+1;
if iters>=maxiters
break;
end
end
x_1=x_kk(1,1)
x_2=x_kk(2,1)
iters
eps
cg-NT.zip_conjugate gradient_flagiyh_共轭梯度法_希尔伯特矩阵_牛顿方法
版权申诉
5星 · 超过95%的资源 84 浏览量
2022-07-13
23:29:15
上传
评论
收藏 1KB ZIP 举报
Kinonoyomeo
- 粉丝: 75
- 资源: 1万+
最新资源
- MyBatis动态SQL是一种强大的特性,它允许我们在SQL语句中根据条件动态地添加或删除某些部分,从而实现更加灵活和高效的数据
- MyBatis动态SQL是一种强大的特性,它允许我们在SQL语句中根据条件动态地添加或删除某些部分,从而实现更加灵活和高效的数据
- 【ERP标准流程-标准流程-进货管理】(DOC 17页).doc
- MyBatis动态SQL是一种强大的特性,它允许我们在SQL语句中根据条件动态地添加或删除某些部分,从而实现更加灵活和高效的数据
- MyBatis动态SQL是一种强大的特性,它允许我们在SQL语句中根据条件动态地添加或删除某些部分,从而实现更加灵活和高效的数据
- Java爬虫项目【项目开发计划】(共12页).docx
- 11111111111
- MyBatis动态SQL是一种强大的特性,它允许我们在SQL语句中根据条件动态地添加或删除某些部分,从而实现更加灵活和高效的数据
- MyBatis动态SQL是一种强大的特性,它允许我们在SQL语句中根据条件动态地添加或删除某些部分,从而实现更加灵活和高效的数据
- MyBatis动态SQL是一种强大的特性,它允许我们在SQL语句中根据条件动态地添加或删除某些部分,从而实现更加灵活和高效的数据
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
评论1