matlab艾克曼螺旋有限深海漂流

1、用matlab可视化艾克曼漂流螺旋和艾克曼螺线clc,clear%艾克曼螺线u10=9;v10=10;fai=21;pa=1.205;ro=1024;Az=0.01;w=7.29*10A(-5);z=li nspace(-50,0,30);zz=-z;a=sqrt(w*si nd(fai)/Az);DO=pi/a;t=pa*sqrt(u10A2+v10A2)*v10;u=t./(sqrt(2).*a.*Az.*ro).*exp(-pi/D0.*zz).*cos(pi./4-pi./D0.*zz);v=t./(sqrt(2).*a.*Az.*ro).*exp(-pi/D0.*zz).*si n(

2、 pi./4-pi./D0.*zz);w1=zeros(1,le ngth(v);x=zeros(1,30);y=zeros(1,30);plot3(x,y,z,'red')hold on%a1=le ngth(v),a2=le ngth(wl)quiver3( x,y, z,u, v, w1,0);%axes('positio n',0.1 0.2 0.3 0.4);%axis equalhold on%z0=-25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25 -25

3、-25;z0=(-50).*o nes(1,30);quiver3( x,y,zO,u, v, w1,0);title('Ekman 螺线')clear,clc%有限深海漂流u10=9;v10=10;fai=21;pa=1.205;u=zeros(1,30);v=zeros(1,30);ro=1024;Az=0.01;w=7.29*10A(-5);a=sqrt(w*si nd(fai)/Az);DO=pi/a;hd=0.1,0.25,0.5,1.25;%hd=h/DOh=DO.*hd;%zz=li nspace(-30,0,30);%z=-zz;%深度为0到-50米z1=li

4、nspace(0,6,30);z2=li nspace(0,15,30);z3=li nspace(0,30,30);z4=li nspace(0,76,30);z=z1;z2;z3;z4;t=pa*sqrt(u10A2+v10A2)*v10;%t 代表 tao;for i=1:4;A=(t./(a.*Az.*ro).*(cosh(a.*h(i).*cos(a.*h(i)+si nh(a.*h(i).*si n(a.*h(i)./. (cosh(2.*a.*h(i)+cos(2.*a.*h(i);B=(t./(a.*Az.*ro).*(cosh(a.*h(i).*cos(a.*h(i)-si

5、nh(a.*h(i).*si n(a.*h(i)./.(cosh(2.*a.*h(i)+cos(2.*a.*h(i);%A0=t./(a.*Az.*ro);A 仁 cosh(a.*h(i).*cos(a.*h(i);%A2=s in h(a.*h(i).*si n(a.*h(i);A3=cosh(2.*a.*h(i)+cos(2.*a.*h(i);%A=A0.*(A1+A2)./A3%B=A0.*(A1-A2)./A3u=A.*s in h(a.*(h(i)-z(i,:).*cos(a.*(h(i)-z(i,:)-B.*cosh(a.*(h(i)-z(i,:).*si n(a.*(h(i)-z(i,:); v=A.*cosh(a.*(h(i)-z(i,:).*si n(a.*(h(i)-z(i,:)+B.*s in h(