用MATLAB编了一个程序，基于下面的公式，小伙伴感觉对这个公式熟悉又陌生；编辑处理的效果通过图，大家可以看出好多，这里我就不去解释了，大家可以看懂；程序代码：%已知参数R=200;%mmpb=180;%MPasigma_fb=2000;%MPaalpha_deg=1:1:60;%角度范围alpha_rad=deg2rad(alpha_deg);%转换为弧度%计算公式h_alpha=(R*pb)./(2*sigma_fb.*cos(alpha_rad).^2);h_theta=(R*pb)./(2*sigma_fb).*(2-tan(alpha_rad).^2);h_sum=h_alpha+h_theta;%%2×2子图布局figure('Position',[1001001000800]);%1.h_alpha与角度subplot(2,2,1);plot(alpha_deg,h_alpha,'b-','LineWidth',2);gridon;xlabel('\alpha(deg)','FontSize',12);ylabel('h_\alpha(mm)','FontSize',12);title('h_\alpha与角度','FontSize',14);set(gca,'FontWeight','bold');%2.h_theta与角度subplot(2,2,2);plot(alpha_deg,h_theta,'r-','LineWidth',2);gridon;xlabel('\alpha(deg)','FontSize',12);ylabel('h_\theta(mm)','FontSize',12);title('h_\theta与角度','FontSize',14);set(gca,'FontWeight','bold');%3.h_alpha和h_theta对比subplot(2,2,3);plot(alpha_deg,h_alpha,'b-','LineWidth',2);holdon;plot(alpha_deg,h_theta,'r--','LineWidth',2);gridon;xlabel('\alpha(deg)','FontSize',12);ylabel('厚度h(mm)','FontSize',12);title('h_\alpha与h_\theta对比','FontSize',14);legend('h_\alpha','h_\theta','Location','best');set(gca,'FontWeight','bold');%4.h_alpha+h_theta与角度subplot(2,2,4);plot(alpha_deg,h_sum,'m-','LineWidth',2);gridon;xlabel('\alpha(deg)','FontSize',12);ylabel('h_\alpha+h_\theta(mm)','FontSize',12);title('总厚度与角度','FontSize',14);set(gca,'FontWeight','bold');sgtitle('厚度与角度关系分析','FontSize',16,'FontWeight','bold');来源：气瓶设计的小工程师
