应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab有限差分法解Couette流出错
31/1返回列表
查看: 981  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fanghd2005

铁虫 (小有名气)

[求助] 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
回复此楼

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

程序

» 猜你喜欢
不会的太多
1楼 2016-08-26 16:30:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +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
赞一下 回复此楼
2楼2016-08-26 17:28:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

512095518

专家顾问 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
我运行了你下你的代码，没发现报错哦。我用的是MATLAB R2016a。不知道是不是版本问题
赞一下 回复此楼
没有极限，只有超越！
3楼2016-08-27 10:15:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fanghd2005 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中科技大学071000，求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 生物学调剂 +4 Surekei 2026-03-21 4/200 2026-03-22 07:18 by ilovexiaobin
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