(1)

clear all;

close all;

clc;

img=double(imread('lena.jpg'));

[m n]=size(img);

w=4;

sigma=1.2;

[x y]=meshgrid(-w:w,-w:w);

%高斯函数对应的二阶偏导,如何求导数请自行脑补

Dxx = 1/(2*pi*sigma^4)*(x.^2/sigma^2-1)*exp(-(x.^2+x.^2)/(2*sigma^2)); Dyy = 1/(2*pi*sigma^4)*(y.^2/sigma^2-1)*exp(-(x.^2+y.^2)/(2*sigma^2)); Dxy = 1/(2*pi*sigma^6)*(x.*y)*exp(-(x.^2+y.^2)/(2*sigma^2));

Ixx=imfilter(img,Dxx,'replicate');

Iyy=imfilter(img,Dyy,'replicate');

Ixy=imfilter(img,Dxy,'replicate');

H=cell(m,n);

for i=1:m

for j=1:n

H{i,j}=[Ixx(i,j) Ixy(i,j);Ixy(i,j) Iyy(i,j)];

end

end

(2)

function H=hessian(f,x,x0)

% H=hessian(f,x) 计算表达式

% H=hessian(f,x,x0) 计算hessian矩阵的值 x0为x的初值

%$copyright by$ LUO sir

switch nargin

case 1

error('please input variables in f(x)')

case 2

H=subhessian(f,x);

case 3

H=subhessian(f,x,x0);

运用matlab求黑塞矩阵,用matlab求黑塞矩阵相关推荐

  1. Matlab实现连通域标记算法求图像连通域

    Matlab实现连通域标记算法求图像连通域 连通域 连通域标记算法 连通域 连通区域(Connected Component)一般是指图像中具有相同像素值且位置相邻的前景像素点组成的图像区域(Regi ...

  2. matlab 曲线拟合求导,基于matlab曲线拟合的数据预测分析

    龙源期刊网 http://www.doczj.com/doc/2f22f716cbaedd3383c4bb4cf7ec4afe05a1b148.html 基于matlab曲线拟合的数据预测分析 作者: ...

  3. 求连通域面积matlab

    matlab有自带函数bwreopen()求连通域,并移除小面积连通域,返回一个mask模板不同连通域用不同的label标志,背景用0标志. function [ mask ] = bwareopen ...

  4. matlab求微分方程精确解,matlab求微分方程精确解及近似解.ppt

    matlab求微分方程精确解及近似解.ppt 还剩 24页未读, 继续阅读 下载文档到电脑,马上远离加班熬夜! 亲,喜欢就下载吧,价低环保! 内容要点: 求微分方程的解q 自牛顿发明微积分以来,微分方 ...

  5. 运用数学软件matlab求无穷积分,matlab积分的计算及其简单应用论文.doc

    积分的计算及其简单应用 摘要:本文简要的概述了MATLAB 在高等数学中积分的计算及应用:利用MATLAB 中符号积分和数值积分的命令,计算定积分和不定积分.同时,也可以通过这些命令来解决一些实际问题 ...

  6. ubuntu14.04扩展屏幕后,打开matlab就死机。求大神帮忙解决!!!万分感谢

    ubuntu14.04扩展屏幕后,打开matlab就死机.求大神帮忙解决!!!万分感谢

  7. matlab怎么求标准化特征向量,matlab特征向量标准化

    % S1: S2: S3: S4: AHP 的 matlab 程序实现构造比较矩阵 求特征根和特征向量 求得最大特征根的特征向量 做一致性检验 clear; A=[1 1.2 1.5 1.5; % . ...

  8. matlab输出的特征向量,关于matlab中的eig函数(求特征值和特征向量)(最新整理)

    <关于matlab中的eig函数(求特征值和特征向量)(最新整理)>由会员分享,可在线阅读,更多相关<关于matlab中的eig函数(求特征值和特征向量)(最新整理)(3页珍藏版)& ...

  9. MATLAB新手简明使用教程(七)——使用matlab建立多项式以及求导,商求导乘积求导等——新手来看,保证看懂。

    前期回顾 上一期中,我们学了下面的知识: 定积分的基本概念和一些简单的几何意义. 使用 int 函数计算不定积分. 使用 int 函数计算定积分. 本期内容 本期我打算给大家介绍一下使用matlab对 ...

  10. matlab 求曲面体积,matlab求两曲面之间的体积

    MATLAB求曲面相交所成空间曲线的图形 放在你程序后也可,单独运行也行:t=-0.1:0.1:2*pi;x=2*cos(t);%交线参数方程z=2*sin(t);y1=sqrt(5)*ones(si ...

最新文章

  1. C/C++各种数据类型转换汇总
  2. 搭建你的Spring.Net+Nhibernate+Asp.Net Mvc 框架
  3. Python技术分享:教你如何爆破隔壁WiFi密码
  4. Prthon编码和解码/is 和 ==
  5. Django RestFramework BaseSerializer
  6. pooled-jms_Hibernate隐藏的宝石:pooled-lo优化器
  7. php 课程节次周次提取,用PHP写一个最简单的解释器Part5(计算器最后一节,下节开始如何写个脚本语言)...
  8. 课程设计——企业网络项目搭建(下)
  9. MYSQL数据库基础概念
  10. 【Codeforces Round #555 (Div. 3) G】Inverse of Rows and Columns【bitset优化暴力...】
  11. 软件项目管理案例教程(第3版)复习内容
  12. 关于虚拟机中IPI中断的思考
  13. java项目编码问题解决
  14. 秘宝 联想正式进入元宇宙 互联网巨头的数字藏品用的哪种技术?
  15. 直流斩波电路在matlab中的建模与仿真,基于MATLAB/Simulink的直流斩波电路分析
  16. 加班申请 ----中间表--系统自动算出---可调休天数
  17. Pandownload 之后,非它莫属
  18. APP绕过模拟器-小白版
  19. 认认真真写博客,踏踏实实编代码
  20. python --- 常见题目 2019.01.03

热门文章

  1. 信息安全概论复习笔记
  2. Autodesk Maya 2009标准培训教材
  3. kali linux 源码免杀,shellcode超级免杀
  4. HTML基本语法总结
  5. 第一章:Ruby 安装 - Windows
  6. linux 网络对时
  7. 嵌入式Linux系统驱动hp1020打印机
  8. 【数据库】MySQL单表查询
  9. java 学习资料 分享
  10. 惠普计算机图标不在桌面,我的电脑图标没了怎么办