function Newton()
x0=[0.1;0.5];
x1=x0-inv(myJacobi(x0))*myfun(x0);
while norm(x1-x0)>1e-3
x0=x1;
x1=x0-inv(myJacobi(x0))*myfun(x0);
end
x1
msgbox('MATLAB编程答疑,请加QQ: 1530497909','MATLAB答疑','help')
web http://url.cn/NSFcAs -browser
end
function f=myfun(x)
syms x1 x2
f1=(15*x1+10*x2)-((40-30*x1-10*x2)^2*(15-15*x1))*5e-4;
f2=(15*x1+10*x2)-((40-30*x1-10*x2)*(10-10*x2))*4e-2;
f=[f1;f2];
x1=x(1);
x2=x(2);
f=eval(f);
end
function J=myJacobi(x)
syms x1 x2
f1=(15*x1+10*x2)-((40-30*x1-10*x2)^2*(15-15*x1))*5e-4;
f2=(15*x1+10*x2)-((40-30*x1-10*x2)*(10-10*x2))*4e-2;
J=[diff(f1,x1) diff(f1,x2)
diff(f2,x1) diff(f2,x2)];
x1=x(1);
x2=x(2);
J=eval(J);
end
[MATLAB项目实例源码]牛顿法求解非线性方程组 源程序代码.zip
134 浏览量
2023-08-27
02:26:45
上传
评论
收藏 1KB ZIP 举报
心兰相随引导者
- 粉丝: 775
- 资源: 5486
最新资源
- 基于python的高性能爬虫程序,使用了多线程+缓存+xpath实现的,这里以彼-岸图库为例,实现,仅用于学习交流
- 中分辨率成像光谱仪(MODIS)烧毁面积产品信息MODIS-C6-BA-User-Guide-1.2.pdf
- Screenshot_20240427_172613_com.huawei.browser.jpg
- 关于学习Python的相关资源网站链接及相关介绍.docx
- (HAL库)基于STM32F103C8T6的温控PID系统[Dht11、ESP8266、无线透传、L298N……]
- VoLTE高丢包优化指导书.xlsx
- Rust资源文件.zip
- 前后端分离实践:使用 React 和 Express 搭建完整登录注册流程
- gradle-publish-to-MavenLocal.zip
- 10份网络优化创新案例.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