janezh
2009-12-24, 16:14
Warning: Explicit integral could not be found.
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit integral could not be found.
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit integral could not be found.
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit integral could not be found.
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit solution could not be found.
> In D:\matlab\toolbox\symbolic\solve.m at line 133
In D:\matlab\toolbox\symbolic\@sym\solve.m at line 49
In D:\matlab\work\NdTransferCoN.m at line 3
function f = test(B)
syms m n x y;
[m,n]=solve(40-30*normcdfa(m,100,50)-20+28*dblquad(normpdfa(x,100,50)*normpdfa(y,100,50),0,m,m+n-x,inf)+28*dblquad(-normpdfa(x,100,50)*normpdfa(y,100,50),m,m+n-y,0,n),...
40-30*normcdfa(n,100,50)-20+28*dblquad(-normpdfa(x,100,50)*normpdfa(y,100,50),0,m,n,m+n-x)+28*dblquad(normpdfa(x,100,50)*normpdfa(y,100,50),m+n-y,inf,0,n))
其中,
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);
function pdf1=normpdfa(t,u,sigma)
t=sym(t);
pdf2=exp(-(t-u)^2/(2*sigma^2));
pdf3=pdf2/(sigma*sqrt(2*pi));
cdf2=int(pdf3,t,0,inf);
pdf1=pdf3/cdf2;
B={[20],[0 500],[100 100],[50 50],[2 2];[40 40],[10 10],[0 0],[0],[12 12]}
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit integral could not be found.
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit integral could not be found.
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit integral could not be found.
> In D:\matlab\toolbox\symbolic\@sym\int.m at line 58
In D:\matlab\work\NdTransferCoN.m at line 3
Warning: Explicit solution could not be found.
> In D:\matlab\toolbox\symbolic\solve.m at line 133
In D:\matlab\toolbox\symbolic\@sym\solve.m at line 49
In D:\matlab\work\NdTransferCoN.m at line 3
function f = test(B)
syms m n x y;
[m,n]=solve(40-30*normcdfa(m,100,50)-20+28*dblquad(normpdfa(x,100,50)*normpdfa(y,100,50),0,m,m+n-x,inf)+28*dblquad(-normpdfa(x,100,50)*normpdfa(y,100,50),m,m+n-y,0,n),...
40-30*normcdfa(n,100,50)-20+28*dblquad(-normpdfa(x,100,50)*normpdfa(y,100,50),0,m,n,m+n-x)+28*dblquad(normpdfa(x,100,50)*normpdfa(y,100,50),m+n-y,inf,0,n))
其中,
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);
function pdf1=normpdfa(t,u,sigma)
t=sym(t);
pdf2=exp(-(t-u)^2/(2*sigma^2));
pdf3=pdf2/(sigma*sqrt(2*pi));
cdf2=int(pdf3,t,0,inf);
pdf1=pdf3/cdf2;
B={[20],[0 500],[100 100],[50 50],[2 2];[40 40],[10 10],[0 0],[0],[12 12]}