matlab最速下降法

function x=fsxsteep(f,e,a,b)

% fsxsteep函数最速下降法

% x=fsxsteep(f,e,a,b)为输入函数f为函数e为允许误差(a,b)为初始点; % fsx

x1=a;x2=b;

Q=fsxhesse(f,x1,x2);

x0=[x1 x2]';

fx1=diff(f,'x1'); %对x1求偏导数

fx2=diff(f,'x2'); %对x2求偏导数

g=[fx1 fx2]'; %梯度

g1=subs(g); %把符号变量转为数值

d=-g1;

while (abs(norm(g1))>=e)

t=(-d)'*d/((-d)'*Q*d); %求搜索方向

x0=x0-t*g1; %搜索到的点

v=x0;

a=[1 0]*x0;

b=[0 1]*x0;

x1=a;

x2=b;

g1=subs(g);

d=-g1;

end;

x=v;

function x=fsxhesse(f,a,b)

% fsxhesse函数求函数的hesse矩阵;

% 本程序仅是简单的求二次函数的hesse矩阵!;

% x=fsxhesse(f)为输入函数f为二次函数x1,x2为自变量;

% fsx

x1=a;x2=b;

fx=diff(f,'x1'); %求f对x1偏导数

fy=diff(f,'x2'); %求f对x2偏导数

fxx=diff(fx,'x1'); %求二阶偏导数对x1再对x1

fxy=diff(fx,'x2'); %求二阶偏导数对x1再对x2

fyx=diff(fy,'x1'); %求二阶偏导数对x2再对x1

fyy=diff(fy,'x2'); %求二阶偏导数对x2再对x2

fxx=subs(fxx); %将符号变量转化为数值

fxy=subs(fxy);

fyx=subs(fyx);

matlab 最速下降,matlab最速下降法相关推荐

  1. c++引用matlab类,matlab调用C++函数浅谈(一)

    由于在下才疏学浅,在网上看各高手指南时亦觉云里雾里,遂决定一切说明从最基础说起,一是方便自己(记性奇差),二是方便似我的小白.以下部分是我从各网站论坛等摘抄.重组.改写过的,以求更加详实明朗,由于参考 ...

  2. 磁盘驱动读取系统MATLAB仿真,matlab读写..doc

    matlab读写. MATLAB二进制数据文件的读写 (2011-06-04 19:44:27) 转载▼ 标签: easleyhux matlab 二进制读写 杂谈分类: MATLAB 所谓二进制格式 ...

  3. 【 MATLAB 】MATLAB 实现模拟信号采样后的重建(三)一阶保持(FOH)内插

    上篇博文采用了零阶保持(ZOH)的方式进行了重构:[ MATLAB ]MATLAB 实现模拟信号采样后的重建(二)零阶保持(ZOH) 这篇博文我们使用一阶保持(FOH)内插来重建信号,采用的案例依然是 ...

  4. 【 MATLAB 】MATLAB 实现模拟信号采样后的重建(二)零阶保持(ZOH)

    上篇博文采样sinc函数内插的方式实现了模拟信号的重建:[ MATLAB ]MATLAB 实现模拟信号采样后的重建(一) 这篇博文我们使用零阶保持器(ZOH)来重建信号,采用的案例依然是上篇博文中的案 ...

  5. Matlab:Matlab中常用的函数、案例详细攻略

    Matlab:Matlab中常用的函数.案例详细攻略 目录 常用函数 1.与文件相关 2.MATLAB GUI不同控件函数间变量传递方法 常用函数 Matlab中的bwmorph函数解释 bwmorp ...

  6. vb调用matlab工具箱,Matlab与VB集成 - 关于VB的经验之谈 - VB爱好者乐园(VBGood) - 关于VB的经验,电子教程,代码,控件,论坛,博客,微博等....

    工业生产中经常遇到复杂的数据信息处理问题,需要大运算量的矩阵计算及对分析结果进行实时.直观的图形化显示.Matlab是集数值分析.矩阵运算.信号处理和图形显示于一体的高性能数学软件,将其强大的计算功能 ...

  7. ccs matlab联调,MATLAB与CCS联调中idelink_ert_tlc和ert_tlc选择

    同问.希望高工给解释一下.我在用idelink_ert_tlc的时候就能通过2014B给DSP直接下载程序.但是用ert_tlc就不行,提示文件路径不对.问题描述如下,请高工帮忙解决一下,谢谢! ## ...

  8. matlab转向语句,MATLAB控制语句

    MATLAB控制语句 目的:研究控制结构(用于,是否,切换,中断,继续,输入/输出功能,读取和存储数据). 如果:如果计算逻辑表达式,并根据表达式的值执行一组语句. If语句的语法 例子 输出: Sw ...

  9. opencv matlab配置,Matlab下运行c++程序的配置(包含opencv的c++程序)

    本文使用的是matlab2012b.VS2010和OpenCV249 1.打开matlab 2.输入mex -setup Welcome to mex -setup.  This utility wi ...

  10. matlab netcdf,Matlab读取Netcdf文件

    NetCDF(网络通用数据格式)文件后缀为.nc,文件中的数据结构包含维(dimensions).变量(variables)和属性(attributes)三种描述类型,每种类型都会被分配一个名字和一个 ...

最新文章

  1. React + Koa 实现服务端渲染(SSR)
  2. 为什么要学习Python编程语言?哪些人适合学习Python?
  3. thrift的中文编码处理
  4. 基于AgileEAS.NET SOA 平台SAAS架构技术的开源分销ERP系统-SmartERP.NET下载配置说明
  5. C# 中几个关键词的使用
  6. 没有事情,错误1503_为什么依靠用户报告错误是您做过的最愚蠢的事情
  7. mysql 数据库 更新_mysql数据库更新
  8. 量化交易很好,但是也存在问题
  9. 试题18 四平方和(枚举法)
  10. echart中线性渐变的应用
  11. Apache 模块 mod_cache应用
  12. 利用VS2010模仿QQ2011登陆界面
  13. JAVA设计模式 — 生成器模式(Builder)
  14. 工具类ConfigTool封装Nacos Config 本地缓存(实战附代码实现)
  15. html 页面数据显示不全,解决了一个Web网页显示不全的BUG
  16. 平安京s9服务器维护,决战平安京S9赛季段位怎么继承_决战平安京S9赛季段位继承详情_素材吧...
  17. spring-data-jpa hiberneate 中 uui的生成策略
  18. boost::string_algo库详解
  19. QQ能联网但是电脑上所有的浏览器都不能上网
  20. 微信小程序子组件向父组件传值的两个方法

热门文章

  1. 竞品分析:得物APP、nice APP、识货APP
  2. R:导入其他样式数据的方法
  3. JAVA经典面试题附答案(持续更新版)
  4. 电视盒子为啥不安装鸿蒙,电视盒子视频播放不了怎么办?只需三步轻松解决!...
  5. navicat mysql 免安装_mysql免安装版配置+navicat测试
  6. jQueryWEUI自定义对话框-带有textarea
  7. mysql 分隔字符串的函数_Mysql 字符串分隔函数
  8. 记一个tcp udp测试工具ethrc
  9. 红外图像、灰度图像、深度图像对比
  10. django 下载返回文件流