绘制随k变化的等熵面积比公式

xmax = 4;
Ma = 0.1:0.05:xmax;
figure;
hold on;
xlim([0,xmax]);
ylim([0,10]);% 预定义k值的向量
k_values = 1.2:0.1:1.4;% 创建一个细胞数组来存储图例标签
legendStrings = cell(1, length(k_values));% 绘制每条曲线并记录图例标签
lines = [];
for idx = 1:length(k_values)k = k_values(idx);y = 1./Ma .* ((1 + (k-1)/2 .* Ma.^2) * 2/(k+1)).^((k+1)/(2*(k-1)));h = plot(Ma,y);lines = [lines, h]; % 存储每条线的对象legendStrings{idx} = sprintf('k=%.1f', k); % 记录图例标签
end% 在循环结束后添加图例
legend(lines, legendStrings, 'Location', 'Best');
xlabel('Ma');
ylabel('A/A_{cr}');%这里识别了Latex代码
title('Relationship between Ma and A/A_{cr}');

效果如图