vinipu
Спасибо большое за разъяснения!

385-NJ=NWT
Жесть - интересно, почему они не могут настроить манипуляции на NWT, а передают NJ. По свертке как раз хорошо видно, что между всеми точками/тире одинаковая пауза.

В ndblist прошла информация по Unid 374 JRM:
Is NDB "JAHROM" on 374 kHz
Coordinates 28°35'07.40"N -- 053°35'02.70"E
GRID LL68TO

A4389/17 NOTAMN
Q) OIIX/QNDCF/IV/BO/AE/000/999/
A) OISJ B) 1712231039 C) PERM
E) REF AIP PAGE AD 2-5 OISJ ITEM 2.19,
NDB/DME JRM CHANGES AS FLW:
NDB/DME JRM FREQ 374 KHZ CH 85X (113.800 MHZ),
AMEND AIP ACCORDINGLY.
CREATED: 23 Dec 2017 10:49:00
SOURCE: OIIIYNYX

Спасибо Alain Désaubry!

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




Ошибки есть у обеих функций. Пирсон чаще ошибается в бОльшую сторону (это уже после доточки её, поначалу почти всегда выбирала 5-10 кратные интервалы). Ошибка наклона у моей с виду чаще - но Пирсон тоже не редко этим страдает. Но

— Хм-м-м, — сказали суровые сибирские мужики и положили рельс.
— КРЯК! — сказала японская бензопила и сломалась.
— Ага-а-а! — сказали суровые сибирские мужики

в плане чувствительности моя функция пока лучше Пирсона. Пример:





Итог отрицательный. Я очень надеялся, что стандартная функция корреляции (пусть и считает дольше) будет давать красивые прямые правильные интервалы. Но нет. Во всяком случае в варианте, когда я считаю матрицами по 4 штуки (просто парами дает вообще плохой результат). Можно увеличивать размер матрицы - но тогда считаться это будет совсем долго (сейчас и так в разы дольше моей функции). Возможно, Пирсон пригодится для автонахождения по известным сверткам.

wazzoo
KarapuZ
Наконец-то и я установил программу. Опробовать решил на записи двухнедельной давности-хороший проход с юга, несколько редких маяков и один "новичок"-иранский ARK-280. Запись тогда изучал тщательно, но после вашей программы обнаружил еще новеньких -пару иранцев (KLH-325, JRM-374),кипрский PHA-328 и греческий THR-307.Замечательная программа, большое спасибо!

vinipu
Спасибо за отзыв!
Рад, что все заработало

vinipu
Замечательная программа
Действительно, мощный инструмент получился для оффлайн анализа и идентификации сигналов NDB. Благодаря стараниям и как оказалось весьма незаурядным способностям Ивана к программированию и обработке сигналов.
Остаётся надеяться, что инструмент этот, примут на вооружение и наши, местные форумчане-снайперы
эфира. Если что - мы поможем в освоении софта.
"Белых пятен" еще много в этом деле, как одно из перспективных направлений - китайское, страна огромна
территориально, а данных по NDB оттуда с "гулькин нос".

KarapuZ
Вадим, спасибо за комплиман :)
Я тоже надеюсь, что активные NDB исследователи и NDB DXеры будут использовать эту прогу. Сам с нетерпением жду, когда доберусь до сетапа, чтобы вернуть Афедри на антенну и попробовать убрать грязь с ДВ - надеюсь дивно поохотиться :)
Заказал еще и усилок, на который ты меня навел (для ред-питайи) - хочу потом повесить его прямо на антенну, тоже поможет в борьбе с помехами (т.к. сейчас усиливает усилок в приемнике - вместе с грязью, прилетевшей через цепи, связывающие с компом).

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

BTW в бытовом плане тоже стало удобно. Я теперь вместо книжки перед сном - беру ноутбук и прощелкиваю файлики - сон навевает очень хорошо, и время использую эффективно :) и уши не вянут

