各位大侠，帮忙看看啊

ccchen · 发表于 2006-4-19 15:52

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

r=rand(22,1); xma=[2 1 0 2 2 3 3 3 2 12 4 3 2 3 3 3 3 3 1 1 1 1]; xmi=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; for i=1:22 x(i)=xmi(i)+xma(i)*unifrnd(0,1); syms t real; a(i)=gamma(r(i)/2); b(i)=t-xmi(i); f(i)=(b(i)^(r(i)/2-1)*exp(-b(i)/2))/(a(i)*2^(r(i)/2)); M=int(f(i),t,xmi(i),x(i)); m(i)=double(M(i)); end 为什么m（i）可以循环求取几个后，就出现问题？ 如果是程序有错误应该一次循环也不行啊？不懂啊。 请知道的大侠帮帮阿 该怎么改？

suffer · 发表于 2006-4-21 09:01

<DIV class=quote>以下是引用ccchen在2006-4-19 15:52:56的发言： r=rand(22,1); xma=[2 1 0 2 2 3 3 3 2 12 4 3 2 3 3 3 3 3 1 1 1 1]; xmi=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; for i=1:22 x(i)=xmi(i)+xma(i)*unifrnd(0,1); syms t real; a(i)=gamma(r(i)/2); b(i)=t-xmi(i); f(i)=(b(i)^(r(i)/2-1)*exp(-b(i)/2))/(a(i)*2^(r(i)/2)); M=int(f(i),t,xmi(i),x(i)); m(i)=double(M(i)); end 为什么m（i）可以循环求取几个后，就出现问题？ 如果是程序有错误应该一次循环也不行啊？不懂啊。 请知道的大侠帮帮阿 该怎么改？</DIV>
t是符号型变量，xmi是数值型变量，当然不能做运算了

账号		自动登录	找回密码
密码			我要加入

各位大侠，帮忙看看啊

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（ccchen）各位大侠，帮忙看看啊