| 查看: 203 | 回复: 0 | ||
[求助]
各位大神看看程序···结果有点问题，找不到原因
|
|
在调试程序时··我分布执行第二个for循环··显示每步的p值··理论上每循环一次r，显示的应该是该r层所有空间点的值·（应该是115个），循环完之后应该是显示200X115个值 但是奇怪的是我单步执行时发现，没循环一次r···就会把所有的值列遍而且也不对····请问下，这是怎么回事··是我对p（x,y）这类的语法理解错误了吗····我是将求解求解域离散成时间层和空间层的网格····程序主要是逐个时间层求解各个值 u=zeros(200,115);p=zeros(200,115); deltat=0.01;%%时间步长 deltax=6;%%57节车，每辆车分成两段，过115个节点，列车管总长度除以节点数，得到空间步长 p(1,1:115)=6;%%所有车，充气前的空气压力初值 u(1,1:115)=0;%%所有车，充气前的空气速度 for y=1:60 p(y,1)=6+(8-6)*(1-exp(-3*deltat*y));%%假定往第一辆车充气，空气压力由指数函数决定 end u(:,115)=0;%%假定终端的速度给定 D=31.75/1000;%%管道直径 Q=0.001;%%管道泄漏量，单位是pa/s，空气压强的单位都是pa c=deltat/deltax; f=0.03;%%管道摩擦系数 for r=1:199 %%在时间层上进行逐层求解 p(r+1,2:115)=p(r,2:115).*(1-c*(u(r,2:115)-u(r,1:114)))-u(r,2:115).*(p(r,2:115)-p(r+1,1:114)).*c-Q*deltat%%问题就在此句时发现的 u(r+1,114:-1:1)=u(r,114:-1:1).*(1-c.*(u(r+1,115:-1:2)-u(r,114:-1:1)))-c./p(r,114:-1:1).*(p(r,114:-1:1)-p(r,115:-1:2))-2*f*abs(u(r,114:-1:1)).*u(r,114:-1:1)/D; end %%求解语句我删掉了···谢谢各位了···· |
回复此楼» 猜你喜欢
一志愿郑州大学材料与化工085600，求调剂
已经有16人回复
-
一志愿北京科技大学材料工程085601，求调剂
已经有14人回复
-
315求调剂
已经有12人回复
-
五邑大学土木工程招调剂生2026
已经有4人回复
-
334求调剂
已经有7人回复
-
求调剂
已经有8人回复
-
285求调剂
已经有8人回复
-
材料295
已经有7人回复
-
295求调剂
已经有4人回复
-
322求调剂
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
求教！哪位大神帮我看看狄拉克方程里的一个问题啊！！
已经有11人回复
- 想不通找不到原因，到底怎么回事，心烦，rna提取问题
已经有15人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- MS 6.0计算弹性弹性常数尽然没有结果文件，大神们看看这是什么问题
已经有4人回复
- 马上要开始做CTAB提乳酸菌DNA了，各位大神帮忙看看实验步骤有木有啥致命问题额
已经有5人回复
- 求助各位大神关于含组氨酸标签蛋白纯化的问题
已经有8人回复
- 各位物理大神们讨论一下滑轮问题吧 呵呵
已经有15人回复
- 各位大神求教，有关沸石分子筛的粒径的问题。
已经有8人回复
- HKU申请提交RP问题，求各位大神指教，不胜感激！
已经有10人回复
- 关于投稿时作者简介的问题，请教各位大神。
已经有5人回复
- 【求助】关于考研的几个问题，请各位大神帮忙解答一下
已经有12人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
