matlab绘制一个标准大气压下海水密度同时随温度和盐度的变化曲线（纵坐标：温度-2～40°C，横坐标：盐度0 ～ 42）

代码如下：

clear,clc;

t=[-2:0.5:40];

S1=0;

S2=7;

S3=14;

S4=21;

S5=28;

S6=35;

S7=42;

 a = 8.24493 * 10 ^ -1 - 4.0899 * 10^ -3 * t + 7.6438 * 10^-5 * t.^2- 8.2467 * 10^ -7 * t.^3 + 5.3875 * 10^-9 * t.^4;

b = -5.72466 * 10^-3 + 1.00227 * 10^ -4 * t - 1.6546 * 10^-6 * t.^ 2;

c = 4.8314 * 10 ^-4;

p = 999.842594 + 6.793952 * 10 ^ -2 * t - 9.095290 * 10^-3 * t.^ 2

    + 1.001685 * 10 ^-4 * t.^3 - 1.120083 * 10 ^ -6 * t.^ 4 + 6.536332 * 10 ^-9 * t .^5;

Y1=p+a*S1+b*S1^(2/3)+c*S1^2;

Y2=p+a*S2+b*S2^(2/3)+c*S2^2;

Y3=p+a*S3+b*S3^(2/3)+c*S3^2;

Y4=p+a*S4+b*S4^(2/3)+c*S4^2;

Y5=p+a*S5+b*S5^(2/3)+c*S5^2;

Y6=p+a*S6+b*S6^(2/3)+c*S6^2;

Y7=p+a*S7+b*S7^(2/3)+c*S7^2;

plot(t,Y1,'b',t,Y2,'g',t,Y3,'r',t,Y4,'c',t,Y5,'m',t,Y6,'y',t,Y7,'k');

legend('s=0','s=7','s=14','s=21','s=28','s=35','s=42');

xlabel('temperature/℃');

ylabel('Seawater Density/kg.m^-3');

结果：