OK cereb,
Let's take it small in steps 
meaning: do not consider this to be the answer.
It is the opening step of a dialog.
Do not worry too much about the the length of the dialog.
Obviously, it has to be finite;
(This is one example of one 'too long' dialog:
http://answers.google.com/answers/threadview?id=500730
)
As you see, before the dialog gets too long, I do usually say that we have to
come to closure in few more installments.
The reason I want to go in small steps is to make sure we are communicating.
The problem with the magnet
______________________________
Waving the magnet around the coil, is a legitimate 'forcing function'
or 'stimulus'. We are able to model electromechanical systems, such as
electric motors, where physical motion and induced EMF interact. The problem
with your original description was that it was too undefined: The induced of EMF
depends on speed, and the speedprofile of the motion, which determines
the shape of voltage induced in L1 was not well defined.
This is illustrated in these two links:
Induced e.m.f. = depends on speed
http://courses.science.fau.edu/~rjordan/rev_notes/28.2.htm
shape of e.m.f. pulse induced in coil by a magnet
http://www.dataharvest.co.uk/datalogging/easysense_fa.html
Also, such a situation may require a nonlinear model if the distance varies
in a wide range. Again, we can solve nonlinear models, but here, in what I
take to be a conceptual stage, it is better to make the model as simple as
possible (but 'not more simple than that :)
http://www.quotationspage.com/quote/3430.html
So, following that adage, let's start really simply  we can add componets later:
Description of the task
_______________________
You do have a stimulus (or forcing function) which is the
'rate of pay as a function of time'.
We consider two rectanglual pulses with same area of 2000 units.
One is narrow and tall (the Xmass bonus), other is wide and low (a raise).
We will call them Rx(t) and Rr(t).
You have 'a worker' (what a control engineer would call 'a plant')
which produces output P(t) (measured in some output units).
The 'worker' is modelled by a linear system K, (kernel)
so that the output (of production) P depends on rate of pay:
P(t) = Sum K(t  p) R(p)
Sum is over p ,ranging from infinite past up to present = t.
meaning: the current rate of production P(t) depends not just on the current
rate of pay R(t), but also on the
last week R(t1),
and week before that R(t2)... etc
and the 'memory of the past' fades at the rate controlled by some
'time constant' tau= R *C .
We know that K can be represented by a circuit (mechanical or electrical)
but whatever circuit we choose, we want to be able, for the given K
and the stimulus R we choose, to obtain a plot of P(t).
In addition, we want to represent K by some simple circuit with small
number of components, so we can tweek the values, obtain a modified K,
and perhaps give some interpretation to their function in terms of the
psychology of the worker.
Questions (RFC)
1) First question: Is this, as described above, the task?
An ideal tool for this would be Scilab (a free simulation paltform) which can
produce plots like this:
http://www.scilab.org/doc/demos_html/node179.html
as described here
http://www.scilab.org/
Commercial programs like
http://www.wolfram.com/
ttp://www.math.ufl.edu/help/matlabtutorial/
or free version of matlab,
http://www.octave.org/
would also do the job.
What job? Job b) using the Kernel K to calculate Response P and plot it
and a) Calculating the K from a diagram of the circuit.
It is a bit more than a plotting program, but not much more.
2) Question two is, do you have some such program installed,
or do you want to install it?
It is not necessary to install one, but if you want to work interactively:
change the component and get a response, interactively, it would
make the task easier.
Alternative would be to start with the plotting applet, and some simple
Kernel, someting like K(d) = e + g* expt( r * d)
and get sample responses to different stimuli ( Rx and Rr ..)
P.S. r = 1/tau can be complex number  in which case we get an oscilatory
response  indication of 'inertia' or inductor in the circuit.
3) Are you familiar with complex numbers, stuff like
http://hyperphysics.phyastr.gsu.edu/Hbase/cmplx.html
http://planetmath.org/encyclopedia/ComplexExponentialFunction.html
http://www.mathworks.com/access/helpdesk/help/toolbox/dspblks/complexexponential.html
Do you want to use or avoid that terminology?
Hedgie 
Request for Answer Clarification by
cerebga
on
18 Mar 2006 12:07 PST
Thorough job, Hedgie, Thank you again.
I made curve generation my principle goal, and without having yet
generated such curve I trust that I can do so with the links etc. you
gave to me. I will download the simulation software and see what I
can do with it to generate a plot shape of the type I wish to
generate, and will then be ready in case more interaction is
indicated.
I knew that the magnet approach would be cumbersome to actually
generate the desired curve on a scope, for example. The "real" issue
for me is this: I have invested 20 years of work related to
human/social/societal behavior. It has been rather novel work,
seminal in parts, and at age 73 I feel the pressure to present it so
it can be understood by many people who search for credible solutions
to seemingly insurmountable obstacles to change the way we manage the
affairs within and among societies.
Since we have only "effective objectivity through collective
subjectivity" to rely on, such as laws, treaties and such, my aim was
to find a way to overcome the constraints of humans being able to be
objective about humans, and I searched for an objective frame of
reference to overcome these constraints. In that I suceeded,
presumptuous as that may sound.
I use phenomenology for explaining the method I employed. Social
scientists could use some of my findings employing expert systems and
AI tools to simulate policy formulation based upon the core model for
individual behavior and its use for a societal core model that uses
the individual core model. But I am more interested in making my work
plausible to many people who have valid concerns for our future, even
make it credible and mainly actionable. So I felt the need to model
much of what I found.
One of the most consequential Findings of my work relates to our
individual as well as collective behavior of "growth" that applies in
economics, in governance and the exponentially rising expansion or
even race for scientific progress. The curve I am after now is but
a start, subsequent graphical illustrations will further illustrate
what drives this obsession with "growth", which in turn creates the
many problems we face everywhere that derive from this obsession.
I can use your approach to model that and other curves we may want to
deal with hereafter, because it may well serve those who want to model
things I describe in my work, and they would demand more rigor in my
reporting. I was hoping merely for something to illustrate to
nonprofessionals something they could use to better relate to these
graphs. Such approach would also more adequately reflect my
inabilities to offer more than "cute" analogies.
Certainly our present task would be completed by using your approach
to create scaleable graphs, and it would be completed upon generating
them. Maybe I will use my LCR/magnet example to just explain verbally
how such curves "could be" generated that way but because of the
modelling requirements for what I am describing, that your approach
was necessary.
This verbiage here may help to explain my desire for the LRC/magnet
model and it may help when we deal with the other curves related to
this "growth" obsession.
I think we are where we wanted to be, and I will now get the
simulation program and see what I can do with all your inputs. Above
all, I want to thank you for having stayed with this issue patiently
and constructively. If you wanted to create two curves to show
scalability using the simulation software link, that would be the
answer that closes this part. I would hope that we can then deal with
the use of that approach to create the other graphs I have not yet
described but that rely upon this one.

