Jacco=Jacco2; Jacco(isb,:)=0; Jacco(:,isb)=0; Jacco(isb,isb)=1; Jacco(isb+nb,:)=0; Jacco(:,isb+nb)=0; Jacco(isb+nb,isb+nb)=1; Jacco(pv+nb,:)=0; Jacco(:,pv+nb)=0; Jacco=Jacco+sparse(nb+pv, nb+pv, ones(npv,1), 2*nb, 2*nb); %----------计算节点电压的半不变量 s1=Jacco\eye(size(Jacco)); gama_v=zeros(2*nb,jieshu); %电压相角、幅值 for i=1:jieshu gama_v(:,i)=(s1.^i)*gama_pq(:,i); end gama_vm=gama_v(nb+pq,:); %pq节点电压幅值的半不变量

这段代码是在计算概率潮流时，对雅可比矩阵进行变换，并计算节点电压的半不变量。首先，将雅可比矩阵Jacco2的行和列对应PV节点和平衡节点的位置的元素都设置为0，将对应位置的元素设置为1，以保证对应节点的电压相角和电压幅值不变。然后，将Jacco2中对应PV节点的行和列都设置为0，将对角线上对应位置的元素设置为1，以保证对应节点的电压幅值不变。最后，将Jacco矩阵的大小扩展为2nb×2nb，并将其对角线上对应位置的元素设置为1，以保证平衡节点的电压相角和电压幅值不变。接着，通过求解线性方程组s1=Jacco\eye(size(Jacco))，得到雅可比矩阵的逆矩阵，并计算电压相角和电压幅值的半不变量。最后，从电压的半不变量中提取出PQ节点的电压幅值的半不变量gama_vm。