G hghhgfhfghghhhggggh
Лабораторная работа № 4
Вариант 5 (по номеру копьютера).
Рекута Анна Александровна ЕКТ-41
vn=30; n=256; t=0:1/vn:10;
x1=sin(2*pi*t/2); x2=sin(2*pi*t/5); x3=sin(2*pi*10*t); x=x1+x2+x3;
q1=[0 0.4/vn 0.4/vn 1]; w1=[1 1 0 0];
q2=[0 5/vn 5/vn 1]; w2=[1 1 0 0];
q3=[0 5/vn 5/vn 1]; w3=[1 1 0 0];
q4=[0 5/vn 5/vn 1]; w4=[0 0 1 1];
%LLF 1/5
N=100; Vs=0.4;
b=FIR1(N,Vs/vn); a=[]; a=[1 a];
y=FILTFILT(b,a,x); h=abs(freqz(b,a,n)); ff=(1/n)*(0:n-1);
%HF 10
N1=50; Vs1=5;
b1=FIR1(N1,Vs1/vn,'high'); a1=[]; a1=[1 a1];
y1=FILTFILT(b1,a1,x); h1=abs(freqz(b1,a1,n));
%LF 1/2 1/5
N2=10; Vs2=5;
b2=FIR1(N2,Vs2/vn); a2=[]; a2=[1 a2];
y2=FILTFILT(b2,a2,x); h2=abs(freqz(b2,a2,n)); ff=(1/n)*(0:n-1);
%LF butter 1/2 1/5
N3=2; Vs3=5;
[b3,a3]=butter(N3,Vs3/vn);
y3=FILTFILT(b3,a3,x); h3=abs(freqz(b3,a3,n)); ff=(1/n)*(0:n-1);
subplot(311), plot(t,x);,xlabel('sin(2*pi*t/2)+sin(2*pi*t/5)+sin(2*pi*t*10)'); title('PRODUCED BY ANNA RECUTA')
subplot(312), semilogx(ff,h,q1,w1);xlabel('N=100');
subplot(313), plot(t,y);xlabel('sin(2*pi*t/2)');
pause
subplot(411),plot(ff,h2,q2,w2);xlabel('N=10');title('PRODUCED BY ANNA RECUTA')
subplot(412),plot(t,y2);xlabel('sin(2*pi*t/2)+sin(2*pi*t/5)');
subplot(413), plot(ff,h3,q3,w3);xlabel('BATTERVORD N=2');
subplot(414), plot(t,y3);xlabel('sin(2*pi*t/2)+sin(2*pi*t/5)');
pause
subplot(211),plot(ff,h1,q4,w4);xlabel('N=50');
subplot(212), plot(t,y1);xlabel('sin(2*pi*t*10)');