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


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

Институт Искусственного ИнтеллектаКафедра Промышленной информатикиПрактическая работа № 1по дисциплине«Передача информации в системах реального времени»«Спектральное представление периодического сигнала»Работу выполнил:Студент группы ______________подпись Работу проверил:к.т.н., доцент Акимов Дмитрий Александрович ______________подписьМосква 2023

Содержание

РТУ МИРЭА 1

Пример 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;endend;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 & xs(1,x)=2*x+0.5;endif x>=N/2 & xs(1,x) = 3*x-1;endend;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-');