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

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有5人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复

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

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

智能机器人

Robot (super robot)

我们都爱小木虫

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

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

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

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

返回列表