矩阵表达式

nkkopd

自定义的函数：
# |# z2 z2 Z9 p3 ^- rfunction f=largejfun(x)
f=(x-1./1:100).^2;
end
9 Y3 D0 T+ i% b/ s$ ^! z+ Y
调用函数的时候：
! d4 V7 P4 O6 l; ~clc
2 \+ ]8 L- R' vx0=10*ones(1,100);
: Z) `4 A" D) u; p% hoptions=optimset('display','iter','TolFun',1e-8);
3 B8 ]8 i( L7 T, z3 O[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)

: D7 L2 p3 f" k1 ^% ?0 q& L
运行结果出错，出错结果显示如下：
: d7 \2 H  m+ P& e( u5 Z& d% D错误使用 fminunc (line 348)
Supplied objective function must return a scalar value.
8 u! @1 p% Y3 W; ~$ D; \8 V
8 W" W0 I* a/ B( E  u5 r: J
2 `- B$ V. U8 e6 C9 [5 T/ a出错 fminu1 (line 4)
4 X# ?4 ~% c: f' m7 q% [[x,fval,exitflag,output]=fminunc(@largejfun,x0,options)


% a8 {- ]% L# [


" N/ T- t( A' B" J( H/ |

! n6 l8 i+ w, M$ E7 Q- E. ^2 _
  E+ N3 [+ i. K/ S+ h; s- ?& Q我将自定义函数改为如下：
function f=largejfun(x)
f=0;
for i=1:100
    f=f+(x(i)-1/i)^2;
+ P! i9 ]) Z3 g- Q: kend

最后运行结果没有任何问题。
那是哪里有问题呢?

kekek

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

Uifhjvv

楼上正解，楼主试试

小小鲁班

nkkopd

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

1 w# {8 T: n% ]4 d谢谢大神指导