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


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



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

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

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

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

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

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

Екі бөлікте.

  1. импульс енінің, спектра тиімді енін тап-  сіздің нұсқаңыз үшін - 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




Нұсқа

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



Достарыңызбен бөлісу:
1   2   3   4   5




©emirsaba.org 2024
әкімшілігінің қараңыз

    Басты бет