%% STFT - okna : Zad.1 - 3_termin_2

Fs=100;
t=0:(1/Fs):20;
% x1 sinc o okresie 2s i amplitudzie dla t=0 równej 5 w czasie t=<0.8)s
x1=5.*sinc((t-5)/2).*(t>=0 & t< 8);
% dla t=<8,20>s
% skladowa harmoniczna o ampl. 2 i liniowo zmiennej czest. od 5 do 29 Hz
x2=2*sin(2*pi*t.*linspace(5,29,length(t)).*(t>=8 & t<=20));
% skladowa harmoniczna  o amp. 0.8 i czest. 7 Hz
x3=0.8*sin(2*pi*t*7.*(t>=8 & t<=20));
x=x1+x2+x3;

% *** STFT  *** -->
okno_A = 512;
okno_B = 32;

N=length(t);
tabA=zeros(okno_A,N-okno_A);
tabB=zeros(okno_B,N-okno_B);

% okno 512
for n=1:(N-okno_A)
    xx=x(1,n:n+okno_A-1);
    WAx=abs(fftshift(fft(xx)));
    tabA(:,n)=WAx';
end
% okno 32
for n=1:(N-okno_B)
    xx=x(1,n:n+okno_B-1);
    WAx=abs(fftshift(fft(xx)));
    tabB(:,n)=WAx;
end
% Wyswietlenie wynikow
figure;
ay = linspace(-fs/2, fs/2, okno_A);
ax = (0:N-okno_A-1)/fs;
subplot(211),imagesc(ax,ay,tabA)
title('512');

ay = linspace(-fs/2, fs/2, okno_B);
ax = (0:N-okno_B-1)/fs;
subplot(212),imagesc(ax,ay,tabB)
title('32');

%% Dekompozycja falkowa na 3 poziomie z wyk. symletu 4 rzedu
% - 2018_termin_1_G2
clear all;clc;close all;

fs=50;
t=-10:(1/fs):10;

% sygnal pi?okszta?tny o ampl. 0.6, okresie 3s i wierzcholkow dla
% t=-10:3:10
x1=(0.6/3)*mod(t-2,3);
% sygnal Gaussa o sredniej t=2s,ampl=0.75,std=2s
x2=0.75*exp(-(t-2).*(t-2)/(2*2*2));
% sygnal harmoniczny o czest. 3Hz i zmiennej liniowo ampl od 3 do 1
x3=linspace(3,1,length(t)).*(sin(2*pi*t*3));
% sygnal harmoniczny o czest. 7 Hz i ampl=0.45
x4=0.45*sin(2*pi*t*7);

x=x1+x2+x3+x4;

% 3 - 3 poziom z wykorzystanem symletu 4 rzedu
[C,S]=wavedec(x,3,'sym4');
C1=C;
C1(1+S(1):S(1)+S(2)+S(3))=0;
C1(1+S(1)+S(2):S(1)+S(2)+S(3)+S(4))=0;
C=waverec(C1,S,'sym4');
subplot(211),plot(t,x,'b',t,C,'k');
subplot(212),plot(t,x1,'r',t,x3,'b');