Please help me someone! How to run for 3D?

1 次查看（过去 30 天）
soe min aung
soe min aung 2020-3-5
编辑： soe min aung 2020-3-9
clc
clear all
syms k k1 k2 x y
v=0.14;
L=100;
h=2;
t=5.95;
g=0.0098;
eta0=2;
w=sqrt(g*k*tanh(k*h));
A = i*k1*(50/pi)^2*(exp(i*100*k1)-1);
B = i*k2*(100/pi)^2*(exp(i*50*k2)+exp(i*150*k2));
C = i*k2*(100/pi)^2*(exp(-i*150*k2)+exp(-i*50*k2));
D = ((sin(w*t))/(w*cosh(k*h)))*((eta0*v)/(2*L));
E = (1-exp(-i*100*k1))/(i*k1);
F = (exp(-i*100*k1))/(1-((50*k1)/pi)^2);
G = (exp(i*150*k2)-exp(i*50*k2))/(i*k2);
H = 1/(1-(100*k2/pi)^2);
I = (4*sin(50*k2))/k2;
J = (exp(-i*50*k2)-exp(-i*150*k2))/(i*k2);
zeta_ = D*(E*A-F*A)*(G*B-H*B+I+J*C-H*C);
IFT1 = ifourier(zeta_,k1,x)
IFT = ifourier(IFT1,k2,y)
%For figure
xvals = 0:100;
yvals = linspace(-200,200,101);
[X,Y] = meshgrid(xvals, yvals);
Zsym = subs(IFT, {x, y}, {X, Y});
Z = double(Zsym);
surf(X,Y,Z)
  13 个评论
显示 11更早的评论
Walter Roberson
Walter Roberson 2020-3-9
subplot(1,3,1)
surf(X, Y, real(Z))
title real
subplot(1,3,2)
surf(X, Y, imag(Z))
title imaginary
subplot(1,3,3)
surf(X, Y, abs(Z))
title abs
soe min aung
soe min aung 2020-3-9
编辑：soe min aung 2020-3-9
By your advice, I have received a figure.Thank you so much for your kindness.

请先登录，再进行评论。

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

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by