- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验6 线性代数方程组的数值解法
实验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外部输入改变1个单位(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
文档评论(0)