Cигнал почему-то получается вверх ногами...

@Bender310 · Регистрация: 04.05.2013

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

А у меня еще вопросик по работе в матлаб: строю спектр сигнала(делаю вручную, т.к. БПФ делает ошибки на НЧ, вроде все нормально, но если обратно собираю сигнал, то он почему-то получается вверх ногами, это видно на картинке, кр-исходный сигнал, зел-собранный). Может кто сталкивался с этим?

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
pz2=load ('pz2.txt');
r=pz2(:,2);
y=pz2(:,3);
h=detrend(r);
T=0:2*pi*0.008/74.735:2*pi;
T=rot90(T);
Signal=h;
h1=0;
h2=0;
for i=1:1:200
    a=(1/pi)*trapz(T,Signal);
 ak(i)=(1/pi)*trapz(T,Signal.*cos(i*T));
 bk(i)=(1/pi)*trapz(T,Signal.*sin(i*T));
t=0:2*pi*0.008/74.735:2*pi;
h1=h1+ak(i).*cos(i*t);
h2=h2+bk(i).*sin(i*t);
Ep=a+(h1+h2);
end
 figure(32);
plot(t,Ep,'g');
figure(55);
plot(t,h,'r');

@R2D2 · 11.05.2013, 23:59

Сообщение от Bender310

pz2=load ('pz2.txt');

Можно повидать вот то самое pz2.txt?

@lomt · 13.05.2013, 10:59

Сообщение от Bender310

строю спектр сигнала(делаю вручную, т.к. БПФ делает ошибки на НЧ,

и вот это можно пояснить, что подразумевается под ошибками БПФ...

@Bender310 · 14.05.2013, 22:16 **[ТС]**

Re: Суммирование семейства кривых в Матлаб
Сообщение от Bender310
строю спектр сигнала(делаю вручную, т.к. БПФ делает ошибки на НЧ,
и вот это можно пояснить, что подразумевается под ошибками БПФ...

смотрела для синусоиды-частоты не совпали

Добавлено через 2 минуты
к сожалению файл загрузить не могу, слишком большой, да это и неважно, для всех сигналов строит перевернутым. В принципе где поставить минус, чтоб сигнал получился таким же-знаю, но с точки зрения физики объяснить не могу.

@vital792 · 14.05.2013, 23:08

Сообщение от Bender310

смотрела для синусоиды-частоты не совпали

а можно взглянуть на файл с которым проводилось тестирование? Функция fft в матлабе не может ошибаться - проверялось тысячу раз. Она работает как обертка в библиотеке fftw которая проверялась много тысяч раз). Уже не первый раз на нашем форуме сомневаются в правильности работы функции fft:
Подскажите пожалуйста, правильно ли работает преобразование Фурье?
Обратное преобразование Фурье
Поверьте, функция работает правильно. А ваш код как то отдаленно напоминает преобразование Фурье, но...

@Bender310 · 15.05.2013, 00:00 **[ТС]**

это и есть преобразование Фурье. попозже скину код для тестового сигнала. просто для fft я еще не очень представляю, как избавиться от некоторых гамоник и собрать обратно сигнал.когда как в моей написанной программе нужные гармоники можно просто обнулить

Зосима · 15.05.2013, 10:41

vital792, заяц пытается сделать ряд Фурье

что-то вроде этого:
http://www.dpva.info/Guide/Gui... ierSeries/

Bender310, без самого сигнала и чёткого задания мы тебе врядли поможем

Расскажи, что тебе нужно сделать (может это можно в одну строчку реализовать), а файлик упакуй в архив и выложи на форум

Добавлено через 1 минуту
PS: у тебя масштабы графиков по Y разные, что также искажает восприятие графической информации

@vital792 · 15.05.2013, 14:06

Сообщение от Зосима

пытается сделать ряд Фурье

я знаю что такое ряды Фурье

В данном коде "заяц" делает и прямое и обратное преобразование в одном цикле, но есть маленькое но... ВременнОе измерения для прямого преобразования зачем то переворачивается, время течет вспять

Matlab M
1
2
T=0:2*pi*0.008/74.735:2*pi;
T=rot90(T);

Предполагаю, это делалось из за того, что исходный сигнал хранился в виде столбца, но в таком случае надо просто транспонировать)

Matlab M
1
2
T=0:2*pi*0.008/74.735:2*pi;
T=T';

В обратном преобразовании время течет как положено:

Matlab M
1
t=0:2*pi*0.008/74.735:2*pi;

Это версия, конечно нужны данные чтобы ответить точно. А вообще я уверен на 100% что fft поможет)

@Bender310 · 16.05.2013, 21:15 **[ТС]**

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

Зосима · 17.05.2013, 09:51

Bender310, выложи, пожалуйста, архив данными!!!

Новые блоги и статьи Все статьи Все блоги /
Модульная разработка через 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-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .

@Bender310 0 / 0 / 0 Регистрация: 04.05.2013 Сообщений: 7
	14.05.2013, 22:16 [ТС]
	Re: Суммирование семейства кривых в Матлаб Сообщение от Bender310 строю спектр сигнала(делаю вручную, т.к. БПФ делает ошибки на НЧ, и вот это можно пояснить, что подразумевается под ошибками БПФ... смотрела для синусоиды-частоты не совпали Добавлено через 2 минуты к сожалению файл загрузить не могу, слишком большой, да это и неважно, для всех сигналов строит перевернутым. В принципе где поставить минус, чтоб сигнал получился таким же-знаю, но с точки зрения физики объяснить не могу. 0

@Bender310 0 / 0 / 0 Регистрация: 04.05.2013 Сообщений: 7
	15.05.2013, 00:00 [ТС]
	это и есть преобразование Фурье. попозже скину код для тестового сигнала. просто для fft я еще не очень представляю, как избавиться от некоторых гамоник и собрать обратно сигнал.когда как в моей написанной программе нужные гармоники можно просто обнулить 0

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	15.05.2013, 10:41
	vital792, заяц пытается сделать ряд Фурье что-то вроде этого: http://www.dpva.info/Guide/Gui... ierSeries/ Bender310, без самого сигнала и чёткого задания мы тебе врядли поможем Расскажи, что тебе нужно сделать (может это можно в одну строчку реализовать), а файлик упакуй в архив и выложи на форум Добавлено через 1 минуту PS: у тебя масштабы графиков по Y разные, что также искажает восприятие графической информации 0

@Bender310 0 / 0 / 0 Регистрация: 04.05.2013 Сообщений: 7
	16.05.2013, 21:15 [ТС]
	задание: есть массив данных, предположительно в нем есть колебательная составляющая-24 часа, ее нужно исключить. для этого я считаю коэффициенты Фурье, строю спектр и вычитаю соответствующую гармонику, далее собираю сигнал обратно. В принципе уже вроде поняла, что ошибка идет не в цикле, а раньше. 0

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	17.05.2013, 09:51
	Bender310, выложи, пожалуйста, архив данными!!! 0