Случайные мысли

[dlinyj]

Сегодня, не смотря на то, что завтра в 7 вставать, таки решил попробовать реализовать коррекцию (нормализацию) кривой распада в прямую. Вывел формулу из периода полураспада, и попробовал нормировать суточный график. Каждая точка среднесуточная.

Было


Стало

На лицо улучшение, но видимо для меня период полураспада немного другой. И графики следует ещё немного корректировать. Хотя на лицо резкий спад до 22 числа, затем "плато", потом снова спад до 7 январа и уже овное плато. Корреляции с космическими объектами не проводил. Но интересно, что другие такие корреляции проводили и веьсма успешно: http://www.chronos.msu.ru/RREPORTS/parkhomov_ritmy/parkhomov_ritmy.htm и вы не поверите! Корреляции ЕСТЬ (хотя быть этого не может!).

Выше приведённый график сделан от отчаинья и нехватки времени на нормальную корректную обработку. Следует делать совершенно иным методом.

Может кто мне скажет, вот существует быстрое преобразование Фурье для функции на плоскости (т.е. графика), а как на счёт преобразования Фурье для поверхности? Беглое гугление не дало ответа на мой вопрос, но думаю это оОооочень интересная штукенция.

=====================================================================================================

Просмотрел все восемь лекций Шноля "От 0 до 80. Симон Шноль" http://rutracker.org/forum/viewtopic.php?t=3747586 - автобиография известнейшего биофизика. Получил грандиозное удовольствие и вынес несколько интереснейших жизненных ценностей для себя.

В первую очередь: говорите со своими родителями! Распрашивайте их об их жизни, о жизни их родителей. Это невероятно интересно! Должна быть преемственность поколений. Родители обязательно уйдут, и мы обязательно уйдём и наши дети. Но нужно, чтобы мы знали и знали наши дети, внуки...

=====================================================================================================

Сегодня специально приехал к отцу и распрашивал его о биографии деда, о том как они познакомились с бабушкой и т.д. Невероятная история, полная невероятных тонкостей. Например, мой дед по бумагам защитил диплом 22 июня 1941. Буквально в течении месяца ушёл добровольцем на фронт. Или невероятные рассказы моей тёти, как они во время эвакуации, когда начиналась бомбёжка поезда, выскакивали из поезда и прятались. И ей, 5-ти летней девочке, очень нравилось сидять и прятаться...

Многое, и многое другое. Не буду здесь рассказывать, хотя без сомнения, стоит многое записать, чтобы не забыть! И попросить батю записать.

Comments

[dlinyj.livejournal.com]

Логично, что существует. Даже n-мерное должно существовать. Другой вопрос: как оно выглядит и как им пользоватся

[eddy-em.livejournal.com]

Ты что, никогда картинки не фильтровал?
Гугли БПФ изображений.

[dlinyj.livejournal.com]

Представляешь, никогда не приходилось. Да и вообще работать с изображениями.

Пока бегло поглядел матчасть и понял, что это как раз то, что нужно. По сути плоскость - это размеры изображения, а высоты - его яркость. За идею поясной поклон!

Единственное, что мне непонятно, это как выглядят гармоники. Например, разложение любой функции на плоскости в сумму синусоид я могу представить. То как выглядит здесь? Ведь, если синусойды (или косинусоиды, как удобнее), идут парралельно какой либо оси, то они никогда не смогут создать сложный рельеф.

Так же, не очень понятен график АЧХ, интуиция подсказывает, что он так же должен быть трёхмерным, но как будет выглядеть опять же не представляю.

И да, быть может подскажите программу работы с графикой, которая может строить такие АЧХ?

[eddy-em.livejournal.com]

> как выглядит здесь?
Все просто: в самом центре у тебя нулевая частота (постоянная составляющая), чем дальше от центра — тем большая пространственная частота данной гармоники. За амплитуду гармоники отвечает абсолютная величина значения в пространстве Фурье; за фазу — фаза этого значения (там величины комплексные).

Почитай Гонсалеса и Вудса. Очень советую.

АЧХ ты не построишь, т.к. этот график четырехмерный, но можешь построить отдельно амплитуду и отдельно фазу.

Для начала советую Octave попробовать, там есть БПФ и куча чего еще, а синтаксис матлабовский — так что примеры из Гонсалеса и Вудса там должны (правда, не все) работать.

[m e (from livejournal.com)]

подсказываю -- это imagemagick

помимо охренненого количества других дел умеет делать бпф

imagemagick.org/Usage/fourier/ (чтобы меня не скринило за ссылки, меня надо добавить в друзья)

ачх (если я правильно понял -- тут речь ведь не идет о линейном процессе?) будет вполне себе трехмерной, и представляться (сюрприз!) такой же двумерной картинкой

гармоники будут произведениями одномерных гармоник, и имеют вид exp(imx)exp(iny) или после разложения cos(mx)cos(ny) -- всего 4 варианта с заменами косинуса на синус

