基于Matlab的PV模型仿真小实验

-Name: Yang , Vorname: Ying Student ID: 359144PV Assignment 1Exercise 1Code:%Exercise 1close allclear all %define input valuesDOY=172;TZ=2;lambda_g=13.2;fi_g=52.3;LT=6:1:20;%call function sundata to get sun enevation and azimuth Am, Ys, SAz = sundata(lambda_g,fi_g,TZ,DOY,LT); plot(SAz,Ys);xlabel(Azimuth in degrees As);ylabel(Altitude Position in degrees Ys);title(Berlin(13.2,52.3) on 21.06.2013(172/365);grid on;Function sundata:function Am, Ys, SAz = sundata(lambda_g,fi_g,TZ,DOY,LT)%UNTITLED Summary of this function goes here% Detailed explanation goes here J=360*DOY/365; TEQ=0.0066+7.3525*cos(J+85.9)*pi/180)+9.9359*cos(2*J+108.9)*pi/180)+0.3387*cos(3*J+105.2)*pi/180);delta=0.3948-23.2559*cos(J+9.1)*pi/180)-0.3915*cos(2*J+5.4)*pi/180)-0.176*cos(3*J+26)*pi/180); % delta units-degree for k=1:1:15 TLT(k)=LT(k)-TZ+ (4*lambda_g+TEQ)/60; W(k)=(12-TLT(k)*15; Ys(k)=asin( cos(W(k)*pi/180)*cos(fi_g*pi/180)*cos(delta*pi/180) +sin(fi_g*pi/180)*sin(delta*pi/180) )*180/pi; SunA(k)=acos( (sin(Ys(k)*pi/180)*sin(fi_g*pi/180)-sin(delta*pi/180) )/(cos(Ys(k)*pi/180)*cos(fi_g*pi/180) )*180/pi; if LT(k)0 Vstep(i+1)=Vstep(i);else Vstep(i+1)=-Vstep(i);end D(i+1)=Vo/(Vo+V(i+1); Vnew(i+1)=V(i+1)+Vstep(i+1);endend Figure 6(G=1000, T=25) would be given together with the figure 7Exercise 7Code:% Exercise 7 close all;clear all;%define the new voltage step as 0.1Vdelta_V=0.1;Vo=30;G=input(Please input the value of Irrad:);T=input(Please input the value of Temperature:); %call function to get D and Vnew D,Vnew = MPP( G,T,delta_V,Vo ); subplot(2,1,1);plot(D);xlabel(Time steps);ylabel(Duty cycle);grid on;subplot(2,1,2);plot(Vnew);xlabel(Time steps);ylabel(Module operating voltage(V);grid on;Figure 6(G=1000, T=25): Figure 7(G=1000, T=25):Compare with the figures of exercise 6 and 7, get the conclusion:The more the Voltage step is, the speed of tracking is faster but the accuracy of tracking is lower. Exercise 8Code:% Exercise 8close all;clear all;Vo=30;G=input(Please input the value of Irrad:);T=input(Please input the value of Temperature:);%call function to get D and Vnew D,Vnew,P = AHCA( G,T,Vo ); subplot(2,1,1);plot(D);xlabel(Time steps);ylabel(Duty cycle);grid on;subplot(2,1,2);plot(Vnew);xlabel(Time steps);ylabel(Module operating voltage(V);grid on;Figure 8(G=1000, T=25): Figure 7(G=1000, T=25):Compare with the figures of exercise 8 and 7, get the conclusion that: the adaptive hill climbing can much faster track the maximum power point and the accuracy is also better.function D, Vnew, P = AHCA( G,T ,Vo)%UNTITLED3 Summary of this function goes here% Detailed explanation goes hereV=zeros(1,400);Vnew=zeros(1,400);Vstep=zeros(1,400);P=zeros(1,400);D=zeros(1,400);dP=zeros(1,400);dV=zeros(1,400);%define the first two stepsV(1)=0;D(1)=Vo/(V(1)+Vo);I=PVmod(V(1),G,T);P(1)=I*V(1);Vstep(1)=0.4;V(2)=V(1)+Vstep(1);I=PVmod(V(2),G,T);P(2)=V(2)*I;Vnew(1)=V(1)+Vstep(1);D(2)=Vo/(V(2)+Vo);%Voltage step is dP/dVfor i=2:1:400 dP(i)=P(i)-P(i-1); dV(i)=V(i)-V(i-1); Vstep(i)=dP(i)/dV(i); Vnew(i)=V(