px4-firmware/matlab/scripts/Inertial Nav EKF/K_VELY.c

117 lines
4.1 KiB
C
Raw Normal View History

t2 = q0*q3*2.0;
t9 = q1*q2*2.0;
t3 = t2-t9;
t4 = q0*q0;
t5 = q1*q1;
t6 = q2*q2;
t7 = q3*q3;
t8 = t4-t5+t6-t7;
t10 = q0*q1*2.0;
t11 = q2*q3*2.0;
t12 = t10+t11;
t13 = q1*vd*2.0;
t14 = q0*ve*2.0;
t26 = q3*vn*2.0;
t15 = t13+t14-t26;
t16 = q0*vd*2.0;
t17 = q2*vn*2.0;
t27 = q1*ve*2.0;
t18 = t16+t17-t27;
t19 = q3*vd*2.0;
t20 = q2*ve*2.0;
t21 = q1*vn*2.0;
t22 = t19+t20+t21;
t23 = q3*ve*2.0;
t24 = q0*vn*2.0;
t28 = q2*vd*2.0;
t25 = t23+t24-t28;
t29 = P[0][0]*t15;
t30 = P[1][1]*t18;
t31 = P[5][4]*t8;
t32 = P[6][4]*t12;
t33 = P[0][4]*t15;
t34 = P[1][4]*t18;
t35 = P[2][4]*t22;
t78 = P[4][4]*t3;
t79 = P[3][4]*t25;
t36 = t31+t32+t33+t34+t35-t78-t79;
t37 = P[5][6]*t8;
t38 = P[6][6]*t12;
t39 = P[0][6]*t15;
t40 = P[1][6]*t18;
t41 = P[2][6]*t22;
t81 = P[4][6]*t3;
t82 = P[3][6]*t25;
t42 = t37+t38+t39+t40+t41-t81-t82;
t43 = t12*t42;
t44 = P[5][0]*t8;
t45 = P[6][0]*t12;
t46 = P[1][0]*t18;
t47 = P[2][0]*t22;
t83 = P[4][0]*t3;
t84 = P[3][0]*t25;
t48 = t29+t44+t45+t46+t47-t83-t84;
t49 = t15*t48;
t50 = P[5][1]*t8;
t51 = P[6][1]*t12;
t52 = P[0][1]*t15;
t53 = P[2][1]*t22;
t85 = P[4][1]*t3;
t86 = P[3][1]*t25;
t54 = t30+t50+t51+t52+t53-t85-t86;
t55 = t18*t54;
t56 = P[5][2]*t8;
t57 = P[6][2]*t12;
t58 = P[0][2]*t15;
t59 = P[1][2]*t18;
t60 = P[2][2]*t22;
t87 = P[4][2]*t3;
t88 = P[3][2]*t25;
t61 = t56+t57+t58+t59+t60-t87-t88;
t62 = t22*t61;
t63 = P[5][3]*t8;
t64 = P[6][3]*t12;
t65 = P[0][3]*t15;
t66 = P[1][3]*t18;
t67 = P[2][3]*t22;
t89 = P[4][3]*t3;
t90 = P[3][3]*t25;
t68 = t63+t64+t65+t66+t67-t89-t90;
t69 = P[5][5]*t8;
t70 = P[6][5]*t12;
t71 = P[0][5]*t15;
t72 = P[1][5]*t18;
t73 = P[2][5]*t22;
t92 = P[4][5]*t3;
t93 = P[3][5]*t25;
t74 = t69+t70+t71+t72+t73-t92-t93;
t75 = t8*t74;
t80 = t3*t36;
t91 = t25*t68;
t76 = R_VEL+t43+t49+t55+t62+t75-t80-t91;
t77 = 1.0/t76;
A0[0][0] = t77*(t29-P[0][4]*t3+P[0][5]*t8+P[0][6]*t12+P[0][1]*t18+P[0][2]*t22-P[0][3]*t25);
A0[1][0] = t77*(t30-P[1][4]*t3+P[1][5]*t8+P[1][0]*t15+P[1][6]*t12+P[1][2]*t22-P[1][3]*t25);
A0[2][0] = t77*(t60-P[2][4]*t3+P[2][5]*t8+P[2][0]*t15+P[2][6]*t12+P[2][1]*t18-P[2][3]*t25);
A0[3][0] = t77*(-t90-P[3][4]*t3+P[3][5]*t8+P[3][0]*t15+P[3][6]*t12+P[3][1]*t18+P[3][2]*t22);
A0[4][0] = t77*(-t78+P[4][5]*t8+P[4][0]*t15+P[4][6]*t12+P[4][1]*t18+P[4][2]*t22-P[4][3]*t25);
A0[5][0] = t77*(t69-P[5][4]*t3+P[5][0]*t15+P[5][6]*t12+P[5][1]*t18+P[5][2]*t22-P[5][3]*t25);
A0[6][0] = t77*(t38-P[6][4]*t3+P[6][5]*t8+P[6][0]*t15+P[6][1]*t18+P[6][2]*t22-P[6][3]*t25);
A0[7][0] = t77*(-P[7][4]*t3+P[7][5]*t8+P[7][0]*t15+P[7][6]*t12+P[7][1]*t18+P[7][2]*t22-P[7][3]*t25);
A0[8][0] = t77*(-P[8][4]*t3+P[8][5]*t8+P[8][0]*t15+P[8][6]*t12+P[8][1]*t18+P[8][2]*t22-P[8][3]*t25);
A0[9][0] = t77*(-P[9][4]*t3+P[9][5]*t8+P[9][0]*t15+P[9][6]*t12+P[9][1]*t18+P[9][2]*t22-P[9][3]*t25);
A0[10][0] = t77*(-P[10][4]*t3+P[10][5]*t8+P[10][0]*t15+P[10][6]*t12+P[10][1]*t18+P[10][2]*t22-P[10][3]*t25);
A0[11][0] = t77*(-P[11][4]*t3+P[11][5]*t8+P[11][0]*t15+P[11][6]*t12+P[11][1]*t18+P[11][2]*t22-P[11][3]*t25);
A0[12][0] = t77*(-P[12][4]*t3+P[12][5]*t8+P[12][0]*t15+P[12][6]*t12+P[12][1]*t18+P[12][2]*t22-P[12][3]*t25);
A0[13][0] = t77*(-P[13][4]*t3+P[13][5]*t8+P[13][0]*t15+P[13][6]*t12+P[13][1]*t18+P[13][2]*t22-P[13][3]*t25);
A0[14][0] = t77*(-P[14][4]*t3+P[14][5]*t8+P[14][0]*t15+P[14][6]*t12+P[14][1]*t18+P[14][2]*t22-P[14][3]*t25);
A0[15][0] = t77*(-P[15][4]*t3+P[15][5]*t8+P[15][0]*t15+P[15][6]*t12+P[15][1]*t18+P[15][2]*t22-P[15][3]*t25);
A0[16][0] = t77*(-P[16][4]*t3+P[16][5]*t8+P[16][0]*t15+P[16][6]*t12+P[16][1]*t18+P[16][2]*t22-P[16][3]*t25);
A0[17][0] = t77*(-P[17][4]*t3+P[17][5]*t8+P[17][0]*t15+P[17][6]*t12+P[17][1]*t18+P[17][2]*t22-P[17][3]*t25);
A0[18][0] = t77*(-P[18][4]*t3+P[18][5]*t8+P[18][0]*t15+P[18][6]*t12+P[18][1]*t18+P[18][2]*t22-P[18][3]*t25);
A0[19][0] = t77*(-P[19][4]*t3+P[19][5]*t8+P[19][0]*t15+P[19][6]*t12+P[19][1]*t18+P[19][2]*t22-P[19][3]*t25);
A0[20][0] = t77*(-P[20][4]*t3+P[20][5]*t8+P[20][0]*t15+P[20][6]*t12+P[20][1]*t18+P[20][2]*t22-P[20][3]*t25);
A0[21][0] = t77*(-P[21][4]*t3+P[21][5]*t8+P[21][0]*t15+P[21][6]*t12+P[21][1]*t18+P[21][2]*t22-P[21][3]*t25);
A0[22][0] = t77*(-P[22][4]*t3+P[22][5]*t8+P[22][0]*t15+P[22][6]*t12+P[22][1]*t18+P[22][2]*t22-P[22][3]*t25);
A0[23][0] = t77*(-P[23][4]*t3+P[23][5]*t8+P[23][0]*t15+P[23][6]*t12+P[23][1]*t18+P[23][2]*t22-P[23][3]*t25);