矩阵表达式

nkkopd

自定义的函数：
function f=largejfun(x)
' @% c) R, T7 @- h3 R% I" k4 Sf=(x-1./1:100).^2;
end

调用函数的时候：
clc
x0=10*ones(1,100);
  @+ e5 G' c6 ^' p! P" \" hoptions=optimset('display','iter','TolFun',1e-8);
0 j: @# |0 L$ }7 o[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)

% w; j5 E# ^3 w! E( u, P
运行结果出错，出错结果显示如下：
错误使用 fminunc (line 348)
' _' v; c7 D& i6 YSupplied objective function must return a scalar value.
  h+ w$ L' D8 {- e. t
' o1 p$ C7 P3 `8 R1 L, N  w% s6 q
出错 fminu1 (line 4)
- W9 V8 l5 C$ [. f3 ~9 }, n, S, U$ g; [[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)




: I- a  V7 y# O+ G0 Z9 v



( O, t8 Z& q$ N1 l7 @" I我将自定义函数改为如下：
function f=largejfun(x)
f=0;
& W: Q) v' x! [9 L/ Y7 H1 e& _0 ufor i=1:100
    f=f+(x(i)-1/i)^2;
. A+ r+ K6 H1 g, p' n1 Aend

最后运行结果没有任何问题。
8 s% k) S8 v$ N5 {3 ~& R那是哪里有问题呢?
) |3 t: }6 r9 J$ D0 o

kekek

很明显啊，f=(x-1./1:100).^2这种写法，甭管你的x是标量还是向量，1./1:100肯定是向量，做差之后还是向量(或矩阵)，再求平方，仍然不会是标量啊。
想转换成标量就用sum或者norm之类的函数，通过求和、求模来转化成标量。

Uifhjvv

楼上正解，楼主试试

小小鲁班

nkkopd

kekek 发表于 2020-8-17 10:48
很明显啊，f=(x-1./1:100).^2这种写法，甭管你的x是标量还是向量，1./1:100肯定是向量，做差之后还是向量( ...

5 E* }$ m. y# q, _% D( P) H谢谢大神指导
/ {; y  a" t6 k9 S) P! W( q