Matlab模拟风电-光伏最大功率点追踪MPPT曲线

wm = [0:0.001:2.0]; R = 80; %radius of wind turbine rou = 1.29; %air density %wind velocity Vw = [7:1:16]; pi = 3.141592653; e = 2.71828; Pm = zeros(9,2001); lambada = zeros(9,2001); lambada1 = zeros(9,2001); Pmppt = zeros(1,2001); for g = 1:1:9 for k = 1:1:2001 lambada(g,k) = wm(k)*R/Vw(g); lambada1(g,k) = 1/(1/lambada(g,k)-0.035); Cp(g,k) = 0.22*(116/lambada1(g,k)-5.0)*e^(-12.5/lambada1(g,k)); Pm(g,k) = 0.5*R*R*pi*rou*Cp(g,k)*Vw(g)^3; end plot(wm,Pm(g,:)); hold on; end for k = 1:1:2001 if (wm(k) > 0.5)&(wm(k) < 1.1) Pmppt(k) = 0.598*1.8e7*wm(k)*wm(k)*wm(k); end if wm(k) >= 1.1 Pmppt(k) = 1.8e7; end end plot(wm,Pmppt); grid on;