矩阵表达式

nkkopd

自定义的函数：
  J) y( ^5 g  m5 Y" z7 ~function f=largejfun(x)
f=(x-1./1:100).^2;
end
# u; Z: r- C1 B" v5 I- |
调用函数的时候：
clc
, z% N  j. n+ ~( t& z1 o) ?8 lx0=10*ones(1,100);
options=optimset('display','iter','TolFun',1e-8);
[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)

4 @: C5 X3 i* m0 h: K. ^9 K
2 o: ~4 W+ N  l3 t+ F+ x运行结果出错，出错结果显示如下：
) n" M0 i3 u2 Y2 `+ {错误使用 fminunc (line 348)
, ?! f8 n, I; q2 U( E3 n) Q5 k9 OSupplied objective function must return a scalar value.


出错 fminu1 (line 4)
) Q4 n9 z& \5 n1 R: ]' g) B$ a[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)



$ V# T5 }1 J- x% d3 ~
0 b" Y) g. L( V+ X
, q- r: F1 M: L' j# o: p
9 w4 i* Z( y/ k/ U) l
6 P% r% M/ D# L( @4 R' a' j
4 {  X! y4 Q9 P  }我将自定义函数改为如下：
function f=largejfun(x)
f=0;
5 ~& Q# \6 m- ]for i=1:100
4 Z' a! g9 i$ I# O( z% ]    f=f+(x(i)-1/i)^2;
end

最后运行结果没有任何问题。
( k% T- _8 h" c那是哪里有问题呢?
, c2 T. D1 ~  H1 F2 S

kekek

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

Uifhjvv

楼上正解，楼主试试

小小鲁班

nkkopd

kekek 发表于 2020-8-17 10:48
+ B  q* E( z" `8 n# |( t( T很明显啊，f=(x-1./1:100).^2这种写法，甭管你的x是标量还是向量，1./1:100肯定是向量，做差之后还是向量( ...

谢谢大神指导