绘制随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}');

效果如图