24小时热门版块排行榜

返回列表

rhyme1990

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 1
在线: 28分钟
虫号: 1957878
注册: 2012-08-27
专业: 计算机科学的基础理论

[求助] 【求助】matlab程序问题

下面粗体代码，我想实现一个Y循环，但是单求一个Y时，第一次很快解出来，多次运行就变慢了···循环更是很难算出来···这是什么问题？是那个求积分吗？我是新手，求帮助啊。

clc,clear
pmm=0.7;poo=0.3;
syms t
Y=30;
v=6.8; %平均旅行速度
k=0.9;
l=225;
ct=1;
for i=6:18 %求i的阶乘
kj(ct)=1;
for j=2:i
      kj(ct)=kj(ct)*j;
end
ct=ct+1;
end
pz1=0;pz2=0;
for i=1:13
p1(i)=5.21^(i+5)/kj(i)*exp(-5.21);
p2(i)=9.26^(i+5)/kj(i)*exp(-9.26);
pz1=pz1+p1(i);
pz2=pz2+p2(i);
end
for i=1:13
pm(i)=p1(i)/pz1;
po(i)=p2(i)/pz2;
end
for i=1:13
p(i)=pmm*pm(i)+poo*po(i);
a(i)=1/(i+5);
sigma(i)=0.1-1/(3*(i+5));
fyi(i)=exp(-(t-a(i))^2/(2*sigma(i)^2))/(sqrt(2*pi)*sigma(i));
end
fy=0;
for i=1:13
fy=fy+p(i)*fyi(i);
end
js=0;
for Y=30:40;
js=js+1;
X=23.8*(Y+1);
for i=1:14
q(i)=int(fy,t,(2*i-1)/(2*(Y+1)),(2*i+1)/(2*(Y+1)));
end
q=subs(q);
temp=1;
for i=1:14
temp=temp*(1-q(i));
end
P(js)=1-temp;
Pi(js)=0;
for i=1:14
for j=1:i-1
      Pi(js)=Pi(js)+q(i)*q(j);  %用Pi表示P'
end
end
T(js)=l*log(1-k)/(v*(Y+1)*log(P(js)));
C(js)=X^2/(2.63*10^5);
end

回复此楼

» 猜你喜欢

依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有9人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有11人回复
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有7人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有8人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有10人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有9人回复
有院领导为了换新车，用横向课题经费买了俩车已经有10人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有6人回复

1楼 2012-08-27 10:55:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 rhyme1990 的主题更新

返回列表