数学建模社区-数学中国

标题: 看看下面程序该如何改正？？谢谢 [打印本页]

作者: 迷走星球 时间: 2009-7-23 10:26
概率密度函数的积分问题，得到的函数好像不能积分不知道怎么改请指教！！！
syms m n i j t
R=200;J=10 ;D=3000;K=1000;u=600;o=196.6292;m=20;n=18;
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;

作者: danli 时间: 2009-7-27 21:29
你程序里面的m,n,j都给出来了，就不要再定义了。字符好像不能积分，我上回算过出错。字符变量算了之后结果还是字符变量。我还不知道怎么转换。

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)