Файл: Практическая работа 1 по дисциплине Передача информации в системах реального времени Спектральное представление периодического сигнала.docx

МИНОБРНАУКИ РОССИИ
Федеральное государственное бюджетное образовательное учреждение высшего образования «МИРЭА – Российский технологический университет» РТУ МИРЭА

Институт Искусственного Интеллекта

Кафедра Промышленной информатики

Практическая работа № 1

по дисциплине

«Передача информации в системах реального времени»

«Спектральное представление периодического сигнала»

Работу выполнил:

Студент группы

______________подпись

Работу проверил:

к.т.н., доцент Акимов Дмитрий Александрович ______________подпись

Москва 2023

Содержание

Пример 1

Результатом выполнения кода (Приложение А) является график (Рисунок 1):



Рисунок 1

Пример 2

Результатом выполнения кода (Приложение Б) является график (Рисунок 2):



Рисунок 2

Пример 3

Результатом выполнения кода (Приложение В) является график (Рисунок 3):



Рисунок 3

Пример 4

Результатом выполнения кода (Приложение Г) является график (Рисунок 4):



Рисунок 4

Задание

Для варианта 3 имеем:

---

a=2, b=0.5, c=3,d=-1,l=1

Был написан код (Приложение Д) результатом выполнения которого является график (Рисунок 5):



Рисунок 5

Приложение А

Листинг исходного кода

N=256;

T=1;

dt = T/(N);

Nyq = N/(2*T);

df = 1/T;

f=20;

w=2*pi*f;

t=linspace(0,T,N);

s=2*sin(w*t./T);

F=fft(s)/N;

subplot(211);

plot(t,s,'b-');

subplot(212);

plot(1:N,abs(F(1:N)),'r-');

Приложение Б

Листинг исходного кода

n=16;

N=64;

t=linspace(0,n,N);

f=exp(-(t-n/2).^2);

F=fft(f);

subplot(411);

plot(t,f);

subplot(412);

plot(1:N,real(F));

subplot(413);

plot(1:N, imag(F));

subplot(414);

plot(1:N, abs(F));

Приложение В

Листинг исходного кода

N=256;

T=1;

dt=T/(N);

Nyq=N/(2*T);

df=1/T;

f=20;

w=2*pi*f;

t=linspace(0,T,N);

s=2*(cos(w*t./T)+i*sin(w*t./T));

F=fft(s)/N;

subplot(211);

plot(t,real(s),'b-');

F1=F(1:N/2+1);

F2=F(N/2+1:N);

F=[F2,F1];

nu=-Nyq+df*(0:N);

subplot(212);

plot(nu(1:N),abs(F(1:N)),'r-');

Приложение Г

Листинг исходного кода

N=256;

T=1;

dt=T/(N);

Nyq = N/(2*T);

df=1/T;

f=50;

w=-2*pi*f;

t=linspace(0,T,N);

for x=1:1:N;

s(1,x)=0;

if x>N/2 -8 & x < N/2+8;

s(1,x)=2;

end

end;

F=fft(s)/N;

subplot(411);

plot(t,s,'b-');

F1=F(1:N/2+1);

F2=F(N/2+1:N);

F=[F2,F1];

nu = -Nyq+df*(0:N);

subplot(412);

plot(nu(1:N),real(F(1:N)));

subplot(413);

plot(nu(1:N),imag(F(1:N)));

subplot(414);

plot(nu(1:N),abs(F(1:N)),'r-');

Приложение Д

Листинг исходного кода

N = 256;

T = 2;

dt = T/(N);

Nyq = N/(2*T);

df = 1/T;

f = 50;

w = -2*pi*f;

t = linspace(0, T, N);

for x =1:1:N;

if x>=N/2-1 & x
s(1,x)=2*x+0.5;

end

if x>=N/2 & x
s(1,x) = 3*x-1;

end

end;

F = fft(s,N);

subplot(411);

plot(t,s,'b-');

F1=F(1:N/2+1);

F2=F(N/2+1:N);

F=[F2, F1];

nu = -Nyq+df*(0:N);

subplot(412);

plot(nu(1:N), real(F(1:N)));

subplot(413);

plot(nu(1:N), imag(F(1:N)));

subplot(414);

plot(nu(1:N), abd(F(1:N)), 'r-');

---