李路稀 李路稀
关注数: 38 粉丝数: 148 发帖数: 2,669 关注贴吧数: 20
求教为何报错? Jacbi=zeros(2*n-2); h=0; k=0; for i=1:n %对PQ节点的处理 if B2(i,2)==1 h=h+1; for j=1:n if j~=isb k=k+1; if i==j %对角元素的处理 Jacbi(2*h-1,2*k-1)=-B(i,j)*B2(i,5)+G(i,j)*B2(i,6)+imag(I(h,1)); Jacbi(2*h-1,2*k)=G(i,j)*B2(i,5)+B(i,j)*B2(i,6)+real(I(h,1)); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1)); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1)); else %非对角元素的处理 Jacbi(2*h-1,2*k-1)=-B(i,j)*B2(i,5)+G(i,j)*B2(i,6); Jacbi(2*h-1,2*k)=G(i,j)*B2(i,5)+B(i,j)*B2(i,6); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1); end if k==(n-1) %将用于内循环的指针置于初始值,以确保雅可比矩阵换行 k=0; end end end end end k=0; for i=1:n %对PV节点的处理 if B2(i,2)==2 h=h+1; for j=1:n if j~=isb k=k+1; if i==j %对角元素的处理 Jacbi(2*h-1,2*k-1)= -B(i,j)*B2(i,5)+G(i,j)*B2(i,6)+imag(I(h,1)); Jacbi(2*h-1,2*k)= G(i,j)*B2(i,5)+B(i,j)*B2(i,6)+real(I(h,1)); Jacbi(2*h,2*k-1)=2*B2(i,6); Jacbi(2*h,2*k)=2*B2(i,5); else %非对角元素的处理 Jacbi(2*h-1,2*k-1)= -B(i,j)*B2(i,5)+G(i,j)*B2(i,6); Jacbi(2*h-1,2*k)= G(i,j)*B2(i,5)+B(i,j)*B2(i,6); Jacbi(2*h,2*k-1)=0; Jacbi(2*h,2*k)=0; end if k==(n-1) %将用于内循环的指针置于初始值,以确保雅可比矩阵换行 k=0; end end end end end %disp('雅克比矩阵为:') ; Jacbi; %列参数跟书上的相反 这里面第一列是对电压相位的偏导 第二列是对幅值的偏导 %求解修正方程,获取节点电压的不平衡量 DetaU=zeros(2*n-2,1); DetaU=inv(Jacbi)*DetaS;
1 下一页