数学建模社区-数学中国

标题: 积分运算的检验出错了 求助?? [打印本页]

作者: 迷走星球    时间: 2009-7-23 08:36
标题: 积分运算的检验出错了 求助??
syms m n i j t
R=200;J=10 ;D=3000;K=1000;u=600;o=196.6292;m=20;n=18;
%f=1/(2*pi)^(1/2)/o*exp(-(x-u)^2/2/o^2)
%C1=(J*m+K)*(1-int(f,x,0,m*n))
%C2=symsum((J*j+D)*int(f,x,(j-1)*n,j*n)+F*symsum(i*int(f,x,(j-1)*n+i-1,(j-1)*n+i),i,1,n),j,1,n)
C1=(J*m+K)*(1-(F(m*n)-F(0)))
C2=symsum(((J*j+D)*(F(j*n)-F((j-1)*n))+R*symsum(i*(F((j-1)*n+i)-F((j-1)*n+i-1)),i,1,n)),j,1,n)
T=int(1-(F(t)-F(0)),t,0,m*n)
y=(C1+C2)/T;
M文件:
function f=F(x)
syms w
p=0.000087719;u=600;o=196.6292;
f1=1/(2*pi)^(1/2)/o*exp(-(w-u)^2/2/o^2);
f=1-(1-int(f1,w,-inf,x))*(1-p)^x;
作者: madio    时间: 2009-8-7 09:02
我用的matlab7.1运算没有报错,只有警告,应该是没有明确的积分形式,结果是
Warning: Explicit integral could not be found.
> In sym.int at 58

T =

int((1-4678340085518391/1208925819614629174706176*pi^(1/2)*5313808070744117^(1/2)*erf(262144/5313808070744117*5313808070744117^(1/2)*t-157286400/5313808070744117*5313808070744117^(1/2))-4678340085518391/1208925819614629174706176*pi^(1/2)*5313808070744117^(1/2))*(9006409152229565/9007199254740992)^t-4678340085518391/1208925819614629174706176*erf(157286400/5313808070744117*5313808070744117^(1/2))*pi^(1/2)*5313808070744117^(1/2)+4678340085518391/1208925819614629174706176*pi^(1/2)*5313808070744117^(1/2),t = 0 .. 360)




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5