реклама на сайте
подробности

 
 
 
Reply to this topicStart new topic
> Статистический анализ в EXCEL
Herz
сообщение Jun 22 2009, 21:03
Сообщение #1


Гуру
******

Группа: Модераторы
Сообщений: 10 983
Регистрация: 23-11-05
Пользователь №: 11 287



Хочу построить в Excel диаграмму-спектр некоторого сигнала. К массиву дискретных отсчётов длиной, кратной степени двойки, применил преобразование Фурье. Получил колонку данных в виде сумм действительных и мнимых слагаемых. Кто знает, что с ними делать дальше? В конце концов, хотел бы смоделировать спектральную фильтрацию: "вырезать" "ненужные" частоты и выполнить обратное преобразование, дабы получить результат в виде, близком к исходному. Но хелп весьма скудный, а книги толковой нету...
Go to the top of the page
 
+Quote Post
EUrry
сообщение Jun 23 2009, 03:52
Сообщение #2


Гуру
******

Группа: Свой
Сообщений: 3 218
Регистрация: 14-11-06
Из: Н. Новгород
Пользователь №: 22 312



Есть книга "Статистика в Excel - Макарова, Трофимец, 2002 г.". Может поможет. PDF, DJVU


--------------------
Все не могут только сеять разумное, доброе, вечное: кому-то надо и пахать!
Природа не терпит пустоты: там, где люди не знают правды, они заполняют пробелы домыслом. © Бернард Шоу
Go to the top of the page
 
+Quote Post
MrYuran
сообщение Jun 23 2009, 05:36
Сообщение #3


Беспросветный оптимист
******

Группа: Свой
Сообщений: 4 640
Регистрация: 26-12-07
Из: Н.Новгород
Пользователь №: 33 646



Так это лучше в маткаде (лабе) делать, а не в экселе извращаться.


--------------------
Программирование делится на системное и бессистемное. ©Моё :)
— а для кого-то БГ — это Bill Gilbert =)
Go to the top of the page
 
+Quote Post
тау
сообщение Jun 23 2009, 08:01
Сообщение #4


.
******

Группа: Участник
Сообщений: 2 424
Регистрация: 25-12-08
Пользователь №: 42 757



Цитата(Herz @ Jun 23 2009, 01:03) *
Кто знает, что с ними делать дальше? В конце концов, хотел бы смоделировать спектральную фильтрацию: "вырезать" "ненужные" частоты и выполнить обратное преобразование, дабы получить результат в виде, близком к исходному.


потренируйтесь следующим образом:
составьте колонку из чисел 1, 2, 3, ....16 (N=16 ) начиная с первой строки столбца А
в соседней колонке наберите формулу =sin(B1 * 2* ПИ()/16) + sin(B1 * 3* 2* ПИ()/16) и размножьте результат на все 16 строк
это первая и третья гармоника (сумма двух синусов) на периоде равном 16 строкам представленная в виде отсчетов .
Примените прямое преобразование Фурье для колонки B . В колонке с результатом увидите спектр. Значения в спектре для 2-й и 16 строки одинаковы кроме знака мнимой части. То же самое и для 4-й и 14 строчки . Таково свойство симметрии ДПФ.
Делаете "фильтрацию" rolleyes.gif , обнуляете в спектре 4 -ю и 14 строчки в результатах прямого преобразования . Запускаете обратное преобразование (галочку "инверсия" в нижнем левом углу формы поставьте) . Получаете отсчеты отфильтрованного синуса первой гармоники .

Сообщение отредактировал тау - Jun 23 2009, 08:03
Go to the top of the page
 
+Quote Post
Herz
сообщение Jun 24 2009, 10:21
Сообщение #5


Гуру
******

Группа: Модераторы
Сообщений: 10 983
Регистрация: 23-11-05
Пользователь №: 11 287



Цитата(EUrry @ Jun 23 2009, 05:52) *
Есть книга "Статистика в Excel - Макарова, Трофимец, 2002 г.". Может поможет. PDF, DJVU

Спасибо, скачал. Надеюсь, поможет.

Цитата(тау @ Jun 23 2009, 10:01) *
в соседней колонке наберите формулу =sin(B1 * 2* ПИ()/16) + sin(B1 * 3* 2* ПИ()/16) и размножьте результат на все 16 строк

Тут, видимо, вместо В имелось в виду А? Эта часть вроде затруднений не вызывает. Только как увидеть этот спектр в графическом виде?
Цитата
Делаете "фильтрацию" rolleyes.gif , обнуляете в спектре 4 -ю и 14 строчки в результатах прямого преобразования . Запускаете обратное преобразование (галочку "инверсия" в нижнем левом углу формы поставьте) . Получаете отсчеты отфильтрованного синуса первой гармоники .

