数学建模社区-数学中国»论坛 › 【数模论坛】资源分类 › 数模问题互助 › 积分运算的检验出错了求助？？

查看: 2767|回复: 1

积分运算的检验出错了求助？？

[复制链接]

字体大小: 正常放大

迷走星球

4 主题	4 听众	62 积分

升级 60%

该用户从未签到

群组: 数学建模

群组: 数学趣味、游戏、IQ等

群组: 电子科大数学建模交流群

群组: 飞虹数模

群组: LINGO

电梯直达

1^#

发表于 2009-7-23 08:36 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

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;

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

madio

3万主题	1311 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

2^#

发表于 2009-8-7 09:02 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

我用的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)