function [r,m]=mulDamp(F,x0,h,u,v,eps)
format long;
if nargin==5
eps=1.0e-6;
end
FI = transpose(F)*F/2;
n = length(x0);
x0 = transpose(x0);
m=1;
tol=1;
while tol>eps
j = 0;
fx = subs(F,findsym(F),x0);
J = zeros(n,n);
for i=1:n
x1 = x0;
x1(i) = x1(i)+h;
afx = subs(F,findsym(F),x1);
J(:,i) = (afx-fx)/h;
end
FIx = subs(FI,findsym(FI),x0);
for i=1:n
x2 = x0;
x2(i) = x2(i)+h;
gradFI(i,1) = (subs(FI,findsym(FI),x2)-FIx)/h;
end
s=0;
while s==0
A = transpose(J)*J+u*eye(n,n);
p = -A\gradFI;
r = x0 + p;
FIr = subs(FI,findsym(FI),r);
if FIr<FIx
if j == 0
u = u/v;
j = 1;
else
s=1;
end
else
u = u*v;
j = 1;
if norm(r-x0)<eps
s=1;
end
end
end
x0 = r;
tol = norm(p);
m=m+1;
if(m>100000) %迭代步数控制
disp('迭代步数太多,可能不收敛!');
return;
end
end
format short;
没有合适的资源?快使用搜索试试~ 我知道了~
MATLAB语言常用算法_非线性方程组求解.rar
共22个文件
m:22个
需积分: 5 0 下载量 163 浏览量
2023-06-06
13:40:23
上传
评论
收藏 11KB RAR 举报
温馨提示
本资料仅供参考学习
资源推荐
资源详情
资源评论
收起资源包目录
MATLAB语言常用算法_非线性方程组求解.rar (22个子文件)
MATLAB语言常用算法_非线性方程组求解
mulGXF2.m 1KB
mulRank1.m 677B
DiffParam1.m 494B
mulNewtonSOR.m 1KB
mulNewton.m 579B
mulConj.m 1KB
mulGXF1.m 863B
mulNewtonStev.m 781B
mulMix.m 1KB
mulStablePoint.m 516B
mulDiscNewton.m 767B
mulDamp.m 1KB
mulBFS.m 688B
SOR.m 849B
DiffParam2.m 688B
mulNumYT.m 895B
mulGSND.m 613B
mulDFP.m 659B
mulDNewton.m 759B
mulSimNewton.m 600B
mulFastDown.m 654B
mulVNewton.m 732B
共 22 条
- 1
资源评论
毕业季zwlzyhzylzyl
- 粉丝: 318
- 资源: 5650
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 瑞萨RL78系列MCU设备文件
- Java电子商务平台案例介绍:实现商品展示、购物车、订单管理和支付功能
- yolov8初始化权重文件
- Nacos控制台简介.pdf
- 基于STM32与ESP8266的智能家居系统,以STM32为核心控制单元,通过传感器监测室内环境,结合微信小程序实现数据传输与远
- ChineseCalendar-1.0.5-SNAPSHOT.jar
- 基于Bagging集成学习方法的情绪分类预测模型研究(数据集+代码).rar
- nacos-server-2.3.0-BETA.zip
- 双11活动数据分析-商品销售数据、促销活动的信息、用户行为数据
- 基于记忆与模型协同过滤的电影推荐系统研究与实践(数据集+代码).rar
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功