Зертханалық жұмыс №10 (4 сағат) Фурье түрленуі. Сигналдардың спектрлік талдауы. Жұмыс мақсаты

T - Фурье түрлендіруінде интегралдау интервалы. dt

жүктеу/скачать 119,62 Kb.

бет	4/5
Дата	07.01.2022
өлшемі	119,62 Kb.
	#20210
түрі	Құрамы

1 2 3 4 5

Байланысты:
Зертханалық жұмыс №10

T - Фурье түрлендіруінде интегралдау интервалы.

dt - дискретизация қадамы.

Ұсынысқа мағынаның талғамының параметр тапсырынды нұсқа үшін биік көрсетілген.
6.Бағдарламаны айтылмыш нұсқаның параметрлерімен қосыңыз. Нәтижені оқытушыға көргізіңіз. Конспекте кестелер салу. Спектр тәртібінің мінезін түсіндіру.Жұптылыққа және тақ екеніне мән беріңіз.
Мысал ретінде берілген программа кодтары
Мысал ретінде тікбұрышы импульс программа кодтары

% тікбұрышты импульс- аналитикалық спектр формулалары
a = 2; % импульс ұзақтығы (с)

T = 20; % жиілігінің графиктағы интервалы (c)

dt = 0.1; %осьтің мәніндегі интервалы (с)
t = -T/2:dt:T/2; % уақыттың дискртетті мәндері (c)
figure;

plot(t, s);

axis([-T/2 T/2 0 1.5]);

xlabel('t (c)');

title('Signal s(t)');

grid on;
fb = 2; % (Hz) жиілігінің интервалы

df = 0.01; % (Hz) графикадағы қадамы
f = -fb:df:fb; % (Hz) жиілігінің интервалы

S = sin(pi*a*f)./f/pi;

Ampl = abs(S); % АЧХ есептейміз

Re = real(S); % спектр дұрыс бөлігін есептейміз

Im = imag(S); % спектр бөлігін есептейміз

MaxSpectr = max(Ampl); % спектр графикасының максимумы

figure;

plot(f, Ampl);

axis([-fb fb -1.2*MaxSpectr 1.2*MaxSpectr]);

xlabel('f (Hz)');

title('Spectrum A(f)');

grid on;
figure;

plot(f, Re);

axis([-fb fb -1.2*MaxSpectr 1.2*MaxSpectr]);

xlabel('f (Hz)');

title('Spectrum Re(S(f))');

grid on;
figure;

plot(f, Im);

axis([-fb fb -1.2*MaxSpectr 1.2*MaxSpectr]);

xlabel('f (Hz)');

title('Spectrum Im(S(f))');

grid on;

% ДПФ тікбұрышты импульсы
a = 2; % импульс ұзақтығы (с)

fb = 4; % (Hz) жиілігінің графиктағы интервалы

T = 20; % Фурье түрлендіруіндегі интегралдау интервалы

dt = 0.1; % дискретизация қадамы (с)

t = -T/2+dt:dt:T/2; %уақыттың дискртетті мәндері

N1 = length(t); %тақ болу үшін тексеру

if mod(N1,2) == 1

N1 = N1 + 1;

T = dt*N1;

t = -T/2+dt:dt:T/2;

end
s = rectpuls(t/a); %сигналдар қорытындысы
figure;

plot(t, s);

axis([-T/2 T/2 0 1.5]);

xlabel('t (c)');

title('Signal s(t)');

grid on;
x=s; % ДПФ сигнал жібереміз

N = length(x); % Фурье түрлендіруінің мүшелер саны

X=fft(x); % Фурье түрлендіруінің тікелей дисктретациясы

k=1:N;

ex=exp(i*2*pi*(k-1)*(N/2-1)/N); %экспоненттер есептейміз

Sp = dt*ex.*X; % тәртіптелмеген спектр есептейміз
S1=Sp(N/2+2:N);

S2=Sp(1:N/2+1);

S = [S1 S2]; % тәртіпті спектр
Ampl = abs(S); %АЧХ есептейміз

Re = real(S); % спектр дұрыс бөлігін есептейміз

Im = imag(S); % спектр бөлігін есептейміз
df=1/T; % жиілік дискретизациясының интервалы

f=(-N/2+1:N/2)*df; % дискретты жиілік жасаймыз

MaxSpectr = max(Ampl); % спектр графикасының максимумы
figure;

plot(f, Ampl);

axis([-fb fb -1.2*MaxSpectr 1.2*MaxSpectr]);

xlabel('f (Hz)');

title('Spectrum A(f)');

grid on;
figure;

plot(f, Re);

axis([-fb fb -1.2*MaxSpectr 1.2*MaxSpectr]);

xlabel('f (Hz)');

title('Spectrum Re(S(f))');

grid on;
figure;

plot(f, Im);

axis([-fb fb -1.2*MaxSpectr 1.2*MaxSpectr]);

xlabel('f (Hz)');

title('Spectrum Im(S(f))');

grid on

Бағдарламаны айтылмыш нұсқаның параметрлерімен орындаңыз. Нәтижені оқытушыға көргізіңіз. Кестелерді конспектілеңіз. Спектр тәртібінің мінезін қараңыз. .Жұптылыққа көзге іл немесе атқаратын қызметімнің емес жұптылығының жасаушы белгі. Спектрлік мінезін және спектр сызығының әртүрлі биіктігін анықта. Найквист жиілігін сендердің нұсқаң үшін табыңыз.

Есеп беру талаптары.

Есеп беру арада алдымда асырауға керекке.

1) жазып ал- тапсырманың параметрлерінің сіздіңнұсқаңыз үшін - жұптылық, f0, a .

2) зарисовать және екі про-грамм ара нәтиже біріншінің жұмысының алған кестелерді - s (t), Re (S (f)), Im (S (f)),АЧХ

3) бағдарламаның жұмысы үшін таңдап алған параметрлерді - T, dt, fb, df жазып ал 4) Найквиста таңдап ал-жиілігін – f Nyquist жазып ал.

Екі бөлікте.

импульс енінің, спектра тиімді енін тап- сіздің нұсқаңыз үшін - t, f, B .

Үш бөлікте.

1) тапсырманың параметрлерін сіздің нұсқаңыз үшін - f1, f2, f3, a1, a2, a3 жазып ал.
2) кестелерді сызыңыз .

3) бағдарламаның жұмысы үшін таңдап алғанпараметрлерді - T, dt, fb, Tb жазып ал.
4) Найквист таңдап ал- жиілігін - fNyquist жазып ал.

Бақылау сұрақтар.

1. Фурье өзгерісінің формуласы.
2. АЧХ , ФЧХ ұғымы.
3. Белгінің базасының ұғымы.
4. Жұптылық және ақиқаттық белгінің спектрлік емес жұптылығы.
5. Астасудың белгісінің спектра дельталық мінезі.
6. Фурье дискретті өзгерісінің формулы.
7. Найквист жиілігінің мағынасы .
8. Белгінің спектрін есептеу алгоритмі мен ДПФ көмегі.
Тапсырманың нұсқаларының зертхананың жұмысына.

Тапсырма 1, 2

Нұсқа

Функция жұптылығы

жүктеу/скачать 119,62 Kb.

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

1 2 3 4 5