close all;
T = [0:0.001:1];       %1000 Hz
h1 = sin(2*pi*2*T); %2Hz
h2 = sin(2*pi*10*T); %10Hz
h3 = sin(2*pi*30*T); %30Hz
s = h1 + h2 + h3;
plot(T,s,'r');
%hold on
%liniowa zmiana fazy
h1f = sin(2*pi*2*T+pi/10); %2Hz
h2f = sin(2*pi*10*T+pi/2); %10Hz
h3f = sin(2*pi*30*T+3*pi/2); %30Hz
sf = h1f + h2f + h3f;
figure, plot(T,sf,'b');
%nieliniowa zmiana fazy
h1n = sin(2*pi*2*T+pi/2); %2Hz
h2n = sin(2*pi*10*T-pi/5); %10Hz
h3n = sin(2*pi*30*T+3*pi/20); %30Hz
sn = h1n + h2n + h3n;
figure,plot(T,sn,'g');
%stałe przesunięcie fazy
h1n = sin(2*pi*2*T+pi/10); %2Hz
h2n = sin(2*pi*10*T+pi/10); %10Hz
h3n = sin(2*pi*30*T+pi/10); %30Hz
sn = h1n + h2n + h3n;
figure,plot(T,sn,'k');

figure,plot(T,h1,'k');
hold on 
plot(T,h2,'r');

figure,plot(T,h1f,'k');
hold on 
plot(T,h2f,'r');

figure,plot(T,h1n,'k');
hold on 
plot(T,h2n,'r');