数学建模社区-数学中国
标题:
关于积分的一个问题
[打印本页]
作者:
在仰望星空
时间:
2012-10-6 20:26
标题:
关于积分的一个问题
clc;clear
$ L* \9 n* F5 b7 W% g3 g2 v
m=40;LSL=38.5;USL=41.5;
3 j, Q9 f/ ?7 a9 s9 X0 m
k1=0.6;k2=0.3;k3=1.0;sig0=0.5;
: F* g/ A* g9 W) c" w0 I
sig1=4.0;
- w- s7 X7 e: ~6 ?/ H
x=[m,sig0];
- X! t3 \8 Z" _7 ~/ _5 M2 f
6 y) `$ @3 Q! a R: ` L: k
6 z* d' p, n% B4 d
7 P% v7 {. R/ r# {$ i% X
EL=k1*quadl(@(y)myfun(x,y,m),LSL,m)+k2*quadl(@(y)myfun(x,y,m),m,USL);%顾客偏好引起的损失
& M4 D4 t4 L, |+ Z: ]) H. T' M
, B$ Q+ A b' ^8 D# Z7 {+ _/ f
" Y3 a1 [! a# L h
function g=myfun(y,x,m)
& I$ F% S# v0 x2 |$ Y/ C3 o9 R
mu=x(1);sig=x(2);
2 ]# Z2 n5 v! V. e! L3 O# c% h
f=(sqrt(2.*pi).*sig).^(-1).*exp((-(y-mu).^2)./(2.*sig^2));
6 Q* T& B Q5 z; ?8 u6 I g$ X
g=(y-m).^2.*f;
* z3 g6 l; P/ T1 p$ b
6 }( B* ]2 Y$ b: f4 M; P
; A% M6 O3 d/ Y S x* \2 G
??? Attempted to access y(13); index out of bounds because numel(y)=2.
0 C4 S* Y' E2 V
: `, V' P) U4 e1 h/ d+ x
Error in ==> quadl at 78
% @* ]# f$ @$ M O& ?
if ~isfinite(y(13))
# ]9 q' C; v1 R) x( l; l: d
' S5 J- D- [2 d5 j; m+ S+ @7 u
/ h3 ~; P8 N" h8 I1 ~! E
这个问题怎么解决啊 求助
未命名.bmp
(769.55 KB, 下载次数: 301)
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
这个发错了
