Question

求助matlab中有限差分法中追赶法的实现

Answer 1

for ii=1:n
B=zeros(M-1,1);
A=zeros(M-2,1);
C=zeros(M-2,1);
S=zeros(M-1,1);
for ii=1:M-2
B(ii)=1+2*r;
A(ii)=-r;
C(ii)=-r;
S(ii)=u(ii+1,1);
end
B(M-1)=1+2*r;
S(M-1)=u(M,1);
u(1,2)=0;
u(M+1,2)=0;

S(1,1)=S(1,1)+r*u(1,2);
S(M-1,1)=S(M-1,1)+r*u(M+1,2);
S(1) = S(1)/B(1);
T = B(1);
k = 2;
while k~=M
B(k-1)=C(k-1)/T;
T=B(k)-A(k-1)/B(k-1);
S(k)=(S(k)-A(k-1)*S(k-1))/T;
k=k+1;
end
k=1;
while k~=M-1
S(M-1-k)=S(M-1-k)-B(M-1-k)*S(M-k);
k=k+1;
end
u(2:M,2)=S;
u(:,1)=u(:,2);
end
如题，我写的一段小程序，实现A*u(j-1,n)+B*u(j,n)+C*u(j+1,n)=u(j,n-1)差分方程的求解,其中u(j,n-1)记为S()， 这是程序中的一小部分，最后的输出结果就是u中的数据基本都是inf无穷大，我查了一下，A和C数组中的数据都正确，但是B数组中的数据变成
-0.056738
-Inf
-0.056738
-Inf
-0.056738
-Inf
1.128
我找了一下午 实在看不到错误在哪里了。。求达人指教指教！ 谢谢