%----------计算节点电压的半不变量 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节点电压幅值的半不变量 %-----------形成支路潮流的灵敏度矩阵 li=shuju.branch(:,1); lj=shuju.branch(:,2); b=size(shuju.branch,1); %线路数 G11=zeros(b,nb); %线路有功对角度的倒数 G12=zeros(b,nb); %线路有功对电压的导数 G21=zeros(b,nb); G22=zeros(b,nb); for i=1:b G11(i,li(i))=-H(li(i),lj(i)); G11(i,lj(i))=H(li(i),lj(i)); G12(i,li(i))=0.01*2*branch(i,14)/bus(li(i),8)-N(li(i),lj(i)); %%化为标幺值 G12(i,lj(i))=N(li(i),lj(i)); G21(i,li(i))=-M(li(i),lj(i)); G21(i,lj(i))=M(li(i),lj(i)); G22(i,li(i))=0.01*2*branch(i,15)/bus(li(i),8)-H(li(i),lj(i)); %%化为标幺值 G22(i,lj(i))=H(li(i),lj(i)); end G=[G11 G12;G21 G22]; t1=G*s1;
时间: 2023-10-30 21:04:51 浏览: 65
这段代码实现了对电力系统进行潮流灵敏度分析的功能。
首先,通过计算节点电压的半不变量,得到电压相角、幅值。计算过程中使用了Jacobi矩阵的逆矩阵s1,按照幂次逐次计算,得到gama_v矩阵。
然后,根据电力系统的支路信息,形成支路潮流的灵敏度矩阵G。灵敏度矩阵G由四部分组成,分别是G11、G12、G21和G22。其中,G11和G21是线路有功对节点相角的导数,G12和G22是线路有功对节点电压幅值的导数。
最后,将灵敏度矩阵G与节点电压的半不变量s1相乘,得到支路潮流的灵敏度矩阵t1。
相关问题
%----------计算节点电压的半不变量 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节点电压幅值的半不变量
这段代码是用来计算节点电压的半不变量,主要分为以下几个步骤:
1. 通过 Jacco 矩阵和单位矩阵求解 s1 矩阵,其中 Jacco 矩阵是电力系统节点导纳矩阵的一部分。
2. 定义一个零矩阵 gama_v,用来存储电压相角和幅值,其中 2*nb 表示节点数的两倍,jieshu 表示计算的最高次数。
3. 使用 for 循环,计算电压相角和幅值。具体来说,对于每个 i,使用 s1 的 i 次方乘以 gama_pq 的第 i 列,得到 gama_v 的第 i 列。
4. 通过切片操作,提取出 pq 节点对应的电压幅值,即 gama_v 的第 nb+pq 行。
这个半不变量在电力系统的节点电压稳定性分析中有重要的作用。
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。
相关推荐
![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)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![m](https://img-home.csdnimg.cn/images/20210720083646.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![xlsx](https://img-home.csdnimg.cn/images/20210720083732.png)
![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)