Dispersion Curves For Love Waves
Dispersion Curves For Love Waves
Dispersion Curves For Love Waves
Satyajeet Behera
Homework-08
In the code used, I have taken the value of β1 = 1, and β2 = 1.8. The dispersion
curves plot is given below:
1
Appendix
The code is in MATLAB, and is given below:
clear all
clc
%% Range o f x and y
x = l i n s p a c e (0 , 60 , 10000);
y = linspace (1 , 2 , 10000);
[ X, Y] = meshgrid ( x , y ) ;
%% D i s p e r s i o n Equation
Z = tan ( (X∗ 1 . 3 ) . /Y . ∗ s q r t (Y. ˆ 2 − 1 ) ) − 3 . ∗ s q r t ( 1 − (Y/ 1 . 8 ) . ˆ 2 ) . / s q r t (Y. ˆ 2 −
% F i n d i n g where t h e e q u a t i o n e q u a l s a p p r o x i m a t e l y z e r o
threshold = 0.01;
i d x = abs ( Z ) < t h r e s h o l d ;
%% P l o t
figure ;
p l o t (X( i d x ) , Y( i d x ) , ’ b . ’ ) ;
x l a b e l ( ’ \ omega ’ ) ;
y l a b e l ( ’ c /\ b e t a 1 ’ ) ;
t i t l e ( ’ P l o t o f D i s p e r s i o n e q u a t i o n f o r Love Waves ’ ) ;
g r i d on ;
s a v e a s ( g c f , ’ p l o t . png ’ ) ;