城市轨道交通供电系统研究（Matlab代码实现）-Toy模板网

这篇具有很好参考价值的文章主要介绍了城市轨道交通供电系统研究（Matlab代码实现）。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

城市轨道交通供电系统是为城市轨道交通运营提供所需电能的系统，不仅为城市轨道交通电动列车提供牵引用电,而且还为城市轨道交通运营服务的其他设施提供电能，如照明、通风、空调、给排水、通信、信号、防灾报警、自动扶梯等，应具备安全可靠、技术先进、功能齐全、调度方便和经济合理等特点。

在城市轨道交通的运营中，供电一旦中断，不仅会造成城市轨道交通运输系统的瘫痪,还会危及乘客生命与财产安全。因此，高度安全可靠而又经济合理的电力供给是城市轨道交通正常运营的重要保证和前提。

📚2 运行结果

城市轨道交通供电系统研究（Matlab代码实现）

部分代码：

clear;
clc;
load('Tram.mat');
pelect=zeros(size(T,2),1);
% Extracting the power demand from the data structure T
for i=1:size(T,2)
pelect(i,1)=T(1,i).pelec;
end
figure;
plot(1:size(T,2),pelect);
title('Power Supply Profile');
% Calculating the maximum power demand of the time period
pmax=max(pelect);
ptotal=0;
for i=1:size(pelect,1)
ptotal=ptotal+pelect(i,1);
end
% ptotal=sum(pelect);
% Calculating the average pwower demand over the time period
pavg=ptotal/size(pelect,1);
disp('Average power: ' );
disp(pavg);
% Calculating the PPH
% If both average and maximum are greater than zero,
% Pavg and Pmax will mostly be greater than zero.
% If it is less than zero, we definitely need storage and hence PPH is set
% to 1.
if(pavg>=0 && pmax>0)
pph=1-(pavg/pmax);
else
pph=1;
end
disp('Potential for hybridization in power: ' );
disp(pph)
energyInBatt=zeros(size(pelect,1)-1,1);
% Pavg is assumed to be provided by the external supplier and the peaks
% shall be managed by the storage, so the energy in battery shall be
% calculated on the basis of the difference between the electrical power
% and the average power.
for i=2:size(pelect,1)
energyInBatt(i,1)=energyInBatt(i-1,1)-(pelect(i-1,1)-pavg);
end
figure;
plot(energyInBatt);
title('Energy in battery (in J)');
% Useful energy of a battery is definaed as the difference between the maximum and
% minimum energy stored in a time cycle
eu=max(energyInBatt)-min(energyInBatt);
disp(strcat('Useful Energy (in Joules): ',num2str(eu)));
if(eu~=0 && pmax>=0)
phe=pmax/eu;
else
phe=Inf;
end
disp('Potential for hybridization in energy: ' );
disp(phe)