I have follwing curve resulting from plotting current in capacitor versus time. (Matlab code is attached).
I want to calculate the charge stored in the capacitor. For that I need to calculate the positive area of the curve. As you can see in the zoomed in image, the current plotted has pulses (switching frequency = 20kHz), and I would like to calculate the charge stored in one fundamental cycle (fundamental frequency = 60Hz).
So theoretically, the capacitor is getting charged when the current flows into it. For the curve it means the value of current that is above zero(positive).
I tried calculating the area of the curve by simplifying the curve into traingles and trapezoids, and then taking the coordinate points values (x,y). But I am pretty sure it is not correct since it does not take into account the presence of pulses.
I also tried the trapz(x,y) in matlab but the calculation does not seem right.
And the expression of the current is not a straightforward function, it is defined in terms of switching functions and duty ratios so I am not sure how to use the integral function in Matlab.