Analysis for Transfer Function
Considering a grid connected system with grid as current harmonic source
Fig. 1: Grid connected System
Considering LCL equivalent as ZLCL;
Fig. 2: Equivalent Circuit
Where,
Z L : Resistive Load
Zg
: Grid Impedance
sC s 3 L2C sL
Z LCL : LCL Equivalent = s 2 LC 1
I Ih : I g I L I LCL
Where:
I L : Load Current
Ig
: Grid Impedance Current
I LCL : LCL filter Current
Now applying KCL we get,
Vo Vo Vo 1 1 1
Ih Vo
Z g Z l Z LCL
Z g Z l Z LCL
Rearranging;
I h Z l Z g Z LCL
Vo Z g Z l Z LCL
Now the transfer function V0/Ih (Output / Input);
Vo Z g Z l Z LCL
I h Z l Z g Z g Z LCL Z g Z LCL
Considering without Grid Impedance the transfer function will be:
Vo ZZ
l LCL
I h Z l Z LCL
MATLAB Code
clear;
clc;
Li = 0.89e-3
Lg = 0.20e-3
C = 30e-6
Rd = 0.8
Rl = 4.88 (Per Phase Load)
Lx = 0.9e-3 (Grid Inductance)
Rg = 0.4 (Grid Resistance)
%***********************System Transfer Function**************************%
s= tf('s')
num1 = [Li*Lg*Lx*C*Rl Rl*C*Li*Lg+Rl*C*Rd*Lx*Li+Rl*C*Rd*Lx*Lg
Rl*Rd*Rg*C*Li+Rl*Rd*Rg*C*Lg+Rl*Lx*Li+Rl*Lx*Lg Rg*Rl*Li+Rg*Rl*Lg 0]
den1 = [Li*Lg*Lx*C Li*Lg*C*Rl+C*Li*Lg*C*Rg+C*Rd*Lx*Li+C*Rd*Lx*Lg
Rd*Rl*C*Li+Rd*Rl*C*Lg+Rd*Rg*C*Li+Rd*Rg*C*Lg+Lx*Li+Lx*Lg
Rl*Li+Rl*Lg+Rg*Li+Rg*Lg+Lx Rg+Rl]
sys1 = tf(num1,den1)
%***************************LPF 4th Order*********************************%
LPF = tf ([0 2*3.14285*10], [1 2*3.14285*10])
LPF_4order = LPF^4
%************************Final Tf*****************************************%
sys2 = tf (LPF_4order*sys1)
%******************PLOTS**************************************************%
subplot (2,1,1)
step (sys2);
title('Step Response')
subplot(2,1,2)
step (sys2/s);
title('Ramp Response')
%*************************************************************************%
Results