内弹道仿真matlab全代码

#内弹道matlab代码

来自丰兄的免费分享

作为需要内弹道仿真的苦逼，在网上找到的代码都需要付费下载，于是我将老师的代码分享给大家供大家免费使用，抵制付费！

副脚本

function dy=ndd_fun(t,y,C)
chi=C(1);lambda=C(2);lambda_s=C(3);chi_s=C(4);Z_s=C(5);mu=C(12);
theta=C(6);B=C(7);V=C(8);Delta=C(9);delta=C(10);alpha=C(11);
Z=y(1);l=y(2);v=y(3);
psi=(Z>=0&Z<1).*(chi*Z.*(1+lambda*Z+mu*Z^2))+...
    (Z>=1&Z<Z_s).*(chi_s.*(1+lambda_s*Z))+...
    (Z>=Z_s)*1;
l_psi=1-(Delta/delta)*(1-psi)-alpha*Delta*psi;
p=(psi-v*v)/(1+l_psi);
dy(1)=sqrt((theta/(2*B))*(p^V)*(Z>=0&Z<=Z_s));
dy(2)=v;
dy(3)=theta*p/2;
dy=[dy(1);dy(2);dy(3)];
end

主脚本

%function ndd
%59-130
A=1.394;
G=33.4;
WO=18.56;
l_g=59.52;
P_0=30000;
fail=1.0