
Re: Dirac Delta Question
Posted:
Jan 17, 2013 5:45 PM


"Sam " <samnuzbrokh@yahoo.com> wrote in message <kd9r7h$38b$1@newscl01ah.mathworks.com>... > Hey y'all, I'm very confused as to how to numerically incorporate dirac delta behaviour into my differential equation. Maybe I'm using a wrong approach? I keep getting error messages to no avail. Any advice would be greatly appreciated! > > function dirac > close all > clc; clear > > %Parameters: > G=1.07; > w=2*pi; > w0=(1.5)*w; > B=(w0)/4; > > % Initial Conditions > y0 = [0,0]; > > % Make a theta vs. time plot > [t,y] = ode23(@f,[0 25],y0,[],G,w,w0,B); > > hold all > plot(t,y(:,1)) > > function dydt = f(t,y,G,w,w0,B) > dydt = [y(2);2*B*y(2)(w0^2)*sin(y(1))+G*(w0^2)*cos(w*t)+dirac(t5)];
 Solve ode on the interval t in [0,5),  take y5 = y(t=5)+[0;1] as initial condition to solve ode on the interval t = (5,25]
Bruno

