数值解析——编程作业.docxVIP

《数值分析》实验报告 第二章： 解线性方程组的直接方法 2、试用MATLAB软件编程实现追赶法求解三对角方程组的算法，并考虑梯形电阻电路问题，电路如下： 其中电路中的各个电流{，，…，}须满足下列线性方程组： 设，，运用求各段电路的电流量。 解： 上述方程组可用矩阵表示为： MatLab程序： %赋初值； a=[0 -2 -2 -2 -2 -2 -2 -2]; b=[2 5 5 5 5 5 5 5]; c=[-2 -2 -2 -2 -2 -2 -2]; d=[8.1481 0 0 0 0 0 0 0]; %三对角方程的追赶法 for i=2:8 %“追”的过程； a(i)=a(i)/b(i-1); b(i)=b(i)-c(i-1)*a(i); d(i)=d(i)-a(i)*d(i-1); end; d(8)=d(8)/b(8); %“赶”的过程； for i=7:-1:1 d(i)=(d(i)-c(i)*d(i+1))/b(i); end; x=d; x 程序运行结果： x = 8.1477 4.