В новой версии будет еще такая вот фишка (давно писал о такой задумке) - фильтр по уровню. Идея в следующем - когда на одной частоте два маяка, один из которых сильнее второго - увидеть слабый на его фоне очень сложно. Свертка при этом работать не будет, т.к. сигнал сильного маяка забивает сигнал слабого. Но - мы можем использовать тот факт, что сигнал сильного маяка лежит заметно выше уровня шума. Теперь берем и все отсчеты, которые выше некоего порога - превращаем.. в шум! Тогда этот сильный маяк уже не будет так сильно мешать свертке слабого маяка. Вот сегодня наткнулся, разбирая старую запись - когда мне это помогло вычислить слабый, при чем почти на уровне шума, маяк на фоне более сильного.

Картинка 1 - это с фильтрацией, обратите внимание на график свертки - я там красненьким подрисовал морзянку. Вообще говоря, порой именно график свертки помогает вычислить последовательность.
Картинка 2 - это без фильтрации. Тут вообще без шансов что-то разглядеть.





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

Еще задумки на будущее.

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

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

Что-то, в ступор меня ввела эта последовательность:), насколько я понял, её можно интерпретировать
как ЮЕ (YuE). Привод в логах NDB list значится как неактивный с 2002 года. Принял его позавчера,
частота 708 kHz:

KarapuZ
В базе есть такой:
708.000 UN 1008 1029 "ID+5""gap " Kurgan
Может, он? По формату похож - тоже 5 сек пауза (точнее 4,5)

wazzoo
Может, он? Да, похож. Но где-же пауза между знаками U и N?)

