How can I solve system of the matrix differential equation within MATLAB?

I have a matrix of second order differential equation as;

X''= A X' + B X + C

where A, B and C are nxn dimensional matrices, and I am seeking for solving this equation via ode functions. I appreciate your feedback

Thanks

Javad

Answer 1

[George Papazafeiropoulos]

eqn = 'D2y - A*Dy - B*y = C';

inits2 = 'y(0)=a, Dy(0)=b';

y=dsolve(eqn,inits2,'x')

where a and b are the initial or boundary conditions

[Javad]

Please note that the A,B and C are n dimensional matrices, as I know dsolve is used for scalar values.

[George Papazafeiropoulos]

You can use the above solution procedure for each relation of your matrix equation separately. This is implemented in the following:

http://www.mathworks.com/matlabcentral/newsreader/view_thread/297700

Except for this, you can adapt the Fint_K_C function of the following submission:

http://www.mathworks.com/matlabcentral/fileexchange/44649-general-single-step-single-solve-integration-algorithm

to solve directly the matrix equation.