[dlinyj.livejournal.com]

Я не знаю как зафрендить гугловский ак :(. Могу лишь сказать, что ссылку в коментах я получу, тоько будет скрыта. А потом раскрою, или не дает оставить коммент с ссылкой?


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

[dlinyj.livejournal.com]

А, ну и главный вопрос, как эти картинки привести к реальным цифрам. У меня типа сигнал ежесекундно представляет график распределения. Я хочу из него делать поверхность и поглядеть периодические составляющие

[m e (from livejournal.com)]

ты лучше расскажи идею полностью (а то она у тебя размазана видимо по постам)

что я понял -- ты хочешь посчитать корреляцию этих данных с положением солнца; если так -- не ясно, зачем тебе двумерное преобразование фурье

тут вроде проще всего посчитать ∫ f(t) exp(-kt) sin(2пt/сутки)dt и ∫ exp(-kt) f(t)cos(2пt/сутки)dt для нескольких десятков или сотен разумных k (<vanga_mode>и увидеть, что на экспоненту можно наплевать</vanga_mode>); и потом прикинуть, какова вероятность получить столь большие интегралы в случае, если распределение случайно

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

[dlinyj.livejournal.com]

На самом деле идея несколько шире, чем поиск космофизических корреляций в случайных процессах. Так получилось, что последнее явилось побочным исследованием, но тоже весьма интересным.

Скажу сразу, я не занимаюсь антинаучным шаманизмом, и стараюсь подходить к своим опытам со всем научным прагматизмом. Колебания температуры и давления никак не сказывается на радиоактивном распаде. И суточные колебания температуры в комнате весьма незначительны, в пределах 5 градусов.

Корреляции существуют. Я в этом посте приводил ссылку на научную публикацию, доказывающую это http://www.chronos.msu.ru/RREPORTS/parkhomov_ritmy/parkhomov_ritmy.htm . Всё делается на приборе, аналогичном моему. Проф. Шноль пошёл дальше, и получил значительно более интересные результаты http://www.chronos.msu.ru/RREPORTS/shnol_ufn.pdf . Говоря кратко, они строят ежеминутные гистограммы распределения радиоактивного распада (по 60-ти точкам), сглаживают его скользящим средним, а затем сравнивают (вручную) полученные гистограммы. Во избежания предвзятости, гистограммы перемешиваются случайным образом. И они находят похожие или почти одинаковые гистограммы.

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

[m e (from livejournal.com)]

> Колебания температуры и давления никак не сказывается на радиоактивном распаде.

колебания температуры сказываются на *регистрирующем* устройстве, в частности, меняют частоту кварца, из-за чего "секунда" становится другой длинны, и количество зарегистрированных распадов в "секунду" тоже меняется

так же может меняться чувствительность аналоговой части и прочее

вполне возможно, ты это дело как-то предусмотрел, но я, т.к. не читал ранние посты, не в курсе

[m e (from livejournal.com)]

> Я в этом посте приводил ссылку на научную публикацию, доказывающую это http://www.chronos.msu.ru/RREPORTS/parkhomov_ritmy/parkhomov_ritmy.htm

открой эту страницу и найди там слова "термостат"

[m e (from livejournal.com)]

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

поподробнее объясни откуда берется поверхность (т.е. не ясно, откуда вторая координата -- первое это время, а вторая?)

[m e (from livejournal.com)]

> На самом деле идея несколько шире, чем поиск космофизических корреляций в случайных процессах.

и как же выглядит идея полностью?

[dlinyj.livejournal.com]

Есть другие средства связи, кроме комментов в ЖЖ? Быть может гталк?

[m e (from livejournal.com)]

послал свой гталк (он же жаббер) тебе приватным сообщением в жж (чтобы спаммерам не светить)

ты в какое время дня/ночи в гталке?

[dlinyj.livejournal.com]

Я обычно в рабочее время в гталке. Но часто бывает работаю работу :))))

[m e (from livejournal.com)]

для луны данных скорее всего мало, потому что меняя k ты видимо сможешь очень сильно менять ∫ f(t) exp(-kt) sin(2пt/месяц) dt, но если тебе повезет, то этот интеграл окажется необычно большим даже при самом худшем k

вообще тема интересная

з.ы. а с друзьями да, меня немного сложно добавлять, т.к. меня зовут модуль 956459 ext_956459, поэтому линк будет www точка livejournal.com/friends/add.bml?user=ext_956459

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

[dlinyj.livejournal.com]

В друзья добавил

[m e (from livejournal.com)]

а еще я припоминаю парочку интересных вещей:

1. слухи о том, что локально около человека время может идти с другой скоростью; емнип там измеряли небольшую разность в частотах (биения наверно?) генератора рядом со спящим человеком и вдалеке

2. некие высокочастоные колебания электормагнитной природы, но весьма нестандарные -- их не гасил слой воды в 10 метров, они детектировались своей антенной, но при этом детектор "жучков" ничего не детектировал