矩阵表达式

nkkopd

自定义的函数：
+ F9 M; C" P* ^8 K# X) c4 Jfunction f=largejfun(x)
f=(x-1./1:100).^2;
0 ?! q9 T6 u5 o# a  Z! d3 Gend

调用函数的时候：
3 A) B8 ?- J  C6 g2 J! bclc
1 E6 s8 e, t" m0 u) J9 b4 ex0=10*ones(1,100);
8 L( O& ?: m: _% ]0 loptions=optimset('display','iter','TolFun',1e-8);
[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)
/ `) R9 s* R3 j7 j7 ?/ f; r, b* U

运行结果出错，出错结果显示如下：
) z- w, l7 `3 a- b$ ^3 N错误使用 fminunc (line 348)
Supplied objective function must return a scalar value.


' |8 E, _( }- o2 y$ r! s( E7 H3 @出错 fminu1 (line 4)
; O; m: [2 c1 H* b! P6 F, ^$ h[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)
8 v8 v7 h  |& g3 S
) I# y- @* D- E* y1 C) c6 {3 M




$ ^* P4 ~7 |+ H" D
4 T- B/ r$ S3 t; I# K. k+ f
' u% \. B( [/ ], A2 \+ M我将自定义函数改为如下：
function f=largejfun(x)
f=0;
* Z# k3 A- z, V: i9 Ofor i=1:100
    f=f+(x(i)-1/i)^2;
end

6 F! }" o. f' \' k; ?# X最后运行结果没有任何问题。
9 r4 j" L; F7 Y' p0 S4 R5 r6 @) z( l那是哪里有问题呢?
/ \* u! n, D- w8 O9 ~

kekek

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

Uifhjvv

楼上正解，楼主试试

小小鲁班

nkkopd

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

' O" \% J2 G/ V$ Y  L谢谢大神指导