clear all;close all;clc
% load handel.mat
% filename = 'handel.wav';
% audiowrite(filename,y,Fs);
duration=8;sampling=8192;minf=100;maxf=2000;
L = sampling*duration;             % Length of signal
time = linspace(0, duration, L);
% signal = chirp(time, minf, time(end), maxf);

f=1000;signal1=sin(2*pi*f*time);
% f=500;
signal2=sin(2*pi*f*time);
% f=1000;
signal3=sin(2*pi*f*time);
signal=signal1+signal2+signal3;

subplot(211)
plot(time, signal, 'b*-');title('Signal in time')
grid on;


subplot(212)
Fs=sampling;
Y = fft(signal);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
plot(f,P1)
% return

filename = 'handel.wav';
audiowrite(filename,signal,Fs);
% return
%%%%%%%%

x=linspace(0,duration,duration*sampling)

[y,Fs] = audioread('handel.wav',[1,duration*sampling]);

tic
sound(y,Fs);
toc