Question: Integral devision by zero

Hello,

For a control systems project I'am working on I need to minimize the actuator effort required to control an input voltage to certain output voltage. As a first experiment I chose a sinusoidal input ug and a sinusoidal output ul. The function to minimize is the funtion Uint. We know that when both the input voltage and output voltage share the frequency, phase and amplitude, that the control effort is zero. Running the animation (see Maple input below) shows that if the input and output voltage frequencies approach eachother the function Uint goes to zero. However, if the frequencies exactly match and I evaluate the function Uint I get a devision by zero notification (this could already be seen from the function of Uint). How can I solve this?

restart

with(plots):

u[G] := sin(omega[G]*t)

u[L] := sin(omega[L]*t)

U[int] := int((u[G]-u[L])^2, t)

omega[G] := 9:

animate(plot, [[u[G], u[L], (1/10)*U[int]], t = 0 .. 5, legend = ["Public grid voltage", "Local grid voltage (control goal)", "Cum. Actuator voltage (scaled)"]], omega[L] = 8 .. 10, gridlines = true, labels = ["Time [s]","Voltage [V]"], labeldirections = ["horizontal", "vertical"], labelfont = ["ARIAL", "bold", 12])

omega[L] := 9; t := 10; evalf(U[int])

NULL

 

Download 20131109_Division_by.mw

 

Thanks

Please Wait...