| 查看: 237 | 回复: 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 %%求解语句我删掉了···谢谢各位了···· |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有35人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有15人回复
-
基金中了
已经有8人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
求教！哪位大神帮我看看狄拉克方程里的一个问题啊！！
已经有11人回复
- 想不通找不到原因，到底怎么回事，心烦，rna提取问题
已经有15人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- MS 6.0计算弹性弹性常数尽然没有结果文件，大神们看看这是什么问题
已经有4人回复
- 马上要开始做CTAB提乳酸菌DNA了，各位大神帮忙看看实验步骤有木有啥致命问题额
已经有5人回复
- 求助各位大神关于含组氨酸标签蛋白纯化的问题
已经有8人回复
- 各位物理大神们讨论一下滑轮问题吧 呵呵
已经有15人回复
- 各位大神求教，有关沸石分子筛的粒径的问题。
已经有8人回复
- HKU申请提交RP问题，求各位大神指教，不胜感激！
已经有10人回复
- 关于投稿时作者简介的问题，请教各位大神。
已经有5人回复
- 【求助】关于考研的几个问题，请各位大神帮忙解答一下
已经有12人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50