数学建模社区-数学中国

标题: 关于积分的一个问题 [打印本页]

作者: 在仰望星空    时间: 2012-10-6 20:26
标题: 关于积分的一个问题
clc;clear
  _. o. K  A; e/ }7 O3 |m=40;LSL=38.5;USL=41.5;) F% e( T5 d$ x' e* @; Z4 v8 J
k1=0.6;k2=0.3;k3=1.0;sig0=0.5;; H5 L" b! O+ ^& Z; O
sig1=4.0;
2 V. U* k" V( J: o" ~. Wx=[m,sig0];4 F% t9 ]1 r# G4 N) w' T4 \
$ ~# H- I/ W2 j  ~

9 m2 o0 n3 _* R) P. ~) R1 G$ P$ e$ E, ^" V3 D4 `
EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
* @9 F) |$ p5 F* T. ~6 L* i$ s$ @8 o/ i) m7 P1 [5 K" N
- i$ w7 J1 F! t5 y, G$ Y4 W' I
function g=myfun(y,x,m)1 g$ J. F& d5 T& e4 m: ]; ?) e
mu=x(1);sig=x(2);
  S; t( R$ [5 g; A0 N9 Z: z7 }9 j5 xf=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
* C: B* u9 E1 R( F* t& A3 G/ xg=(y-m).^2.*f;
- o" a) r1 R; j0 M$ r. q; u/ L' |( @8 c

9 r. g7 _$ E, k- P' m??? Attempted to access y(13); index out of bounds because numel(y)=2.
% M" S  ?5 h( N, z2 I
: J+ N0 b2 \% G5 _Error in ==> quadl at 78
& z& e5 ]% I/ g; @( H- x/ M8 D+ Lif ~isfinite(y(13))
! S  u2 r! c2 D2 ^$ z2 ]/ p! u$ p ) k$ Q! i* Z8 \, ]
! [8 T/ Y+ R; m0 r7 f% u; S. v
这个问题怎么解决啊  求助

未命名.bmp (769.55 KB, 下载次数: 259)

未命名.bmp


作者: 在仰望星空    时间: 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