clear?all;
close?all;
clc;
load（‘yuanshifengdian.mat‘）;
????pin=1/6000;
????tao=1./（2*pi*pin）;
????P=zeros（1length（yuanshi））;%儲(chǔ)能吸收功率
????Tc=10;%數(shù)據(jù)采樣間隔（s）
????E=zeros（11000）;%儲(chǔ)能容量
????E（11）=0;%儲(chǔ)能初始容量
????lvbo=zeros（1length（yuanshi））;%風(fēng)電經(jīng)儲(chǔ)能濾波后功率
????lvbo（11）=yuanshi（11）;
????Pmax=10000;
????Emax=10000;
????for?i=2:length（yuanshi）
????????lvbo（1i）=（tao/（tao+Tc））*lvbo（1i-1）+（Tc/（Tc+tao））*yuanshi（i1）;
????????P（1i）=（tao/（tao+Tc））*（yuanshi（i）-lvbo（1i-1））;
????????E（1i）=E（1i-1）+P（1i）*Tc/3600;
????????if?E（1i）<0
????????????P（1i）=-3600*E（1i-1）/Tc;
????????????E（1i）=0;
????????end
????end
????figure（1）
plot（P‘b‘）;hold?on;
plot（E‘r‘）;
legend（‘儲(chǔ)能充電功率‘‘儲(chǔ)能能量變化‘）;
figure（2）
plot（yuanshi‘b‘）;hold?on;
plot（lvbo‘r‘）;
legend（‘平抑前風(fēng)電‘‘平抑后風(fēng)電‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????887??2013-06-16?19:59??風(fēng)電濾波\ditonglvbo.m

?????文件??????26259??2013-04-23?19:42??風(fēng)電濾波\yuanshifengdian.mat

?????目錄??????????0??2013-10-24?22:05??風(fēng)電濾波

-----------?---------??----------?-----??----

????????????????27146????????????????????3