实验6 线性代数方程组的数值解法.docVIP

实验6 线性代数方程组的数值解法 [实验目的] 1.???? 学会用MATLAB软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作初步分析； 2.???? 通过实例学习用线性代数方程组解决简化的实际问题。 [实验内容] 5-5 输电网络：一种大型输电网络可简化为图5.5（见书）所示电路，其中R1，R2，…，Rn表示负载电阻，r1，r2，…，rn表示线路内阻，I1，I2，…，In表示负载上的电流。设电源电压为V。 (1)列出求各负载电阻R1，R2，…，Rn的方程； (2)设I1=I2=…=In=I，r1=r2=…=rn=r，在r=1，I=0.5，V=18，n=10的情况下求R1，R2，…，Rn及总电阻R0。 [问题分析、模型建立及求解] (1) 设电源负极为电势为0，电阻R1上对应节点电压为V1，对于任意节点，根据KCL定律列出方程： 而，可得： k=2，3，……，n-1； k=1时 ，为与上式形式一致，化为 k=m（）时 k=n时 设以上方程组的矩阵形式为： 则 ? (2)代入参数： ，，V=18，n=10， ? 在命令窗口输入MATLAB程序如下： clear all; n=10; %由题目要求设定 A11=sparse(1:n-1,1:n-1,-1,n,n); %定义A的对角元素，除（n,n) A12