Three Phase Shunt Active Power Filter using pq theory
- 17 files
- 557 downloads
It solves the problem of harmonic current in power system. It also compensate reactive power and balances load on three phase(if load is unbalanced). It uses Clarke transformation to calculate value of p0,p and q, filters out constant part of p and thus finding the compensating part of p0,p and q. Then, it calculates compensating current by inverse-Clarke transformation. The compensating current is injected to the network via three phase inverter using hysteresis band current control(HBCC) method. In this way, harmonic and reactive compnent of load is supplied by this shunt active power filter.Also, compensating p0 means balancing load on three phase.
Cite As
Niranjan Bhujel (2019). Three Phase Shunt Active Power Filter using pq theory (https://www.mathworks.com/matlabcentral/fileexchange/61215-three-phase-shunt-active-power-filter-using-pq-theory), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (47)
Updates
| 2.0.7.0 | Exported to run on older version. |
|
| 2.0.6.0 | Exported for previous version (R2013a or latest). |
|
| 2.0.5.0 | Some bug fixed. |
|
| 2.0.4.0 | . |
|
| 2.0.3.0 | . |
|
| 2.0.2.0 | Made more informative. |
|
| 2.0.1.0 | Problem on PI controller has been fixed. |
|
| 2.0.0.0 | Problems of unbalanced load fixed. |
|
| 1.10.0.0 | Change in description. |
|
| 1.1.0.0 | Preview image is added. |
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Krishna Perka (view profile)
tetuko (view profile)
mr Niranjan Bhujel please repaired your model..thank u
tetuko (view profile)
any error after i run this model : some notification are :
Clear
05:56 PM Elapsed: 7 sec
Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.667089e-16. [4 similar]
The current directory 'c:\program files\matlab\r2018a\bin' is reserved for MATLAB files.
Please change your current directory to a writable directory outside of the MATLAB installation area.
Component:Stateflow | Category:General error
The current directory 'c:\program files\matlab\r2018a\bin' is reserved for MATLAB files.
Please change your current directory to a writable directory outside of the MATLAB installation area.
Component:Stateflow | Category:General error
The current directory 'c:\program files\matlab\r2018a\bin' is reserved for MATLAB files.
Please change your current directory to a writable directory outside of the MATLAB installation area.
Component:Stateflow | Category:General error
The current directory 'c:\program files\matlab\r2018a\bin' is reserved for MATLAB files.
Please change your current directory to a writable directory outside of the MATLAB installation area.
Component:Stateflow | Category:General error
The current directory 'c:\program files\matlab\r2018a\bin' is reserved for MATLAB files.
Please change your current directory to a writable directory outside of the MATLAB installation area.
Component:Stateflow | Category:General error
The current directory 'c:\program files\matlab\r2018a\bin' is reserved for MATLAB files.
Please change your current directory to a writable directory outside of the MATLAB installation area.
Mr. Niranjan , can you fix your ,model ?
Kaf khan (view profile)
can anyone provide 2016a model file please at > kafeel.khan23@gmail.com
tetuko (view profile)
do you have a email address ?
Muza Iqbal (view profile)
i need the analysis of this model (Report) , pleas at > muzammiliqbal351@gmail.com
moaz sardina (view profile)
i need the analysis of this model (Report) , pleas at > moaz.ahmed@pua.edu.eg
Chinmay Gadre (view profile)
Hey,really good simulation,being of great help. I had a few queries regarding PI control in conjunction with SPWM controller and a few more. If you could revert back at 'chinmaygadre1@gmail.com' it would be of great help. Thanks in advance.
fadhil zaidi (view profile)
how to insert it in matlab?
ASHOK KUMAR (view profile)
For Latest Simulation Projects,Please Contact http://asokatechnologies.in/
Madhura (view profile)
Can u pls send me model in .mdl matlab 2010R version at mail madhura.t@rediffmail.com
MOHD MARWAN NOOR DIN (view profile)
Binay Paudyal (view profile)
Eslavath Lakpathi (view profile)
Sir i wants this model in 2016b. can you please provide a link or send to my mail:lakpathieslavath@gmail.com
SINGARAVELAN S (view profile)
Sir i also wants this model in 2016a. can you please provide a link or send to my mail:cingameee@gmail.com
Shubh Yadav (view profile)
Shubh Yadav (view profile)
sir please send me 2015 or 2013 model at........8392shubh@gmail.com
MOHD MARWAN NOOR DIN (view profile)
Sir Niranjan Bhujel. sorry sir. can I know your email sir ? Or this is my Email sir mohdmarwan5217@gmail.com . I want ask you something urgent sir. I hope sir.
Niranjan Bhujel (view profile)
Only 6 mosfets are used (6 IGBTs can also be used). Fourth leg of inverter is connected to neutral of the system at the point of common coupling. G1, G2 and G3 are gate signals for the mosfets (or IGBTs) of the inverter.
MOHD MARWAN NOOR DIN (view profile)
fourth leg inverter sir (8 IGBT) and what is the function of G sir. For G1, G2,G3 in the inverter. TQ for answering my questions sir.
Niranjan Bhujel (view profile)
Do you mean complete inverter or fourth leg of inverter?
MOHD MARWAN NOOR DIN (view profile)
Done run sir Niranjan Bhujel. Btw, where can I put the voltage source inverter (VSI) with fourth leg ?
mahmood (view profile)
Could anyone send me the model in R2015a.
m.almahari@gmail.com
Niranjan Bhujel (view profile)
R2016b or latest.
MOHD MARWAN NOOR DIN (view profile)
what version of matlab can I use to run this simulation sir ?
khechiba kamel (view profile)
is not working in matlab 2013a it givs this Msg "object[id]!=NULL" how can i solve this pb thak's so much
muhammad fateh (view profile)
hai sir, would you please send this to me but in r2016a . my email @mhmmdfateh94@gmail.com
RENJITH R SHENOI (view profile)
RENJITH R SHENOI (view profile)
Could you please update how we are selecting the values of P and I for the PI controller for the calculation of Ploss ( Loss component) to keep the DC link voltage constant
RENJITH R SHENOI (view profile)
Venkat Montu.... I have forwarded 2016a file
venkat jagilinki (view profile)
thanks a lot sir...its working fine
montu vermah (view profile)
sir , please send it to my email....i cn't download it.
venkat jagilinki (view profile)
Sir, I have implemented the model but facing issue with respect to controller part. Could you please provide me the model in 2016a Version..it will be very helpful for me.
my mail id: venkat.jagilinki@gmail.com
Niranjan Bhujel (view profile)
Which file??? This file can be downloaded by clicking download button.
montu vermah (view profile)
sir ,pliz mail me this file..... mantukoiri1990@gmail.com
RENJITH R SHENOI (view profile)
Could you please update how we are selecting the values of P and I for the PI controller for the calculation of Ploss ( Loss component) to keep the DC link voltage constant.
RENJITH R SHENOI (view profile)
Niranjan :Could you please forward the same. Thank you
RENJITH R SHENOI (view profile)
Niranjan
Thank you . Mail id is renjith.shenoi@gmail.com
RENJITH R SHENOI (view profile)
Niranjan
Thank you . Mail id is renjith.shenoi@gmail.com
Niranjan Bhujel (view profile)
The fundamental freq of oscilation is twice the supply frequency. Supply freq is 50 Hz. So, it is 100 Hz. Check your mail for R2016a version of this file.
RENJITH R SHENOI (view profile)
Niranjan Bhujel.
Could you please pass the link of model in 2016a Version?
1) The fundamental freq of oscillation of instantaneous power is 100 HZ. Is that fixed?
Niranjan Bhujel (view profile)
RENJITH R SHENOI,
1) Total q needs to be compensated because we are compensating harmonics as well as reactive component of current.
2) There is no gain of -1.
3) The fundamental freq of oscillation of instantaneous power is 100 HZ. There is no component with frequency between dc(0 Hz) and 100 Hz. So, cut off frequency of 50 Hz works fine.
4) We need to calculate oscillating part but LPF gives constant part. So, constant part is subtracted from total part of p to obtain oscillating part. Ploss has been subtracted because it is to be drawn.
RENJITH R SHENOI (view profile)
Could you please provide a link to open the model 2016a version.
RENJITH R SHENOI (view profile)
Thank you, I have few more queries on this.
Vabc converted to VAlpha and Vbeta then P and Q calculated.P = Paverage + Pharmonic. Similarlly Q=Q average+ Qharmonic. In this model LPF is applied only to separate P so that P harmonic need to compensate.Iwould like to get some clarifications from the model.
So,
1) No need to filter Q? Total Q needs to compensate and directly taken for calculation of Ic1 and Ic2.
2) Why a gain of -1 is placed after Q?
2) Why Cut off frequency of LPF is 50? to separate the dc component it should be low like 5 Hz or 10 Hz. Comment please
3) After Low pass Filtering of P , DC component (Paverage ) is added with Loss component (Ploss )and subtracted from P to get AC component of real power. (Pharmonic) . So here it is coming negative right. Please comment on this.
mail : renjith.shenoi@gmail.com
Niranjan Bhujel (view profile)
Yes, q can be negative. It is instantaneous imaginary power. It is different than conventional reactive power. We are fitering p because only oscillating part of p is to be supplied by SAPF. So total p0,oscillating part of p and total q is used to calculate reference current.
RENJITH R SHENOI (view profile)
Hello, Voltage and current transformed to Alpha Beta co-ordinates and corresponding P and Q calculated. Will Q is negative?
Why we are filtering P and Q anf after filter out what value we will use for computation of reference current? Thank you
Vasu The King Maker (view profile)
How to control voltage and frequency for parallel operated generators
Mail *: thekingmakervasu@gmail.com