Технология цифровой связи

Проектирование цифровой линии

Разработка программы

Основной модуль - содержит основной текст программы, вызовы необходимых функций (вычисления прямого и обратного преобразования Фурье, вычисление интегралов методом прямоугольника), а также вычисление динамики постоянной составляющей и максимального отклонения от постоянной составляющей в сторону меньших значений.

clear;

n=10^ (-5); %коэффициент

dt=0.01. *n; %шаг приращения времени

t1=-3. *n: dt: 0; %Задание интервала времени для отрицательного промежутка

t2=0: dt: 3. *n; %Задание интервала времени для положительного промежутка

t=-3. *n: dt: 3. *n; %Задание общего интервала времени

dw=0.5001. * (1. /n); %Задание шага приращения частоты

w=- (100. * (1. /n)): dw: 100. * (1. /n); %Задание интервала частоты

dwc= (10^7-10^5). /49; wc=10^5: dwc: 10^7; % Задание интервала частоты

n_v_gr= ; %Номера выводимых графиков

n_str_g=size (n_v_gr,2) +1; n_g=1; sch=1; %Вспомогательные переменные для построения графиков

u1t=10. *exp (25. *t1. / (2. *pi. *n)). *sin (35. *t1. /n); %Задание исходного сигнала для промежутка t1

u2t=10. *exp (-25. *t2. / (2. *pi. *n)). *sin (35. *t2. /n); %Задание исходного сигнала для промежутка t2

for i=1: size (t1,2) ut (i) =u1t (i); end %Занесение в единый массив для импульса составляющей промежутка t1

for i=1: size (t2,2) ut (size (t1,2) +i-1) =u2t (i); end %Занесение в единый массив для импульса составляющей промежутка t2

ut_sr=integral (ut,t). / (t (size (t,2)) - t (1));

fw=priam_pr_fur (w,t,ut); %Вызов функции прямого преобразования Фурье

subplot (n_str_g,2,n_g); n_g=n_g+1; plot (t,ut); grid on; %Построение исходного сигнала

subplot (n_str_g,2,n_g); n_g=n_g+2; plot (w,abs (fw)); grid on; %Построение преобразованного сигнала

for i=1: size (wc,2) %Цикл с количеством итераций равным количеству частот

R=relax _otr (w, wc (i)); %Вызов функции задающей функцию отражения

fw_n=fw. *R; %Перемножение спектра и функции отражения

ut_n=obr_pr_fur (t,w,fw_n); %Обратное преобразование Фурье

if n_v_gr (sch) ==i; %Проверка необходимости построения графика

subplot (n_str_g,3,n_g); n_g=n_g+1; plot (w,abs (R)); grid on; %Построение функции отражения

subplot (n_str_g,3,n_g); n_g=n_g+1; plot (w,abs (fw_n)); grid on; %Построение спектра умноженного на функцию отражения

subplot (n_str_g,3,n_g); n_g=n_g+1; plot (t,ut_n); grid on; %Построение импульса после обратного преобразования Фурье

sch=sch+1;_post (i) = integral (ut_n,t). / (t (size (t,2)) - t (1)); % Вычисление постоянной составляющей

ut_otcl_min (i) =min (real (ut_n-ut_sr)); % Вычисление максимального отклонения в сторону меньших значений

end

subplot (n_str_g,1,n_str_g. *2-1); plot (wc,ut_post); grid on; %Построение динамики постоянной составляющей

subplot (n_str_g,1,n_str_g. *2); plot (wc,ut _otcl_min); grid on; %Построение динамики максимального отклонения от постоянной составляющей в сторону меньших значений

Модуль integral. Вычисление интеграла методом прямоугольника. Содержит цикл для последовательного сложения вычисленных значений функции на определённом промежутке.

%Вспомогательная функция вычисления интеграла переданной функции

%методом прямоугольника

function vozvr=integral (y,x)

dx= (x (size (x,2)) - x (1)). / (size (x,2) - 1);=0;i=1: size (x,2)=sum+y (i). *dx;=sum;

Модуль priam_pr_fur. Вычисление прямого преобразования Фурье. Содержит два цикла для вычисления преобразования Фурье при заданной частоте. В функцию передаётся интервал частоты, интервал времени и сама функция, зависящая от времени.

%Функция прямого преобразования Фурье

function vozvr=priam_pr_fur (w,t,ut)

dt= (t (size (t,2)) - t (1)). / (size (t,2) - 1);

for i1=1: size (w,2)=0;i2=1: size (t,2)=sum+ut (i2). *cos (w (i1). *t (i2)) - j. *ut (i2). *sin (w (i1). *t (i2));(i1) =sum. *dt;=fw;

Модуль obr_pr_fur. Вычисление обратного преобразования Фурье. Содержит 2 цикла для вычисления обратного преобразования Фурье. В функцию передаётся интервал частоты, интервал времени и сама функция, зависящая от частоты.

%Функция обратного преобразования Фурье

function vozvr=obr_pr_fur (t,w,fw) Перейти на страницу: 1 2


Другое по теме:

Расчет технико-экономических показателей системы контроля В соответствии с ГОСТ 16504 система контроля - это совокупность средств контроля (СК), ОК и оператора, взаимодействующих по правилам, установленным нормативно-техническими документами (НТД). Получаемая с помощью СК информация об исправности ...