2 Зертханалық жұмыс

жүктеу/скачать 1,44 Mb.

бет	1/3
Дата	19.10.2023
өлшемі	1,44 Mb.
	#186846

1 2 3

Байланысты:
2 лабка ОЦОС АСПАН

2 Зертханалық жұмыс
Жұмыс мақсаты

Қазақстан Республикасы ғылым және жоғарғы білім министрлігі
Коммерциялық емес акционерлік қоғам
«ҒҰМАРБЕК ДӘУКЕЕВ АТЫНДАҒЫ
АЛМАТЫ ЭНЕРГЕТИКА ЖӘНЕ БАЙЛАНЫС УНИВЕРСИТЕТІ»
Институт: Коммуникациялық және ғарыштық инженерия
Кафедра: Телекоммуникациялық инженерия

№2 Зертханалық жұмыс

Пән: «Сигналдарды цифрлық өндеу негіздері»

Тақырып: «Сандық сигналдарды және олардың спектрін зерттеу»
Мамандық: 6В06201-«Радиотехника, электроника және телекоммуникация»
Орындаған: Сапархан Аспандияр Тобы: РЭТк 21-5
Тексерген: Жумай Инкар Нурлыбековна

Алматы 2023

Жұмыс мақсаты: Сандық сигналдарды және олардың спектрін зерттеу, сондай-ақ оларды MATLAB-да енгізу. Фурье түрлендіруін пайдалана отырып аналогтық сигналдың аналитикалық спектрін алу. Дискретті Фурье түрлендіру (ДФТ) арқылы сандық сигнал спектрін есептеу.
1 -нұсқа
Зертханалық жұмысты орындау тәртібі мен тапсырмалар.
Кесте

№	А₁	А₂	А₃	₁	₂	₃	t₁	t₂	t₃
1	1.0	0.5	1.0	2	3	3	0	2	4

Осы жұмыстың тапсырмалары:
1) Фурье түрлендіруінің көмегімен аналогтық сигналдың (t) аналитикалық спектрін алу. Бұл жағдайда Фурье түрленуінің сызықтығын және кешіктірілген сигналдың Фурье түрлендіру сипатын қолдануға болады (бір тікбұрышты импульстің аналитикалық спектрі белгілі). Амплитудасы мен фазалық спектрін есептеңіз.
2) Сигналды құрайтын жеке импульстардың спектрін көрсету, аналогты және цифрлық сигналдардың спектрін графиктер түрінде көрсету және оларды салыстыру.
A1 = 1,0;
tau1 = 2/32;
t1 = 0/32;
f = (0:0.01:28);
Fp1 = (A1 * sin(pi * f * tau1) ./ (pi * f)) .* exp(-j * pi * f * (2 * t1 + tau1));
Ap1 = abs(Fp1);
figure;
plot(f, Ap1);
grid on;
title('Amplitude spectrum p1');

Pp2 = atan(imag(Fp1) ./ real(Fp1));

figure;
plot(f, Pp2);
grid on;
title('Phase spectrum p1');

1 сурет - p1(t) импульсінің амплитудалық спектрі

2 сурет - p1(t) импульсінің фазалық спектрі
A2 = 0.5;
tau2 = 3/32;
t2 = 2/32;
f = (0:0.01:28);
Fp2 = (A2 * sin(pi * f * tau2) ./ (pi * f)) .* exp(-j * pi * f * (2 * t2 + tau2));
Ap2 = abs(Fp2);
figure;
plot(f, Ap2);
grid on;
title('Amplitude spectrum p2');

Pp2 = atan(imag(Fp2) ./ real(Fp2));

figure;
plot(f, Pp2);
grid on;
title('Phase spectrum p2');

3 сурет – p2(t) импульсінің амплитудалық спектрі

4 сурет – p2(t) импульсінің фазалық спектрі

A3 = 1.0;

tau3 = 3/32;
t3 = 4/32;
f = (0:0.01:28);
Fp3 = (A3 * sin(pi * f * tau3) ./ (pi * f)) .* exp(-j * pi * f * (2 * t3 + tau3));
Ap3 = abs(Fp3);
plot(f, Ap3);
grid on;
title('Amplitude spectrum p3');

Pp3 = atan(imag(Fp3) ./ real(Fp3));

figure;
plot(f, Pp3);
grid on;
title('Phase spectrum p3');

5 сурет – p3(t) импульсінің амплитудалық спектрі

6 сурет – p3(t) импульсінің фазалық спектрі
A1 = 1,0;
A2 = 0.5;
A3 = 1.0;
tau1 = 2/32;
tau2 = 3/32;
tau3 = 3/32;
t1 = 0/32;
t2 = 2/32;
t3 = 4/32;
f = 0:0.01:128;

Fp1 = (A1 * sin(pi * f * tau1) ./ (pi * f)) .* exp(-1i * pi * f * (2 * t1 + tau1));

Fp2 = (A2 * sin(pi * f * tau2) ./ (pi * f)) .* exp(-1i * pi * f * (2 * t2 + tau2));
Fp3 = (A3 * sin(pi * f * tau3) ./ (pi * f)) .* exp(-1i * pi * f * (2 * t3 + tau3));

Fs = Fp1 + Fp2 + Fp3;

figure;
plot(f, abs(Fs));

grid on;

7 сурет - s(t) аналогтық сигналдың амплитудалық спектрі
N=128;
A1=1,0; A2=0.5; A3=1.0;
tau1=2/32; tau2=3/32; tau3=3/32;
t1=0/32; t2=2/32; t3=4/32;
f=0:0.01:N;
for m=1:N
S(m)=0;
end
for m=(t1*N+1):(t1+tau1)*N
S(m)=A1;
end
for m=(t2*N+1):(t2+tau2)*N
S(m)=A2;
end
for m=(t3*N+1):(t3+tau3)*N
S(m)=A3;
end
n=(0:N-1);
stem(n,S,'.')
title('Сандық S(n)')

8-сурет Сандық S(n)
N = 128;
A1 = 1;0; A2 = 0.5; A3 = 1.0;
tau1 = 2/32; tau2 = 3/32; tau3 = 3/32;
t1 = 0/32; t2 = 2/32; t3 = 4/32;
f = 0:0.01:N;
S = zeros(1, length(f));

for m = round(t1 * length(f)):(round(t1 * length(f)) + round(tau1 * length(f)))

S(m) = A1;
end

for m = round(t2 * length(f)):(round(t2 * length(f)) + round(tau2 * length(f)))

S(m) = A2;
end

for m = round(t3 * length(f)):(round(t3 * length(f)) + round(tau3 * length(f)))

S(m) = A3;
end

F = zeros(1, N);

for k = 0:N-1

W = exp(-1j * 2 * pi * f * k / N);
F(k+1) = sum(S .* W);
end

stem((0:N-1), abs(F), '.'), grid

title('Сандық спектрдің амплитудасы S(n)')

9- сурет Сандық спектрдің амплитудасы S(n)

жүктеу/скачать 1,44 Mb.

Достарыңызбен бөлісу:

1 2 3