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节点电压幅值的半不变量
时间: 2024-02-13 14:03:35 浏览: 28
这段代码是在计算概率潮流时,对雅可比矩阵进行变换,并计算节点电压的半不变量。首先,将雅可比矩阵Jacco2的行和列对应PV节点和平衡节点的位置的元素都设置为0,将对应位置的元素设置为1,以保证对应节点的电压相角和电压幅值不变。然后,将Jacco2中对应PV节点的行和列都设置为0,将对角线上对应位置的元素设置为1,以保证对应节点的电压幅值不变。最后,将Jacco矩阵的大小扩展为2nb×2nb,并将其对角线上对应位置的元素设置为1,以保证平衡节点的电压相角和电压幅值不变。接着,通过求解线性方程组s1=Jacco\eye(size(Jacco)),得到雅可比矩阵的逆矩阵,并计算电压相角和电压幅值的半不变量。最后,从电压的半不变量中提取出PQ节点的电压幅值的半不变量gama_vm。
相关问题
%-修正平衡节点和pv节点对应的系数 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);
这段代码主要是修正雅可比矩阵,以便进行牛顿-拉夫逊潮流计算。具体解释如下:
1. Jacco=Jacco2;:将计算得到的雅可比矩阵Jacco2赋值给Jacco。
2. Jacco(isb,:)=0; Jacco(:,isb)=0; Jacco(isb,isb)=1;:将平衡节点对应的行和列全部置为0,同时将平衡节点对应的对角线元素置为1。这是因为平衡节点的电压相角和电压幅值是固定的,不参与计算。
3. Jacco(isb+nb,:)=0; Jacco(:,isb+nb)=0; Jacco(isb+nb,isb+nb)=1;:同理,将平衡节点对应的虚部行和列全部置为0,同时将平衡节点对应的虚部对角线元素置为1。
4. Jacco(pv+nb,:)=0; Jacco(:,pv+nb)=0;:将PV节点对应的行和列全部置为0,因为PV节点的无功功率不参与计算。
5. Jacco(pv+nb,pv+nb)=1;:将PV节点对应的无功功率的对角线元素置为1,因为PV节点的有功功率是由其他节点的电压幅值来调节的。
6. Jacco=Jacco+sparse(nb+pv, nb+pv, ones(npv,1), 2*nb, 2*nb);:将Jacco中的所有元素都加上一个稀疏矩阵,其大小为(nb+pv)×(nb+pv),其中稀疏矩阵的对角线元素全部为1,其余元素全部为0。这相当于向Jacco的对角线元素中添加了一些额外的1,以保证Jacco的行列式不为0,从而避免牛顿-拉夫逊潮流计算中的奇异性问题。
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节点电压幅值的半不变量
这是一段 MATLAB 代码,代码中定义了一个名为 s1 的矩阵,使用了 size 函数根据 Jacco 矩阵的大小创建了一个同样大小的单位矩阵。同时定义了一个名为 gama_v 的零矩阵,并使用了 for 循环计算了 gama_v 的值。最后,通过切片操作获取了 gama_vm 矩阵。具体的算法和应用需要结合上下文才能理解。
相关推荐
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![application/x-zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![application/x-rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)