Отводов на все деньги) А если конкретнее, 1024. Код буквально тот же, только входной вектор транспонирован. Т.е. это не ПВОС, а частотный режектор после ПОСа, вторым прогоном
А Вы можете показать на этой выборке больший коэффициент подавления? Всем будет интересно посмотреть.
Особо не подбирал параметры, выборка та же, её объем тот же, но добавил time-domain обработку.
У SNRа где-то 2 дБ прирост, спектр и осциллограммы веселее:
1. Наследуем Марка Твена: Есть ложь, наглая ложь, статистика и величина бина над шумом после БПФа. =) Этот размах зависит от настроек спектроанализатора (RBW у аппаратного, длина выборки у БПФа) и в отрыве от них не показателен.
Спасибо за статью, возомню себя рецензентом и дам пару замечаний:
1. «От максимума сигнала до шума примерно 60 дБ» — по осциллограмме не похоже. Там весь динамический диапазон 72 дБ, а шум далеко не пару разрядов.
2." А в широкой полосе маленькие остатки помех «тонут» в шуме. Ниже шума подавить помеху у нас не получилось. Физика. Хотя во временной области гармоника остается видимой." — физика или нет, это интерпретации для самоуспокоения) Если даже человек видит признаки, по которым можно отделить помеху от полезных процессов, то математика тем более найдет за что зацепиться при правильной постановке задачи. Это конкретный используемый алгоритм при конкретных настройках не смог справиться.
3. Символами R, r всё-таки обычно обозначают саму корреляционную/ковариационную матрицу, а не результат её инверсии
А нет ли информации, как отправка пакета с борта привязана ко времени? Борт отправляет когда вздумается или, например, строго в начале миллисекунды? Если второе, то миллисекунды по каким часам?
Гениально) Только получается, что мы таким образом вынесли GPS из ADS-B приемника в самолет. И если сломают GPS на самолетах, то сломается и MLAT. Но, Flight Radar вряд ли это беспокоит. Вероятно, всё ради удешевления ADS-B приемника? При размещении приемника его координаты определяются единожды сторонним устройством, например, телефоном?
Почему станции FlightRadar должны стоять не ближе 15 км?
А информации по их синхронизации нет? Ведь чтобы за один шаг сравнить время прилета сигнала, они должны жить в одной шкале времени. Если она берется из GPS, то глушение GPS убьет и MLAT.
Спасибо за статью! Чем отличается система в Пулково от MLAT FlightRadar? Как они осуществляют синхронизацию наземных приемников? Где можно почитать про Пулковский вариант? И какая точность определения координат получается?
Отличная статья, спасибо! А необходимость иметь корреляционные матрицы не ставит крест на алгоритме, они вполне себе оцениваются при обработке реальных сигналов
Сам документ 12 01 или 12 02, это электронный документ
Имеется в виду какая-та специальная электронная система документооборота? Или просто вордовский документ 12 01 с парой страниц: титульник, лист регистрации изменений и лист с табличкой?
В ЕСПД используется смущающее многих программистов обозначение исходных текстов программы и результата сборки “документами”. Документ “текст программы”, согласно 19.101-77, имеет обозначение 12. Дальше принято, что исходники обозначаются как 12 01 — т.е. 01(первый) документ вида 12, а бинарники — как 12 02 — т.е. второй документ вида 12.
12 01 при этом оформляет в виде бумажки или в виде диска с исходниками? И вообще, что из себя физически представляет, скажем, 12 01 — записанный CD с обложкой, на котором написан номер документа?
Судя по фразе «The low orbit Satellite Based Augmentation System» из этой заметки, имеется в виду функциональное дополнение к существующей Beidou. То есть они со спутника будут передавать дифпоправки, с помощью которых можно будет компенсировать часть ошибок в существующей системе. Иначе говоря, это будет не навигационный сигнал, а просто система передачи данных. Но так до 10 см в реальном времени погрешность не опустить.
Ждал инсайдиков про Сферу в комментах, но так и не нашел. В каком виде они собираются навигацию прикручивать? Тем более с такой точностью на низких орбитах?
6 ГГц при полосе за 50 МГц, вход внешней опоры и синхроимпульс — явно для объединения нескольких таких в массив с последующей пространственной обработкой. Похоже на пеленгатор, короче)
Мы сейчас начнем спорить о терминах, что малоприятно, но моя позиция такая: шум (noise) не может быть синусоидой, синусоида — это уже помеха (interference).
Под шумом обычно понимают слабокоррелированный (широкополосный) случайный процесс, как правило тепловой природы. В контексте этой задачи это дает то, что в отдельном бине этого шума мало мощности.
Теорема Котельникова говорит об условиях, при которых аналоговый сигнал может быть взаимнооднозначно представлен дискретной выборкой. В приложении к данной работе под «аналоговым сигналом» следует понимать шум+синусоиду. Если вы хотите точно воспроизвести шум+синусоиду, деваться некуда, берите 2 МГц частоту семплирования.
Но дальше происходит подмена понятий. Полезной нарекается только синусоида. А для её отображения вообще достаточна нулевая частота дискретизации (запишите на бумажке амплитуду и фазу и не мучайте АЦП). Да, когда мы вместо 2 МГц взяли частоту 10 кГц, шумы многократно наложились, мы не выжали потенциальный SNR. Но при расчете отношения сигнал/шум надо рассматривать шум не в полосе 5 кГц, а только в бине самой синусоиды. А там SNR может быть уже вполне достаточный (не максимальный, а достаточный!).
Иными словами, если ваш полезный сигнал — синусоида, то поставьте фильтр на выходе процедуры дискретизации, а не на входе.
Такие мысли возникли у меня при прочтении статьи по диагонали. Кто выиграет $1000, все деньги можете оставить себе, я не претендую ;-)
Человек из правильных вводных данных сделал неправильные выводы, бывает.
Эффекты ОТО действительно заметны, но жить совершенно не мешают. Для удобства на спутниках ГЛОНАСС частота опорного генератора немного понижена, что-то в духе 4.999сколько-то-там-девяток99 МГц вместо 5.0 МГц, насколько помню. Но это только для удобства.
Часы спутника уходят, это принципиальный факт и его нужно учитывать. А дальше, пусть они хоть на сутки за сутки уходят, главное, чтобы мы знали как они уходят. В случае с ОТО всё просто — линейный, элементарно компенсируемый тренд. Жить не помешал бы, даже если бы привентивных мер не принимали бы. Было бы в модели псевдодальности ни десять, а одиннадцать слагаемых. Ну и ЧСВ разработчиков аппаратуры повысилось: эвона какие высокие материи учитывать в навигационной задаче приходится!
Я не знаю, ни разу P код не принимал. ИКД, естественно, тоже не видел =) 500кГц — уже сложнее, потеряем примерно 13 дБ. Но всё равно реально, особенно при условиях, в которых обычно работает геодезический приемник.
«если мы берем немодулированную гармонику» т.е. немодулированный гармонический сигнал то количество бит ЦАП модулятора = бесконечности, ибо это аналоговый сигнал, тогда надо использовать бесконечно-битный АЦП, что собственно все и делают если используют аналог.
Мы с вами сошлись на мысли, что гармоника (та самая, которую модулируют) — аналоговый сигнал. Модулятор её дёргает по фазе на 180 градусов. По свойствам тригонометрических функций это эквивалентно умножению на плюс или минус единицу. Так что же в этой процедуре отнимает у сигнала его аналоговость и ставит на новую ступень абстракции, требующую однобитного и обязательно однобитного АЦП?
Подождите, вы меня запутали. Вы продекларировали настолько красивое и изящное правило, что я захотел стать вашим последователем: по вашим словам, количество бит ЦАП модулятора = количество бит АЦП приемника. Почему же для немодулированного сигнала вы ломаете эту стройную картину мира и предлагаете 1-битный АЦП?