主要内容

数字系统数据

数字系统数据概述

航空航天工具箱产品可以使美国空军(USAF)数字系统文件到MATLAB®环境使用datcomimport函数。有关更多信息,请参见datcomimport函数引用页面。本节解释如何从美国空军数字系统文件导入数据。

下列主题中的示例使用作为一个可用航空航天工具箱的例子。你可以通过输入运行示例astimportddatcom在MATLAB命令窗口。

美国空军数字系统文件

下面是一个示例输入文件为美国空军数字系统wing-body-horizontal tail-vertical尾配置运行在5α,两个马赫数,计算静态和动态两个海拔和衍生品。您还可以查看这个文件通过输入类型astdatcom.in在MATLAB命令窗口。

美元FLTCON NMACH = 2.0马赫(1)= 0.1,0.2 $ $ FLTCON NALT = 2.0, ALT (1) = 5000.0, 8000.0 $ $ FLTCON NALPHA = 5, ALSCHD (1) = -2.0, 0.0, 2.0, ALSCHD(4) = 4.0, 8.0,循环= 2.0 $ $ OPTINS SREF = 225.8, CBARR = 5.75, BLREF = 41.15 $ $合成器XCG = 7.08, ZCG = 0.0, XW = 6.1, ZW = -1.4, ALIW = 1.1, XH = 20.2, ZH型= 0.4,ALIH = 0.0,十五= 21.3,ZV = 0.0, VERTUP = .TRUE。$ $身体NX = 10.0 X (1) = -4.9, 0.0, 3.0, 6.1, 9.1, 13.3, 20.2, 23.5, 25.9, R (1) = 0.0, 1.0, 1.75, 2.6, 2.6, 2.6, 2.0, 1.0, 0.0 $ $ WGPLNF CHRDTP = 4.0, SSPNE = 18.7, SSPN = 20.6, CHRDR = 7.2, SAVSI = 0.0, CHSTAT = 0.25, TWISTA = -1.1, SSPNDD = 0.0, DHDADI = 3.0, DHDADO = 3.0, = 1.0型naca - w - 6 - 64美元美元a412 HTPLNF CHRDTP = 2.3, SSPNE = 5.7, SSPN = 6.625, CHRDR = 0.25, SAVSI = 11.0, CHSTAT = 1.0, TWISTA = 0.0, = 1.0型naca - h - 4 - 0012美元VTPLNF CHRDTP = 2.7, SSPNE = 5.0, SSPN = 5.2, CHRDR = 5.3, SAVSI = 31.3, CHSTAT = 0.25, TWISTA = 0.0, = 1.0美元naca - v - 4 - 0012 CASEID SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL尾配置湿下情况

美国空军数字系统生成的输出文件为同一wing-body-horizontal tail-vertical尾配置运行在5α,两个马赫数,两个高度可以被进入类型astdatcom.out在MATLAB命令窗口。

数据从系统文件

使用datcomimport函数把数字系统数据到MATLAB环境中。

alldata = datcomimport (astdatcom。的,真的,0);

导入系统的数据

datcomimport函数创建一个单元阵列结构包含数据的数字系统输出文件。

data = = {1} alldata中的数据结构体字段:案例:“SKYHOGG BODY-WING-HORIZONTAL TAIL-VERTICAL尾配置“马赫:[0.1000 - 0.2000]alt:α(5000 8000):[2 0 2 4 8]nmach: 2 nalt: 2 nalpha: 5 rnnub:[]大卖场:0循环:2 sref: 225.8000 cbar: 5.7500 blref: 41.1500暗淡:“英尺”引出:“度”stmach: 0.6000 tsmach: 1.4000保存:0药栓:[]修剪:0潮湿:1建立:1部分:0 highsym: 0 highasy: 0 highcon: 0 tjet: 0 hypeff: 0磅:0压水式反应堆:0接地:0 wsspn: 18.7000 hsspn: 5.7000 ndelta: 0三角洲:[]deltal: [] deltar:[]已:0 grndht:[]配置:cd [1 x1 struct]: [5 x2x2双重的]cl: [5 x2x2双重的]厘米:[5 x2x2双]的cn: [5 x2x2双]的ca: [5 x2x2双]的xcp: [5 x2x2双重的]cla: [5 x2x2双重的]cma: [5 x2x2双)的地方:[5 x2x2双重的]cnb: [5 x2x2双重的]clb: [5 x2x2双重的]qqinf: [5 x2x2双)的每股收益:[5 x2x2双重的]depsdalp: [5 x2x2双重的]clq: [5 x2x2双]cmq:, [5 x2x2双)的复合:[5 x2x2双重的]cmad: [5 x2x2双]的clp: [5 x2x2双重的]cyp: [5 x2x2双重的]cnp: [5 x2x2双重的]中国北车:[5 x2x2双]的clr: [5 x2x2双)

