关于积分的一个问题

在仰望星空

clc;clear
m=40;LSL=38.5;USL=41.5;
+ K% r+ D4 Z9 ~8 rk1=0.6;k2=0.3;k3=1.0;sig0=0.5;
6 R1 b/ p( {3 ^+ e" asig1=4.0;
: B3 b) v- s* w9 \. p0 Hx=[m,sig0];


* ?7 {: w( m0 s$ |/ \( X7 P
2 s; N' g. _# R* a8 c+ K& v2 `( REL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
& _5 y) E) U- n* Q% j0 @/ l0 f

2 E3 q# o/ a( P* \8 mfunction g=myfun(y,x,m)
, `5 t+ b5 t" C" ?mu=x(1);sig=x(2);
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
g=(y-m).^2.*f;
7 ]" m2 g6 u" A5 W3 T

2 M0 `+ I+ K! s4 U, r3 v  r??? Attempted to access y(13); index out of bounds because numel(y)=2.

; E/ x7 }  N5 z9 Y6 O; i) aError in ==> quadl at 78
if ~isfinite(y(13))

* z3 y3 g. O! u: F
- P/ ^1 K: ], }这个问题怎么解决啊  求助

在仰望星空

帖子怎么删除  这个发错了

hehainan86

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

hehainan86

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