文档介绍：弹道仿真C源程序/*铅垂平面内飞行器弹道仿真C源程序*/#include""#include""#include""#[5][6],ACY[5][6],AMZAF[5][6],AMZWZ[5][2],AXG[11][2],AJZ[3][2],AP[2],AMC[2],AGC[2],ANDM[2],ANDAF[2],Y[8],B[3],L,S,SONIC,RHO;/*ACX[5][6]:阻力系数数组;ACY[5][6]:升力系数数组;*//*AMZAF[5][6]:某一质心位置下的静稳定性导数数组;AMZWF[5][2]:阻尼力矩系数导数数组;*//*AXG[10][2]:质心位置Xg随时间的变化规律;AJZ[3][2]:转动惯量Jz随时间的变化规律;*//*AP[2]:起飞,续航发动机的推力值;AMC[2]:起飞,续航发动机燃料质量秒流量值;*//*AGC[2]:质心位置变化的始末值;ANDM[2]:气动数据插值所需要的Ma数的最大最小值;*//*ANDAF[2]:气动数据插值所需要的攻角的最大最小值;Y[8]存放积分结果的数组;*//*B[3]:存放三个时间值的数组;L:特征长度;S:特征面积;SONIC:声速;RHO:大气密度;*/doubleMA,ABC_ALPHA,CX,CY,MZAF,MZWZ,XG,JZ,ALPHA;/*MA:马赫数;ABC_ALPHA:攻角的绝对值;CX:阻力系数;CY:升力系数;MZAF:静稳定性导数值;*//*MZWZ:阻尼力矩导数值;XG:质心位置;JZ:转动惯量;ALPHA:攻角;*/FILE*fi,*fo;/*主函数*/voidmain(){voidrk();/*调用rk函数*/voidresult();/*调用result函数*/inti,j;doubleh;/*h为步长*/fi=fopen("","r");for(i=0;i<5;i++){for(j=0;j<6;j++)fscanf(fi,"%lf",&ACX[i][j]);}for(i=0;i<5;i++){for(j=0;j<6;j++)fscanf(fi,"%lf",&ACY[i][j]);}for(i=0;i<5;i++){for(j=0;j<6;j++)fscanf(fi,"%lf",&AMZAF[i][j]);}for(i=0;i<5;i++){for(j=0;j<2;j++)fscanf(fi,"%lf",&AMZWZ[i][j]);}for(i=0;i<10;i++){for(j=0;j<2;j++)fscanf(fi,"%lf",&AXG[i][j]);}for(i=0;i<3;i++){for(j=0;j<2;j++)fscanf(fi,"%lf",&AJZ[i][j]);}for(i=0;i<