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

1 次查看（过去 30 天）
soe min aung
soe min aung 2020-7-9
回答： KSSV 2020-7-9
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))

请先登录，再回答此问题。

采纳的回答

KSSV
KSSV 2020-7-9
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))

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by