Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: Solution help

  1. Apr 30, 2010 #1
    Could someone help me to solve the equation below?

    cos(wt).exp(jθ)

    I want to find something like

    cos(wt+θ)

    thanks from now on
     
  2. jcsd
  3. Apr 30, 2010 #2

    berkeman

    User Avatar

    Staff: Mentor

    Welcome to the PF. That is not an equation, since there is no equal sign "=".

    Do you mean you have a function like that, and wand want to change the form? Is there a particular question associated with this? And exponents with j in them generally mean a complex number. Are you wanting to find the conversion from polar to rectangular form for complex numbers?
     
  4. May 1, 2010 #3
    I mean i want to provide a form like cos(wt+θ). to get rid of exp(jθ) part.
     
  5. May 1, 2010 #4
    exp(jθ) = cos(jθ) + i sin(jθ)
     
  6. May 1, 2010 #5
    No I know but I want to learn what cos(wt).exp(jθ) equals in form of sin or cosine fonk.
     
  7. May 1, 2010 #6
    I mean cos(wt) multiplied by exp(jθ)
     
  8. May 1, 2010 #7

    berkeman

    User Avatar

    Staff: Mentor

    Well, it's still pretty hard to understand what you are looking to do, but how about this:

    [tex]e^{j\theta}[/tex] is a complex number, with Real and Imaginary components as suggested above. You can picture a complex number in a 2-D plane, with the Real axis pointing to the right, and the Imaginary axis pointing up (kind of like x,y axes).

    The function you are wanting to multiply that complex number by is real (as far as I can see), so it would only scale the complex number. Since it varies in time, the scaling of the complex number would vary with time. If theta is constant, then the complex number is constant, so the complex vector (from the origin of the Real and Imaginary axes to the point exp(jθ) ) would just grow and shrink with cos(wt). If theta varies with time as well somehow, then the complex vector would rotate, at the same time cos(wt) modulates its amplitude.

    You would need to definethe relationship between theta and omega further to get a more complete answer.
     
  9. May 2, 2010 #8

    collinsmark

    User Avatar
    Homework Helper
    Gold Member

    Here is something that might help.

    [tex] cos(\theta) = \frac{e^{i \theta} + e^{-i \theta}}{2} [/tex]

    [tex] sin(\theta) = \frac{e^{i \theta} - e^{-i \theta}}{2i} [/tex]

    (1) First, convert your sin() or cos() function to the appropriate exponential above.

    (2) Perform your multiplication with the existing exponential, noting that

    [tex] e^a e^b = e^{a + b} [/tex]

    (3) convert back to sin() or cos().

    This should give you what you are looking for.

    [Edit: Hint: After you do the above once, you'll likely find that you can convert something of the form [itex] cos(\omega t)e^{j \theta} [/itex] to the form [itex] cos(\omega t + \theta) [/itex] in your head forever after. It turns out to be a rather simple conversion. 'Trivial in fact. :wink:]
     
    Last edited: May 2, 2010
  10. May 2, 2010 #9
    Let’s make it easier if you like

    cos(θ1).exp(iθ2)

    where θ1 and θ2 are constants in rad.

    As you recommended I found in the following expression

    [Exp(i(θ1+ θ2)) + Exp(i(θ2- θ1))] /2

    İf I try to turn back to the cosine form I can not provide something like

    cos(θ1+ θ2) or cos(θ1- θ2) or cos(θ1- θ2+ PI/2)....

    OR

    sin(θ1+ θ2) or sin(θ1- θ2) or sin(θ1- θ2+ PI/2)....

    This is from topic of electric physics .

    You have a voltage source as cos(wt) or cos(θ1) w=2.PI.f t: time

    I have a low pass filter that is described as

    exp(iθ2) (say transfer fonction in sinusoidal continuous form)

    which is to make θ2 rad (or degrees) shift in phase with respect to input signal say cos(wt)

    What am I going to see on scope at output ?=How the output will be seen in time domain ?

    I am interested in reel part of the output

    If I use MATLAB with simulink I can see a output voltage wave form that is shifted by θ2 with respect to input signal. (Lets ignore change in amplitude)

    With best regards
     
  11. May 5, 2010 #10

    collinsmark

    User Avatar
    Homework Helper
    Gold Member

    Sorry for the late reply, but I've been giving this thread some thought.

    The best I can surmise, is that you might be incorrectly mixing terms in the time domain and the frequency domain. Terms such as [itex] cos(w_0 t) [/itex] typically are found in the time domain. Terms like [itex] e^{j \theta} [/itex] are typically found in the frequency domain.

    Unfortunately, I didn't realize that until after I replied, and I think I might have just made the confusion worse. Sorry about that. Let me start over.

    Suppose you have a time domain signal,

    [tex] v_s(t) = Acos(\omega_0 t) [/tex]

    Call its frequency domain representation [itex] V(\omega _0) [/itex]. And suppose that in the frequency domain (sometimes called phasor representation), you determined the frequency response of some other signal to be [itex] V_c(\omega)/V_s(\omega) = H(\omega) = e^{j \theta} [/itex]. In other words, suppose that you have determined that

    [tex] V_c(\omega) = V_s e^{j \theta}. [/tex]

    When converting this to the time domain, the result is simply,

    [tex] v_c(t) = Acos(\omega_0 t + \theta), [/tex]

    when working with sinusoidal steady-state analysis. That's what I meant by trivial in my last post. :wink:

    In general, when performing sinusoidal steady-state analysis (ignoring transients like when you close a switch or something), if you have

    [tex] Y(\omega) = X(\omega)H(\omega) [/tex]

    then

    [tex] |Y| = |X||H| [/tex]

    and

    [tex]\angle Y = \angle X + \angle H[/tex]

    Perhaps this is a good time to show a concrete example. Suppose you have an AC voltage source connected to a series RC circuit, and you want to measure the voltage across the capacitor.

    You can draw the circuit the normal way, but it makes it a little easier to draw the circuit in the phasor representation, since that makes it automatically ready for frequency domain analysis.

    attachment.php?attachmentid=25619&d=1273113361.gif

    Finding the sinusoidal steady-state current is easy enough. We just divide the voltage by the total impedance,

    [tex] I(\omega) = \frac{V_s}{R + \frac{1}{j \omega C}} [/tex]

    Multiply the current by the impedance of the capacitor to find the voltage we are looking for,

    [tex] V_c(\omega) = \frac{V_s \frac{1}{j \omega C}}{R + \frac{1}{j \omega C}} = \frac{V_s \frac{1}{j \omega C}}{\frac{1+ j \omega RC}{j \omega C}} = \frac{V_s}{1 + j \omega RC} [/tex]

    The transfer function is then,

    [tex] H(\omega) = \frac{V_c(\omega)}{V_s(\omega)} = \frac{1}{1+ j \omega RC} [/tex]

    The denominator is a complex number. It can be written in polar form.

    [tex] H(\omega) = \frac{1}{\sqrt{1 + (\omega RC)^2}e^{j \ atan(\frac{\omega RC}{1})}} = \frac{1}{\sqrt{1 + (\omega RC)^2}}e^{-j \ atan(\omega RC)} [/tex]

    The magnitude and phase of H are,

    [tex] |H| = \frac{1}{\sqrt{1 + (\omega RC)^2}} [/tex]

    [tex] \angle H = \theta_H = - atan(\omega RC) [/tex]

    Earlier I contended that if we want to find the sinusoidal steady-state analysis vc(t), knowing that [itex] v_s(t) = Acos(\omega_0 t) [/itex], it is simply

    [tex] v_c(t) = A|H|cos(\omega_0 t + \theta_H) [/tex]

    And I think that's what your question is about. Correct me if I'm wrong, but I think you were originally asking, 'where does the phase shift come from -- mathematically speaking?' To answer that, we need to convert [itex] v_s (t) = cos(\omega_0 t) [/itex] to the frequency domain -- not just its symbol, [itex] V_s(\omega) [/itex] -- but the whole thing. And we can do that using the Fourier transform. This is where it gets fun (or if you wish, phun :biggrin:). I suggest you work through this at least once. You may never have to do it again for sinusoidal steady-state analysis, but going though it will answer your questions.

    We define the Fourier transform and inverse Fourier transform as,

    [tex] X(\omega) = F\{x(t)\} = \int _{-\infty} ^{\infty} x(t) e^{-j \omega t} dt [/tex]

    [tex] x(t) = F^{-1} \{X(\omega)\} = \frac{1}{2 \pi}\int _{-\infty} ^{\infty} X(\omega) e^{j \omega t} d\omega [/tex]

    Please note that these are not the only way to define the Fourier transform and inverse! In another way to define them, instead of the [itex] 1/2 \pi [/itex] in the inverse (like I used), a [tex] 1/\sqrt{2 \pi} [/tex] is used in both the forward and inverse. In yet another way to define them, no terms are to the left of the integral, and a [itex] 2 \pi [/itex] is in the exponential terms. The way I defined them here is very common in electrical engineering courses. But it is not absolute. I highly suggest using the definition that your textbook/coursework uses, whatever that may be.

    The Fourier trasform of [itex] v_s(t) [/tex] in our circuit is,

    [tex] V_s(\omega) = \int _{-\infty} ^{\infty} A \ cos(\omega_0 t)e^{-j \omega t} = \pi A \left( \delta(\omega + \omega_0) + \delta(\omega - \omega_0) \right) [/tex]

    where [itex] \delta() [/itex] is the Dirac delta function (sometimes called the impulse function, and sometimes simply the delta function). Now you might be asking, "where in the world did he get that?" Honestly, the above integral is a little tricky to do explicitly, and it involves some reformulation of the integral, then taking some limits and such. But if you want a weak proof, take the inverse Fourier transform of [itex] \pi A \left( \delta(\omega + \omega_0) + \delta(\omega - \omega_0) \right) [/itex] and you'll find it easily produces [itex] A \ cos(\omega_0 t) [/itex].

    Recall

    [tex] V_c(\omega) = H(\omega)V_s(\omega) [/tex]

    [tex] = \frac{1}{\sqrt{1 + (\omega RC)^2}}e^{-atan(\omega RC)} \pi A \left( \delta(\omega + \omega_0) + \delta(\omega - \omega_0) \right) [/tex]

    Thus,

    [tex] V_c(\omega) = \frac{\pi A}{\sqrt{1 + (\omega RC)^2}} \left( \delta(\omega + \omega_0) + \delta(\omega - \omega_0) \right) e^{-j \ atan(\omega RC)}[/tex]

    Now let's take the inverse Fourier transform.

    [tex] v_c(t) = \frac{1}{2 \pi} \int _{-\infty} ^{\infty} \frac{\pi A}{\sqrt{1 + (\omega RC)^2}} \left( \delta(\omega + \omega_0) + \delta(\omega - \omega_0) \right) e^{-j \ atan(\omega RC)}e^{j \omega t} d \omega[/tex]

    [tex] = \frac{A}{2} \int _{-\infty} ^{\infty} \frac{1}{\sqrt{1 + (\omega RC)^2}} \left( \delta(\omega + \omega_0) + \delta(\omega - \omega_0) \right) e^{j \left( \omega t - atan(\omega RC) \right)} d \omega[/tex]

    Noting that [itex] \int f(x)\delta(x - x_0)dx = f(x_0) [/itex], the integral evaluates to

    [tex] v_c(t) = \frac{A}{2 \sqrt{1 + (\omega_0 RC)^2}} \left(e^{j(\omega_0 t - atan(\omega_0 RC))} + e^{-j(\omega_0 t - atan(\omega_0 RC))} \right) [/tex]

    [tex] = \frac{A}{\sqrt{1 + (\omega_0 RC)^2}} cos \left( \omega_0 t - atan(\omega_0 RC) \right) [/tex]

    thus for sinusoidal steady-state analysis,

    [tex] v_c(t) = A|H|cos(\omega_0 t + \theta_H) [/tex]

    And that's where the phase shift comes from! :biggrin:
     

    Attached Files:

Share this great discussion with others via Reddit, Google+, Twitter, or Facebook