实验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=sparse(n,n,-0.5,n,n); %定义(n,n) A1=A11+A12; %对角元素 A2=sparse(1:n-1,2:n,0.5,n,n); %输入A的上次对角元素 A3=sparse(2:n,1:n-1,0.5,n,n); %输入A的下次对角元素 A=A1+A2+A3; b1=0.5*ones(n,1); %b的除第一项元素 b2=sparse(1,1,18,n,1); %b的第一项元素 b=b1-b2; R=A\b 输出结果如下： R = ? 26.0000 17.0000 9.0000 2.0000 -4.0000 -9.0000 -13.0000 -16.0000 -18.0000 -19.0000 所以各阻值为 (R1，R2，…，R10)=(26，17，9，2，-4，-9，-13，-16，-18，-19) 总电阻R0（即输入等效电阻）为，又 得到 ? 5-6　有5个反应器连接如图5.6（见书），各个Q表示外部输入、输出及反应器间的流量（m3/min），各个c表示外部输入及反应器内某物质的浓度 (mg/m3)。假定反应器内的浓度是均匀的，利用质量守恒准则建立模型，求出各反应器内的浓度c1~c5，并讨论反应器j外部输入改变１个单位(mg/min)所引起的反应器i浓度的变化。 [问题分析、模型建立及求解] 　　当反应容器中反应物浓度稳定时，输入物质质量与输出物质质量平衡，即输入物质质量等于输出物质质量。由此分析，分别对1~5容器列质量平衡方程，以输入物质为“+”，输出物质为“－” ，代入数据整理为 由矩阵表示： 　， ? 编写MATLAB程序如下： A=[-6,0,1,0,0 3,-3,0,0,0 0,1,-9,0,0 0,1,8,-11,2 3,1,0,0,-4]; b=[-50,0,-160,0,0]; c=A\b %计算并输出c pause; dQc01=[1,0,0,0,0]; dQc02=[0,0,1,0,0]; c11=A\(b+dQc01); %01的输入减少1个单位 dc11=c-c2; c12=A\(b-dQc01); %01的输入增加1个单位 dc12=c-c1; c21=A\(b+dQc02); %03的输入减少1个单位 dc21=c-c4; c22=A\(b-dQc02) ; %03的输入增加1个单位 dc22=c-c3; [c,c11,dc11,c12,dc12,c21,dc21,c22,dc22]%外部输入改变引起的 %反应器浓度变化列表比较 MATLAB给出结果如下： c = ? 11.5094 11.5094 19.0566 16.9983 11.5094 所以改变前各反应器浓度为 (c1，c2，c3，c4，c5)=( 11.5094，11.5094，19.056