Subplot displacement against time
3 views (last 30 days)
Show older comments
%02/05/2016
clear
%Define the total time for simulation
T_total=10;
%Set the computation increment
dt=0.1;
%calculate the number of steps
ntstep=T_total/dt+1;
%Define Input Waves
alphaa=1; %velocity of wave a
alphab=2; %velocity of wave a
fa=1; %frequency of wave a
fb=0.5; %frequency of wave b
ampa=6; %amplitude of wave a in cm
ampb=3; %amplitude of wave b in cm
%set the source depth
depth=1;
%Calc wave parameters using input
Ta=1/fa; %period of wave a
Tb=1/fb; %period of wave a
wa=2*pi*fa; %angular frequency of wave a
wb=2*pi*fb; %angular frequency of wave b
la=Ta*alphaa;%wavelength of wave a
lb=Tb*alphab;%wavelength of wave b
ka=2*pi/la; %wavenumber of wave a
kb=2*pi/lb; %wavenumber of wave b
x=(-10:1:10); %x
y=(-10:1:10); %y
z=(-20:1:0); %z
t=zeros(length(x),length(y), length(z)); %start time
% setup the figure
figure(1);
%label axes
xlabel('Distance (m)');
ylabel('Distance (m)');
zlabel('Depth (m)');
hold on
[X,Y,Z]=meshgrid(x,y,z);
R=sqrt(X.^2+Y.^2+(Z+depth).^2);
%loop over nsamp
for n=1:ntstep
t=t+dt; %time
arga=wa*t-ka*R; %argument of wave a
argb=wb*t-kb*R; %argument of wave b
a=ampa*sin(arga)./max(1,R);
a(R>t.*alphaa)=0; % causality condition
a(R<t.*alphaa-la/2)=0; % limit lenght of input wave to one half cycle
b=ampb*sin(argb)./max(1,R);
b(R>t.*alphab)=0; % causality condition
b(R<t.*alphab-lb/2)=0; % limit lenght of input wave to one half cycle
figure(1);
hold off;
wf=a+b;
slice(X,Y,Z,wf,x,y,z);
axis square;
shading INTERP
caxis([0,1]);
colormap([[0:0.1:1]',zeros(length([0:0.1:1]'),1),zeros(length([0:0.1:1]'),1)]);
col=colorbar;
alpha(0.05);
hold on;
%label axes
xlabel('Distance (km)');
ylabel('Distance (km)');
zlabel ('Depth (km)');
ylabel(col,'Amplitude (cm)'); % label the colour bar
contour(X(:,:,21),Y(:,:,21),wf(:,:,21)); %contour plots at the surface only
figure(2);
subplot(2,2,1);
plot(t(:,:,21), a(:,:,21));
axis square;
title(strcat('time = ', num2str(t(1)), ' s'));
pause (dt);
end
I want to create two subplots (instead of figure 1 & figure 2), with figure 2 being the displacement (the wave amplitude) against time at the surface at any chosen point [x,y,0]
I want there to be a dot or 'marker' and it just sort of draws the graph as time progresses if that makes sense?
Any help appreciated.
0 Comments
Answers (1)
CS Researcher
on 3 May 2016
You should do this:
figure(1);
subplot(2,1,1);
Change figure(2) to figure(1) and change subplot(2,2,1) to
subplot(2,1,2);
2 Comments
CS Researcher
on 4 May 2016
if x is the vector for displacement and t is the vector for time, you can just use:
plot(t,x);
See Also
Categories
Find more on Line Plots in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!