数学建模社区-数学中国

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

---

作者: 在仰望星空    时间: 2012-10-6 20:26
标题: 关于积分的一个问题
clc;clear
m=40;LSL=38.5;USL=41.5;
k1=0.6;k2=0.3;k3=1.0;sig0=0.5;
sig1=4.0;
" P, r& T. o; J* O- O  ~( b/ b" C7 Bx=[m,sig0];



EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失

+ s  c4 ], Q- P' ~$ w0 A" E% O
function g=myfun(y,x,m)
' Q3 s& h/ |* f+ i, p& ^7 Lmu=x(1);sig=x(2);
1 S5 ^8 G6 e2 z, qf=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
g=(y-m).^2.*f;
$ @6 U5 T) K9 B( C  d" Q
6 o2 z, u! @1 w. e8 W6 O
) i* P& D/ G( }: m??? Attempted to access y(13); index out of bounds because numel(y)=2.

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


- f  s8 T5 P# Z0 O这个问题怎么解决啊  求助

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

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