- /
-
Ocean Simulator
on 29 Nov 2023
- 35
- 79
- 0
- 4
- 1211
drawframe(1);
Write your drawframe function below
function drawframe(f)
rng(13); % setting seed for random numbers
i=meshgrid(1:129);
x=linspace(-10,10,129);
X=meshgrid(x);
a=meshgrid(linspace(-1,1,129)*pi*129/200) ;
W= sqrt(hypot(a,a')*10);
K=a.^2+a'.^2;
k=sqrt(K);
w=a./k*cosd(30)+a'./k*sind(30);
P=1e-7./K.^2.*exp(-1e-6./K).*w.^2;
P(K==0|w<0)=0 ;
h=sqrt(P/2).*(randn(129)+1i*randn(129));
w=exp(1i*W*f/2);
Z=real(ifft2(h.*w+conj(rot90(h,2)).*conj(w)).*(2*(mod(i+i',2)==0)-1));
surf(X,X',Z); hold on;
m=-2e-4;
patch([10 -10 x],-10+[0 0 0*x],[m m Z(1,:)],[m m Z(1,:)]);
patch([0 0 0*x]-10,[10 -10 x],[m;m;Z(:,1)],[m;m;Z(:,1)]);
colormap(linspace(0.4,1,25)'*[0 0 1]);
shading interp;
axis square off;
axis([-10 10 -10 10 m -m]);
lightangle(-45,30);
hold off;rng(13); % setting seed for random numbers
i=meshgrid(1:129);
x=linspace(-10,10,129);
X=meshgrid(x);
a=meshgrid(linspace(-1,1,129)*pi*129/200) ;
W= sqrt(hypot(a,a')*10);
K=a.^2+a'.^2;
k=sqrt(K);
w=a./k*cosd(30)+a'./k*sind(30);
P=1e-7./K.^2.*exp(-1e-6./K).*w.^2;
P(K==0|w<0)=0 ;
h=sqrt(P/2).*(randn(129)+1i*randn(129));
w=exp(1i*W*f/2);
Z=real(ifft2(h.*w+conj(rot90(h,2)).*conj(w)).*(2*(mod(i+i',2)==0)-1));
surf(X,X',Z); hold on;
m=-2e-4;
patch([10 -10 x],-10+[0 0 0*x],[m m Z(1,:)],[m m Z(1,:)]);
patch([0 0 0*x]-10,[10 -10 x],[m;m;Z(:,1)],[m;m;Z(:,1)]);
colormap(linspace(0.5,1,25)'*[0 0.2 1]);
shading interp;
axis square off;
axis([-10 10 -10 10 m -m]);
lightangle(-45,30);
hold off;
end