Anyone can help me. Could you please check my code with attach file for fig 8 page 74 by using equation 2.30 page 67.

Question

soe min aung 2021-3-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/774267-anyone-can-help-me-could-you-please-check-my-code-with-attach-file-for-fig-8-page-74-by-using-equat

评论： darova 2021-3-17

Modeling of tsunami generation and propagation by a spreading curvilinear seismic faulting in linearized shallow-water wave theory -edit1.pdf

clc
clear all
close all
syms k1 k2
eta0 = 1;     % Maximum displacement(m)
L = 100;      % Propagate length (km)
W = 100;      % Source width (km)
v = 0.14;      % Ruptuer velocity to obtain maximum surface amplitude(km/s)
h = 2;        % Water depth (km)
g = 0.0098;   % Acceleration due to gravity (km/s^2)
t1 = 50/v;    % (s)
t = (t1);     % (s)
k = sqrt(k1.^2+k2.^2);
w = sqrt(g.*k.*tanh(k.*h));
% A = ((sin(w*t))./(w.*cosh(k*h))).*((eta0*v)./(2*L));
A = ((cos(w*t))./(cosh(k*h))).*((eta0*v)./(2*L));
B = (1-exp(-1i*100*k1))./(1i*k1);
C = ((exp(-1i*100*k1))./(1-((50*k1)./pi).^2)).*(1i*k1).*((50./pi).^2).*((exp(1i*100*k1))-1);
D = (exp(1i*150*k2)-exp(1i*50*k2))./(1i*k2);
E = (1./(1-(((100*k2)./pi).^2))).*(1i*k2).*((100./pi).^2).*(exp(1i*50*k2)+exp(1i*150*k2));
F = (4.*sin(50*k2))./k2;
G = (exp(-1i*50*k2)-exp(-1i*150*k2))./(1i*k2);
H = (1./(1-(((100*k2)./pi).^2))).*(1i*k2).*((100./pi).^2).*(exp(-1i*150*k2)+exp(-1i*50*k2)); 
K1 = 0:100;
K2 = -150:150;
for x1 = 1:101;
    for x2 = 1:301;
        k1 = K1(x1);
        k2 = K2(x2);
        if (k==0)
             eta_(x1,x2)= 200*eta0*v*t;
        else
             eta_(x1,x2)= (A.*(B-C)).*((D-E)+F+(G+H)); 
        end
    end
end
Zsym = subs(eta_);
Z = double(Zsym);
ETA = ifft2(Z)
mesh(abs(ETA))
figure(2)
mesh(real(ETA))
figure(3)
mesh(imag(ETA))
% M = 101;
% N = 301;
% Zsym = subs(eta_);
% Z = double(Zsym);
% ETA = abs(ifft2(Z));
% I = -M/2:1:M/2-1; 
% J = -N/2:1:N/2-1; 
% mesh(I,J,ETA)
% % X = abs(ifft2(Z,K1,K2));
% % X = ifftshift(X);
% % % F = [-K1/2:K1/2-1]/K1; 
% % plot(F,X)