matlab如何形成节点导纳矩阵,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程...
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
求大神帮忙看看这个程序,解释一下!
%形成节点导纳矩阵
%考虑参考节点
%考虑互感支路
function [d_g,d_b,u_g,u_b]=calyy(itoi,noels_hgln,noelsnum_hgln,bs_noels,bs_cp,gb_hgln,gb_noels,gb_cp,lv_bs,lv_hgln,lv_noels,lv_cp,type_bs,iy,jy,remv_noels)
%lv_bs
d_g=zeros(1,lv_bs);
d_b=zeros(1,lv_bs);
u_g=zeros(1,length(jy));
u_b=zeros(1,length(jy));
type_noels=zeros(1,lv_noels);
for i=1:lv_hgln
%%支路自阻抗
lnnum=noelsnum_hgln(i);
lnnum;
gb_noels;
for ii1=1:lnnum
iln=noels_hgln(i,ii1);
[d_g,d_b,u_g,u_b]=getlnselfgb(iln,bs_noels,type_bs,itoi,gb_noels,iy,jy,d_g,d_b,u_g,u_b,lv_bs);
type_noels(iln)=1;
end
for ii1=1:lnnum
iln=noels_hgln(i,ii1);
for ii2=ii1+1:lnnum
zln=noels_hgln(i,ii2);
b=gb_hgln(i,ii1,ii2);
[d_g,d_b,u_g,u_b]=getlnhugb(iln,zln,bs_noels,type_bs,itoi,0,b,iy,jy,d_g,d_b,u_g,u_b,lv_bs);
end
end
end
for i=1:lv_noels
if remv_noels(i)==1%被除去线路
continue;
end
ibs=bs_noels(i,1);
zbs=bs_noels(i,2);
if type_noels(i)==1
continue;
end
type_noels(i)=1;
iibs=itoi(ibs);
zzbs=itoi(zbs);
g=gb_noels(i,1);
b=gb_noels(i,2);
j=0;
if iibs~=zzbs
j=getlocal(iy,jy,iibs,zzbs,lv_bs,1);
end
bch=gb_noels(i,3);
%分接头所在端的自导纳要+=(g+jb)/t^2
if bch<0
bch=-bch;
if gb_noels(i,4)==1%变比位于末端
ibs=iibs;
zbs=zzbs;
else %变比位于首段
ibs=zzbs;
zbs=iibs;
end
if ibs<=lv_bs
d_b(ibs)=d_b(ibs)+b;
d_g(ibs)=d_g(ibs)+g;
end
if zbs<=lv_bs
d_b(zbs)=d_b(zbs)+b/bch/bch;
d_g(zbs)=d_g(zbs)+g/bch/bch;
end
if j~=0
u_g(j)=u_g(j)-g/bch;
u_b(j)=u_b(j)-b/bch;
end
continue;
end
%disp('线路支路');
if iibs<=lv_bs
d_b(iibs)=d_b(iibs)+b;
d_b(iibs)=d_b(iibs)+bch;
d_g(iibs)=d_g(iibs)+g;
end
if zzbs<=lv_bs
d_b(zzbs)=d_b(zzbs)+b;
d_b(zzbs)=d_b(zzbs)+bch;
d_g(zzbs)=d_g(zzbs)+g;
end
if j~=0
u_g(j)=u_g(j)-g;
u_b(j)=u_b(j)-b;
end
end
for i=1:lv_cp
ibs=bs_cp(i,1);
if ibs==0
continue;
end
headbs=itoi(ibs);
if headbs>lv_bs
continue;
end
d_b(headbs)=d_b(headbs)+gb_cp(i,1);
end
k=iy(lv_bs);
u_g=u_g(1:k-1);
u_b=u_b(1:k-1);
function [d_g,d_b,u_g,u_b]=getlnselfgb(iln,bs_noels,type_bs,itoi,gb_noels,iy,jy,d_g,d_b,u_g,u_b,lv_bs)
g=gb_noels(iln,1);
b=gb_noels(iln,2);
bch=gb_noels(iln,3);
if bch<0
disp('error!!!');
return;
end
ibs=bs_noels(iln,1);
zbs=bs_noels(iln,2);
iibs=itoi(ibs);
zzbs=itoi(zbs);
if iibs~=lv_bs+1
d_b(iibs)=d_b(iibs)+b+bch;
d_g(iibs)=d_g(iibs)+g;
end
if zzbs~=lv_bs+1
d_b(zzbs)=d_b(zzbs)+b+bch;
d_g(zzbs)=d_g(zzbs)+g;
end
j=0;
if iibs~=zzbs
j=getlocal(iy,jy,iibs,zzbs,lv_bs,1);
end
if j==0
disp('找不到支路位置');
%exit;
return;
end
u_g(j)=u_g(j)-g;
u_b(j)=u_b(j)-b;
function [d_g,d_b,u_g,u_b]=getlnhugb(iln,zln,bs_noels,type_bs,itoi,g,b,iy,jy,d_g,d_b,u_g,u_b,lv_bs)
ibsi=bs_noels(iln,1);%首
zbsi=bs_noels(iln,2);%末
ibsz=bs_noels(zln,1);%首
zbsz=bs_noels(zln,2);%末
first_head=itoi(ibsi);
first_tail=itoi(zbsi);
second_head=itoi(ibsz);
second_tail=itoi(zbsz);
if first_head==second_head&&first_head~=lv_bs+1
disp('add');
d_b(first_head)=d_b(first_head)+2*b;
d_g(first_head)=d_g(first_head)+2*g;
end
if first_tail==second_tail&&first_tail~=lv_bs+1
disp('add');
d_b(first_tail)=d_b(first_tail)+2*b;
d_g(first_tail)=d_g(first_tail)+2*g;
end
if first_head==second_tail&&first_head~=lv_bs+1
disp('add');
d_b(first_head)=d_b(first_head)-2*b;
d_g(first_head)=d_g(first_head)-2*g;
end
if first_tail==second_head&&first_tail~=lv_bs+1
disp('add');
d_b(first_tail)=d_b(first_tail)-2*b;
d_g(first_tail)=d_g(first_tail)-2*g;
end
[d_g,d_b,u_g,u_b]=addlnval(first_head,second_head,g,b,iy,jy,d_g,d_b,u_g,u_b,lv_bs);
[d_g,d_b,u_g,u_b]=addlnval(first_head,second_tail,-g,-b,iy,jy,d_g,d_b,u_g,u_b,lv_bs);
[d_g,d_b,u_g,u_b]=addlnval(second_head,first_tail,-g,-b,iy,jy,d_g,d_b,u_g,u_b,lv_bs);
[d_g,d_b,u_g,u_b]=addlnval(second_tail,first_tail,g,b,iy,jy,d_g,d_b,u_g,u_b,lv_bs);
function [d_g,d_b,u_g,u_b]=addlnval(first_head,second_head,g,b,iy,jy,d_g,d_b,u_g,u_b,lv_bs)
if first_head==second_head
return;
end
j=getlocal(iy,jy,first_head,second_head,lv_bs,1);
if j==0
return;
end
u_b(j)=u_b(j)+b;
u_g(j)=u_g(j)+g;
matlab如何形成节点导纳矩阵,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程...相关推荐
- matlab求节点导纳矩阵,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 求大神帮忙看看这个程序,解释一下! %形成节点导纳矩阵 %考虑参考节点 %考虑互感支路 function [d_g,d_b,u_g,u_b]=calyy( ...
- matlab节点导纳阵求逆,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 求大神帮忙看看这个程序,解释一下! %形成节点导纳矩阵 %考虑参考节点 %考虑互感支路 function [d_g,d_b,u_g,u_b]=calyy( ...
- 导纳矩阵 matlab,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 求大神帮忙看看这个程序,解释一下! %形成节点导纳矩阵 %考虑参考节点 %考虑互感支路 function [d_g,d_b,u_g,u_b]=calyy( ...
- matlab如何建立矩阵导纳,关于利用矩阵稀疏技术求解节点导纳矩阵的MATLAB编程
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 求大神帮忙看看这个程序,解释一下! %形成节点导纳矩阵 %考虑参考节点 %考虑互感支路 function [d_g,d_b,u_g,u_b]=calyy( ...
- MATLAB 电力系统潮流计算(采用稀疏矩阵,可计算1000节点)
MATLAB电力系统潮流计算 用MATLAB利用牛顿-拉夫逊法计算电科院22节点系统(使用稀疏技术),计算出各节点的电压相角.各支路的功率等.使用了稀疏计算技术,能计算1000节点以上的系统,且计算时 ...
- matlab找多项式最值,如何利用matlab计算最小多项式
7.1.1 分段线性插值 所谓分段线性插值就是通过插值点用折线段连接起来逼近原曲线,这也是计算机绘制图形的基本原理.实现分段线性插值不需编制函数程序,MATLAB自身提供了内部函数interp1其主要 ...
- matlab 求矩阵秩,求矩阵秩的两种方法及MATLAB的应用
摘 要: 高等代数是一门逻辑思维比较强和理论知识比较深的学科, 它具有丰富的数学知识, 涉及许多重要的数学思想, 其在数学领域的应用很广泛, 如行列式.矩阵的相关计算和求解线性方程组的解方面的应 ...
- 利用对位相乘法计算线性卷积-附Matlab代码
目录 1.线性卷积原理 2.利用对位相乘法计算线性卷积 3.Matlab实操 3.1源代码 3.2仿真结果 线性卷积原理 公示(1)为两个离散非周期序列进行线性卷积 一般我们计算线性卷积时使用图解法最 ...
- matlab将表格数据做成图,入excel表格数据做图-如何向matlab中导入excel表格 利用表格数据做出图......
如何向matlab中导入excel表格 利用表格数据做出图... 1.将待导入的矩阵结构的数据录入Excel中,录入时注意行列要跟原矩阵一一对应 2.录入完以后保存数据,为了后续步骤使用方便,命名时我 ...
最新文章
- Ubuntu 系统 下 AndroidStudio 工具 使用ctrl+z 回退上一步的方法
- GT Transceiver中的重要时钟及其关系(5)QPLL的工作原理介绍
- Mybatis是如何实现SQL语句复用功能的?
- 面试:说说你对“零拷贝”的理解?
- xgboost调参指南
- 哄骗DB2look 重新创立优化器访问筹划(3)
- 纹理对象纹理单元纹理目标_网页设计理论:纹理
- JS调试的时候遇到无限debugger怎么办?
- 主流视频客户端核心代码的实现
- SQL Server 2012入门T-SQL基础篇:(7)Where子句与Having子句的区别
- java通用排序工具类
- 证件照如何弄蓝白渐变色背景
- HC32F460 浮点运算开启
- linux floating ip,neutron分配浮动IP(floating ip)的方法
- 搜索插件像百度那样的智能感知效果
- 刻意练习(重点摘抄)
- android微信7,微信7.0安卓版之初体验
- 视频文件的原声如何批量更换
- 时光穿梭机(删除文件)
- 30 款 IDEA 宝贝插件
热门文章
- cad中lisp文件给恶作剧_[软件心得]我用的CAD,LISP加载文件
- 校园人到企业人的转变
- 小程序中使用echarts
- 关于马云的“996”
- Java找到休息日问题
- html鼠标滑过图片透明,鼠标滑过图片透明度发生改变的特效 - YangJunwei
- The type okio.ByteString cannot be resolved. It is indirectly referenced from required .class files
- Postman:传递的参数是List类型时 传参格式的写法
- N-Gram文件格式介绍 - ARPA
- Google signin 引发的崩溃(人)