Fancy3.mws


vector stuff

op

electric field plot

envelope detector

driven oscillator


E&M Packet Problem ABW solution Potential and field plots of three point charges Learning objectives: Translation of standard 1/r potential to off-origin points Physical intuition regarding potential Identify equilibrium points Relate potential plot and electric field vectors MAPLE: plot3d, grad, fieldplot commands

arbitrary units are used throughout in Coulomb's law

Set plot range; use non-even values to avoid singular points in grids

3d countour plot of potential function

Get electric field from negative gradient of potential

(Note: "-grad(U)" doesn't seem to work, so using "grad(-U)"

Create new vector field of unit vectors E/|E|

Compress dynamic range of E-vector lengths by taking ln(1+|E|)
log is taken twice to put all lengths within range of fieldplot command


Plot the phase and amplitude response of a driven, damped oscillator as a function of frequency.

This problem illustrates dealing with complex numbers, and the use of the 'op' command to parse an output expression.

It is interesting that very heavy damping (b=10, k=10, m=1) crushes the response away from zero frequency.

Phase of the response as a function of frequency. The conventional definition is the phase of the driver minus the phase of the driven system. The relative amplitude response will be taken as exp(-I alph).

The driver is in phase with the driven system at very low frequencies, and pulls ahead until it is Pi/2 ahead at resonance, and goes on to be Pi ahead at very high frequencies.

Amplitude response as a function of frequency


In AM (amplitude-modulated) radio, the amplitude of a high-frequency signal is modulated by the lower-frequency audio information. AM radios use 'envelope detectors' involving RC circuits in order to remove the high-frequency (the 'carrier wave') and leave the audio signal.

This problem gives you a simple modulated signal. It also gives you the 'rectified' signal we would see after passing through a diode (which lets current flow in one direction only.).

You are to pass the signal through a series RC circuit, first through a resistor R and then a capacitor C to ground. The 'envelope' signal is to be the voltage across the capacitor.

Letting v be the incoming signal, the equation to be solved is v -IR -q/C = 0 .

You must replace the current I by a function of q, the charge on the capacitor. This (call it eq1) will be the generic equation to be solved for q the charge on the capacitor as a function of time. By substituting values of R and C into eq1, you get eq2. This new equation is to be solved numerically for q (the charge on the capacitor) as a function of time.

a) Get the rectified signal first, by using the Heaviside function. The input signal is

b) Set up the differential equation, using q(t) for the charge on the capacitor C . C is connected to ground and to the resistor R . R has the input voltage on one side and C on the other. [The current will be dq/dt]

c) Substitute values into the DE to make a new DE which will be solved numerically. Start off by setting RC equal to 1/10 the period of the 'carrier frequency'.

d) Solve the new DE numerically. Plot (use odeplot) and see what this gives you for an envelope function Compare this to the plot of the modulated function.

From here on down, the students should be doing the work.

RC is 0.45 ms. carrier pd= 0.17 ms. modulation pd .= 20 ms.

Note that we don't get all of the signal back, nor is the shape a perfect fit to the envelope.


End of Fancy3.mws