| 查看: 383 | 回复: 0 | ||
[求助]
【求助】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 |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有10人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
找到一些相关的精华帖子，希望有用哦~
Matlab中save问题
已经有5人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 关于Matlab作图的问题
已经有14人回复
- matlab调试问题
已经有3人回复
- matlab非线性参数拟合问题
已经有7人回复
- MATLAB 画图小问题
已经有11人回复
- 同样的MATLAB程序，为何在别人的电脑上运行不了
已经有21人回复
- matlab 程序问题
已经有4人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
