75611

РАЗЛОЖЕНИЕ ФУНКЦИЙ В КОМПЛЕКСНЫЙ РЯД ФУРЬЕ

Лекция

Коммуникация, связь, радиоэлектроника и цифровые приборы

Это и есть разложение в комплексный ряд Фурье. Коэффициенты Сk называются комплексными коэффициентами Фурье и, подобно действительным коэффициентам Фурье, вычисляются как скалярные произведения

Русский

2015-01-15

60.5 KB

1 чел.

ОС. Лекция 5.

РАЗЛОЖЕНИЕ ФУНКЦИЙ В КОМПЛЕКСНЫЙ РЯД ФУРЬЕ

Система функций на отрезке  образует ортонормированную систему функций. Значит, произвольная функция f(t) может быть представлена по этой системе следующим образом:

Это и есть разложение в комплексный ряд Фурье. Коэффициенты Сk называются комплексными коэффициентами Фурье и, подобно действительным коэффициентам Фурье, вычисляются как скалярные произведения f(t) и ejkt:

Если период функции не равен , а, например, равен Т, то получим следующее общее выражение для комплексных коэффициентов:

Коэффициенты Фурье являются комплексными числами, но f(t) является действительной функцией, а значит правая часть  последнего выражения должна быть действительной. Так оно и есть на самом деле, потому что коэффициенты Ck и C-k являются сопряженными. Если взяты целые положительные значения k, то функцию f(t) можно записать в виде:

Но, учитывая то, что Ck и C-k являются сопряженными, получим:

Ниже приведена программа разложения дискретизированной функции y=x2 содержащего N значений в комплексный ряд Фурье на интервале [-T,T] с М членами разложения, M<N, и  последующего восстановления. Для сравнения приведены результаты, полученные с помощью стандартных функций fft (БПФ) и ifft (ОБПФ) MATLAB.

%Разложение функции t^2 в комплексный ряд Фурье

%в дискретизированном виде на интервале [0,T]

%Восстановление функции производится по формуле

% fв(i)=y(i)=sum(ck*exp(j*2*pi*0*i/N)), k=[1,M],

% i= [0,N-1]

%Чем больше М, тем точнее восстановление

 

T=4;%Значение T (произвольное)

N=128;%количество значений функции на интервале [0,T](произвольное)

M=8;

for i=1:N

   f(i)=2*T*(i-1)^2/N;  %исходная дискретизированная функция      

end

for k=1:M

  C(k)=0;

for i=1:N

   C(k)=C(k)+f(i)*exp(-j*2*pi*k*(i-1)/N);   

end

C(k)=C(k)*(1/N);

end

 

for i=1:N

   y(i)=0;

   f3(i)=0;

   for k=1:M    

   y(i)=y(i)+C(k)*exp(j*2*pi*k*(i-1)/N);   

   end

end

for k=1:M

   koef(k)=C(k);

   koef2(k)= exp(j*2*pi*k*(i-1)/N);

end    

   f3(1)=C(1);

   f3(2)=C(2);

   f3(3)=C(3);

   f3(4)=C(4);

   f3(5)=C(5);

   f3(5)=C(6);

 

   

   i=1:N;

   f3=f3/max(f3);

 

%2. ФУНКЦИОНАЛЬНОЕ ПРЕОБРАЗОВАНИЕ (БПФ)

i=1:N;

bpfy=fft(f,N);%БПФ

bpf=(bpfy.*conj(bpfy));%БПФ

bpf=bpf/max(bpf);%%%%%%%

f2=ifft(bpfy);

figure

hold on;

plot(i(1:10),bpf(1:10));

plot(i(1:10),f3(1:10));

axis tight;

title('Frequency domain')

xlabel('Количество периодов')

 

hold off;

figure

plot(i,f2);

axis tight;

%нахождение макс. знач. функции БПФ для массива Y

C=max(bpf);

for i=1:N %поиск количества периодов, соответствующих максимуму БПФ

   if (bpf(i)==C)          

       kpbpf=(i-1);         

       break

   end

end

kp_bpf=kpbpf

i=1:N;

figure

%plot(i(1:N),y(1:N));%отображение графика y линией красного цвета

plot(i,y,'r-');

axis tight;

title('Time domain. Восст. Y  и f2')

xlabel('Points number')

%hold on;

figure

plot(i,f2);

axis tight;

%hold off;

pause;

