matlab - matlab -- fminsearch 与矩阵的使用

Question

如果我有一个功能，例如：

k=1:100
func=@(s) sum(c(k)-exp((-z(k).^2./s)))

其中 c 和 z 是具有相同大小的矩阵（例如 1x100），有没有办法使用 fminsearch 来查找“s”值？

Answer 1

我猜你想找到argmin你的符号函数， 在数组中使用最大值和最小值的索引

或者

Marco Cococcioni 的 ARGMAX/ARGMIN：

function I = argmax(X, DIM)
%ARGMAX    Argument of the maximum
%   For vectors, ARGMAX(X) is the indix of the smallest element in X. For matrices,
%   MAX(X) is a row vector containing the indices of the smallest elements from each
%   column. This function is not supported for N-D arrays with N > 2.
%
%   It is an efficient replacement to the use of [Y,I] = MAX(X,[],DIM);
%   See ARGMAX_DEMO for a speed comparison.
%
%   I = ARGMAX(X,DIM) operates along the dimension DIM (DIM can be 
%   either 1 or 2).
%
%   When complex, the magnitude ABS(X) is used, and the angle
%   ANGLE(X) is ignored. This function cannot handle NaN's.
%
%   Example:
%       clc
%       disp('If X = [2 8 4; 7 3 9]');
%       disp('then argmax(X,1) should be [2 1 2]')
%       disp('while argmax(X,2) should be [2 3]''. Now we check it:')
%       disp(' ');
%       X = [2 8 4; 
%            7 3 9]
%       argmax(X,1)
%       argmax(X,2)
%
%   See also ARGMIN, ARGMAXMIN_MEX, ARGMAX_DEMO, MIN, MAX, MEDIAN, MEAN, SORT.

%   Copyright Marco Cococcioni, 2009.
%   $Revision: 1.0 $  $Date: 2009/02/16 19:24:01$

if nargin < 2,
    DIM = 1;
end

if length(size(X)) > 2,
    error('Function not provided for N-D arrays when N > 2.');
end

if (DIM ~=1 && DIM ~= 2),
    error('DIM has to be either 1 or 2');
end

if any(isnan(X(:))),
    error('Cannot handle NaN''s.');    
end

if not(isreal(X)),
    X = abs(X);
end

max_NOT_MIN = 1; % computes argmax
I = argmaxmin_mex(X, DIM, max_NOT_MIN);

Answer 2

fminsearch在第二个参数中需要一个初始条件，而不是边界条件（尽管某些选项可能支持边界）。

打电话

 fminsearch(func,-0.5)

您在向量中看到的示例是跨多个系数的多维搜索，并且向量是每个系数的初始值。不限制搜索空间。

你也可以使用

fminbnd(func, -0.5, 1);

它执行约束最小化。

但我认为你应该最小化误差的范数，而不是总和（最小化总和会导致很大的误差幅度——非常非常负）。

如果您有优化工具箱，那么lsqnonlin可能会很有用。