用MATLAB编写程序：求解下列三角线性方程组。算法必须有效地利用系数矩阵的稀疏性。 4m1+m2=3 m1+4m2+m3=3 m2+4m3+m4=3 m3+4m4+m5=3 …… m48+4m49+m50=3 m49+4m50=3

该三角线性方程组可以用矩阵形式表示为Ax=b，其中

A=[4 1 0 0 ... 0; 1 4 1 0 ... 0; 0 1 4 1 ... 0; ... ; 0 0 0 ... 4 1; 0 0 0 ... 1 4];

b=[3;3;3;...;3;3];

由于系数矩阵A是一个三角矩阵，因此可以使用前向回代法求解该方程组。

步骤如下：

1. 初始化向量x=zeros(50,1)。

2. 对于i=1:50，计算：

x(i)=(b(i)-A(i,1:i-1)*x(1:i-1))/A(i,i);

3. 输出向量x即为方程组的解。

MATLAB代码如下：

A=zeros(50); A(1,1:2)=[4,1]; for i=2:49 A(i,i-1:i+1)=[1,4,1]; end A(50,49:50)=[4,1]; b=3*ones(50,1); x=zeros(50,1); for i=1:50 x(i)=(b(i)-A(i,1:i-1)*x(1:i-1))/A(i,i); end disp(x);