ode45 matlab 出错,Matlab中ode45求解微分方程组出错。

ode45解微分方程组，结果数量级居然为10e304，请各位帮忙看看啊。

问题为求解一个7自由度系统(两个移动5个转动)在初始位移激励下个自由度的的位移和加速度。

目标函数

function dq=Est(t,q)

global  H I1 I4 I5 tha10 tha20 tha30 tha40 tha50 q0 dq0 m1 m2 m3 m4 m5 I2 I3 r1 r2 r3 rc4 k1f c1f k1v c1v k4f c4f k4v c4v kt1 ct1 kt2 ct2 kt3 ct3 kt4 ct4 g;

%global value,accounting for different mass ,height and sitting postures

stha10=sind(tha10);stha20=sind(tha20);stha30=sind(tha30);stha40=sind(tha40);stha50=sind(tha50);

ctha10=cosd(tha10);ctha20=cosd(tha20);ctha30=cosd(tha30);ctha40=cosd(tha40);ctha50=cosd(tha50);

ctha1020=cosd(tha10-tha20);ctha1030=cosd(tha10-tha30);ctha2030=cosd(tha20-tha30);ctha4050=cosd(tha40-tha50);

l1=0.06*H;l2=0.06*H;l4=0.2*H;l5=0.21*H;r4=0.4278*l4;r5=0.4284*l5;

%q=[q1,q2,q3,q4,q5,q6,q7,q1dot,q2dot,q3dot,q4dot,q5dot,q6dot,q7dot];