矩阵表达式

nkkopd

自定义的函数：
( J) R1 O4 H. k; Ufunction f=largejfun(x)
: B! h3 e5 g- J7 E6 Ef=(x-1./1:100).^2;
) V. Q/ E2 B. Aend

调用函数的时候：
clc
x0=10*ones(1,100);
$ e6 C/ w3 _1 h$ m+ Goptions=optimset('display','iter','TolFun',1e-8);
[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)


7 ^8 q( H: t: G! ^9 }运行结果出错，出错结果显示如下：
/ m( j3 h' L0 f, o3 E+ Z" u错误使用 fminunc (line 348)
" J# _! W" s6 r+ J) fSupplied objective function must return a scalar value.

* e; V1 B. y& Q' O7 H
出错 fminu1 (line 4)
[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)
6 }+ x0 j+ s) y  i4 o1 r
5 j. E. }! g7 ?% I2 E- K
  f- A, H- F' \8 L' Q- q: j


; Y( d; o( P$ r; R( \: s! c
- O' s9 M: [5 W+ O  g

我将自定义函数改为如下：
function f=largejfun(x)
9 M% a- {/ B2 L) uf=0;
for i=1:100
    f=f+(x(i)-1/i)^2;
end

5 }' ]+ m: K) q5 ]! a7 ~8 e1 b$ q最后运行结果没有任何问题。
5 @( j3 a' D4 ?那是哪里有问题呢?
' w& Z+ E1 b2 I# ^0 X/ b6 U7 o( z

kekek

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

Uifhjvv

楼上正解，楼主试试

小小鲁班

nkkopd

kekek 发表于 2020-8-17 10:48
$ t  h, t( y5 S1 J% g" s2 Y很明显啊，f=(x-1./1:100).^2这种写法，甭管你的x是标量还是向量，1./1:100肯定是向量，做差之后还是向量( ...

谢谢大神指导