18
\$\begingroup\$

The picture below shows a RLC circuit. A RLC circuit is an electrical circuit consisting of a resistor (R), an inductor (L), and a capacitor (C), connected in series or in parallel. (1)

enter image description here

In order to simplify computations, it's common to work in the frequency (Laplace) domain instead of the time domain.

Your task is:

Take the values R, L and C as input, and return the voltages VR, VL and VC

The conversion to the Laplace domain is as follows:

R = R
XL = j*w*L      // OK, XL = w*L, and ZL = j*XL, but don't mind this here.  
XC = 1/(j*w*C)  // I haven't ruined physics, it's only a minor terminology tweak

where j = sqrt(-1), and w = 2*pi*50 (The frequency is 50 Hz).

The combined impedance, when the components are in series is Z = R + XL + XC. You might remember U = R*I from high school physics lectures. It's almost the same, but a bit more complex now: VS = Z*I. The current is calculated by dividing the voltage VS by the total impedance Z. To find the voltage over a single component, you need to know the current, then multiply it by the impedance. For simplicity, the voltage is assumed to be VS = 1+0*j.

Equations you might need are:

XL = j*w*L
XC = 1/(j*w*C)
Z = R + XL + XC   // The combined impedance of the circuit
I = VS / Z         // The current I (Voltage divided by impedance)
VR = I * R        // Voltage over resistance (Current times resistance)
VL = I * XL       // Voltage over inductor (Current times impedance)
VC = I * XC       // Voltage over capacitor (Current times impedance)

The input is from either STDIN or as function arguments. The output/result must be three complex numbers, in a list, string or whatever is most practical in your language. It's not necessary to include names (ex VR = ...), as long as the results are in the same order as below. The precision has to be at least 3 decimal points for both the real and imaginary part. The input and output/results can be in scientific notation if that's default in your language.

R and L are >= 0, and C > 0. R, L, C <= inf (or the highest possible number in your language).

A simple test case:

R = 1, L = 1, C = 0.00001

VR = 0.0549 + 0.2277i
VL = -71.5372 +17.2353i
VC = 72.4824 -17.4630i

For the results above, this could be one (of many) valid ouput format:

(0.0549 + 0.2277i, -71.5372 +17.2353i, 72.4824 -17.4630i)

Some valid ouput formats for one voltage value are:

1.234+i1.234,   1.23456+1.23456i,   1.2345+i*1.2345,   1.234e001+j*1.234e001.

This list is not exclusive, so other variants can be used, as long as the imaginary part is indicated by an i or a j (common in electrical engineering as i is used for current).

To verify the result for other values of R,L and C, the following must be true for all results: VR + VL + VC = 1.

The shortest code in bytes win!

By the way: Yes, it's voltage over a component, and current through a component. A voltage has never gone through anything. =)

\$\endgroup\$
  • 2
    \$\begingroup\$ Actually, reactances are real numbers, so XL = omega*L. The impedance of the inductor is Z = jXL. (This does not affect the problem, it's only a correction) \$\endgroup\$ – Voitcus Sep 1 '15 at 10:09
  • \$\begingroup\$ @Voitcus, true... I simplified it a bit, to not make the question too confusing. I included the j in the XL/XC terms, when going to the frequency domain. I never said that the reactance was complex though (although I called it X, and not jX) =) But I agree with you! I actually called it impedance too. \$\endgroup\$ – Stewie Griffin Sep 1 '15 at 10:18
  • \$\begingroup\$ Can I take a list of 3 numbers as a function input, or does it have to be 3 separate arguments? \$\endgroup\$ – Martin Ender Sep 1 '15 at 10:37
  • \$\begingroup\$ @MartinBüttner, list is OK. \$\endgroup\$ – Stewie Griffin Sep 1 '15 at 10:57
7
\$\begingroup\$

Pyth, 30 29 28 bytes

L*vw*100.l_1)K[Qy0c1y1)cRsKK

Try it online.

\$\endgroup\$
9
\$\begingroup\$

Mathematica, 33 bytes

So close to Pyth...

