弹道仿真C源程序

《弹道仿真C源程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、弹道仿真C源程序/*铅垂平面内飞行器弹道仿真C源程序*/#include"math.h"#include"stdio.h"#include"stdlib.h"#defineRAD57.3doubleACX[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

2、][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]：气动数据插值所需要的攻角

3、的最大最小值；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();/*

4、调用rk函数*/voidresult();/*调用result函数*/inti,j;doubleh;/*h为步长*/fi=fopen("input.dat","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]

5、[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

6、=0;i<3;i++){for(j=0;j<2;j++)fscanf(fi,"%lf",&AJZ[i][j]);}for(i=0;i<2;i++)fscanf(fi,"%lf",&AP[i]);for(i=0;i<2;i++)fscanf(fi,"%lf",&AMC[i]);for(i=0;i<2;i++)fscanf(fi,"%lf",&AGC[i]);for(i=0;i<

7、;2;i++)fscanf(fi,"%lf",&ANDM[i]);for(i=0;i<2;i++)fscanf(fi,"%lf",&ANDAF[i]);for(i=0;i<8;i++)fscanf(fi,"%lf",&Y[i]);for(i=0;i<3;i++)fscanf(fi,"%lf",&B[i]);fscanf(fi,"%lf%lf%lf%lf%lf",&

8、L,&S,&SONIC,&RHO,&h);fclose(fi);fo=fopen("output.dat","w");fprintf(fo,"%s","TIMEALPHAVXY");fprintf(fo,"%s