矩阵表达式

nkkopd

自定义的函数：
8 q' _, B" ^2 b- k& j& Zfunction f=largejfun(x)
( E* j! g. H0 `f=(x-1./1:100).^2;
$ N+ \; C9 @2 i9 nend
. V! ^2 N% j" y( `& v! C+ t
调用函数的时候：
clc
' Y- H) q7 v; F/ sx0=10*ones(1,100);
options=optimset('display','iter','TolFun',1e-8);
1 A" `% {2 M# T( p[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)


/ q& {2 ^" q2 F4 C# t2 [# c/ }; g运行结果出错，出错结果显示如下：
+ y& B8 R* @2 b( y# r# Z6 U错误使用 fminunc (line 348)
Supplied objective function must return a scalar value.
5 C- O( {% w* o% W

& @2 L  \% [: R1 p+ K: b( M出错 fminu1 (line 4)
[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)


& L4 t( J. r% Y


; t2 P' P4 \- R, |7 J" E
0 t9 V& o4 F; w4 o0 Q) O& T
/ c: o. l$ i( q3 }, n' h6 t4 p
我将自定义函数改为如下：
. I; z7 ~2 N' P9 D4 Bfunction f=largejfun(x)
4 u+ T5 T- S& v  g4 ?f=0;
% s; J9 X5 O; G- R+ k4 ]% Nfor i=1:100
7 y. c9 J- x2 X  M) r$ L    f=f+(x(i)-1/i)^2;
7 @  t* l$ J( w& i* S  @- lend
" |+ X1 Z4 U; P( d% e' U  y
最后运行结果没有任何问题。
那是哪里有问题呢?

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肯定是向量，做差之后还是向量( ...

, W0 H5 @" R- v" X* P; O! K谢谢大神指导
/ l- @; M  ^! Z9 C! B& G" r/ p