MATLAB: How to plot this equation to obtain the figure

3dhomework

Best Answer

N = 50;   %subdivisions per dimension
xmin = -1; xmax = 101;
ymin = -160; ymax = 160;
tmin = 0; tmax = 3600;
xvec = linspace(xmin, xmax, N);
yvec = linspace(ymin, ymax, N);
tvec = linspace(tmin, tmax, N);
[X, Y, T] = ndgrid(xvec, yvec, tvec);
maskx = 0 <= X & X <= 100;
masky1 = -150 <= Y & Y < -50;
masky2 = -50 <= Y & Y < 50;
masky3 = 50 <= Y & Y <= 150;
xi = zeros(size(X));
mask1 = maskx & masky1;
mask2 = maskx & masky2;
xi(mask1) = xi0 * v*T(mask1)/(2 * L) .* (1 - cos(pi/50*X(mask1))) .* (1 - cos(pi/100*(Y(mask1) + 150)));
xi(mask2) = xi0 * V*T(mask2) / L .* (1-cos(pi/50*X(mask2)));
xi(mask3) = xi0 * v*T(mask3)/(2 * L) .* (1 - cos(pi/50*X(mask3))) .* (1 - cos(pi/100*(Y(mask3) - 150)));
random_time_idx = randi(length(tvec));
random_time = tvec(random_time_idx);
x_for_t = X(:,:,random_time_idx);
y_for_t = Y(:,:,random_time_idx);
zi_for_t = xi(:,:,random_time_idx) / xi0;
surf(x_for_t, y_for_t, zi_for_t)
xlabel('x (km)');
ylabel('y (km)');
zlabel('$\frac{\zi(x,y,t1)}{\zi_0}', 'interpreter', 'latex')
title( sprintf('time = %.2f', random_time) );

Related Solutions

MATLAB: How to plot

YOu need to proceed like this:

m = 300 ; n = 100 ; 
x = linspace(0,100,m) ; 
eta0 = 2 ; 
L = 100 ; 
W = 100 ; 
v = 0.14 ; 
t = 5.95 ; 
y1 = linspace(-150,-50,n) ; 
[X,Y1] = meshgrid(x,y1) ; 
T1 = eta0*v*t/(2*L)*(1-cos(pi/50*X)).*(1-cos(pi/100*(Y1+150))) ; 
y2 = linspace(-50,50,n) ; 
[X,Y2] = meshgrid(x,y2) ; 
T2 = eta0*v*t/L*(1-cos(pi/50*X)) ; 
y3 = linspace(50,150,n) ; 
[X,Y3] = meshgrid(x,y3) ; 
T3 = eta0*v*t/(2*L)*(1-cos(pi/50*X)).*(1+cos(pi/100*(Y3-150))) ; 
X = [X ; X; X] ; 
Y = [Y1 ; Y2 ;Y3] ; 
eta = [T1 ; T2 ; T3] ; 
surf(X,Y,eta) ;

You need to make some changes in the above code...which you should figure out.

This is the intitial deformation for tusnami.By the way I am interested to know the reference of this document; can you share the document.

MATLAB: How to write piecewise function in matlab for the following function. Please give me matlab code.

Creating an anonymous function for ‘T’ is much easier and faster than using piecwise and the Symbolic Math Toolbox for this:

zeta0 = 0.5
L = 100;
W = 100;
h = 2;
g = 0.0098;
vt=sqrt(g*h);
c=1;
v=c*vt;
t1 = L/v;
t = t1;
T1 = @(x,y) zeta0*v*t/(2*L)*(1-cos(pi/50*x)).*(1-cos(pi/100*(y+150))) ; 
T2 = @(x,y) zeta0*v*t/L*(1-cos(pi/50*x)) ;
T3 = @(x,y) zeta0*v*t/(2*L)*(1-cos(pi/50*x)).*(1-cos(pi/100*(y-150))) ; 
T = @(x,y) ((0<=x) & (x<=100) & (-150<=y) & (y<=-50)).*T1(x,y) + ((0<=x) & (x<=100) & (-50<=y) & (y<=50)).*T2(x,y) + ((0<=x) & (x<=100) & (50<=y) & (y<=150)).*T3(x,y);
[X,Y] = ndgrid(-100:5:200);
figure
surf(X,Y,T(X,Y))
grid on

See Anonymous Functions for details, if you are not familiar with them

Best Answer

Related Solutions

MATLAB: How to plot

MATLAB: How to write piecewise function in matlab for the following function. Please give me matlab code.

Related Question