Задание 4 к практическим занятиям (ММ РУиС) - История изменений

Korogodin: /* Задача 4.1 */

2015-03-25T12:44:34Z

‎Задача 4.1

Korogodin: /* Задача 4.1 */

2015-03-25T12:43:55Z

‎Задача 4.1

Korogodin: /* Задача 4.1 */

2015-03-25T12:39:58Z

‎Задача 4.1

Korogodin в 19:00, 1 апреля 2014

2014-04-01T19:00:31Z

Korogodin: /* Задача 4.1 */

2014-04-01T18:53:16Z

‎Задача 4.1

Korogodin в 18:52, 1 апреля 2014

2014-04-01T18:52:28Z

Dneprov в 17:22, 27 марта 2014

2014-03-27T17:22:49Z

Dneprov: /* Задача 4.1 */

2014-03-27T17:22:37Z

‎Задача 4.1

Korogodin: /* Задача 4.1 */

2014-03-24T11:41:21Z

‎Задача 4.1

Korogodin: /* Задача 4.1 */

2014-03-18T17:50:46Z

‎Задача 4.1

@@ Строка 21: / Строка 21: @@
 t = 0:Td:Tmod; % [sec]
-S = cos(2*pi* (500 + 0 * rand(1,1)) * t);
+df = 100;
 figure;

@@ Строка 29: / Строка 29: @@
 </source>
-Подсказка: <code>Tmod = T - Td</code>
+Попробуйте применить оконную обработку для улучшения результатов.
 == Задача 4.2 ==

← Предыдущая		Версия 12:39, 25 марта 2015
Строка 8:		Строка 8:

	== Задача 4.1 ==		== Задача 4.1 ==
		+
		+	'''Постановка задачи''':
		+
		+	Угадайте частоту загаданного сигнала и выведите значение в Гц в Command Window:
		+
		+	<source lang = matlab>
		+	clear all; close all; clc;
		+
		+	Fd = 40e3; % [Hz]
		+	Td = 1/Fd; % [sec]
		+	Tmod = 0.3; % [sec]
		+	t = 0:Td:Tmod; % [sec]
		+
		+	S = cos(2pi (500 + 0 * rand(1,1)) * t);
		+
		+	figure;
		+	plot(t, S)
		+
		+	% freq of S - ?
		+	</source>
		+
		+	Подсказка: <code>Tmod = T - Td</code>
		+
		+
		+	== Задача 4.2 ==

	'''Постановка задачи''':		'''Постановка задачи''':

@@ Строка 6: / Строка 6: @@
 * рассмотреть пример использования ДПФ в алгоритмах обработки сигналов;
 * научиться применять весовые окна.
 == Задача 4.1 ==

@@ Строка 111: / Строка 111: @@
 subplot(5,2,9)
-stem(t*1e6, abs(y_rej));
+stem(t*1e6, real(y_rej));
 xlabel('t, \mus');
 ylabel('y after rejector');

← Предыдущая		Версия 18:52, 1 апреля 2014
Строка 120:		Строка 120:
	ylabel('fft(y after rejector)');		ylabel('fft(y after rejector)');
	</source>		</source>
		+
		+
		+	[[Категория: ММ РУиС (дисциплина)]]

← Предыдущая		Версия 17:22, 27 марта 2014
Строка 120:		Строка 120:
	ylabel('fft(y after rejector)');		ylabel('fft(y after rejector)');
	</source>		</source>
−	\|}

← Предыдущая		Версия 17:50, 18 марта 2014
Строка 23:		Строка 23:

	Требуется построить осциллограммы и спектры сигналов на входе и выходе режектора для различных значений <math>f_0</math>, с использования оконной обработки (окно Парзена) и без неё.		Требуется построить осциллограммы и спектры сигналов на входе и выходе режектора для различных значений <math>f_0</math>, с использования оконной обработки (окно Парзена) и без неё.
		+
		+	'''Пример листинга программы''':
		+	<source lang = matlab>
		+	clear all; clc; close all;
		+
		+	use_window = 0;
		+
		+	Fd = 44.1e6; % Hz
		+	N = 64;
		+	T = 1/Fd * N;
		+	t = ( (1:N) - 1 ) / Fd;
		+	f = ( (1:N) - 1 ) / T;
		+	stdS = 8; % RMS of signal
		+	JtoS_dB = 30; %dB, jammer-to-signal
		+	JtoS_amp = 10^(JtoS_dB/20);
		+
		+	S = randn(1, N) * stdS; % Signal
		+
		+	f0 = Fd / 8; % Jammer intermediate frequency
		+	A_jam = stdS * JtoS_amp; % Amplitude of jam
		+	Jam = A_jam * sin(2pif0*t); % Jam
		+
		+	y = S + Jam; % ADC output
		+
		+	% Graphics
		+	figure(1);
		+	subplot(5,2,1)
		+	stem(t*1e6, S)
		+	xlabel('t, \mus');
		+	ylabel('S');
		+
		+	subplot(5,2,2)
		+	stem(f/1e6, abs(fft(S)))
		+	xlabel('f, MHz');
		+	ylabel('fft(S)');
		+
		+	subplot(5,2,3)
		+	stem(t*1e6, Jam)
		+	xlabel('t, \mus');
		+	ylabel('Jam');
		+
		+	subplot(5,2,4)
		+	stem(f/1e6, abs(fft(Jam)))
		+	xlabel('f, MHz');
		+	ylabel('fft(Jam)');
		+
		+	subplot(5,2,5)
		+	stem(t*1e6, y)
		+	xlabel('t, \mus');
		+	ylabel('y');
		+
		+	yf = fft(y);
		+	subplot(5,2,6)
		+	stem(f/1e6, abs(yf))
		+	xlabel('f, MHz');
		+	ylabel('fft(y)');
		+
		+	h = 400;
		+	signal_clear = 0;
		+
		+	if use_window
		+	yw = y .* parzenwin(N)';
		+	else
		+	yw = y;
		+	end
		+	yfw = fft(yw);
		+
		+	subplot(5,2,7)
		+	stem(t*1e6, yw);
		+	xlabel('t, \mus');
		+	ylabel('y after window');
		+
		+	subplot(5,2,8)
		+	stem(f/1e6, abs(yfw));
		+	xlabel('f, MHz');
		+	ylabel('fft(y after window)');
		+
		+	yf_rej = yfw;
		+	for k = 1:N
		+	if abs(yf_rej(k)) > h
		+	yf_rej(k) = 0;
		+	else
		+	yf_rej(k) = yf_rej(k);
		+	end
		+	end
		+	y_rej = ifft(yf_rej);
		+
		+	subplot(5,2,9)
		+	stem(t*1e6, y_rej);
		+	xlabel('t, \mus');
		+	ylabel('y after rejector');
		+
		+	subplot(5,2,10)
		+	stem(f/1e6, abs(fft(y_rej)));
		+	xlabel('f, MHz');
		+	ylabel('fft(y after rejector)');
		+	</source>