%--------------用半不变量求Gram-charlie展开系数 g_vm=zeros(size(gama_vm)); g_xianlu_p=zeros(size(gama_xianlu_p)); g_xianlu_q=zeros(size(gama_xianlu_q)); for i=1:jieshu g_vm(:,i)=gama_vm(:,i)./(gama_vm(:,2).^(i/2)); g_xianlu_p(:,i)=gama_xianlu_p(:,i)./(gama_xianlu_p(:,2).^(i/2)); g_xianlu_q(:,i)=gama_xianlu_q(:,i)./(gama_xianlu_q(:,2).^(i/2)); end
时间: 2024-03-31 14:37:37 浏览: 60
这段代码是用半不变量求 Gram-Charlie 展开系数,具体包括:
- `g_vm`、`g_xianlu_p` 和 `g_xianlu_q` 分别表示节点电压幅值、支路有功和无功的 Gram-Charlie 展开系数,它们都是大小为 `nb x jieshu` 或 `b x jieshu` 的矩阵;
- `gama_vm`、`gama_xianlu_p` 和 `gama_xianlu_q` 分别表示节点电压幅值、支路有功和无功的半不变量,它们都是大小为 `nb x jieshu` 或 `b x jieshu` 的矩阵;
- `i` 表示需要计算的展开系数的阶数;
- `gama_vm(:,2)`、`gama_xianlu_p(:,2)` 和 `gama_xianlu_q(:,2)` 分别表示节点电压幅值、支路有功和无功的半不变量的二阶半不变量;
- `gama_vm(:,i)./(gama_vm(:,2).^(i/2))`、`gama_xianlu_p(:,i)./(gama_xianlu_p(:,2).^(i/2))` 和 `gama_xianlu_q(:,i)./(gama_xianlu_q(:,2).^(i/2))` 分别表示将第 `i` 阶半不变量除以对应的二阶半不变量的 `i/2` 次方,得到 Gram-Charlie 展开系数。
具体来说,Gram-Charlie 展开系数是用于描述电力系统的无穷小扰动响应的一种数学工具。它们可以通过半不变量求解,具体方法是将半不变量表示为 Gram-Charlie 展开式,然后对展开式求导,得到展开系数。
阅读全文