Дискретизация аналоговых сигналов

ГУАП

КАФЕДРА №22

ОТЧЕТ
ЗАЩИЩЕН С ОЦЕНКОЙ

ПРЕПОДАВАТЕЛЬ

РАБОТУ ВЫПОЛНИЛ

Санкт-Петербург 20__

1. Цель работы: изучение основных свойств спектральной функции дискретных сигналов и метода ее расчета на ЦВМ.
2. Выполнение лабораторной работы

Задание №11

[pic 1]

Рисунок 1 – Функция сигнала S(t) и ее параметры

Листинг программы:

function [s,s0,S,S0]=mysignal(t,f)

% Функция расчета:

% s - дискретного сигнала

% s0 - аналогового сигнала

% S - спектра дискретного сигнала

% S0 - спектра аналогового сигнала

% Аналоговый сигнал - прямоугольный видеоимпульс.

% Входные параметры функции:

% t - массив значений времени, t=(0:dt:T-dt)',

% f - массив значений частоты, f=(-0.5/dt:df:0.5/dt-df)',

%Глобальные переменные

global T N Ti dt df tau w

% T  - интервал наблюдения сигнала

% N  - количество отсчетов сигнала на интервале наблюдения

% Ti - длительность импульса

% dt - период дискретизации

% df - шаг расчета спектра по оси частот, df=1/(N*dt) 

K=64; %Количество подынтервалов, на которое разбивается период дискретизации при численном расчете S0

M=N*K; %Общее Количество подынтервалов при численном расчете S0 %Количество точек при микродискретизации

Dt=T/M; %Длительность подынтервала %% Шаг микродискретизации

tnew=(0:M-1)'*Dt; %Вектор отсчетов времени при численном расчете S0  % вектор времени при микродискретизации

Df=1/(M*Dt); %Соответствующий шаг по частоте при численном расчете S0

fnew=(-0.5*M:0.5*M-1)'*Df; %Вектор отсчетов частоты при численном расчете S0

T=2;

tau=0.5*T;

w=2*pi*10;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Далее следует расчет спектров аналогового и дискретного сигналов в случае

% прямоугольного видеоимпульса

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

s0=zeros(M,1); %Обнуляем массив значений аналогового сигнала

Ind=find((0<=tnew)&(tnew<=tau)); %Определение номеров отсчетов, удовлетворяющих условию 0<=t<=Ti

s0(Ind,1)=(1-(abs((2*tnew(Ind)-tau)/tau))).*exp(1i*w*tnew(Ind));

s=s0(1:K:end);%Вычисление дискретного сигнала