24小时热门版块排行榜    

查看: 3599  |  回复: 0

半导体jlu

新虫 (初入文坛)

[交流] 隐式差分的编程问题

最近刚学有限差分,遇到了这样的一个问题,在附图里面。这是我编程的程序。
clc
clear all
k=1;
nx=6
dx=1/(nx-1);
dt=0.02;
n=100;
s=k*dt/(dx^2);

A = sparse(nx,nx);
for i=2:nx-1
A(i,i-1) = -s;
A(i,i ) = (1+2*s);
A(i,i+1) = -s;
end

A(1 ,1 ) = 1;
A(nx,nx) = 1;

rhs = zeros(nx,1);
rhs(2:nx-1) = Told(2:nx-1);
rhs(1) = 0;
rhs(nx) = 100;

但是不知道时间上的T(n+1)和T(n)该如何区分和表示,请教诸位大侠,该如何处理呢?

隐式差分的编程问题
7a899e510fb30f243688b745cd95d143ad4b0361.jpg
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 半导体jlu 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见