untitled.m

能根据温盐资料计算海洋的浮性频率能根据温盐资料计算海洋的浮性频率能根据温盐资料计算海洋的浮性频率

[Pa,Ta,S,t,sb,se,sr,yu,ui]=textread('data.dat','%f %f %f %f %f %f %f %f %f');
clear sb
clear se
clear sr
clear yu
clear ui
P(1)=1.01325*10^5;
pi=3.1415926;
h(1)=0;
for i=2:1:size(Pa)+1
    h(i)=Pa(i-1);
end
clear Pa
g=9.780318*(1+0.0053024*(sin(-25.83*2*pi/180))^2);
n=10^(-5);%压力匹配系数
T=1.00024*Ta;%t90换算到t68
for i=2:size(S)+1;
        denw(i)=999.842594+6.793952*10^(-2)*T(i-1)-9.059290*10^(-3)*(T(i-1))^2+...
            1.001685*10^(-4)*(T(i-1))^3-1.120083*10^(-6)*(T(i-1))^4+...
            6.536332*10^(-9)*(T(i-1))^5;
        A(i)=8.24493*10^(-1)-4.0899*10^(-3)*T(i-1)+7.6438*10^(-5)*(T(i-1))^2-...
            8.2467*10^(-7)*(T(i-1))^3+5.3875*10^(-9)*(T(i-1))^4;
        B(i)=-5.72466*10^(-3)+1.0227*10^(-4)*T(i-1)-1.6546*10^(-6)*(T(i-1))^2;
        C=4.8314*10^(-4);
        den(i)=denw(i)+A(i)*S(i-1)+B(i)*(S(i-1))^(3/2)+C*(S(i-1))^2;%标准大气压下求得的密度
        KW(i)=19652.21+148.4206*T(i-1)-2.327105*(T(i-1))^2+...
            1.360477*10^(-2)*(T(i-1))^3-5.155288*10^(-5)*(T(i-1))^4;
        AW(i)=3.239908+1.43713*10^(-3)*T(i-1)+1.16092*10^(-4)*(T(i-1))^2-...
            5.77905*10^(-7)*(T(i-1))^3;
        BW(i)=8.50935*10^(-5)-6.12293*10^(-6)*T(i-1)+5.2787*10^(-8)*(T(i-1))^2;
        K0(i)=KW(i)+(54.6746-0.603459*T(i-1)+1.09987*10^(-2)*(T(i-1))^2-...
            6.167*10^(-5)*(T(i-1))^3)*S(i-1)+(7.944*10^(-2)+1.6483*10^(-2)*T(i-1)-...
            5.3009*10^(-4)*(T(i-1))^2)*(S(i-1))^(3/2);
        A1(i)=AW(i)+(2.2838*10^(-3)-1.0981*10^(-5)*T(i-1)-1.6078*10^(-6)*(T(i-1))^2)*S(i-1)+...
            1.91075*10^(-4)*(S(i-1))^(3/2);
        B1(i)=BW(i)+(-9.9348*10^(-7)+2.0816*10^(-8)*T(i-1)+9.1697*10^(-10)*(T(i-1))^2)*S(i-1);
        P(i)=den(i)*g*(h(i)-h(i-1))+P(i-1);
        K(i)=K0(i)+A1(i)*(n*P(i))+B1(i)*(n*P(i))^2;
        denth(i)=den(i)*(1-n*P(i)/K(i))^(-1);% 高压下求得的密度
end
for i=1:size(S)
    d(i)=den(i+1);
    de(i)=denth(i+1);
end
dz=2;
d=reshape(d,86,1);
de=reshape(de,86,1);
j=1;
for i=1:size(d)-1
    N(j)=real((-g/d(i)*(d(i)-d(i+1))/dz)^0.5);
    N(j+1)=N(j);
    j=j+2;
end
save 'e:\homework\p.mat' N
clear