Someone help me. How to plot for this code to obtain figure?

1 view (last 30 days)
clc
clear all
close all
syms x y
zeta0 = 1;
h = 2;
g = 0.0098;
vt=sqrt(g*h);
c=1;
v=c*vt
t1 = 50/v;
t2=200/v;
t=t1;
zeta=(1/(2*pi)^2)*200*zeta0*v*t*((exp(i*100*x)/(i*x))-(1/(i*x)))*((exp(-i*50*y)-exp(-i*150*y))/(i*y))
x = 0:100;
y = linspace(-150,-50,101);
[X1,Y1] = meshgrid(x,y);
Zsym = subs(zeta, {x,y},{X1,Y1});
Z = double(Zsym);
surf(X1,Y1,abs(Z))

Accepted Answer

KSSV
KSSV on 9 Jul 2020
clc
clear all
close all
zeta0 = 1;
h = 2;
g = 0.0098;
vt=sqrt(g*h);
c=1;
v=c*vt
t1 = 50/v;
t2=200/v;
t=t1;
x = 0:100;
y = linspace(-150,-50,101);
[X,Y] = meshgrid(x,y);
zeta=(1/(2*pi)^2)*200*zeta0*v*t*((exp(i*100*X)./(i*X))-(1./(i*X))).*((exp(-i*50*Y)-exp(-i*150*Y))./(i*Y)) ;
surf(X,Y,abs(zeta))

More Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!