高斯牛顿法 matlab代码实现
高斯牛顿法的解释请参见
http://blog.csdn.net/dsbatigol/article/details/12448627
%% Gauss-Newton method
clc;
clear;
format long;
syms x1 x2 x3 x4;
error=1e-5; %停机门限
f = [x1+10*x2; 5^(1/2)*(x3-x4); (x2-2*x3)^2; 10^(1/2)*(x1-x4)^2];
v=[x1 x2 x3 x4];
j=jacobian(f,v); %求jacobian行列式
x=[3;-1;0;1]; %初始迭代点
J=subs(j,v,x'); %初点数值带入表达式
F=subs(f,v,x');
k=0; %迭代次数
tic
while (sum((J'*F).^2))^(1/2)>error %判断停机与否, 这里相当于求梯度的模d=-inv(J'*J)*J'*F; %搜索方向x=x+d; %新的迭代点J=subs(j,v,x'); %新的迭代点数值带入表达式 F=subs(f,v,x');k=k+1; %迭代次数加1
end
disp('Gauss-Newton algorithm');
toc
k
x
sigma=(sum((J'*F).^2))^(1/2)
F=(sum(F.^2))^(1/2) %显示迭代次数,变量取值,停机表达式值,目标函数值
%%
高斯牛顿法 matlab代码实现相关推荐
- 【最优化理论】牛顿法+Matlab代码实现
文章目录 1 牛顿法简介 2 牛顿法原理 3 牛顿法推导 4 Matlab代码实现 5 低版本Matlab报错 1 牛顿法简介 牛顿迭代法(Newton's method)又称为牛顿-拉夫逊(拉弗森) ...
- 高斯牛顿法----MATLAB实现
一次作业:参考于下面这个博客:https://blog.csdn.net/zhubaohua_bupt/article/details/74973347 感觉这个博客写的是真的不错. function ...
- MATLAB图像处理--高斯低通滤波、高斯高通滤波(代码及示例)
本文目录 常用的滤波器 高斯低通滤波器 matlab代码 高斯高通滤波器 Butterworth 低通滤波器 常用的滤波器 常用的滤波器有以下几种: 频域滤波器 频域平滑滤波器 理想低通滤波器 巴特沃 ...
- Ceres 库:基础使用,以手写高斯-牛顿法为例
Ceres 库 简介 Ceres库为Google开发的开源C++非线性优化库,被广泛使用于求解最小二乘问题. Ceres库的Github主页如下: 安装 首先,下载Cere的源码: git clone ...
- 高斯牛顿算法matlab代码,matlab实现高斯牛顿法、Levenberg–Marquardt方法
高斯牛顿法: function [ x_ans ] = GaussNewton( xi, yi, ri) % input : x = the x vector of 3 points % y = th ...
- 基于BP-Adaboost神经网络的高斯函数模型预测研究-含Matlab代码
目录 一.引言 二. 预测建模及原理 2.1 BP 神经网络模型 2.2 BP-Adaboost 模型 三.模型预测结果对比 3.1 数据模拟 3.2 预测结果对比 四.参考文献 五.Matlab代码 ...
- 【slam十四讲第二版】【课本例题代码向】【第七讲~视觉里程计Ⅱ】【使用LK光流(cv)】【高斯牛顿法实现单层光流和多层光流】【实现单层直接法和多层直接法】
[slam十四讲第二版][课本例题代码向][第七讲~视觉里程计Ⅱ][使用LK光流(cv)][高斯牛顿法实现单层光流和多层光流][实现单层直接法和多层直接法] 0 前言 1 使用LK光流(cv) 1.1 ...
- 【角点检测】 基于各向异性高斯方向导数滤波器实现图像角点检测附matlab代码
1 内容介绍 为了改进噪声鲁棒性和定位准确性,利用各向异性高斯方向导数滤波器,提出多方向角点检测算法.该算法利用一组各向异性高斯方向导数滤波器对输入图像进行卷积处理得到各个方向的滤波器响应.对于每个像 ...
- 【图像去噪】基于二维双边高斯滤波实现图像去噪附matlab代码
1 简介 图像是生活中重要的信息来源,处理图像有助于理解信息的基本信息.但图像本身可能存在一些被干扰的信息或者噪声.研究了基于高斯滤波和双边滤波算法的数字图像处理技术用于对图像的噪声进行消除.通过对图 ...
最新文章
- python导入外部包_您会喜欢的10个外部Python软件包
- 模仿并超越人类围棋手,KL正则化搜索让AI下棋更像人类,MetaCMU出品
- 实验室中的机械臂-资料汇总
- [luoguP1640] [SCOI2010]连续攻击游戏(二分图最大匹配)
- in作为介词的用法_英文论文写作:极易混淆的短语用法汇总如下
- linux-04-磁盘命令+进程命令
- mysql call procedure into_mysql存储过程之创建(CREATE PROCEDURE)和调用(CALL)及变量...
- 英语学习过程中的几点体会(2)
- 上市前夕再陷版权纠纷,快手面临的风险却不止于此
- 大数据Hadoop基本概念介绍
- cs 5 下载地址 、视频教程、安装教程
- 智能问答机器人python_帮帮智能问答机器人中TaskBot任务对话算法实践
- 计算机网络 ——物理层
- codewars练习(javascript)-2021/3/16
- Codeforces140CNew Year Snowmen
- 【猿说VUE】Vue列表渲染
- 2021年秋招面试真题以及面试技巧分享
- 纯CSS实现三角形图标
- 陶瓷电容NP0、C0G、X7R、X5R、Y5V、Z5U的区别
- Dubbo 支持哪些通信协议?
热门文章
- 机器视觉:光源专业词汇中英文详解
- C++中回调函数(CALLBACK)初探
- 外文翻译 matlab与矩阵处理,外文文献及译 MATLAB 介绍.doc
- matlab里的scalar,scalar_scalar与vector之间的区别
- Solr 中 Schema 结构说明
- pfSense添加子网的几种方式
- 基于Office 365 无代码工作流分析-需求基本分析!
- 真机开包!国产至强5600服务器35张图赏
- 代码理解(2009.11.20)
- 一周只指定一个主要目标,当主要目标完成时,就已经完成了任务的80%