24小时热门版块排行榜

返回列表

fanghd2005

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 722.3
散金: 20
红花: 1
帖子: 60
在线: 33小时
虫号: 3065504
注册: 2014-03-18
性别: GG
专业: 塑料

[求助] matlab有限差分法解Couette流出错已有2人参与

各位前辈，我用有限差分法解一维couette流时计算结果与教材中的结果不一样，
我用matlab写的程序，我对比了一下中间计算结果，与教材中的结果相同，矩阵与线性方程右端系数都相同，只有求解之后的结果不同。
请各位前辈帮我看看我的程序有什么问题，请指正。
%定义模型
Ylength=1;                %平板长度
detaY=0.05;
My=Ylength/detaY;
detat=12.5;                %时间步长
t=15000;                   %结晶计算时间
Nt=t/detat;
%定义初始条件
u=zeros(My+1,1);          %各点速度
u(My+1)=1;
%计算过程
for AC=1:Nt
%流动场,Crank-Nicolson算法
ReD=5000;
A=-detat/(2*ReD*detaY^2);       %差分系数
B=1+detat/(ReD*detaY^2);
for j=1:My-1                   %循环矩阵
      P(j,j)=B;
      if j>1
         P(j-1,j)=A;P(j,j-1)=A;
      end
end
for i=2:My                      %右端常数项
      f(i-1,AC)=(1-detat/(detaY^2*ReD))*u(i,AC)+detat/(2*detaY^2*ReD)*(u(i+1,AC)+u(i-1,AC));
end
u(2:My,AC+1)=P\f(:,AC); %此处出错
u(1,AC+1)=0;
u(My+1,AC+1)=1;
end

回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢

基元I理论下三大核心空间现象精准推导与细节解析已经有0人回复
基于基元 I 统一理论的反重力理论推导已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有58人回复
基于基元I统一理论的量子力学本源推导已经有1人回复
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费已经有20人回复
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）已经有0人回复
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）已经有0人回复

不会的太多

1楼 2016-08-26 16:30:36

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

FMStation

至尊木虫 (知名作家)

应助: 591 (博士)
贵宾: 0.03
金币: 18247.2
红花: 97
帖子: 8891
在线: 1487.9小时
虫号: 2400059
注册: 2013-04-04
专业: 计算机应用技术

【答案】应助回帖

感谢参与，应助指数 +1

u'

ans =

  Columns 1 through 13

   0    0    0    0    0    0    0    0    0    0    0    0    0

  Columns 14 through 21

   0    0    0    0    0    0    0    1