Выделить частоты, входящие в сигнал.

28.03.2013, 23:01. **Показов** 2798. **Ответов** 1

Студворк — интернет-сервис помощи студентам

Всем доброго времени суток!
У меня есть спектр какого то сигнала и мне нужно выделить из него частоты входящие в исходный сигнал. Вот код:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
fs = 1000;
t = (1:fs)/fs;
%--------------------------------------------------------------------------
A1 = 1;
w1 = 18;
f1 = 0;
A2 = 0.7;
w2 = 30;
f2 = -60;
A3 = 0.45;
w3 = 60;
f3 = 40;
A4 = 0.2;
w4 = 120;
f4 = 100;
S1 = A1*sin(2*pi*w1*t + f1);
S2 = A2*sin(2*pi*w2*t + f2);
S3 = A3*sin(2*pi*w3*t + f3);
S4 = A4*sin(2*pi*w4*t + f4);
S = (S1 + S2 + S3).*S4;
subplot(3, 4, 1);
plot(t, S1);
title('S1');
subplot(3, 4, 2);
plot(t, S2);
title('S2');
subplot(3, 4, 3);
plot(t, S3);
title('S3');
subplot(3, 4, 4);
plot(t, S4);
title('S4');
subplot(3, 3, [4 5 6]);
plot(t, S);
title('S = (S1 + S2 + S3)*S4');
%--------------------------------------------------------------------------
k = 16;
NFFT = 2^nextpow2(fs*k);
s = fft(S, NFFT)/NFFT*10;
w = fs/2*linspace(0, 1, NFFT/2 + 1);
subplot(3, 3, [7 8]);
plot(w, k*2*abs(s(1:NFFT/2 + 1)));
xlim([0 200]);
xlabel('Frequency (Hz)')
ylabel('|Y(f)|');

как теперь автоматически выделить исходные частоты? Если кто знает подскажите, заранее спасибо.

@vital792 · 29.03.2013, 13:08

могу предложить такой способ:
1) Рассчитываете спектральную плотность мощности сигнала. Удобнее использовать не фурье спектр, а более сглаженную спектральную оценку, например метод уэлча(функция pwelch). Можно конечно использовать и fft, но тогда результат придется сгладить нч фильтром, чтобы получить огибающую спектра.
2) ищите перегибы огибающей, соответствующие максимумам.
3) сортируете полученные максимумы по убыванию. Положения первых n из них будут определять искомые частоты(подразумевается, что Вам известно число гармоник в сигнале).

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .

@vital792 2014 / 1286 / 61 Регистрация: 05.06.2010 Сообщений: 2,213
	29.03.2013, 13:08
	могу предложить такой способ: 1) Рассчитываете спектральную плотность мощности сигнала. Удобнее использовать не фурье спектр, а более сглаженную спектральную оценку, например метод уэлча(функция pwelch). Можно конечно использовать и fft, но тогда результат придется сгладить нч фильтром, чтобы получить огибающую спектра. 2) ищите перегибы огибающей, соответствующие максимумам. 3) сортируете полученные максимумы по убыванию. Положения первых n из них будут определять искомые частоты(подразумевается, что Вам известно число гармоник в сигнале). 0