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----矩阵求逆的123!
方阵求逆 奇异矩阵:行列式为0的矩阵 非奇异矩阵:行列式不为0的矩阵,即可逆矩阵 对于方阵A,如果为非奇异方阵,则存在逆矩阵inv(A) 对于奇异矩阵或者非方阵,并不存在逆矩阵,但可以使用pinv(A ...
- matlab 分数计算器,Matlab节点导纳矩阵计算.doc
第一章 导纳矩阵的计算简介 1.1变压器的∏型等值电路 在电力系统潮流计算中,往往要计算节点导纳矩阵,而我们计算节点导纳矩阵采用节点电压法来实现,如在变压统 图1双绕组变压电 如果采用标计参数应归时标 ...
- 用不同的姿势求逆序对(复习篇)
用不同的姿势求逆序对(复习篇) 文章目录 用不同的姿势求逆序对(复习篇) 前言 讲解 归并排序 树状数组 线段树 题目 思路 代码 归并排序求逆序对 树状数组求逆序对 线段树求逆序对 历届试题 小朋友 ...
- 标量 向量 标量求导链式法则_矩阵求导术(上)
矩阵求导的技术,在统计学.控制论.机器学习等领域有广泛的应用.鉴于我看过的一些资料或言之不详.或繁乱无绪,本文来做个科普,分作两篇,上篇讲标量对矩阵的求导术,下篇讲矩阵对矩阵的求导术.本文使用小写字母 ...
- 全选主元matlab,全选主元高斯-约当法实矩阵求逆及传统方法求逆
/* 算法没有做过什么改动,加上了内存释放的部分. 算法介绍 矩阵求逆在3D程序中很常见,主要应用于求Billboard矩阵.按照定义的计算方法乘法运算,严重影响了性能.在需要大量Billboard矩 ...
最新文章
- R语言dim函数获取dataframe、matrix的维度、shape实战
- android-技术教程-调试程序的基础,在控制台上打印出想打印的东西 转载
- 基于图文界面的蓝牙扫描工具btscanner
- Tomcat 处理 HTTP 请求源码分析(上)【转】
- 合法的python变量名import_python 环境变量和import模块导入方法(详解)
- xml vs db.properties
- eclipse加载maven工程提示pom.xml无法解析org.apache.maven.plugins:maven-resources-plugin:2.4.3解决方案...
- 【STM32】按键---有关GPIO输入的那些事
- js 导出pdf上传至oss_前端上传图片到oss,压缩图片后上传至oss(补充图片文件旋转90度问题)...
- 试着当个“刺头”,不要被客户“牵着”鼻子走
- 去除PreferenceScreen过期提示
- 使用xlsx.core.min.js进行导入(表格)导出(Excel)功能
- java输战舰_什么是最好的战舰AI?
- 通过路由器端口映射,对局域网内服务器进行ssh登录
- html制作好看的五子棋(源码)
- 泛化、过拟合与欠拟合
- 志愿者积分兑换小程序开发制作
- 山还是山 水还是水
- html提取excel指定单元格数据,如何利用Excel函数公式快速提取出单元格中指定的数据...
- dell 服务器面板显示屏,Dell服务器面板错误码详解
热门文章
- 曦智科技发布最新光子计算处理器PACE
- 运维质变育新机,华为云能否引领政企运维破局?
- 震惊!程序员要放弃 Python 了!?发生了啥?
- 15年大厂经历!大佬总结:0基础如何学习Python?
- 【不吹不黑】详解容器技术架构、网络和生态
- 微服务落地,我们在考虑什么?| 技术头条
- Spark精华问答 | Spark的三种运行模式有何区别?
- linux grub 下载,GRUB 2.04发布下载,附新功能介绍
- mysql驱动profilesql_java连接mysql的驱动配置问题
- 工作流实战_21_flowable 加签 任务向前加签 向后加签