实验二源码.docx
- 文档编号:12862111
- 上传时间:2023-06-08
- 格式:DOCX
- 页数:12
- 大小:15.32KB
实验二源码.docx
《实验二源码.docx》由会员分享,可在线阅读,更多相关《实验二源码.docx(12页珍藏版)》请在冰点文库上搜索。
实验二源码
实验二:
第1小题源码如下:
clc
clear
A=[100;-1/1023/150-1/50;-1/15-1/5018/150];
B=[1000];
[n,m]=size(A);
fori=1:
n
A(i,i)=1/A(i,i);
forj=i+1:
n
A(i,j)=A(i,j)*A(i,i);
end
fork=i+1:
n
forj=i+1:
n
A(k,j)=A(k,j)-A(k,i)*A(i,j);
end
end
end
disp('矩阵A的因子表为:
');
disp(A)
fori=1:
n
B(i)=B(i)*A(i,i);
forj=i+1:
n
B(j)=B(j)-A(j,i)*B(i);
end
end
fori=n-1:
-1:
1
forj=i+1:
-1:
2
B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);
end
end
disp('在因子表的基础上求解线性方程组的解为:
x=');
disp(B)
第2小题源码如下:
2V
clc
clear
A=[100;-1/1023/150-1/50;-1/15-1/5018/150];
B=[200];
[n,m]=size(A);
fori=1:
n
A(i,i)=1/A(i,i);
forj=i+1:
n
A(i,j)=A(i,j)*A(i,i);
end
fork=i+1:
n
forj=i+1:
n
A(k,j)=A(k,j)-A(k,i)*A(i,j);
end
end
end
disp('矩阵A的因子表为:
');
disp(A)
fori=1:
n
B(i)=B(i)*A(i,i);
forj=i+1:
n
B(j)=B(j)-A(j,i)*B(i);
end
end
fori=n-1:
-1:
1
forj=i+1:
-1:
2
B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);
end
end
disp('在因子表的基础上求解线性方程组的解为:
x=');
disp(B)
4V
clc
clear
A=[100;-1/1023/150-1/50;-1/15-1/5018/150];
B=[400];
[n,m]=size(A);
fori=1:
n
A(i,i)=1/A(i,i);
forj=i+1:
n
A(i,j)=A(i,j)*A(i,i);
end
fork=i+1:
n
forj=i+1:
n
A(k,j)=A(k,j)-A(k,i)*A(i,j);
end
end
end
disp('矩阵A的因子表为:
');
disp(A)
fori=1:
n
B(i)=B(i)*A(i,i);
forj=i+1:
n
B(j)=B(j)-A(j,i)*B(i);
end
end
fori=n-1:
-1:
1
forj=i+1:
-1:
2
B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);
end
end
disp('在因子表的基础上求解线性方程组的解为:
x=');
disp(B)
6V
clc
clear
A=[100;-1/1023/150-1/50;-1/15-1/5018/150];
B=[600];
[n,m]=size(A);
fori=1:
n
A(i,i)=1/A(i,i);
forj=i+1:
n
A(i,j)=A(i,j)*A(i,i);
end
fork=i+1:
n
forj=i+1:
n
A(k,j)=A(k,j)-A(k,i)*A(i,j);
end
end
end
disp('矩阵A的因子表为:
');
disp(A)
fori=1:
n
B(i)=B(i)*A(i,i);
forj=i+1:
n
B(j)=B(j)-A(j,i)*B(i);
end
end
fori=n-1:
-1:
1
forj=i+1:
-1:
2
B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);
end
end
disp('在因子表的基础上求解线性方程组的解为:
x=');
disp(B)
8V
clc
clear
A=[100;-1/1023/150-1/50;-1/15-1/5018/150];
B=[800];
[n,m]=size(A);
fori=1:
n
A(i,i)=1/A(i,i);
forj=i+1:
n
A(i,j)=A(i,j)*A(i,i);
end
fork=i+1:
n
forj=i+1:
n
A(k,j)=A(k,j)-A(k,i)*A(i,j);
end
end
end
disp('矩阵A的因子表为:
');
disp(A)
fori=1:
n
B(i)=B(i)*A(i,i);
forj=i+1:
n
B(j)=B(j)-A(j,i)*B(i);
end
end
fori=n-1:
-1:
1
forj=i+1:
-1:
2
B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);
end
end
disp('在因子表的基础上求解线性方程组的解为:
x=');
disp(B)
第3小题例17.6源码:
clear
n=5
nl=5
B=[120.03i01.050;230.08+0.3i0.5i10;340.015i01.051;250.1+0.35i010;350.04+0.25i0.5i10];
X=[10;20;30;40;50];
Y=zeros(n);
fori=1:
n
ifX(i,2)~=0;
p=X(i,1);
Y(p,p)=1./X(i,2);
end
end
fori=1:
nl
ifB(i,6)==0
p=B(i,1);q=B(i,2);
else
p=B(i,2);q=B(i,1);
end
Y(p,q)=Y(p,q)-1./(B(i,3)*B(i,5));
Y(q,p)=Y(p,q);
Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)^2)+B(i,4)./2;
Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;
end
disp('导纳矩阵Y=');
disp(Y)
例17.7源码:
clc
clear
n=3;
nl=6;
B=[012i010;024i010;122i010;0320i010;238i010;135i010];
m=0;Z=zeros(n);
fork1=1:
nl
p=B(k1,1);q=B(k1,2);
ifB(k1,6)==0
l=1./B(k1,5);
elsel=B(k1,5);
end
ifp==0
ifq>m%追加接地树支
Z(q,q)=B(k1,3);m=m+1;
else%追加接地连支
fork=1:
m,
Z(k,m+1)=-Z(k,q);Z(m+1,k)=-Z(q,k);
end
Z(m+1,m+1)=Z(q,q)+B(k1,3);
forl1=1:
m
fork=1:
m
Z(l1,k)=Z(l1,k)-Z(l1,m+1)*Z(m+1,k)./Z(m+1,m+1);
end
Z(l1,m+1)=0;
end
fork=1:
m+1
Z(m+1,k)=0;
end
end
elseifq>m%追加不接地树支
fork=1:
m
Z(k,q)=Z(k,p)*l;Z(q,k)=Z(p,k)*l;
end
Z(q,q)=l^2*Z(p,p)+l^2*B(k1,3);
m=m+1;
else
fork=1:
m%追加不接地连支
Z(k,m+1)=l*Z(k,p)-Z(k,q);
Z(m+1,k)=l*Z(p,k)-Z(q,k);
end
Z(m+1,m+1)=l^2*Z(p,p)+Z(q,q)-2*l*Z(p,q)+l^2*B(k1,3);
forl1=1:
m
fork=1:
m
Z(l1,k)=Z(l1,k)-Z(l1,m+1)*Z(m+1,k)./Z(m+1,m+1);
end
Z(l1,m+1)=0;
end
fork=1:
m+1
Z(m+1,k)=0;
end
end
end
end
Z1=zeros(n);
fori=1:
n
forj=1:
n
Z1(i,j)=Z(i,j)
end
end
disp('阻抗矩阵Z1=');
disp(Z1)
第4小题源码如下:
n=input('请输入节点数:
n=');
nI=input('请输入支路数:
nI=');
B=input('请输入由支路参数形成的矩阵:
B=');
X=input('请输入由节点号及其对地阻抗形成的矩阵:
X=');
m=0;Y=zeros(n);Q=zeros(n);V=zeros(n)
fori=1:
n
ifX(i,2)~=0
p=X(i,1);
Y(p,p)=1./X(i,2);
end
end
fori=1:
nI
ifB(i,6)==0
p=B(i,1);q=B(i,2);
else
p=B(i,2);q=B(i,1);
end
Y(p,q)=Y(p,q)-1./(B(i,3)*B(i,5));
Y(q,p)=Y(p,q);
Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)^2)+B(i,4)./2;
Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;
end
disp('导纳矩阵:
Y=');
disp(Y)
D=input('请输入由于系统变更形成的支路参数矩阵:
D=');
V=input('请输入由于系统变更形成的对地阻抗的矩阵;V=');
ifD
(2)>n
ifV
(2)~=0
p=V
(1);
Q(p,p)=1./V
(2);
end
Y(n,nI)=Y(n,nI)+Q(p,p)+1./D(3)+D(4)./2;
Y(n+1,nI+1)=1./D(3)+D(4)./2;
Y(D
(1),nI+1)=-1./D(3);
Y(n+1,D
(1))=Y(D
(1),nI+1);
elseifD(5)~=0&D(5)~=1
Y(D
(1),D
(1))=Y(D
(1),D
(1))+1./(D(3)*(D(5)^2-B(i,5)^2));
Y(D
(2),D
(2))=Y(D
(2),D
(2));
Y(D
(1),D
(2))=Y(D
(1),D
(2))-1./(D(3)*(D(5)-B(i,5)));
Y(D
(2),D
(1))=Y(D
(1),D
(2));
elseD(5)==1
Y(D
(1),D
(1))=Y(D
(1),D
(1))+1./D(3)+D(4)./2;
Y(D
(2),D
(2))=Y(D
(2),D
(2))+1./D(3)+D(4)./2;
Y(D
(1),D
(2))=Y(D
(1),D
(2))-1./D(3);
Y(D
(2),D
(1))=Y(D
(1),D
(2));
end
end
disp('新形成的导纳矩阵:
Y=');
disp(Y)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 源码