内弹道及枪膛合力Matlab程序

clear;

close all;

format long

d=0.0127;

S=0.82*0.0127^2;

V0=2.04e-5;

l_0=V0/S;

lg=0.924;

f=1000000;

alpha=0.001;

w=0.017;

rou=1600;

theta=0.2;

phi=1.45;

chi=0.79825;

lamda=0.1387;

mu=-0.043956;

e1=0.00052/2;

u1=7.5991e-10;

Is=e1/u1;

chi_s=1.2645;

lamda_s=-0.31322;

zk=1.4434;

%Ik=447000;

m=0.048;

p0=30e6;

delta=800;

psi0=(1/delta-1/rou)/(f/p0+alpha-1/rou); sigma0=sqrt(1+4*lamda*psi0/chi);

z0=2*psi0/chi/(sigma0+1);%(sigma0-1)/2/lamda; %====赋予初值====%

v(1)=0;

l(1)=0;

p(1)=p0;

z(1)=z0;

psi(1)=psi0;