Clarification of Answer by
hedgiega
on
19 Mar 2006 22:52 PST
cerebga
http://answers.google.com/answers/threadview?id=705898
I = V / ( R1 + R2)
http://www.pa.uky.edu/~phy211/graph_applets/plot_graph.html
current
http://answers.google.com/answers/threadview?id=708116
I made curve generation my principle goal, and without having yet
generated such curve I trust that I can do so with the links etc. you
gave to me.
OK.Let's now move to the curve generation. Scilab software may become
very useful, particularly if you have experimental data to fit and will
need to introduce nonlinear components. However, it is a handful and
there is a steep learning curve.
So, for now, we will plot simple solutions using the plotting applet.
We will use the simple kernel which we found in the previous question,
a solution which has a a) constant and b) one exponential with time
constant 1/l1= R1 *C1
Asymptotic regime

The constant is a0= 1/(R1 +R2) and it defines 'asymptotic behavior' aka
steady state: when the 'transients die away' and new (dynamic) equilibrium
is reached your circuit is
P(oo)= a0 * R(oo)
where oo is symbol for 'infinity' = long time after last change of R = rate of pay,
=long time after time 1/l1.
It says that 'rate of production' P is proportional to 'rate of pay'  which is
somewhat unrealistic, since workers productivity is limited, but not surprising
in a linear model. (Ur can be fixed by saying that R=log(pay) or so, but that is a
different issue).
Transient regime

In the time interval between change in pay R, and equilibrium, the
production P will oscillate around new level P(oo) like this
P(delta) = (a1/l1) * exp ( delta * l1) * sin( delta * omega1)
parameters of the curve a1 l1 omega1 comes from the kernel.
Actually, 'solving the circuit' means to express these parameters of the kernel
in terms of the values of components ( L1, R1, R2, C1 ...).
Let's plot the transient P(delta) or P(t) , which is P(t) shifted in time
(t is time, delta is time counted from the
moment t0 in which the change of pay was made) (we will set it to zero)
Let call the applet
http://www.pa.uky.edu/~phy211/graph_applets/plot_graph.html
and enter
((5*exp( .03 *x) * sin (.1 *x ) )) +5 into function box
min x to be zero
max x to be 100
and press [plot the function] button
we get a transient, starting at P=5 , swinging up to P=8, down to 4, ...
to see settling we increase x.max
max x to be 100 and press [plot the function] button
To get further, we meed to find a 'step function' s(t,t0) = 0 if t<t0
and s(t,t0)=1 if t>t0
unfortunately, clicking on the Special Functions link produces error.
Meaning, we have to find a better plotting program.
Please, let me know if you have a preference
(but it has to be multiplatform, run not only on windows)
Hedgie