丢失的系统数据

默认情况下,丢失的数据点设置为99999,数据点设置为南,不存在系统方法或方法不适用的地方。

可以看出在数字系统输出文件并检查导入的数据 C Y β , C n β , C l , C 有数据只在第一个alpha值。这是导入的数据值。

数据。中青文ans (:: 1) = 1.0 e + 004 * -0.0000 -0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 ans (:,: 2) = 1.0 e + 004 * -0.0000 -0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 data.cnb ans (:,: 1) = 1.0 e + 004 * 0.0000 0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 ans (:,: 2) = 1.0 e + 004 * 0.0000 0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999数据。clq ans (:: 1) = 1.0 e + 004 * 0.0000 0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 ans (:,: 2) = 1.0 e + 004 * 0.0000 0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999数据。cmq ans (:,: 1), e = 1.0 + 004 * -0.0000 -0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 ans (:,: 2) = 1.0 e + 004 * -0.0000 -0.0000 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999 9.9999

丢失的数据点会充满第一α的值,因为这些数据点是用于所有α值。

aerotab ={'的地方' ' cnb ' ' clq ' ' cmq的,};k = 1:长度(aerotab) m = 1:数据。nmach h = 1:数据。nalt data.(aerotab{k})(:,m,h) = data.(aerotab{k})(1,m,h); end end end

下面是更新的进口数据值。

数据。中青文ans (:,: 1) = -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 ans (:,: 2) = -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 -0.0035 data.cnb ans (:,: 1) = 1.0 e - 003 * 0.9142 0.8781 0.9142 0.8781 0.9142 0.8781 0.9142 0.8781 0.9142 0.8781 ans (:,: 2) = 1.0 e - 003 * 0.9190 0.8829 0.9190 0.8829 0.9190 0.8829 0.9190 0.8829 0.9190 0.8829数据。clq ans (:,: 1) = 0.0974 0.0984 0.0974 0.0984 0.0974 0.0984 0.0974 0.0984 0.0974 0.0984 ans(:,: 2) = 0.0974 0.0984 0.0974 0.0984 0.0974 0.0984 0.0974 0.0984 0.0974 0.0984数据。cmq ans, (:,: 1) = -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 ans (:,: 2) = -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899 -0.0892 -0.0899

空气动力系数

你现在可以画出空气动力系数:

策划升程曲线的时刻

h1 =图;figtitle ={升程曲线”};k = 1:2次要情节(2,1,k)情节(data.alpha,排列(data.cl (k,::),[1 3 2]))网格ylabel(['升力系数(马赫= ' num2str (data.mach (k))')'])标题(figtitle {k});包含(迎角(度))

策划拖极地的时刻

h2 =图;figtitle ={拖极地”};k = 1:2次要情节(2,1,k)情节(排列(data.cd (k,::),[1 3 2]),排列(data.cl (k,::),[1 3 2]))网格ylabel(['升力系数(马赫= ' num2str (data.mach (k))')'])标题(figtitle {k})端包含(阻力系数)

策划俯仰力矩

h3 =图;figtitle ={俯仰力矩”};k = 1:2次要情节(2,1,k)情节(排列(data.cm (k,::),[1 3 2]),排列(data.cl (k,::),[1 3 2]))网格ylabel(['升力系数(马赫= ' num2str (data.mach (k))')'])标题(figtitle {k})端包含(俯仰力矩系数)

另请参阅