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

Name: Yang , Vorname: Ying Student ID: 359144

PV Assignment 1

Exercise 1

Code:

%Exercise 1

close all

clear all

%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)

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;

)/(cos(Ys(k)*pi/180)*cos(fi_g*pi/180)) )*180/pi;

if LT(k)<=12

SAz(k)=180-SunA(k);

else

SAz(k)=180+SunA(k);

Figure 1:

Exercise 2

Code:

close all;

clear all;

%define V,G,T to call function

V=0:0.01:44.8;

subplot(2,1,1);

plot(V,I);

xlabel('Voltage(V)');

ylabel('Current(A)');

grid on;

for i=1:1:4481

P(i)=I(i)*V(i);

end

subplot(2,1,2);plot(V,P);

xlabel('Voltage(V)');

Function PVmod:

% Detailed explanation goes here

%Definations of constants:

k=1.38*10^(-23);

q=1.60*10^(-19);

Vg=1.12;

Voc=44.8;

Isc=5.5;

Tstc=25+273;

Tk=T+273;