Request for Answer Clarification by
cerebga
on
20 Mar 2006 13:51 PST
The function to plot worked as you described. I made a typo in
entering the function by typing .3 instead of .03 and was rather
pleased to see the very shape of a curve I hoped to generate. It
basically looks like a damped version of the function you offered.
I agree, the Scilab software calls for an amount of learning I am not
prepared to get involved with, although I think it is an incredibly
powerful general engineering tool, and I am glad you steered me toward
it.
You referred back to previous question where you suggested a Laplace
transform as well as to your "description of the task". Under the
title "transient regime" you offered the production equation stating
that the parameters of that equation come from the circuit components
and their values. And that is what I would like to see  the
correlation between curves and circuit elements and their values. So
far so good. The correlation itself I do not yet understand as yet.
You clearly understand what I am searching for, and I will therefore
follow your lead. I trust that in the end I will be able to follow
the steps you have already outlined, and others that will be needed,
to generate the curves. If you have a suitable plotting program, just
advise what it is, how it was generated, whether I have such program
or not, and I will accept that as the answer. Thank you for your
patience thus far.
Re other plotting programs  I have no preference. You suggested the
free version of matlab.

Clarification of Answer by
hedgiega
on
21 Mar 2006 19:20 PST
Octave  a free version of matlab is good, but like scilab, has more
complexity then we need. After looking around a bitm I picked
gnuplot. If the graphing software has the same capabilties, you can use that
of course. I will use gnuplot to visualise the response function,
and post in the next RFC. In the meantime, you may install it and work
through a tutorial.
Description
http://www128.ibm.com/developerworks/library/lgnuplot/
tutorials
http://www.gnuplot.info/help.html
examples
http://t16web.lanl.gov/Kawano/gnuplot/misc1e.html
homepage  download Current officially released version is gnuplot 4.0
http://www.gnuplot.info/

Clarification of Answer by
hedgiega
on
22 Mar 2006 05:20 PST
So,
I did entered the equations into gnuplot.
It seems to do what you want  a response to rectangular pulses.
I mean in a qualitative way. I used arbitrary numerical values.
I can post the pictures of plot, if you need that, but it
will be much better if you creeate them on your version of gnuplot,
so you can change the values and tweek the axes and labels.
What follows is a gnuplot script.
Meaning:
1) you put the text between lines \\ //
into a file called RCL.dem
and type
gnuplot RCL.dem
into a shell. I believe MS people call that 'command console' .
Then, at each prompt 'Hit return to continue'
you 'hit return' (it is, as you see, easy)
and script waks you thrue creation of stimulus pulses and responses
2) then you open the file RCL.dem in the editor and
either, modify values
a0, a1, lambda1 omega1, which describe the linear system
and
start end height , start2 height2, end2
which represent stimuli
(batch mode)
or
lift the command from the *.dem file and paste them into fnuplot console.
That should complete the portion
'how to plot response of the circuit'
of the question.
We will have to start thinking of the closure,
but please, do ask if these instrucyions are not clear.
The script
\\\
#
# $Id: RCL response of a simple linear system to a step function
#
# Requires data files "pulse.dat", kernel.dat from this directory,
# so change current working directory to this directory before running.
# gnuplot> set term <termtype>
# gnuplot> load 'RCL.dem'
#
set key left box
set title " Test availability of functions"
set samples 150
set zeroaxis
set xlabel "Date\nTime"
set yrange [5:5]
set xrange [1:50]
plot [1:10] cos(x), 3*sgn(x), exp(.5*x)
pause 1 "Hit return to continue"
set key right nobox
set title " Plot rectangular pulses : tall (bonus)"
set samples 100
plot 1+sgn(x)
# select pulse 1 == rise
start=3
end=33
height=1.5
#define pulse
pulse(x) =( sgn(xstart)  sgn(xend))* height * .5
#plot it
set yrange [.5:20]
set xrange [0:50]
plot pulse(x)
pause 1 "Hit return to continue"
set title " Plot both stimuli a) rise (pulse2) and b) bonus (pulse) "
set key left box
set samples 200
start2=3
end2=33
height2=5
pulse2(x)=( sgn(xstart2)  sgn(xend2))* height2 *.5
start=3
end=5
height=10
set yrange [.5:200]
set xrange [0:50]
set zeroaxis
plot pulse(x), pulse2(x)
pause 1 "Hit return to continue"
set title " Rise stimulus and Response to the 'rise' "
set zeroaxis
set yrange [2:100]
R=0.5
a1=2
a0= 1/R
lambda1=.2
omega1=.9
kernel(x) = a0 + a1*exp(lambda1*x)*cos(omega1*x)
response(x) = (x<0) ? 0 : kernel(x)
prod(x)= height*(1/R+response(xstart)  response(xend))
plot prod(x)
pause 1 "Hit return to continue"
set title " Response to the 'bonus timulus' "
set yrange [.5:50]
plot height2*(1/R+response(xstart2)  response(xend2))
pause 1 "Hit return to exit"
reset
///
Hedgie