А вот с обратным что-то не получается... Здесь, по-моему, уже должен был получить вещественные числа.
[attachment=33933:Clipboard01.jpg]
Go to the top of the page
 
+Quote Post
тау
сообщение Jun 24 2009, 12:04
Сообщение #6


.
******

Группа: Участник
Сообщений: 2 424
Регистрация: 25-12-08
Пользователь №: 42 757



Цитата(Herz @ Jun 24 2009, 14:21) *
Тут, видимо, вместо В имелось в виду А? Эта часть вроде затруднений не вызывает.

Да, конечно, сылка на те ячейки где номера 1...16, если номера в А столбце , то и соответсвенно Ахх

файл экселя попытаюсь вложить




на левом графике типа сумма двух синусов постороенная из отсчетов.
На правом - график спектра , первую строчку в диапазон брать не надо при выводе диаграммы спектра - там находится постоянная составляющая сигнала.
Эскизы прикрепленных изображений
Прикрепленное изображение
 

Прикрепленные файлы
Прикрепленный файл  FFT_XLS.rar ( 2.74 килобайт ) Кол-во скачиваний: 60
 
Go to the top of the page
 
+Quote Post
тау
сообщение Jun 24 2009, 20:00
Сообщение #7


.
******

Группа: Участник
Сообщений: 2 424
Регистрация: 25-12-08
Пользователь №: 42 757



Цитата(Herz @ Jun 24 2009, 14:21) *
А вот с обратным что-то не получается... Здесь, по-моему, уже должен был получить вещественные числа.


Вы убрали (обнулили ) результаты FFT в 4 и 16 строке , а надо было в 4 и 14-й. Поэтому случилось нарушение условия симметрии действительных данных в отсчетах и они после обратного преобразования у Вас получились комплексными. Если данные на входе прямого преобразования Фурье являются чиста "реальными" то в спектре выполняется условие X(k)=X(N-k) где N - число отсчетов , а к номер гармоники от 0 до N-1.
Можно брать во внимание только N/2 гармоник спектра , если входные отсчеты не комплексные.
Go to the top of the page
 
+Quote Post
Herz
сообщение Jun 25 2009, 09:32
Сообщение #8


Гуру
******

Группа: Модераторы
Сообщений: 10 983
Регистрация: 23-11-05
Пользователь №: 11 287



Теперь получилось. Спасибо за помощь. Кстати, в Вашем примере результаты Фурье-преобразований сохранены как текст, кроме того, десятичный разделитель - запятая. Как это могло получиться? У меня на такой тип данных Excel ругается. Вы не применяли дополнительно никакого форматирования?
А для построения диаграммы, выходит, не хватало только функции МНИМ.АВS. rolleyes.gif

P.S. Excel всё-таки, довольно мощный инструмент. Не имею достаточного опыта для всестороннего сравнения, однако он вполне достоин конкурировать с тем же Матлабом в сборе, обработке и визуализации больших массивов данных, ИМХО. Извращение - это когда не умеешь пользоваться. biggrin.gif
Go to the top of the page
 
+Quote Post
тау
сообщение Jun 25 2009, 10:14
Сообщение #9


.
******

Группа: Участник
Сообщений: 2 424
Регистрация: 25-12-08
Пользователь №: 42 757



Цитата(Herz @ Jun 25 2009, 13:32) *
Кстати, в Вашем примере результаты Фурье-преобразований сохранены как текст, кроме того, десятичный разделитель - запятая. Как это могло получиться? У меня на такой тип данных Excel ругается. Вы не применяли дополнительно никакого форматирования?

Я делал на соседнем дочкином компе, там надстройка для анализа стоит, в моем нету (жлобство при установке офиса). Формат ячеек не как "текст" а скорее всего "общий" типа универсальный. Поэтому при вводе чиста цифр может меняться на числовой формат. Я внимания этому не придавал.
Насчет запятых -разделителей, ну это как повезет, меняется в настройках экселя где-то там в параметрах. Кажется по умолчанию для русских именно запятые.
Дополнительно ничего не форматировал, все было как по умолчанию. С запятыми - точками - часто распространенное непонимание между разными экселями на разных машинах.
Насчет мощности экселя - так там еще и VBA есть, можно вообще все автоматизировать.
Удачи!
Go to the top of the page
 
+Quote Post

Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 30th June 2025 - 22:58
Рейтинг@Mail.ru


Страница сгенерированна за 0.01979 секунд с 7
ELECTRONIX ©2004-2016