【小白本白】使用MATLAB计算公式，并把结果输出到Excel

科技前沿 • 2025-02-18 13:31 • 阅读 57

大家好，我是讯享网，很高兴认识大家。

【前情提要】小白从未使用过MATLAB软件，刚刚上手，尝试简化一些日常的计算工作。因为我知道自己写的代码非常冗长也很呆，并且自己对软件使用十分不熟练。本帖仅是为了记录一下学习的过程，完全不能作为教学贴使用。希望各位大佬看到后轻喷，也希望能够有大佬能留下更常用更便捷的方法。

【正文部分】想要获得气温、空气密度与海拔的关系，并计算温度、空气密度随海拔上升的数值。

根据空气动力学书籍，我们能够得到温度与海拔的关系式：

讯享网

大气密度与海拔的关系式：

0-11000m

11000-20000m

20000-32000m

接下来编写MATLAB程序：

为了之后输出到Excel表格方便查看，先为Excel表格编写表头：

%制作表头 rowName=cell(1,3); rowName{1,1}='海拔高度h(m)'; rowName{1,2}='温度T(℃)'; rowName{1,3}='空气密度ρ(kg/m^3)';

讯享网

在桌面创建表格，并把表头写入

讯享网xlswrite('C:\Users\ChenMX\Desktop\T_Rho.xlsx',rowName(1,1),'sheet1','A1'); xlswrite('C:\Users\ChenMX\Desktop\T_Rho.xlsx',rowName(1,2),'sheet1','B1'); xlswrite('C:\Users\ChenMX\Desktop\T_Rho.xlsx',rowName(1,3),'sheet1','C1'); %将表头分配到Excel表格的第一行各个位置。

给定初始条件

%已知的初始条件，海平面上的标准值，密度rho，温度T0，压强P0 rho0=1.225;T0=288.15;P0=;

根据上述的公式，将代码分为三段，第一段为0-11000m段。

设置海拔高度从0-11000m，每隔500m计算一次，h为一个一维行向量，元素个数为（11000/500+1=23个），之后T与Rho的计算都需要h的参与，所以都是一维行向量。那么A就是一个3行23列的矩阵。希望Excel表能够竖着输出，所以将A转置后成为B。因为表格的第一行已经写了表头，所以写入表格的时候选择sheet1，并且从A2开始写。

讯享网%海拔为0-11000m阶段 Pa=P0; h=0:500:11000; T=T0-0.0065*h; Rho=rho0*(T/T0).^4.25588; A=[h;T;Rho]; B=A'; xlswrite('C:\Users\ChenMX\Desktop\T_Rho.xlsx',B,'sheet1','A2')

同理写出在11000-20000m和20000-m时候的代码，其实只要复制前面的代码稍微改改就可以。注意在11000-20000m平流层段的时候温度不变，那么需要单独设一个元素都相同的一维向量，从11500-20000m中间间隔500共有（20000-11500）/500+1=18个，使用repmat语句创建。

%海拔为11000m-20000m阶段 h=11500:500:20000; T=repmat(216.65,1,18); Rho=0.36392*exp((-h+11000)/6341.62); A=[h;T;Rho]; B=A'; xlswrite('C:\Users\ChenMX\Desktop\T_Rho.xlsx',B,'sheet1','A25') %海拔为20000m-32000m阶段 h=20500:500:32000; T=216.65+0.001*(h-20000); Rho=0.088035*(T/216.65).^-35.1632; A=[h;T;Rho]; B=A'; xlswrite('C:\Users\ChenMX\Desktop\T_Rho.xlsx',B,'sheet1','A43')

最后输出的结果如图所示：

【小白本白】使用MATLAB计算公式，并把结果输出到Excel

相关推荐