filip31Print

%część 1:

N=4000;
fp=1000;
t=0:1/fp:(N-1)/fp;

%sygnał - sinusoida o zmiennej częstotliwości
x=chirp(t, 100, 4.0, 400);
subplot(3,2,1);
plot(t,x);
xlabel('czas[s]');
ylabel('sygnal');




%moduł widma
Nf=2^12;
N21=Nf/2+1;
v=fft(x,Nf);
wx=abs(v);
f=linspace(0,fp/2,N21);
subplot(3,2,2);
plot(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');






%odpowiedź impulsowa
M=201;
h=fir1(M-1,0.5);
th=0:1/fp:(M-1)/fp;
subplot(3,2,3);
plot(th,h);
xlabel('czas[s]');
ylabel('odp. impulsowa');



















%moduł transmitancji
v=fft(h,Nf);
wh=abs(v);
subplot(3,2,4);
plot(f,wh(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul transmitancji');










%sygnał po filtracji
y=filter(h,1,x);
subplot(3,2,5);
plot(t,y);
xlabel('czas[s]');
ylabel('sygnal po filtracji');












%widmo amplitudy po filtracji
v=fft(y,Nf);
wy=abs(v);
subplot(3,2,6);
plot(f,wy(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('widmo aplitudy po filtracji');