矩阵表达式

nkkopd

自定义的函数：
function f=largejfun(x)
+ p4 ~5 F2 ^& Q. l: m8 z# ]9 vf=(x-1./1:100).^2;
. b* b- }1 ^' @# R) qend

调用函数的时候：
clc
x0=10*ones(1,100);
0 J9 V2 O7 U9 v/ x) Moptions=optimset('display','iter','TolFun',1e-8);
& c1 j- x2 U& U" z[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)

( r7 {4 @7 h5 U7 p* s' ~! I0 H& j
运行结果出错，出错结果显示如下：
错误使用 fminunc (line 348)
Supplied objective function must return a scalar value.
- k! I' J; ~5 [5 \; h

+ j% ]+ J: e( @. E2 X出错 fminu1 (line 4)
[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)








  T4 p& c$ U$ H* }2 U$ S/ h我将自定义函数改为如下：
function f=largejfun(x)
/ }  j7 {4 E' X3 df=0;
" `5 X3 i. q9 T& jfor i=1:100
    f=f+(x(i)-1/i)^2;
end

最后运行结果没有任何问题。
那是哪里有问题呢?
  L; o: Q3 a. [: n/ s2 u

kekek

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

Uifhjvv

楼上正解，楼主试试

小小鲁班

nkkopd

kekek 发表于 2020-8-17 10:48
: A+ o# ]! ~+ a3 c; [' K很明显啊，f=(x-1./1:100).^2这种写法，甭管你的x是标量还是向量，1./1:100肯定是向量，做差之后还是向量( ...

谢谢大神指导
) w1 i$ P* T) A3 ?& J/ y" h