www.gusucode.com > 太阳能三种mppt算法包括干扰观察法，电导增益法，变步长电导增益法 > mppt/PandOt.m

    sun=SUN('radio.txt');
T=zeros(size(sun));
T(:)=25;
v=0.05;
P_max=MAXPOW(sun,v,T);
P_track=PandO(sun,v,T);
err=P_max-P_track;
kpm=KPM(P_track,P_max);
figure(1)

subplot(2,2,1)
plot(P_max,'b'),hold on
grid on
plot(P_track,'k:'),hold on
plot(err,'r--')
title('采用干扰观察算法')
xlabel('t(s)');
ylabel('P(W)');
kpm

P_max=MAXPOW(sun,v,T);
P_track=INC1(sun,v,T);
err=P_max-P_track;
kpm=KPM(P_track,P_max);
subplot(2,2,2)
plot(P_max,'b'),hold on
grid on
plot(P_track,'k:'),hold on
plot(err,'r--');
title('采用INCMPPT算法')
xlabel('t(s)');
ylabel('P(W)');
kpm


P_max=MAXPOW(sun,v,T);
P_track=INC2(sun,v,T);
err=P_max-P_track;
kpm=KPM(P_track,P_max);
subplot(2,2,3)
plot(P_max,'b'),hold on
grid on
plot(P_track,'k:'),hold on
plot(err,'r--')
title('采用变步长MPPT算法')
xlabel('t(s)');
ylabel('P(W)');
kpm
legend('Pmax','Ptrack','err',4)