% script stat

%% Sätt upp en del parameterar

f = 1000;                   % Signalens frekvens
omega = 2*pi*f;             % Signalens vinkelfrekvens
fs = 44100;                 % Sampelfrekvens
Ts = 1/fs;                  % Sampelintervall
k = 0:1:(2*fs);             % Diskret tidsvariabel över 2 sek
s1 = cos(omega*k*Ts);       % Signalen skapas
figure(1);plot(s1,'o-');    % Signalen ritas upp

%% Spela upp signalen s1

sound(s1,fs);

%% Skapa två nya signaler och spela upp

s2=sin(omega*k*Ts);
s3=cos(2*omega*k*Ts);
sound(s1,fs);
sound(s2,fs);
sound(s3,fs);

%% Rita upp nåpgra signaler

figure(1);subplot(3,1,1);plot(s1+s3,'o-');
figure(1);subplot(3,1,2);plot(s2+s3,'o-');

%% Spela upp samma signaler

sound(s1+s3,fs);pause;
sound(s2+s3,fs);pause;
sound(s3,fs);

%% Rita upp och spela upp ännu en variant

figure(1);subplot(3,1,3);plot(s1+s2+s3,'o-');
sound(s1+s3,fs);sound(s1+s2+s3,fs);

%% Undersökning av vikningsdistortion

fs2=10000;
Ts2=1/fs2;
k=0:1:(2*fs2);
for frek=1000:1000:10000,
  pause;
  signal=cos(2*pi*frek*k*Ts2);
  figure(2);plot(signal,'o-');
  title(sprintf('Signalens frekvens = %d Hz', frek));
  sound(signal,fs2);
end,
