Subplot displacement against time
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
%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 commentaires
Réponses (1)
CS Researcher
le 3 Mai 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 commentaires
CS Researcher
le 4 Mai 2016
if x is the vector for displacement and t is the vector for time, you can just use:
plot(t,x);
Voir également
Catégories
En savoir plus sur Line Plots dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!