数学建模社区-数学中国
标题:
关于积分的一个问题
[打印本页]
作者:
在仰望星空
时间:
2012-10-6 20:26
标题:
关于积分的一个问题
clc;clear
5 B8 x* K0 ~) c# h, l3 x& @: `
m=40;LSL=38.5;USL=41.5;
7 Z5 k7 j3 {. M K) ?, R2 _" R
k1=0.6;k2=0.3;k3=1.0;sig0=0.5;
3 I# ^. e# H; C5 k/ B
sig1=4.0;
% y' Q: i9 N f+ d9 X! }
x=[m,sig0];
6 R6 |- C9 @! }
! y' F! `! Z. I2 A' }
2 z8 \- @" U( K9 s- P8 G( n7 c
) s( E) H& b$ I+ J4 v6 Y2 U
EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
9 y( Z5 o) k9 C! q0 Z, t2 x
, T5 l3 x z- V% h5 d& P" T
' L2 F* g1 z1 B6 j r
function g=myfun(y,x,m)
; h! v! `: N2 J
mu=x(1);sig=x(2);
5 N1 V& X9 j0 R, @+ c
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
2 {$ j5 I4 P9 r5 n2 y" P9 {
g=(y-m).^2.*f;
# Q( z6 ]) K- R" Z: M( }, A* o5 g
, p0 |2 W6 p# ?5 e
9 @ i/ N3 \7 k: n8 ~: ^# k- U0 R
??? Attempted to access y(13); index out of bounds because numel(y)=2.
# ~* v/ ?8 i8 H2 y% s% P7 ~4 ?
* l. y8 u. w3 s" [) g( h( W
Error in ==> quadl at 78
! F7 P0 p8 E( b2 c
if ~isfinite(y(13))
9 z6 k8 H/ s% U+ R
' m- u% Y$ e3 w& n
6 |: g7 `2 V$ G, o, r2 t; f
这个问题怎么解决啊 求助
未命名.bmp
(769.55 KB, 下载次数: 282)
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
这个发错了