KarapuZ
Меня всегда такие маяки удивляют :)
Тут еще ладно - последовательность получается "непонятная". А вот как быть с 100% м NJ ( http://www.radioscanner.ru/uploader/2019/385.000_l402_nj_unid.png ) - который типа NWT ))

Тут еще ладно - последовательность получается "непонятная". А вот как быть с 100% м NJ ( http://www.radioscanner.ru/uploader/2019/385.000_l402_nj_unid.png ) - который типа NWT ))
Да, тяжёлый случай случай тоже однозначно:)
Кстати, опробовал активный мини-диполь Швабра за городом на выезде на предмет приёма маяков
в участке 200-1000 kHz. На "низах" 200-500 почему-то "не пошло", шумовая полка высокая вышла. От 500
и до 1000 kHz вроде как сносно, этот кусок сейчас и терзаю в Pskov NDB.

KarapuZ
Эх, жаль что не пошло. Хотя на выезде - если есть место - я бы бросил веревку подлиннее, если связка ноут+СДР помех не генерит.

Я сегодня готовил сетап на установку к себе - через неделю поеду:
1. Трансформаторный БП на 5Вольт на 1,5А (такой же как для Максуса, Сергей тогда посоветовал, до сих пор их делают).
- фильтры Murata BNX016 на выходах (2 выхода - на Афедри и на медиаконвертер)
- фильтр -60дБ на сетевом входе (нашел еще более качественные на Али).
Все это собрал в корпуса, промерил по NWT фильтры
2. Опторазвязка на медиаконвертерах - пашет как часы.
3. Однозвенный фильтр на антенну - а точнее просто конденсатор керамический (3,3 и 6,8 нФ - на месте выберу), померил в NWT на всяк случай - прямая наклонная характеристика - то что надо, по идее позволит выкрутить усиление на максимум как раз на ДВ участке для маяков - при этом не пострадают ни СВ ни тропики.

KarapuZ

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

wazzoo
AOR
Да, бухта полёвки была конечно с собой метров 100, но случился форс-мажор в виде проливного дождя,
Швабра оказалась весьма к стати для таких условий. В итоге хоть что-то удалось "поймать" и сделать соотв.
вывод - лишний инструмент не помешает. Пока, судя по данным из лога за эту сессию, радиус по приёму NDB
в среднем порядка 2000 km вышел, что для Швабры я считаю нормальным резалтом, ну и повод потестить её
за городом предоставился лишний раз.

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

AOR
Вот кстати я удивлен, но почитав ndblist наблюдаю весьма приличные DX на обычный мини-вип в авторском исполнении (Jay Law там). 3-5 тысячи км берет вообще не напрягаясь. На фоне скандинавов с их длинными антеннами смотрится очень прилично - не конкурент, но подпирает.

wazzoo
наблюдаю весьма приличные DX на обычный мини-вип в авторском исполнении
Резалты Roelof Bakker, PAoRDT на его мини-вип в кусте боярышника ещё более впечатляют... Сопоставимы
по дальности с уловами финнов.

Вот бы ему минивип от Хайо, было бы еще лучше...

AOR
Вот бы ему минивип от Хайо, было бы еще лучше... Швабра на ДВ и СВ работает значительно хуже чем оригинальный мини-вип

Собственно, вот лог, который снят со Швабры в участке 500-1000. Могло быть и лучше, но развернуться
толком помешал дождь, мини-диполь просто расположил между веток дерева, на высоте ~170 см, верт. поляризация:
Receiver QTH: KN95bg
Receiver: SDR Maxus
Antenna: Active mini-dipole Shvabra

Total loggings: 65
Total NDBs: 65
Total ITUs: 11
ITUs: KAL KAZ RSE RUS TKM TUR UKR UNID UZB XUE XUU
Max distance: 3016 km
________________________________________________
| |
| Distance, km NDBs ITUs |
|________________________________________________|

unknown distance 7 RUS UKR UNID XUE XUU
0-999 26 RUS TUR UKR
1000-1999 27 KAL KAZ RUS
2000-2999 4 RUS TKM UZB
3000-3999 1 RSE

__________________________________________________________________________________________________________
| |
| DD.MM.YYYY HHmm KHZ ID LSB USB Sec. Power Km Azi ITU QTH |
|__________________________________________________________________________________________________________|

01.03.2019 2108 434.0 ER -19 10,05 643 155 TUR Erzurum / Mudurge
01.03.2019 2108 438.0 OR -958 14,98 1634 23 RUS Kirov / Pobedilovo
01.03.2019 2108 439.0 MR -1018 4,28 2149 105 TKM Mary
01.03.2019 2108 564.0 OK -1016 1020 20,00 1290 20 RUS Lyskovo
01.03.2019 2108 565.0 KS -1014 14,99 1179 357 RUS Opalikha
01.03.2019 2108 570.0 DO -411 394 4,98 XUE UNID
01.03.2019 2108 570.0 FE -1017 1019 20,00 998 3 RUS Oktyabrskiy
01.03.2019 2108 581.0 DM -1012 10,05 1149 360 RUS Moscow - Ramenskoe
01.03.2019 2108 591.0 RN -397 399 20,00 1005 27 RUS Penza
01.03.2019 2108 595.0 A -387 395 15,09 733 6 RUS Voronezh / Chertovitskoye
01.03.2019 2108 600.0 US -1019 1022 3,88 UNID UNID
01.03.2019 2108 610.0 D -986 3,35 XUU UNID
01.03.2019 2108 642.0 KN -1019 1022 30,01 1228 359 RUS Moscow / Kostino
01.03.2019 2108 643.0 RZ -1018 20,00 1447 21 RUS Sharanga
01.03.2019 2108 644.0 SL -1050 1053 15,07 646 55 RUS Solodniki
01.03.2019 2108 659.0 DO -1015 1021 30,00 1122 359 RUS Moscow / Domodedovo
01.03.2019 2108 662.0 SM -1021 1030 15,05 76 136 RUS Smolenskaya
01.03.2019 2108 672.0 RS -1012 1019 30,08 636 332 UKR Rashivka
01.03.2019 2108 676.0 BO -1016 1030 10,05 RUS Nizhny Novgorod / Strigino
01.03.2019 2108 680.0 BP -1030 1035 30,03 1182 359 RUS Chelobityevo
01.03.2019 2108 680.0 M -1011 1026 20,00 1230 32 RUS Ulyanovsk / Vostochny
01.03.2019 2108 682.0 RD -1013 1033 10,05 RUS Rostov-na-Donu / Platov
01.03.2019 2108 687.0 R -992 14,75 412 327 UKR Dnipropetrovsk
01.03.2019 2108 687.0 P 992 14,75 416 327 UKR Dnipropetrovsk
01.03.2019 2108 690.0 DM -1011 1015 30,03 242 277 UKR Dmytrivka
01.03.2019 2108 694.0 AL -1015 1025 20,00 971 14 RUS Algasovo
01.03.2019 2108 696.0 BT -998 1004 30,56 1218 26 RUS Alatyr
01.03.2019 2108 700.0 MR -1007 1014 8,00 1192 357 RUS Moscow / Sheremetievo
01.03.2019 2108 705.0 BK 401 15,20 1857 35 RUS Perm / Bolshoye Savino
01.03.2019 2108 705.0 PX -400 15,19 1868 35 RUS Perm / Bolshoye Savino
01.03.2019 2108 708.0 UN -1020 1049 7,68 2228 49 RUS Kurgan
01.03.2019 2108 715.0 D -977 987 12,23 589 285 UKR Odessa / Central
01.03.2019 2108 718.0 NF -1016 1025 20,00 485 112 RUS Nalchik
01.03.2019 2108 718.0 BX -392 405 30,00 604 350 RUS Belgorod
01.03.2019 2108 722.0 O -1020 1025 20,01 911 37 RUS Saratov / Tsentralny
01.03.2019 2108 730.0 KT 960 5,76 321 92 RUS Stavropol / Shopakovskoye
01.03.2019 2108 732.0 AO -1016 1025 20,00 1098 1 RUS Aksinyino
01.03.2019 2108 733.0 ZK -1015 1026 10,11 UKR Kyiv / Zhuliany
01.03.2019 2108 735.0 FF -988 14,95 980 13 RUS Tambov / Donskoye
01.03.2019 2108 740.0 H -1010 1031 10,05 807 106 RUS Makhachkala / Uytash
01.03.2019 2108 745.0 BG -1018 1023 20,00 1199 350 RUS Karmanovo
01.03.2019 2108 749.0 K -1017 1024 20,00 1509 36 RUS Nizhnekamsk / Begishevo
01.03.2019 2108 752.0 BB -1018 1023 20,00 2658 29 RUS Beryozovo
01.03.2019 2108 755.0 T -1023 10,05 1071 74 KAZ Atyrau
01.03.2019 2108 760.0 R 1020 10,05 1140 0 RUS Ramenskoe
01.03.2019 2108 770.0 N -1003 1019 8,00 1192 358 RUS Sheremetyevo
01.03.2019 2108 777.0 ST -1008 1014 15,16 613 41 RUS Sirotinskaya
01.03.2019 2108 795.0 UH -1012 1024 15,01 2231 96 UZB Bukhara
01.03.2019 2108 815.0 DD -1004 994 5,76 831 7 RUS Lipetsk
01.03.2019 2108 825.0 KB -976 1017 15,24 786 319 UKR Kiev / Boryspil
01.03.2019 2108 825.0 KE 1022 14,76 776 318 UKR Kiev / Boryspil
01.03.2019 2108 840.0 KR -1008 1033 10,05 1633 317 KAL Kaliningrad / Khrabovo
01.03.2019 2108 876.0 WG -1007 1039 10,05 608 48 RUS Volgograd / Gumrak
01.03.2019 2108 900.0 NE -1017 1029 20,01 1307 360 RUS Nerl
01.03.2019 2108 905.0 UD -1015 1026 20,00 1242 20 RUS Buturlino
01.03.2019 2108 905.0 ZG -1053 964 29,70 821 79 RUS Zelenga (Astrakhan province)
01.03.2019 2108 910.0 C 1004 15,14 UNID UNID
01.03.2019 2108 932.0 KR -1026 1039 15,00 850 38 RUS Krasnoarmeysk
01.03.2019 2108 949.0 OE -1008 1019 10,05 1146 357 RUS Moscow / Vnukovo
01.03.2019 2108 960.0 RP -1017 1029 20,00 1224 15 RUS Sosnovkoye
01.03.2019 2108 965.0 SZ -400 411 29,85 3016 54 RSE Barabinsk
01.03.2019 2108 970.0 US -1022 1038 60,02 947 37 RUS Yelshanka
01.03.2019 2108 975.0 GD -1006 1016 30,55 926 356 RUS Maloye-Skuratovo
01.03.2019 2108 985.0 FK -1016 1025 20,00 1165 350 RUS Gagarin
01.03.2019 2108 995.0 DW -907 919 15,06 1259 35 RUS Nikolskoye

на СДВ и ДВ Швабра отстает от оригинального мини-випа, на СВ Швабра ведет себя посредственно относительно мини-випа, на КВ практически одинаково. Но так как расположение антенн у каждого индивидуально и уникально, то результаты могут отличатся в тут или другу стороны, хотя думаю что мини-вип все же будет выигрыше на НЧ. Правда Швабра менее подвержена техногенным помехам, но не сказать чтобы сильно выделялась этим.

KarapuZ
Насколько сложно было анализировать среди вещалок?

wazzoo
Насколько сложно было анализировать среди вещалок?
Оказалось - не сложно, даже много быстрее, чем "внизу". Анализировал по шаблонам и при наличии
несущих-кандидаток.

Сию Piaggio Avanti облётывает аэропорт Гагарин в Саратове.
Если есть у кого в Поволжье возможность, послушайте пожалуйста эфир на предмет появления новых приводов и их частот.

Ну что, готов конвуляционный фильтр. Еще более глубокий анализ в условиях "несколько маяков на частоте". Фильтры по уровню отдыхают.
Конвуляционный фильтр позволяет достать совсем слабые маяки под сильными или средней силы маяками - при чем в полуавтоматическом режиме - т.е. фильтруется уже не уровень, а конкретно сигнал мешающего маяка (его отсчеты). После чего жмакается кнопка повторного авто-поиска интервала и находится следующий более слабый маяк. Теперь фильтруются уже его сигналы. Снова жмакаем авто-поиск интервала - и уже третий, совсем на уровне шума маяк выплывает сам.

Ниже даю скриншоты - видео не стал писать, по сути это просто рутинные кликанья от оператора (потом запишу). Картинки лучше показывают процесс и результат.

1. Встаем на частоту


2. Настраиваем фильтр (красным подсвечиваются фильтруемые области свертки)


3. Применяем фильтр, вот что получаем на свертке в старом интервале


4. Жмакаем авто-поиск интервала. Получаем следующий маяк.


5. Фильтруем теперь его


6. Получаем в его интервале такую картину


7. Жмакаем автоинтервал - находим последний, 3-й маяк, с трудом (но читается!) RS


8. Теперь ради интереса - зная интервал 3-го маяка, посмотрим, как выглядит свертка в этом интервале без фильтра:

Для иллюстрации - видео. Записал на другой частоте - суть та же. Три маяка на частоте. Фильтр+автоинтервал позволяют увидеть все три маяка довольно быстро.

Применение фильтра по свертке

Интерфейс в стадии экспериментов - потом чуть поменяется.

wazzoo
Ну что, готов конвуляционный фильтр. Еще более глубокий анализ в условиях "несколько маяков на частоте". Фильтры по уровню отдыхают.
Конвуляционный фильтр позволяет достать совсем слабые маяки под сильными или средней силы маяками - при чем в полуавтоматическом режиме - т.е. фильтруется уже не уровень, а конкретно сигнал мешающего маяка (его отсчеты). После чего жмакается кнопка повторного авто-поиска интервала и находится следующий более слабый маяк. Теперь фильтруются уже его сигналы. Снова жмакаем авто-поиск интервала - и уже третий, совсем на уровне шума маяк выплывает сам.
Иван, на мой взгляд - мега круто! Фильтрация, да ещё ступенчатая временнЫх отсчётов - в анализе маяков
такого ещё не было. Чешуться руки попробовать и вывернуть спектр наизнанку:)
На какой длительности файла записи пробовал? Мне кажется, на коротких файлах будут проблемы с накоплением "полезной свертки" - при наложении последовательностей по времени.