关于积分的一个问题

在仰望星空

clc;clear
. T1 r: A2 v/ A; b2 b# Um=40;LSL=38.5;USL=41.5;
2 X3 L# d, s% \5 T7 kk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
, }$ U( P0 H+ ^( ]: {! b, {sig1=4.0;
! L6 A1 X( O5 Z4 s7 ex=[m,sig0];


  p# T# ~0 Z; a9 B9 v
EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失


, l7 t( q6 q0 S0 ~function g=myfun(y,x,m)
mu=x(1);sig=x(2);
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
g=(y-m).^2.*f;

; y; L9 z+ N, F% W
; J- q& _1 c2 L- |# w6 M??? Attempted to access y(13); index out of bounds because numel(y)=2.

Error in ==> quadl at 78
if ~isfinite(y(13))


" w* z, g! O) I8 b. n& G% \这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

你的函数myfun（x,y,m）在前面的一个式子出现两次里跟后面的定义的不一样，x与y调换了，是故意为之？？

hehainan86

错误提示的意思是：你访问的y已经超出了边界，y的维数为2。