matlab绘制绘制海水密度同时随温度、盐度和压力的变化图（xyz坐标）_matlab中使用gsw库计算海洋密度

代码：

[t,p,s]=meshgrid([-2:0.5:40],[0:10^8/84:10^8],[0:0.5: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;

pw = 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;

Y=pw+a.*s+b.*s.^(2/3)+c.*s.^2;

Kw=19652.21+148.4206.*t-2.327105.*t.^2+1.360477*10^-2.*t.^3-5.155288*10^-5.*t.^4;

Aw=3.239908+1.43713*10^-3.*t+1.16092*10^-4.*t.^2-5.77905*10^-7.*t.^3;

Bw=8.50935*10^-5-6.12293*10^-6.*t+5.2787*10^-8.*t.^2;

A=Aw+(2.283*10^-3-1.0981*10^-5.*t-1.6078*10^-6.*t.^2).*s+1.91075*10^-4.*s.^(2/3);

B=Bw+(-9.9348*10^-7+2.0816*10^-8.*t+9.1697*10^-10.*t.^2).*s;

k=Kw+(54.6746-0.603459.*t+1.09987*10^-2.*t.^2-6.1670*10^-5.*t.^3).*s+(7.944*10^-2+1.6483*10^-2.*t-5.3009*10^-4.*t.^2).*s.^(2/3);

K=k+A.*p.*10^-5+B.*(p.*10^-5).^2;

C=Y./(1-(p.*10^-5)./K);
xslice = [-2:6:40];   
yslice = 0;
zslice = 0;
slice(t,p,s,C,xslice,yslice,zslice)
colorbar

xlabel('温度/℃');

ylabel('压强/Pa');

zlabel('盐度');

结果：