习题4-12的Matlab程序
%代数方程形式AiTi=CiTi+1+BiTi-1+Di
mdim=10;%计算的节点数
x=linspace(1,3,mdim);%生成A、C、B、T数据的基数；
A=cos(x);%TDMA的主对角元素
B=sin(x);%TDMA的下对角线元素
C=cos(x)+exp(x); %TDMA的上对角线元素
T=exp(x).*cos(x); %温度数据
%由A、B、C构成TDMA
coematrix=eye(mdim,mdim);
for n=1:mdim
    coematrix(n,n)=A(1,n);
    if n>=2
        coematrix(n,n-1)=-1*B(1,n);
    end
    if n<mdim
        coematrix(n,n+1)=-1*C(1,n);
    end
end
%计算D矢量
D=(coematrix*T')';
%由已知的A、B、C、D用TDMA方法求解T
 

1.下载地址 205.30 KB (需下载币0个)