l/Tr[l={#,#2(x=100Pi*I),1/x/#3}]&

This is an unnamed function, which takes R, L and C as its three arguments and returns a list of complex numbers as the result (in the required order VR, VL, VC). Example usage:

l/Tr[l={#,#2(x=100Pi*I),1/x/#3}]&[1, 1, 0.00001]
(* {0.0548617 + 0.22771 I, -71.5372 + 17.2353 I, 72.4824 - 17.463 I} *)
\$\endgroup\$
3
\$\begingroup\$

Octave / Matlab, 53 51 bytes

function f(R,L,C)
k=-.01j/pi;Z=[R L/k k/C];Z/sum(Z)

Try it online

Thanks to @StewieGriffin for removing two bytes.

\$\endgroup\$
  • \$\begingroup\$ @StewieGriffin 100j?! So many years using Matlab and I didn't know that could be done! :-) (I did know 1j, but I thought it was just that). Thanks! \$\endgroup\$ – Luis Mendo Sep 1 '15 at 16:06
  • \$\begingroup\$ ... aaaand: Apparently I know more about you than you do! Because, you do/did know that it is possible! =) \$\endgroup\$ – Stewie Griffin Sep 1 '15 at 16:10
  • \$\begingroup\$ @StewieGriffin Ooooh. It happened to me again. Bad memory!!:-D (I never really use that notation) \$\endgroup\$ – Luis Mendo Sep 1 '15 at 16:11
  • \$\begingroup\$ You can save another byte if you start with the inverse of k, like this: k=-.01j/pi;Z=[R,L/k,k/C];Z/sum(Z), or k=-.01j/pi;[R L/k k/C]/(R+L/k+k/C). =) \$\endgroup\$ – Stewie Griffin Sep 2 '15 at 9:02
  • \$\begingroup\$ @StewieGriffin Good idea! Edited \$\endgroup\$ – Luis Mendo Sep 2 '15 at 9:34
3
\$\begingroup\$

APL (Dyalog Unicode), 27 24 bytesSBCS

Full program. Prompts for C, L, R in that order.

(⊢÷+/)(⎕,⎕∘÷,÷∘⎕)÷○0J100

Try it online!

0J100 100 i

 π times that

÷ reciprocal of that

() apply the following tacit function:

÷∘⎕ divide the argument by input (C)

⎕∘÷, prepend input (L) divided by the argument

⎕, prepend input (R)

() apply the following tacit function:

+/ sum the arguments

⊢÷ divide the arguments by that

\$\endgroup\$
  • \$\begingroup\$ @StewieGriffin I am note sure what you mean. the "high minus" ¯ is APL's negative number prefix, to distinguish from the function (i.e math operator) negate -. Anyway, it wouldn't be fair to not count APL chars as single bytes, it is just a matter of encoding, and there are many APL systems that use single bytes to store APL code. E.g. Dyalog has both Unicode and Classic (single-byte) versions of their interpreter. \$\endgroup\$ – Adám Sep 3 '15 at 7:18
  • 1
    \$\begingroup\$ I agree, if you have used encoding where every character is a single byte, then the number of chars should equal the number of bytes. Can you verify that it's the case (i'm not too familiar with apl an different encoding systems). Also, i wasn't familiar with the high minus sign. My bad... \$\endgroup\$ – Stewie Griffin Sep 3 '15 at 8:17
  • \$\begingroup\$ I just pasted the code into a "count number of bytes" box and got back 31. If it's not correct, then of course you'll have a score of 28 :-) Although, ..,49J¯17.4.. would mean the first part is imaginary and the second is real in any other language (or in mathematical notation in general), so it might violate the rule "as long as the imaginary part is indicated by an i or a j". Have a +1 for teaching me about "high minus", and a nice answer, but I'm not sure I can pick it as the accepted answer, when that day comes. \$\endgroup\$ – Stewie Griffin Sep 3 '15 at 8:28
  • 1
    \$\begingroup\$ @StewieGriffin Ninja'd ;) \$\endgroup\$ – Beta Decay Sep 3 '15 at 8:31
  • \$\begingroup\$ @StewieGriffin See help.dyalog.com/14.1/Content/Language/System%20Functions/… \$\endgroup\$ – Adám Sep 3 '15 at 8:40
2
\$\begingroup\$

Octave, 41 bytes

@(R,L,C)(Z=[R L/(k=-.01j/pi) k/C])/sum(Z)

1/(100*j*pi) can be shortened to -.01j/pi which is a lot shorter. By assigning it to the variable k inline, the variable can be used twice. Assigning the entire vector to the variable Z costs 4 bytes, but allows us to divide by sum(Z), which is 5 bytes shorter than (R+L/k+k/C).

\$\endgroup\$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.