Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
 
 
主题工具 显示模式
旧 2009-12-24, 10:06   #1
janezh
初级会员
 
注册日期: 2009-12-17
帖子: 2
声望力: 0
janezh 正向着好的方向发展
默认 为什么说我的数组超出范围啊?!!请大侠指点啊

function NdTransferCoN(B)
syms x y m n u v;
[m,n]=solve(B{2,1}(1)+B{2,3}(1)-(B{2,1}(1)+B{2,3}(1)-B{2,2}(1))*normcdfa(m,B{1,3}(1),B{1,4}(1))-B{1,1}(1)...
+(B{2,1}(2)+B{2,3}(2)-B{2,2}(1)-B{2,4}(1))*(int(int(normpdfa(x,B{1,3}(1),B{1,4}(1))*normpdfa(y,B{1,3}(2),B{1,4}(2)),y,m+n-x,inf),x,0,m))...
+(B{2,1}(1)+B{2,3}(1)-B{2,2}(2)-B{2,4}(2))*(int(int(-normpdfa(x,B{1,3}(1),B{1,4}(1))*normpdfa(y,B{1,3}(2),B{1,4}(2)),x,m,m+n-y),y,0,n)),...
B{2,1}(2)+B{2,3}(2)-(B{2,1}(2)+B{2,3}(2)-B{2,2}(2))*normcdfa(n,B{1,3}(2),B{1,4}(2))-B{1,1}(1)...
+(B{2,1}(2)+B{2,3}(2)-B{2,2}(1)-B{2,4}(1))*(int(int(-normpdfa(x,B{1,3}(1),B{1,4}(1))*normpdfa(y,B{1,3}(2),B{1,4}(2)),y,n,m+n-x),x,0,m))...
+(B{2,1}(1)+B{2,3}(1)-B{2,2}(2)-B{2,4}(2))*(int(int(normpdfa(x,B{1,3}(1),B{1,4}(1))*normpdfa(y,B{1,3}(2),B{1,4}(2)),x,m+n-y,inf),y,0,n)));
m=double(m)
m=double(n)



其中,
function pdf1=normpdfa(x,u,sigma)
syms t
x=sym(x);
pdf2=exp(-(t-u)^2/(2*sigma^2));
pdf3=pdf2/(sigma*sqrt(2*pi));
cdf2=int(pdf3,t,0,inf);
pdf1=pdf3/cdf2;

function cdf1=normcdfa(x,u,sigma)
syms t
x=sym(x);
pdf1=exp(-(t-u)^2/(2*sigma^2));
pdf2=pdf1/(sigma*sqrt(2*pi));
cdf2=int(pdf2,t,0,inf);
pdf3=pdf2/cdf2;
cdf1=int(pdf3,t,0,x);



B={[20],[0 500],[100 100],[50 50],[2 2];[40 40],[10 10],[0 0],[0],[12 12]}
janezh 当前离线   回复时引用此帖
 


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 04:14


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.