数学建模社区-数学中国
标题:
关于积分的一个问题
[打印本页]
作者:
在仰望星空
时间:
2012-10-6 20:26
标题:
关于积分的一个问题
clc;clear
! o5 r- F" R; ?, L7 i
m=40;LSL=38.5;USL=41.5;
6 Y& _, Q/ |/ F" F" h& V
k1=0.6;k2=0.3;k3=1.0;sig0=0.5;
* D2 R$ @$ c/ o, w# _1 |
sig1=4.0;
q* D+ ^) ]" [* t
x=[m,sig0];
3 P+ K( j' Q' x. f2 s
) O0 m$ e w. i0 e! G5 l
) a- w2 W* K0 ]/ g2 q' e
& T7 s7 U+ {8 c2 v, w/ ` y
EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
6 q0 i+ W9 b R4 a8 L1 _6 j# y
* s) K" E' K7 ~+ }5 N4 I
& p8 v/ l& A: k' ~6 o7 E4 a5 r% O
function g=myfun(y,x,m)
4 N4 T' I9 j+ U# h4 ^( p( N7 |
mu=x(1);sig=x(2);
: i. W0 N' }2 P# Q' Y" o) D N
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
4 n* s( {# P% Y e2 b% f
g=(y-m).^2.*f;
; W/ k6 e$ b( c/ s0 W6 Y
+ |7 J" [; J! w
+ A1 F, u& u0 E& m. {! r
??? Attempted to access y(13); index out of bounds because numel(y)=2.
' ?8 C- B; T4 \. x& W9 Q5 R7 m9 V+ f
) c! V, ]# P+ C7 _8 [: J5 K- Q
Error in ==> quadl at 78
9 M' c; e5 b0 U) t$ c
if ~isfinite(y(13))
9 a b8 `$ i# q Q; i( @
- f1 D/ z# T6 x0 s/ ^
( u/ L- Y% i) C
这个问题怎么解决啊 求助
未命名.bmp
(769.55 KB, 下载次数: 327)
2012-10-6 20:24 上传
点击文件名下载附件
作者:
在仰望星空
时间:
2012-10-6 20:32
帖子怎么删除
这个发错了
作者:
hehainan86
时间:
2012-10-8 17:28
你的函数myfun（x,y,m）在前面的一个式子出现两次里跟后面的定义的不一样，x与y调换了，是故意为之？？
作者:
hehainan86
时间:
2012-10-8 17:31
错误提示的意思是：你访问的y已经超出了边界，y的维数为2。
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
这个发错了
