cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Community Tip - Did you know you can set a signature that will be added to all your posts? Set it here! X

Translate the entire conversation x

Total Harmonic Distortion of AC current using Fourier Analysis

cadtelsim
10-Marble

Total Harmonic Distortion of AC current using Fourier Analysis

So, I've been able to work through all the values so far but seems like I have a problem with the summation notation since the THD of alpha does not converge to a solution:

 

VM_10380528_0-1681678396577.png

 

VM_10380528_1-1681678415655.png

 

VM_10380528_2-1681678434105.png

 

 

ACCEPTED SOLUTION

Accepted Solutions
Werner_E
25-Diamond I
(To:cadtelsim)

Not sure why, but the problem is "only" n=2 and function I.a

 

Werner_E_0-1681692024874.png

The function you integrate does not look malicious

Werner_E_1-1681692585902.png

Some additional observations:

Werner_E_6-1681693559643.png

 

 

EDIT:

Solution - Set the system variable TOL to a lower value. Either write TOL:=10^-8 at the top of the sheet in front of the function definitions or use the menu.

Doing so makes the integral converge

Werner_E_3-1681692804476.png

 

The plot below with variable values of alpha takes very long. I stopped it and used a larger step width to speed it up.

It looks like THD is independent from alpha

Werner_E_4-1681693179472.pngWerner_E_5-1681693497267.png

 

View solution in original post

34 REPLIES 34
ppal
17-Peridot
(To:cadtelsim)

File?

cadtelsim
10-Marble
(To:ppal)

I get a notification saying I can't post the file using the email, is there any way to attach the file here?

ppal
17-Peridot
(To:cadtelsim)

ppal_0-1681682107484.png

Drag and drop.

cadtelsim
10-Marble
(To:ppal)

 
Werner_E
25-Diamond I
(To:cadtelsim)

Not sure why, but the problem is "only" n=2 and function I.a

 

Werner_E_0-1681692024874.png

The function you integrate does not look malicious

Werner_E_1-1681692585902.png

Some additional observations:

Werner_E_6-1681693559643.png

 

 

EDIT:

Solution - Set the system variable TOL to a lower value. Either write TOL:=10^-8 at the top of the sheet in front of the function definitions or use the menu.

Doing so makes the integral converge

Werner_E_3-1681692804476.png

 

The plot below with variable values of alpha takes very long. I stopped it and used a larger step width to speed it up.

It looks like THD is independent from alpha

Werner_E_4-1681693179472.pngWerner_E_5-1681693497267.png

 

Great Werner_E! Well done! Congratulation

walter

Its not that difficult to track down errors like this using "try..on error".

Werner_E_1-1681737653957.png

It was surprising that just the single value n=9 threw the error and all others worked fine.

The usual way would be to increase the value of CTOL and see if convergence could be achieved that way. Not sure why I had the idea to decrease TOL 😉

 

Thank you, Werner for the tips and solutions. It wouldn't have occurred to me. I'm a new mathcad user and despite having serious issues with it every day I'm more interested in it than MATLAB. Any suggestions for improvement?


Vina Mofidi

Sch Elect Comptr & Energy Engr

Arizona State University


email: Vmofidi@asu.edu



Thank you for the tip. Its interesting to find out why value of n=9 creates an error. I appreciate your help.

Werner_E
25-Diamond I
(To:cadtelsim)


@cadtelsim wrote:

Thank you for the tip. Its interesting to find out why value of n=9 creates an error. I appreciate your help.


Yeah, if you find out, let us know.

ttokoro
21-Topaz I
(To:cadtelsim)

image.png

ttokoro_0-1681790832948.png

 

image.pngimage.png

ttokoro_1-1681790882387.png

 

image.pngimage.png

ttokoro_2-1681790920966.png

 

image.pngimage.png

Thank you for the detailed analysis. This complete analysis will be instrumental in learning features of Mathcad along with the limitations of THD of ac current related to the delay angle. I'll fully study this document.

image.png

ttokoro
21-Topaz I
(To:ttokoro)

To calculate correctry, TOL is sensitive and integral from -L to L gives bug error at 80 to 100 degree.

Divide it from -L to 0 and from 0 to L makes good answer. Prime 9.

 

image.pngimage.pngimage.pngimage.pngimage.pngimage.png

Thank you for the detailed response. I will do my best to replicate this to understand Mathcad features better. 

MP_10884228
6-Contributor
(To:ttokoro)

ttokoro, could you tell me please what is the meaning of the signal On(t, alfa) and FC(f,N,L) ?

On(t,alpha) is unit step function to make rectifired signal(t) at angle alpha. It is pink plot below to make blue voltage signal.

FC(f,N,L) is a function in Mathcad help or tutorial of Fourier transform function. f is waveform, N is evaluated harmonics number and L is a half of fundamental period T.

 

image.png

MP_10884228
6-Contributor
(To:ttokoro)

I am using Prime 9, and I looked under the FFT examples provided under Help in Resources menu, but I couldn't find it.

sorry, I found it eventually.

MP_10884228
6-Contributor
(To:ttokoro)

Hi ttokoro,

 

Could you tell me please, how did you fill in the area under the On(t,alfa) graph with color?

 

Just change the x-y plot format to bar type.

MP_10884228
6-Contributor
(To:ttokoro)

Hi ttokoro,

 

I would like to know one more thing. Could you tell me please, why is it needed to multiply the input ac voltage with the step function and then to divide it to the Rload for calculating the ac supply current, and not just dividing the ac input voltage to the Rload?

 

Thank you,

Ofcourse if you using R only and not includes L or C, the shape of i(t) is as same as v(t). Then the  Harmonic Distortion of AC current is as same as Voltage waveform. But to find the Total Harmonic Distortion of AC current, it must use Fourier Analysis. And it is very good question to calcilate RMS value of wave signals with using both 

image.pngimage.pngimage.pngimage.png

MP_10884228
6-Contributor
(To:ttokoro)

Could you tell me please where is it defined phi (t) (greek name)l signal that is used for calculating the signal On(t,alfa)?

Could you share please the definition of  phi(t)?

 

MP_10884228_0-1740573306986.png

Thank you.

LucMeekes
23-Emerald III
(To:MP_10884228)

The functionLucMeekes_0-1740573324234.pngis the Heaviside function, or unit step function. It is predefined in Mathcad and Prime. It results in 0 for t<0, and 1 for t>0. The value for t=0 is debatable, that is, the result differs depending on who you ask:

Numerics in Prime:LucMeekes_1-1740573400203.png

Numerics in Mathcad 11:LucMeekes_3-1740573524132.png

I think the symbolic processor in Mathcad 11 gives the best value:LucMeekes_4-1740573593040.png

 

Success!
Luc

 

 

 

 

 

MP_10884228
6-Contributor
(To:LucMeekes)

Hello Luc,

 

For the Heaviside Step Function, it is written that to type Φ, press F, Ctrl+G, which I followed, and the Mathcad tool is saying that Φ variable is undefined.

Could you tell me please what I did wrong?

 

Thank you

LucMeekes
23-Emerald III
(To:MP_10884228)

You need capital phi

Type an F, then [CTRL-G]

 

Success!
Luc

MP_10884228
6-Contributor
(To:LucMeekes)

it works, thank you!

MP_10884228
6-Contributor
(To:LucMeekes)

I have one more question, related to the definition of the variables using index notation.

Could you tell me please, how could be done this without predefining the index values, like for Vrms, Tcycle, and so on?

Announcements

Top Tags