clc
clear
g = 9.81;
h = 15;
l = 500;
n = 1;
L1 = 2*l/n;
disp('At node = 1');
disp(['Wavelength L = ' num2str(L1)]);
T1 = (2*l)/(n*sqrt(g*h));
disp(['Wave Period at T = ' num2str(T1)]);
g = 9.81;
h = 15;
l = 500;
a = 2;
L2 = 2*l/a;
disp('At node = 2');
disp(['Wavelength L = ' num2str(L2)]);
T2 = (2*l)/(a*sqrt(g*h));
disp(['Wave Period at T = ' num2str(T2)]);
g = 9.81;
h = 15;
l = 500;
b = 3;
L3 = 2*l/b;
disp('At node = 3');
disp(['Wavelength L = ' num2str(L3)]);
T3 = (2*l)/(b*sqrt(g*h));
disp(['Wave Period at T = ' num2str(T3)]);
t=0;
H = 1;
x = 0:l;
k1 = 2*pi/L1;
k2 = 2*pi/L2;
k3 = 2*pi/L3;
sigma1 = 2*pi/T1;
sigma2 = 2*pi/T2;
sigma3 = 2*pi/T3;
N1 = H/2*cos(k1*x)*cos(sigma1*t);
N2 = H/2*cos(k1*x)*cos(sigma2*t);
N3 = H/2*cos(k1*x)*cos(sigma3*t);
subplot(3,1,1);
N1 = H/2*cos(k1*x)*cos(sigma1*t);
plot(x,N1, 'r')
hold on
yline(0);
ylim([-1 1])
title ('mode n = 1, L = l/2')
xlabel ('x(m)')
ylabel ('eta (m)')
legend('H/2*cos(k1*x)*cos(sigma2*t)')
subplot(3,1,2)
N2 = H/2*cos(k2*x)*cos(sigma2*t);
plot (x,N2,'r');
hold on
yline(0);
ylim([-1 1])
title ('mode n = 2, L = l')
xlabel ('x(m)')
ylabel ('eta (m)')
legend('H/2*cos(k2*x)*cos(sigma2*t)')
subplot(3,1,3)
N3 = H/2*cos(k3*x)*cos(sigma3*t);
plot(x, N3, 'r')
hold on
yline(0);
ylim([-1 1])
title ('mode n = 3, L= 2l/3')
xlabel ('x(m)')
ylabel ('eta (m)')
legend('H/2*cos(k3*x)*cos(sigma2*t)')
0 Comments
Sign in to comment.