clear all
clc

% Time array
t=(0:0.2:1800)';

% Monochromatic example
primary=1.2*sin(2*pi*t/10);
mono=[t primary];
save exp1a_input_waves.dat mono -ascii

% Modulated example
secondary=1.0*sin(2*pi*t/9.5);
prim_sec=primary+secondary;
modu=[t prim_sec];
save exp1b_input_waves.dat modu -ascii

figure(1)
plot(t,primary)

figure(2)
plot(t,prim_sec)