Основной модуль - содержит основной текст программы, вызовы необходимых функций (вычисления прямого и обратного преобразования Фурье, вычисление интегралов методом прямоугольника), а также вычисление динамики постоянной составляющей и максимального отклонения от постоянной составляющей в сторону меньших значений.
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 система контроля - это совокупность средств контроля (СК), ОК и оператора, взаимодействующих по правилам, установленным нормативно-техническими документами (НТД). Получаемая с помощью СК информация об исправности ...