PPS kolokwium

%1. (17)Wczytanie sygnału i decymacja:
[x,fp]=wavread('S1.wav');
y=decimate(x,5);
subplot(111);
plot(y);

%2. (18)Policzenie widma amplitudowego sygnału:
Nf=1024;
v=fft(x,Nf);
w=abs(v);
f=linspace(0,fp,Nf);
subplot(212);
plot(f,w);
ylabel('czestotliwosc[Hz');
%xlabe('modul widma');

Nf=2^11;
N21=Nf/2+1;
w=fft(x,Nf) %x to sygnał
v=abs(w);
Nf=length(v);
tf=linspace(0,Fs/2,N21); %Fs czestotliwosc probkowania
plot(tf,v(1:N21));
xlabel('czestotliwosc[Hz]')
ylabel('modul widma')


%3. (19)Filtr fir o zadanej długości i f granicznej
h=fir1( 51,[1400]);
plot(h);
freqz(h,1);
y=filter(h,1,x);


%4. (20)Filtracja i wyznaczenie modułu widma sygnału po filtracji

h=fir1(21,[300/500]);
plot(abs(fft(h(512))),'r')
y=filter(h,1,x);
plot(x)
hold on





%5. (21)Generacja sygnału o zadanej długości i częstotliwości
fp=150; %próbkowanie
t1=0:1/fp:2;
f1=5;
y1=sin(t1*2*pi*f1);
t2=2:1/fp:4;
f2=3.5;
fi2=pi/4;
y2=sin(t2*2*pi*f2+fi2);
t3=4:1/fp:8;
A3=1;
y3=2*rand(1,601)-A3; %rozkład normalny szumu
t4=8:1/fp:10;
f4=5;
fi4=pi/2;
y4=sin(2*pi*f4*t4+fi4);
y=[y1 y2 y3 y4];
t=[t1 t2 t3 t4];
plot(t,y);

%6. (22)Policzenie autokorelacji sygnału z 5tki z danym tau_max
[r,T]=xcorr(x,400);
plot(T,r);

%7. (23)Korelacja sygnału z 5tki z zadanym sygnalem sinusoidalnym
f7=5;
fi7=pi;
t7=0:1/fp:10;
y7=sin(2*pi*f7+fi7);
[r,T]=xcorr(y,y7,400);
plot(T,r);


%Kod z zajęć
% WIDMO AMPLITUDOWE SYGNAŁU W ZAKRESIE 0-FP
N=200;
fp=1000;
t=0:1/fp:(N-1)/fp;
%suma sinusow i szumu
x=6*sin(2*pi*250*t)+3*sin(2*pi*350*t)+2*randn(1,N);
subplot(211);
plot(t,x)
xlabel('czas[s]');
ylabel('sygnał');
%widmo
Nf=1024;
v=fft(x,Nf);%widmo zespolone
w=abs(v);%;widmo amplitudowe
f=linspace(0,fp,Nf); %generowanie osi czestotliwosci
subplot(212);
plot(f,w)
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

%generowanie szumu
x=rand(1,1000);
plot(x);
hold on
x=randn(1,1000);
plot(x,'r');
x=rand(1,1000);
%sygnał o zmiennej częstotliwości
x=chirp(t,fo,t1,f1)

%FFT
y=fft(x);
y=abs(abs(y);
y=unurap(angle(y));
y=real(y);
y=imag(y);

%Filtr FIR- skonczona odpowiedz impulsowa
h=fir1(K,fg);
h=fir(K,[fd fg]);
h=fir1(K,fd,'high');
%filtracja
y=filter(h,1,x);


%decymacja -zmiejszenie czestliwosci próbkowania
[x,y]=sinus(1,2,200,0,2)
plot(x,y);
y=decimate(x,y)
ny=...;
plot(ny,y);

%Histogram-jakie prawdopodobienstwo wystapienia sygnalu o danej
%czestotliwosci
sz=randn(1,1000);
hist(sz,51)%51 liczba przedziałów(liczone ile probek wpada do hist)
hist(x,101);

%Autokorelacja mówi nam o zależności w obrębie jednego sygnału
r=xcorr(x,tmax)
[x,n]=sinus(1,2,2000,0,2);
[r,t]=xcor(x);
plot(t,r);
[r,t]=xcorr(x,100);
plot(t,r);

%Korelacja wzajemna
[r,t]=xcorr(x,y,tmax);
[r,t]=xcorr(x,sz,200);
plot(t,r);

%Wave read
[x,fs]=wavread('mbi04popis.wav');
%N - nowa zmienna pomocnicza
y=fft(x,fs);
w=abs(y);
f=linspace(0,fs/2,N);
subplot(211)
plot(f,w(1:N));grid on; zoom on;
xlabel('Czestotliwość [Hz]');
ylabel('Moduł widma');
title('Wykres widma amplitudowego');