电力系统分析1.doc

用程序实现节点导纳矩阵的形成：主程序：%节点导纳矩阵的形成%网络的节点数为nn=input(网络的节点数);A=input(节点相关性矩阵);%1表示两节点相连，0表示两节点不连Y=zeros(n,n);for i=n:-1:1 for j=i-1:-1:1 if (A(i,j) i j Y(i,j)=input(i,j支路之间的导纳值); Y(j,i)=Y(i,j); end end Y(i,i)=sum(Y(i,:);endY%对节点导纳矩阵的修改%0表示不修改；%1表示给原有网络增加节点，且节点处的阻抗为Zij;%2表示在节点i,j之间增加一条阻抗为Zij支路；%3表示在原有网络i,j之间切除一条阻抗为Zij支路；%4表示原有网络i,j之间支路阻抗由Zij变成Z1ijm=1;for m=1:100m=input(导纳矩阵修改的情况分类);switch(m) case(0) break case(1) i=input(在第i点增加阻抗); Zij=input(增加节点的阻抗); Y(n+1,:)=zeros(1,n); Y(:,n+1)=zeros(n+1,1); Y(n+1,n+1)= Y(i,i); Y(i,n+1)=-1/Zij; Y(n+1,i)= Y(i,n+1); Y(i,i)= Y(i,i)+ Y(i,i); Y case(2) i=input(请输入节点i); j=input(请输入节点j); Zij=input(请输入在节点i,j间增加的阻抗); Y(i,i)= Y(i,i)+1/Zij; Y(j,j)=Y(j,j)+1/Zij; Y(i,j)=Y(j,j)-1/Zij; Y(j,i)=Y(i,j); Y case(3) i=input(请输入节点i); j=input(请输入节点j); Zij=input(请输入在节点i,j间增加的阻抗); Y(i,i)= Y(i,i)-1/Zij; Y(j,j)=Y(j,j)-1/Zij; Y(i,j)=Y(j,j)+1/Zij; Y(j,i)=Y(i,j); Y case(4) i=input(请输入节点i); j=input(请输入节点j); Z1ij=input(请输入在节点i,j间变换的阻抗); Y(i,i)= Y(i,i)-1/Y(i,j)+1/Z1ij; Y(j,j)=Y(j,j)-1/Y(i,j)+1/Z1ij; Y(i,j)=Y(j,j)+1/Y(i,j)-1/Z1ij; Y(j,i)=Y(i,j); Y otherwise disp(输入数据有误，请重新输入数据);endend运行结果：网络的节点数3节点相关性矩阵1,1,1;1,1,0;1,0,1i = 3j = 1i,j支路之间的导纳值5i = 2j = 1i,j支路之间的导纳值2Y = 7 2 5 2 2 0 5 0 5导纳矩阵修改的情况分类1在第i点增加阻抗2增加节点的阻抗4Y = 7.0000 2.0000 5.0000 0 2.0000 4.0000 0 -0.2500 5.0000 0 5.0000 0 0 -0.2500 0 2.0000导纳矩阵修改的情况分类2请输入节点i3请输入节点j1请输入在节点i,j间增加的阻抗8Y = 7.1250 2.0000 7.0000 0 2.0000 4.0000 0 -0.2500 7.0000 0 5.1250 0 0 -0.2500 0 2.0000导纳矩阵修改的情况分类3请输入节点i2请输入节点j3请输入在节点i,j间增加的阻抗1Y = 7.1250 2.0000 7.0000 0 2.0000 3.0000 5.1250 -0.2500 7.0000 5.1250 4.1250 0 0 -0.2500 0 2.0000导纳矩阵修改的情况分类4请输入节点i3请输入节点j3请输入在节点i,j间变换的阻抗6Y = 7.1250 2.0000 7.0000 0 2.0000 3.0000 5.1250 -0.2500 7.0000 5.1250 4.0542 0 0 -0.2500 0 2.0000导纳矩阵修改的情况分类7输入数据有误，请重新输入数据导纳矩阵修改的情况分类0三角分解法求解线性方程组：主程序：%用三角分解法解线性组方程AX=bn=input(系数矩阵的阶数);A=input(线性方程组的系数矩阵);b=input(线性方程组的常数项);L,U,m=zhjLU(A,n);y=zeros(1,4);if (m=1) x=b/L/U; x end调用程序：function L,U,m=zhjLU(A,n) RA=rank(A); if RA=ndisp(请注意：因为A的n阶行列式hl等于零，所以A不能进行LU分解.A的秩RA如下:), RA,hl=det(A);m=0;returnendif RA=n for p=1:nh(p)=det(A(1:p, 1:p);endhl=h(1:n);for i=1:nif h(1,i)=0disp(请注意：因为A的r阶主子式等于零，所以A不能进行LU分解.A的秩RA和各阶顺序主子式值hl依次如下：), hl;RAm=0;returnendend if h(1,i)=0 m=1;disp(请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解.A的秩RA和各阶顺序主子式值hl依次如下：)for j=1:nU(1,j)=A(1,j);endfor k=2:nfor i=2:n for j=2:n L(1,1)=1;L(i,i)=1; if ijL(1,1)=1;L(2,1)=A(2,1)/U(1,1); L(i,1)=A(i,1)/U(1,1);L(i,k)=(A(i,k)- L(i,1:k-1)*U(1:k-1,k)/U(k,k);elseU(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j);endendendendhl;RA,U,Lendend运行结果：系数矩阵的阶数3线性方程组的系数矩阵1,2,3;1,7,9;3,6,8线性方程组的常数项1