Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лабы / L 4 / L4_1 / aaa

.doc
Скачиваний:
34
Добавлен:
16.04.2013
Размер:
254.46 Кб
Скачать

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)');

Соседние файлы в папке L4_1