close all;


 

А также другие работы, которые могут Вас заинтересовать

41607. Розрахунок площ адміністративних та побутових приміщень 80 KB
  Визначаємо очікувану кількість чоловіків і жінок друкарні, використавши для цього що характерне співвідношення для поліграфічних підприємств: чоловіків - 45%, жінок - 55%. Тоді очікувана кількість чоловіків і жінок відповідно становитиме
41608. МЕРЫ ПО УЛУЧШЕНИЮ ФИНАНСОВОГО СОСТОЯНИЯ ПРЕДПРИЯТИЯ ОАО «ММК-МЕТИЗ» 362.58 KB
  Рассмотреть теоретические аспекты анализа финансового состояния предприятия; Дать общую характеристику предприятия и проанализировать его финансовое состояние; Выявить проблемы финансового состояния предприятия; Разработать рекомендации по улучшению финансового состояния предприятия
41609. Решение системы линейных уравнений методом простых итераций и методом Чебышева 45.92 KB
  Требуется написать программу реализующая 2 метода решение системы линейных уравнений: 1методом простых итераций; 2методом Чебышева. Теория: 1Метод простых итераций Требуется решить систему уравнений 1 где – симметрическая положительно определенная матрица. Метод простых итераций имеет вид...
41610. Первинні засоби пожежогасіння. Вибір типу та визначення необхідної кількості первинних засобів пожежогасіння 309.89 KB
  Головним критерієм вибору виду вогнегасників є величина можливого осередку пожежі. Визначаємо рекомендовані типи вогнегасників. Користуючись рекомендаціями таблиці Д5 щодо порошкових вогнегасників визначаємо що для захисту промислових обєктів рекомендованими є такі типи переносних порошкових вогнегасників: ВП5 ВП6 ВП9 ВП12 записуємо в табл. Визначаємо кількість вогнегасників.
41611. Диференціальні рівняння в частинних похідних 44.88 KB
  resizeN1; forint i = 0; i N1; i u[0][i] = conditionih; forint j = 1; j NT; j { file T: htj endl; forint i = 1; i N; i f[i] = u[j1][i1] 2u[j1][i] u[j1][i1] 2hh u[j1][i]1 ht Q 2; l[2] = c; b[2] = f[1] c; u[j][0] = 0; u[j][N] = 0; forint i = 2; i N; i { l[i1] = c l[i]; b[i1] = f[i] b[i] c l[i]; } forint i = N1; i 0; i u[j][i] = l[i1]u[j][i1] b[i1]; int emx = 0; for int i = 0; i N; i { file x: ih ...
41612. Створення нової бази даних в середовищі MySQL 138.93 KB
  Створити нову базу даних та заповнити її даними. Короткі теоретичні відомості: Основи роботи з phpMydmin При установці Denwer також встановлюється на комп'ютер phpMydmin за допомогою якого можна керувати базою даних MySQL через вебінтерфейс. У цьому полі латинськими буквами записується назва бази даних наприклад exmple і натискається створити.
41613. Приближенное вычисление интеграла методом Симпсона и методом Гаусса 92.3 KB
  Требуется вычислить интеграл: Требуется использовать: метод Симпсона метод Гаусса Теория: 1 Метод Симпсона Для приближённого вычисления интеграла чаще всего подынтегральную функцию заменяют близкой ей вспомогательной функцией интеграла от которой вычисляется аналитически. В частности если при вычислении подынтегральную функцию заменить интерполяционным многочленом второй степени построенным по значениям функции в трёх...
41614. Состояние дерматовенерологических больных в Винницкой области 354.5 KB
  Проблема совершенствования лекарственного обеспечения населения регионов Украины остается актуальной. Особое значение в её решении имеет региональный подход к изучению фармацевтического рынка, его насыщенности и рациональному использованию лекарственных средств. С этой целью широко используются метод фармакоэкономического анализа
41615. Решение уравнения f(x)=0 методами простых итераций и Ньютона 134.65 KB
  Если же то вычисления заканчивают и за приближённое значение корня принимают величину . Абсциссы вершин этой ломанной представляют собой последовательные приближения корня . Из рисунков видно что если на отрезке то последовательные приближения колеблются около корня если же производная положительна то последовательные приближения сходятся к корню монотонно. Если через точку с координатами провести касательную то абсцисса точки пересечения этой касательной с осью и есть очередное приближение корня уравнения .