How to use event solver to compute numerical solution to impulsively forced van der pol oscillator?
1 次查看(过去 30 天)
显示 更早的评论
I want to use the event solver to find numerical solution to the van der pol equation with impulsive forcing on the right hand side.
I was told that I could use event solver to incorporate impulsive force on the system by first integrating using ode45 and by stopping the integration at time t when i want impulse force to exert and looping it over for multiple n.
How should I write a code on matlab if I want it to compute this for me? Should I start with unforced van der pol system and add the integration step then use the event solver? Or is there any other way to incorporate dirac delta function so that there won't be any problems in the integration process?
0 个评论
回答(1 个)
Vaibhav Awale
2016-1-13
Hi,
One method I can think of is to approximate the impulse using a "heaviside" function and using that as an input to ode45 function. I am attaching a sample 'odeuse.m" script and "vdp10.m" function. I have approximated the impulse function using the below command:
>> u = 10*(heaviside(ut-10) - heaviside(ut-10.1));
You can adjust the peak of impulse and time period appropriately such that the integral over [-Inf, Inf] is 1. In this case the approximated impulse occurs at t = 10s.
Regards,
Vaibhav
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!