fm_sae2.m

一个较好的MATLAB潮流程序

            /2+t41*t1)*t65+t80*t61/t88*(2.0*t71+2.0*t72*t1*t42);

        % C(diff(Qa,vb),optimized);
        t1 = va*va;
        t2 = xa0*xa0;
        t4 = t1/t2;
        t5 = vb*vb;
        t6 = xb0*xb0;
        t7 = 1/t6;
        t8 = t5*t7;
        t9 = va*vb;
        t10 = delta-theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t4+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t54 = sin(t10);
        t60 = va*t12+vb*t43*t11;
        t63 = t12*(t37*vb*t43*t54+t40*t60);
        t66 = t4+t8+2.0*t9*t13*t43;
        t67 = 1/t66;
        t68 = vb*t7;
        t69 = va*t11;
        t72 = 2.0*t68+2.0*t69*t12*t14;
        t77 = va*t40*t51;
        t88 = t66*t66;
        DAE.J22(ha,hb) = DAE.J22(ha,hb) -va*t41*t51*t63*t67*t72/2-t77*t12*(t37*t43*t54+t41*t60*t72/2+t40* ...
            t43*t11)*t67+t77*t63/t88*(2.0*t68+2.0*t69*t12*t43);

        % C(diff(Qa,delta),optimized);
        t1 = va*va;
        t2 = xa0*xa0;
        t3 = 1/t2;
        t4 = t1*t3;
        t5 = vb*vb;
        t6 = xb0*xb0;
        t8 = t5/t6;
        t9 = va*vb;
        t10 = delta-theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t4+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t52 = t51*t3;
        t54 = t37*vb;
        t55 = sin(t10);
        t56 = t43*t55;
        t61 = va*t12+vb*t43*t11;
        t63 = t54*t56+t40*t61;
        t66 = t4+t8+2.0*t9*t13*t43;
        t67 = 1/t66;
        t69 = vb*t55;
        t90 = t66*t66;
        DAE.J21(ha,ha) = DAE.J21(ha,ha) + t1*t41*t52*t63*t67*t69*t14-va*t40*t51*t12*(t54*t43*t11-t41*t61*va* ...
            t69*t12*t14-t40*vb*t56)*t67-2.0*t1*t40*t52*t63/t90*t69*t43;

        % C(diff(Qa,theta),optimized);
        t1 = va*va;
        t2 = xa0*xa0;
        t3 = 1/t2;
        t4 = t1*t3;
        t5 = vb*vb;
        t6 = xb0*xb0;
        t8 = t5/t6;
        t9 = va*vb;
        t10 = delta-theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t4+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t52 = t51*t3;
        t54 = t37*vb;
        t55 = sin(t10);
        t56 = t43*t55;
        t61 = va*t12+vb*t43*t11;
        t63 = t54*t56+t40*t61;
        t66 = t4+t8+2.0*t9*t13*t43;
        t67 = 1/t66;
        t69 = vb*t55;
        t90 = t66*t66;
        DAE.J21(ha,hb) = DAE.J21(ha,hb) -t1*t41*t52*t63*t67*t69*t14-va*t40*t51*t12*(-t54*t43*t11+t41*t61*va ...
            *t69*t12*t14+t40*vb*t56)*t67+2.0*t1*t40*t52*t63/t90*t69*t43;

        % C(diff(Pb,va),optimized);
        t1 = va*va;
        t2 = xa0*xa0;
        t3 = 1/t2;
        t4 = t1*t3;
        t5 = vb*vb;
        t6 = xb0*xb0;
        t8 = t5/t6;
        t9 = va*vb;
        t10 = -delta+theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t4+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t59 = sin(t10);
        t60 = t12*t59;
        t63 = t43*(t37*(-vb*t43-va*t12*t11)+t40*va*t60);
        t66 = t4+t8+2.0*t9*t13*t43;
        t67 = 1/t66;
        t68 = va*t3;
        t69 = vb*t11;
        t72 = 2.0*t68+2.0*t69*t12*t14;
        t77 = vb*t40*t51;
        t89 = t66*t66;
        DAE.J12(hb,ha) = DAE.J12(hb,ha) + vb*t41*t51*t63*t67*t72/2+t77*t43*(-t37*t12*t11+t41*va*t60*t72/2+t40 ...
            *t12*t59)*t67-t77*t63/t89*(2.0*t68+2.0*t69*t12*t43);

        % C(diff(Pb,vb),optimized);
        t1 = va*va;
        t2 = xa0*xa0;
        t4 = t1/t2;
        t5 = vb*vb;
        t6 = xb0*xb0;
        t7 = 1/t6;
        t8 = t5*t7;
        t9 = va*vb;
        t10 = -delta+theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t4+t8+2.0*t9*t13*t14-t38);
        t41 = 1/xb0;
        t42 = xt1*xt1;
        t45 = xt2*xt2;
        t49 = 1/(t19+t29+t12+t41-1/t42*t24-1/t45*t34);
        t58 = t12*sin(t10);
        t61 = t41*(t37*(-vb*t41-va*t12*t11)+t40*va*t58);
        t64 = t4+t8+2.0*t9*t13*t41;
        t65 = 1/t64;
        t68 = 1/t40;
        t71 = vb*t7;
        t72 = va*t11;
        t75 = 2.0*t71+2.0*t72*t12*t14;
        t80 = vb*t40*t49;
        t89 = t64*t64;
        DAE.J12(hb,hb) = DAE.J12(hb,hb) + t40*t49*t61*t65+vb*t68*t49*t61*t65*t75/2+t80*t41*(-t37*t41+t68*va* ...
            t58*t75/2)*t65-t80*t61/t89*(2.0*t71+2.0*t72*t12*t41);

        % C(diff(Pb,delta),optimized);
        t1 = vb*vb;
        t2 = va*va;
        t3 = xa0*xa0;
        t4 = 1/t3;
        t5 = t2*t4;
        t6 = xb0*xb0;
        t7 = 1/t6;
        t8 = t1*t7;
        t9 = va*vb;
        t10 = -delta+theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t5+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t58 = t40*va;
        t59 = sin(t10);
        t60 = t12*t59;
        t62 = t37*(-vb*t43-va*t12*t11)+t58*t60;
        t67 = t5+t8+2.0*t9*t13*t43;
        t68 = 1/t67;
        t79 = t59*t59;
        t91 = t67*t67;
        DAE.J11(hb,ha) = DAE.J11(hb,ha) + t1*t41*t51*t43*t62*t68*va*t60*t14+vb*t40*t51*t43*(-t37*va*t60+t41* ...
            t2*t4*t79*vb*t14-t58*t13)*t68-2.0*t1*t40*t51*t7*t62/t91*va*t59*t12;

        % C(diff(Pb,theta),optimized);
        t1 = vb*vb;
        t2 = va*va;
        t3 = xa0*xa0;
        t4 = 1/t3;
        t5 = t2*t4;
        t6 = xb0*xb0;
        t7 = 1/t6;
        t8 = t1*t7;
        t9 = va*vb;
        t10 = -delta+theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t5+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t58 = t40*va;
        t59 = sin(t10);
        t60 = t12*t59;
        t62 = t37*(-vb*t43-va*t12*t11)+t58*t60;
        t67 = t5+t8+2.0*t9*t13*t43;
        t68 = 1/t67;
        t79 = t59*t59;
        t91 = t67*t67;
        DAE.J11(hb,hb) = DAE.J11(hb,hb) -t1*t41*t51*t43*t62*t68*va*t60*t14+vb*t40*t51*t43*(t37*va*t60-t41* ...
            t2*t4*t79*vb*t14+t58*t13)*t68+2.0*t1*t40*t51*t7*t62/t91*va*t59*t12;

        % C(diff(Qb,va),optimized);
        t1 = va*va;
        t2 = xa0*xa0;
        t3 = 1/t2;
        t4 = t1*t3;
        t5 = vb*vb;
        t6 = xb0*xb0;
        t8 = t5/t6;
        t9 = va*vb;
        t10 = delta-theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t4+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t54 = sin(t10);
        t60 = vb*t43+va*t12*t11;
        t63 = t43*(-t37*va*t12*t54+t40*t60);
        t66 = t4+t8+2.0*t9*t13*t43;
        t67 = 1/t66;
        t68 = va*t3;
        t69 = vb*t11;
        t72 = 2.0*t68+2.0*t69*t12*t14;
        t77 = vb*t40*t51;
        t88 = t66*t66;
        DAE.J22(hb,ha) = DAE.J22(hb,ha) -vb*t41*t51*t63*t67*t72/2-t77*t43*(-t37*t12*t54+t41*t60*t72/2+t40* ...
            t12*t11)*t67+t77*t63/t88*(2.0*t68+2.0*t69*t12*t43);

        % C(diff(Qb,vb),optimized);
        t1 = 1/xb0;
        t2 = vb*t1;
        t3 = va*va;
        t4 = xa0*xa0;
        t6 = t3/t4;
        t7 = vb*vb;
        t8 = xb0*xb0;
        t9 = 1/t8;
        t10 = t7*t9;
        t11 = va*vb;
        t12 = delta-theta;
        t13 = cos(t12);
        t14 = 1/xa0;
        t15 = t13*t14;
        t16 = 1/ab0;
        t21 = 1/xt1;
        t22 = m1*m1;
        t26 = 1/(t21+a1/t22);
        t31 = 1/xt2;
        t32 = m2*m2;
        t36 = 1/(t31+a2/t32);
        t39 = -b1/m1*t21*t26-b2/m2*t31*t36;
        t40 = t39*t39;
        t42 = sqrt(t6+t10+2.0*t11*t15*t16-t40);
        t43 = xt1*xt1;
        t46 = xt2*xt2;
        t50 = 1/(t21+t31+t14+t1-1/t43*t26-1/t46*t36);
        t58 = t2+va*t14*t13;
        t61 = t1*(-t39*va*t14*sin(t12)+t42*t58);
        t64 = t6+t10+2.0*t11*t15*t1;
        t65 = 1/t64;
        t68 = 1/t42;
        t71 = vb*t9;
        t72 = va*t13;
        t75 = 2.0*t71+2.0*t72*t14*t16;
        t80 = vb*t42*t50;
        t88 = t64*t64;
        DAE.J22(hb,hb) = DAE.J22(hb,hb)+ 2.0*t2-t42*t50*t61*t65-vb*t68*t50*t61*t65*t75/2-t80*t1*(t68*t58*t75 ...
            /2+t42*t1)*t65+t80*t61/t88*(2.0*t71+2.0*t72*t14*t1);

        % C(diff(Qb,delta),optimized);
        t1 = vb*vb;
        t2 = va*va;
        t3 = xa0*xa0;
        t5 = t2/t3;
        t6 = xb0*xb0;
        t7 = 1/t6;
        t8 = t1*t7;
        t9 = va*vb;
        t10 = delta-theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;
        t30 = m2*m2;
        t34 = 1/(t29+a2/t30);
        t37 = -b1/m1*t19*t24-b2/m2*t29*t34;
        t38 = t37*t37;
        t40 = sqrt(t5+t8+2.0*t9*t13*t14-t38);
        t41 = 1/t40;
        t43 = 1/xb0;
        t44 = xt1*xt1;
        t47 = xt2*xt2;
        t51 = 1/(t19+t29+t12+t43-1/t44*t24-1/t47*t34);
        t53 = t37*va;
        t54 = sin(t10);
        t55 = t12*t54;
        t60 = vb*t43+va*t12*t11;
        t62 = -t53*t55+t40*t60;
        t67 = t5+t8+2.0*t9*t13*t43;
        t68 = 1/t67;
        t91 = t67*t67;
        DAE.J21(hb,ha) = DAE.J21(hb,ha) + t1*t41*t51*t43*t62*t68*va*t55*t14-vb*t40*t51*t43*(-t53*t13-t41*t60* ...
            va*vb*t54*t12*t14-t40*va*t55)*t68-2.0*t1*t40*t51*t7*t62/t91*va*t54*t12;

        % C(diff(Qb,theta),optimized);
        t1 = vb*vb;
        t2 = va*va;
        t3 = xa0*xa0;
        t5 = t2/t3;
        t6 = xb0*xb0;
        t7 = 1/t6;
        t8 = t1*t7;
        t9 = va*vb;
        t10 = delta-theta;
        t11 = cos(t10);
        t12 = 1/xa0;
        t13 = t11*t12;
        t14 = 1/ab0;
        t19 = 1/xt1;
        t20 = m1*m1;
        t24 = 1/(t19+a1/t20);
        t29 = 1/xt2;