数学建模2021年B题代码

function zhuangliaobi3 x1=[250 275 300 350 400]; y1=[1.4 3.4 6.7 19.3 43.6 ]; y2=[6.32 8.25 12.28 25.97 41.08]; y3=[2.8 4.4 6.2 16.2 45.1 ]; y4=[3.26 4.97 9.32 22.88 38.7]; x2=[250 275 300 325 350 400]; y5=[0.4 0.6 1.1 3.3 6.0 21.1 ]; y6=[2.85 5.35 7.61 7.74 13.81 21.21]; y7=[0.5 1.1 3.0 6.1 9.6 33.5 ]; y8=[6.62 6.62 5.05 8.33 13.1 21.45]; y9=[2.8 7.5 12.6 15.9 27.0 63.2 ]; y10=[4.5 4.79 8.77 16.06 22.41 30.48]; p1=polyfit(x1,y1,4) p2=polyfit(x1,y2,4) p3=polyfit(x1,y3,4) p4=polyfit(x1,y4,4) p5=polyfit(x2,y5,5) p6=polyfit(x2,y6,5) p7=polyfit(x2,y7,5) p8=polyfit(x2,y8,5) p9=polyfit(x2,y9,5) p10=polyfit(x2,y10,5) x=250:0.1:400; y11=polyval(p1,x); y22=polyval(p2,x); y33=polyval(p3,x); y44=polyval(p4,x); y55=polyval(p5,x); y66=polyval(p6,x); y77=polyval(p7,x); y88=polyval(p8,x); y99=polyval(p9,x); y101=polyval(p10,x); hold on plot(x1,y1,'o',x,y11,'k'); plot(x1,y2,'*',x,y22,'r'); plot(x1,y3,'+',x,y33,'b'); plot(x1,y4,'d',x,y44,'g'); plot(x2,y5,'x',x,y55,'c'); plot(x2,y6,'^',x,y66,'k--'); plot(x2,y7,'p',x,y77,'r--'); plot(x2,y8,'s',x,y88,'b--'); plot(x2,y9,'v',x,y99,'g--'); plot(x2,y10,'h',x,y101,'c--'); xlabel('反应的温度'),ylabel('产物的选择性'); axis([250,400,-5,70]) legend('','b1乙醇转化率','','b1C4烯烃选择性','','B2乙醇转化率','','B2C4烯烃选择性','','B3乙醇转化率','','B3C4烯烃选择性','','B4乙醇转化率','','B4C4烯烃选择性','','B6乙醇转化率','','B6C4烯烃选择性','Location','northwest');