楼主: QINGTIAN607
1665 3

[其他] 求助,matlab,hessian语句修改。急求急求急求!!! [推广有奖]

  • 1关注
  • 1粉丝

已卖:74份资源

教授

74%

还不是VIP/贵宾

-

威望
0
论坛币
8957 个
通用积分
50.2815
学术水平
4 点
热心指数
3 点
信用等级
3 点
经验
94226 点
帖子
1078
精华
0
在线时间
890 小时
注册时间
2013-8-24
最后登录
2025-8-3

楼主
QINGTIAN607 发表于 2014-6-2 21:19:17 |AI写论文
10论坛币
语句如下:
H = hessian(f,x,varargin);
% PURPOSE: Computes finite difference Hessian
% -------------------------------------------------------
% Usage:  H = hessian(func,x,varargin)
% Where: func = function name, fval = func(x,varargin)
%           x = vector of parameters (n x 1)
%    varargin = optional arguments passed to the function
% -------------------------------------------------------
% RETURNS:
%           H = finite differnce hessian
% -------------------------------------------------------

% Code from:
% COMPECON toolbox [www4.ncsu.edu/~pfackler]
% documentation modified to fit the format of the Ecoometrics Toolbox
% by James P. LeSage, Dept of Economics
% University of Toledo
% 2801 W. Bancroft St,
% Toledo, OH 43606
% jlesage@spatial-econometrics.com

eps = 1e-6;

n = size(x,1);
fx = feval(f,x,varargin{:});

% Compute the stepsize (h)
h = eps.^(1/3)*max(abs(x),1e-2);
xh = x+h;
h = xh-x;   
ee = sparse(1:n,1:n,h,n,n);

% Compute forward step
g = zeros(n,1);
for i=1:n
  g(i) = feval(f,x+ee(:,i),varargin{:});
end

H=h*h';
% Compute "double" forward step
for i=1:n
for j=i:n
  H(i,j) = (feval(f,x+ee(:,i)+ee(:,j),varargin{:})-g(i)-g(j)+fx)/H(i,j);
  H(j,i) = H(i,j);
end
end
错误:Undefined function or variable 'f'。小弟是新手,急求!!!

最佳答案

klldmofashi 查看完整内容

这是个MATLAB函数,你需要在外面调用,H = hessian(f,x,varargin) 这个里面 f x 和 varargin 是输入的参数,H是返回给你的hessian矩阵输出
关键词:Hessian MATLAB matla atlab HESS matlab

沙发
klldmofashi 发表于 2014-6-2 21:19:18
这是个MATLAB函数,你需要在外面调用,H = hessian(f,x,varargin) 这个里面 f x 和 varargin 是输入的参数,H是返回给你的hessian矩阵输出

藤椅
Kīssˊ貝! 发表于 2014-6-2 21:25:28
哈哈哈哈

板凳
QINGTIAN607 发表于 2014-6-3 09:41:05
klldmofashi 发表于 2014-6-3 05:53
这是个MATLAB函数,你需要在外面调用,H = hessian(f,x,varargin) 这个里面 f x 和 varargin 是输入的参数, ...
谢谢您的回复,我是新手,请问一下在外面怎么调用,具体语句是什么,谢谢!!!

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
jg-xs1
拉您进交流群
GMT+8, 2026-1-1 02:59