Numerical solution of the controlled Duffing oscillator by hybrid functions

A numerical method for solving the controlled Duffing oscillator is presented. The method is based upon hybrid functions approximation. The properties of hybrid functions which consists of block-pulse functions plus Legendre polynomials are given. The method is based upon expanding various time functions in the system as their truncated hybrid functions. The operational matrix of product is introduced. This matrix together with the operational matrix of integration are utilized to reduce the solution of controlled Duffing oscillator to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique.