[交流] 帮忙用matlab联立求解这九个方程已有1人参与

Vp=(Pp-12)/(-3.75*rho_f*Apipe);                             %(6)

Pu(t)=Pp*(t-Tu)+Z0*Vp*(t-Tu)-Z0*Vu(t);                   %(7)

Pp(t)=Pu*(t-Tu)-Z0*Vu*(t-Tu)+Z0*Vp(t);                   %(8)

Pd(t)=2*(PHTS+Koutlet*rho_f*Vd(t)^2/2)-Pd*(t-2*Td)-Z0*Vd*(t-2*Td)+Z0*Vd(t);%(9)

Vv=2.402*10^(-5)*Cv*((Pu-Pd)/rho_f)^0.5/Apipe;          %(10)

Vu=Vv+Aseat/Apipe*y(2);                           %(11)

Vd=Vv-(Aplug-Aseat)/Apipe*y(2)-(Acage-Aplug)*Vb;        %(12)

Vb=-(Acage-Abellow)/(Acage-Aplug)*y(2);                 %(13)

Pb=Pd-8*3.14*miu*L*Vb/((Aplug+Acage)-2*(Acage-Aplug)/log(Acage/Aplug));%(14)

        rho_f=1105;

        Apipe=1.91*10^(-3);

        Z0=1.56*10^6

        Tu=6.67*10*(-4);

        Td=1.65*10^(-2);

        PHTS=0.01*10^6; Abellow=1.22*10^(-3);        Acage=1.44*10^(-3);        Ad=1.42*10^(-1);        Apipe=1.91*10^(-3);        Aplug=1.43*10^(-3);        Aseat=2.85*10^(-4);

        Cv=0.12*y(1);              k=9.25*10^5;        Koutlet=17;        L=0.0254;

        miu=4.6*10^(-4);