RU2494478C1 - Передискретизация в комбинированном банке фильтров транспозитора - Google Patents
Передискретизация в комбинированном банке фильтров транспозитора Download PDFInfo
- Publication number
- RU2494478C1 RU2494478C1 RU2012116009/08A RU2012116009A RU2494478C1 RU 2494478 C1 RU2494478 C1 RU 2494478C1 RU 2012116009/08 A RU2012116009/08 A RU 2012116009/08A RU 2012116009 A RU2012116009 A RU 2012116009A RU 2494478 C1 RU2494478 C1 RU 2494478C1
- Authority
- RU
- Russia
- Prior art keywords
- analysis
- synthesis
- frequency
- signal
- window
- Prior art date
Links
- 238000004458 analytical method Methods 0.000 claims abstract description 295
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 273
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 273
- 230000017105 transposition Effects 0.000 claims abstract description 139
- 238000012545 processing Methods 0.000 claims abstract description 51
- 230000009466 transformation Effects 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 47
- 238000005070 sampling Methods 0.000 claims description 46
- 238000006243 chemical reaction Methods 0.000 claims description 31
- 230000008859 change Effects 0.000 claims description 23
- 230000014509 gene expression Effects 0.000 claims description 6
- 230000010363 phase shift Effects 0.000 claims description 6
- 238000000844 transformation Methods 0.000 claims 1
- 230000005236 sound signal Effects 0.000 abstract description 11
- 239000000126 substance Substances 0.000 abstract 1
- 230000006870 function Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 5
- 238000013507 mapping Methods 0.000 description 4
- 230000000737 periodic effect Effects 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000003775 Density Functional Theory Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- RVRCFVVLDHTFFA-UHFFFAOYSA-N heptasodium;tungsten;nonatriacontahydrate Chemical compound O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W].[W] RVRCFVVLDHTFFA-UHFFFAOYSA-N 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000009527 percussion Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
- G10L19/02—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders
- G10L19/0204—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders using subband decomposition
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
- G10L19/02—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders
- G10L19/022—Blocking, i.e. grouping of samples in time; Choice of analysis windows; Overlap factoring
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
- G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
- G10L19/26—Pre-filtering or post-filtering
- G10L19/265—Pre-filtering, e.g. high frequency emphasis prior to encoding
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/038—Speech enhancement, e.g. noise reduction or echo cancellation using band spreading techniques
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/04—Time compression or expansion
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/118—Prognosis of disease development
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Quality & Reliability (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Genetics & Genomics (AREA)
- Analytical Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Molecular Biology (AREA)
- Microbiology (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Networks Using Active Elements (AREA)
- Static Random-Access Memory (AREA)
Abstract
Изобретение относится к кодированию звуковых сигналов, а именно к способам восстановления высоких частот, включая гармонический транспозитор частотной области. Техническим результатом является повышение качества звучания посредством восстановления высоких частот звукового сигнала усовершенствованной гармонической транспозицией при низкой дополнительной вычислительной сложности. Система содержит банк фильтров анализа (501), содержащий модуль преобразования анализа (601), имеющий коэффициент разрешения по частоте Δf; и окно анализа (611), имеющее длительность DA; причем банк фильтров анализа (501) сконфигурирован с возможностью создания набора сигналов субполос анализа из низкочастотной компоненты сигнала; модуль нелинейной обработки (502, 650), сконфигурирован с возможностью определять набор сигналов субполос синтеза по части набора сигналов субполос анализа, причем часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т; и банк фильтров синтеза (504), содержащий модуль преобразования синтеза (602), имеющий разрешение по частоте QΔf, и окно синтеза (612), имеющее длительность DS; причем банк фильтров синтеза (504) сконфигурирован с возможностью создавать высокочастотную компоненту сигнала по набору сигналов субполос синтеза; где Q - коэффициент разрешения по частоте, такой, что Q≥1 и меньше порядка транспозиции Т; и при этом значение произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа выбирается, исходя из коэффициента разрешения по частоте Q. 6 н. и 19 з.п. ф-лы, 12 ил.
Description
Область техники
Настоящее изобретение относится к кодированию звуковых сигналов, а именно к способам восстановления высоких частот, включая гармонический транспозитор частотной области.
Уровень техники
Технологии восстановления высоких частот (HFR), например технология восстановления спектральной полосы (SBR), позволяют значительно повысить эффективность кодирования традиционных перцептуальных аудио кодеков. В сочетании с форматом ААС стандарта MPEG-4, технологии HFR позволяют создавать весьма эффективные аудио кодеки, которые в настоящее время используются в составе систем спутникового радиовещания ХМ (ХМ Satellite Radio) и цифрового радиовещания Mondiale (Digital Radio Mondiale), а также стандартизированы в рамках проектов 3GPP, DVD Forum (Международная ассоциация производителей оборудования, носителей, программ и фильмов, и других пользователей DVD) и других организаций. Комбинация ААС и SBR получила название aacPlus. В качестве составной части она внесена в стандарт MPEG-4 под названием «Высокоэффективный метод кодирования ААС» (НЕ-ААС). В целом, технология HFR может быть использована в сочетании с любым перцептуальным кодеком, обеспечивая совместимость в прямом и обратном направлении, предоставляя таким образом возможность для модернизации существующих систем вещания, например MPEG Layer-2, применяемой в системе вещания Eureka DAB. Способы транспозиции HFR могут использоваться также в сочетании с речевыми кодеками для получения более широкой полосы частот речевых сигналов при сверхнизкой скорости передачи данных.
Основной принцип, лежащий в основе HRF, состоит в том, что между характеристиками высокочастотной полосы спектра сигнала и характеристиками низкочастотной полосы спектра того же сигнала обычно наблюдается сильная корреляция. Таким образом, при помощи транспозиции сигнала из низкочастотной области в высокочастотную область может быть достигнуто хорошее приближение для восстановления исходной высокочастотной области сигнала.
Данный принцип транспозиции был обоснован в WO 98/57436, включенной в настоящий документ как неотделимая часть посредством ссылки, как способ для восстановления полосы высоких частот из полосы более низких частот звукового сигнала. При использовании данного принципа возможно получение значительной экономии в скорости передачи данных при кодировании звуковых и речевых сигналов. В дальнейшем речь будет идти о кодировании звуковых сигналов, но необходимо отметить, что описываемые способы и системы в равной степени пригодны для кодирования речевых сигналов и комплексного кодирования речевых и звуковых сигналов (USAC).
В системе кодирования звуковых сигналов, основанной на принципе HFR, сигнал с малой шириной полосы частот поступает на вход кодировщика основной формы сигнала для кодирования, а более высокие частоты восстанавливаются в декодере, который использует транспозицию сигнала с малой шириной полосы частот и дополнительную вспомогательную информацию, как правило, закодированную с очень низкой скоростью передачи данных и описывающую требуемую форму спектра. Для малых скоростей передачи данных, когда ширина полосы основного кодированного сигнала мала, восстановление или синтезирование полосы высоких частот имеет еще большее значение, т.к. высокочастотные составляющие звукового сигнала придают ему положительные для восприятия характеристики.
Одной из основных проблем, характерных для алгоритмов гармонического HFR, являются несовместимые ограничения, связанные с достижением требуемого высокого разрешения по частоте для получения высококачественной транспозиции для стационарных звуковых сигналов и с переходной характеристикой системы для импульсных звуковых сигналов или звуковых сигналов ударных инструментов. То есть, в то время как использование высокого разрешения по частоте оказывает положительное влияние на транспозицию стационарных сигналов, такое высокое разрешения по частоте, как правило, требует большей величины окна, что ухудшает характеристики системы при работе с импульсными элементами сигнала. Один из подходов к решению этой проблемы может заключаться в адаптивном изменении окон транспозитора, например, с помощью изменения величины окна как функции характеристик входного сигнала. Как правило, для достижения высокого разрешения по частоте для стационарных элементов сигнала используются более длинные окна, тогда как короткие окна используются для импульсных элементов сигнала с целью достижения хорошей переходной характеристики, т.е. высокого временного разрешения транспозитора. Тем не менее, недостаток такого подхода заключается в том, что операции по анализу сигнала, например определение наличия импульсов и т.п., должны быть встроены в систему транспозиции. Такие операции по анализу сигнала часто содержат этапы принятия решения, например решения о наличии импульса, которые инициируют изменения в процессе обработки сигнала. Более того, такие действия, как правило, влияют на надежность системы и могут вносить искажения сигнала при изменении процесса обработки сигнала, например при переключении между окнами различной величины.
Для повышения качества звучания и для синтезирования необходимой ширины полосы высоких частот сигнала алгоритмы гармонического HFR, как правило, используют несколько порядков транспозиции. Для выполнения совокупности транспозиций различных порядков транспозиции в устройствах известного уровня техники требовалось использование совокупности банков фильтров либо на этапе анализа, либо на этапе синтеза, или же на обоих этапах. Как правило, для каждого отдельного порядка транспозиции требовался отдельный банк фильтров. Более того, в ситуациях, когда кодер основной формы сигнала работает на меньшей частоте дискретизации, чем частота дискретизации конечного выходного сигнала, как правило, возникает дополнительная необходимость в преобразовании основного сигнала до частоты дискретизации выходного сигнала, и эта повышающая дискретизация основного сигнала обычно достигается при помощи еще одного банка фильтров. В итоге, с ростом числа различных порядков транспозиции значительно возрастает вычислительная сложность.
Настоящий документ предлагает решения вышеописанных проблем в отношении переходной характеристики гармонической транспозиции и в отношении вычислительной сложности. В результате при низкой дополнительной сложности достигается усовершенствование гармонической транспозиции.
Краткое описание изобретения
Согласно одному из признаков изобретения, предлагается система, сконфигурированная для создания высокочастотной компоненты сигнала из низкочастотной компоненты этого сигнала. Система может содержать банк фильтров анализа, содержащий модуль преобразования анализа с разрешением по частоте Δf. Модуль преобразования анализа может быть сконфигурирован производить, например, преобразование Фурье, быстрое преобразование Фурье, дискретное преобразование Фурье или вейвлет-преобразование. Кроме того, банк фильтров анализа может включать в себя окно анализа длительностью DA. Окно анализа может иметь вид, например, окна Гаусса, косинус-окна, окна Хэмминга, окна Ханна, прямоугольного окна, окна Бартлетта или окна Блэкмана. Банк фильтров анализа может быть сконфигурирован с возможностью создавать набор сигналов субполос анализа из низкочастотной компоненты сигнала.
Система может содержать модуль нелинейной обработки, сконфигурированный с возможностью определять набор сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т. В частности, сигналы субполос могут содержать комплексные величины и фазовый сдвиг может быть выполнен путем умножения фазы комплексных величин субполос на порядок Т.
Система может содержать банк фильтров синтеза, включающй модуль преобразования синтеза с разрешением по частоте QΔf. Модуль преобразования синтеза может быть сконфигурирован с возможностью производить соответствующее преобразование, обратное тому преобразованию, которое произведено модулем преобразования анализа. Более того, банк фильтров синтеза может включать в себя окно синтеза, имеющее длительность DS и имеющее любую из вышеперечисленных форм. Q является коэффициентом разрешения по частоте, при этом Q≥1 и меньше, чем порядок транспозиции Т. В одном из конкретных вариантов воплощения, коэффициент разрешения по частоте выбирается, как Q>1. Банк фильтров синтеза может быть сконфигурирован с возможностью создавать высокочастотную компоненту сигнала по набору сигналов субполос синтеза.
Как правило, значение произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа выбирается, исходя из коэффициента разрешения по частоте Q. В частности, произведение ΔfDA может быть пропорционально
. В одном из вариантов воплощения, значение произведения ΔfDA меньше или равно
. Более того, произведение ΔfDA может быть больше, чем
. Значение произведения ΔfDA банка фильтров анализа может быть равно значению произведения QΔfDS банка фильтров синтеза. Выбирая банк фильтров анализа и/или синтеза в соответствии с вышеописанными соотношениями, можно уменьшить либо полностью устранить искажения, обусловленные гармонической транспозицией импульсов сигнала, при сохранении невысокой вычислительной сложности гармонического транспозитора.
Кроме того, система может содержать второй модуль нелинейной обработки, сконфигурированный с возможностью определять второй набор сигналов субполос синтеза по набору сигналов субполос анализа, используя второй порядок транспозиции T2; причем второй набор сигналов субполос синтеза определяется исходя из части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на второй порядок транспозиции T2. Порядок транспозиции T и второй порядок транспозиции T2 могут быть различны. Кроме того, система может содержать суммирующий модуль, сконфигурированный с возможностью производить суммирование набора сигналов субполос синтеза и второго набора сигналов субполос синтеза, получая таким образом суммарный набор сигналов субполос синтеза, поступающий на вход банка фильтров синтеза. Суммирующий модуль может быть сконфигурированным с возможностью производить сложение или усреднение соответствующих сигналов субполос из первого набора сигналов субполос синтеза и второго набора сигналов субполос синтеза. Другими словами, суммирующий модуль может быть сконфигурирован с возможностью совместить сигналы субполос синтеза первого набора сигналов субполос синтеза и второго набора сигналов субполос синтеза в соответствии с перекрывающимися интервалами частот.
В одном из вариантов воплощения, банк фильтров анализа может содержать некоторое число KA субполос анализа, где KA>1, и k является индексом субполосы анализа, принимающим значения k=0,…,KA-1. Банк фильтров синтеза может содержать некоторое число NS субполос синтеза, где NS>0, и n является индексом субполосы синтеза, принимающем значения n=0,…,NS-1. В подобных случаях модуль нелинейной обработки может быть сконфигурирован с возможностью производить определение сигнала n-ой субполосы синтеза из набора сигналов субполос синтеза, исходя из сигнала k-ой субполосы анализа и сигнала (k+1)-ой субполосы анализа из набора сигналов субполос анализа. В частности, модуль нелинейной обработки может быть сконфигурирован с возможностью производить определение фазы сигнала n-ой субполосы синтеза в виде суммы фазового сдвига сигнала k-ой субполосы анализа и фазового сдвига сигнала (k+1)-ой субполосы анализа. Более того, модуль нелинейной обработки может быть сконфигурирован с возможностью производить определение величины сигнала n-ой субполосы синтеза в виде произведения возведенной в степень величины сигнала k-ой субполосы анализа и возведенной в степень величины сигнала (k+1)-ой субполосы анализа.
Индекс субполосы анализа k для сигнала субполосы анализа, вносящего вклад в субполосу синтеза с индексом субполосы синтеза n, может выражаться целым числом, полученным в результате выделения целой части выражения
; при этом остаток r может выражаться как
. В таких случаях модуль нелинейной обработки может быть сконфигурирован с возможностью производить определение фазы сигнала n-ой субполосы синтеза в виде суммы фазы сигнала k-ой субполосы анализа, умноженной на Т(1-r), и фазы сигнала (k+1)-ой субполосы анализа, умноженной на T(r), т.е. путем выполнения линейной интерполяции фазы. Более того, модуль нелинейной обработки может быть сконфигурирован с возможностью производить определение величины сигнала n-ой субполосы синтеза в виде произведения величины сигнала k-ой полосы анализа, возведенной в (1-r) степень, и величины сигнала (k+1)-ой субполосы анализа, возведенной в r степень, т.е. путем определения среднего геометрического значения этих величин.
Банк фильтров анализа и банк фильтров синтеза могут быть четным образом набраны так, чтобы центральная частота субполосы анализа соответствовала kΔf, a центральная частота субполосы синтеза соответствовала nQΔf. В альтернативном варианте воплощения, банк фильтров анализа и банк фильтров синтеза могут быть нечетным образом набраны так, чтобы центральная частота субполосы анализа соответствовала
, а центральная частота субполосы синтеза соответствовала
; и разность между порядком транспозиции Т и коэффициентом разрешения Q была четной.
Частота дискретизации для низкочастотной компоненты может быть равна fA. Модуль преобразования анализа может выполнять дискретное М точечное преобразование. Длина окна анализа может составлять LA отсчетов, и/или сдвиг окна анализа может составлять ΔsA отсчетов по низкочастотной компоненте. В таких случаях разрешение по частоте может иметь вид
, длительность окна анализа может иметь вид
, и/или физический временной шаг банка фильтров анализа может иметь вид
.
Частота дискретизации для высокочастотной компоненты может иметь вид fS=QfA. Модуль преобразования синтеза может выполнять дискретное М точечное преобразование, в частности, он может выполнять соответствующее преобразование, обратное тому преобразованию, которое произведено модулем преобразования анализа. Длина окна синтеза может составлять LS отсчетов, и/или окно синтеза может быть сдвинуто на величину сдвига синтеза, длительностью ΔsS отсчетов по высокочастотной компоненте. В таких случаях разрешение по частоте может иметь вид
, длительность может иметь вид
, и/или физический временной шаг банка фильтров синтеза может иметь вид
.
В соответствии со следующим признаком изобретения, предлагается система для создания выходного сигнала, содержащего высокочастотную компоненту, по входному сигналу, содержащему низкочастотную компоненту, при помощи транспозиции порядка Т. Данная система может содержать модуль окна анализа, сконфигурированный с возможностью применять окно анализа длиной LA отсчетов, получая таким образом фрейм входного сигнала. Система может содержать модуль преобразования анализа порядка М и имеющий разрешение по частоте Δf, сконфигурированный с возможностью производить преобразование LA отсчетов в М комплексные коэффициенты. Система может содержать модуль нелинейной обработки, сконфигурированный с возможностью производить изменение фазы комплексных коэффициентов, используя порядок транспозиции Т. Изменение фазы может заключаться в сдвиге фазы комплексных коэффициентов в порядке, изложенном в настоящем документе. Система может содержать модуль преобразования синтеза порядка М, имеющий разрешение по частоте QΔf, сконфигурированный с возможностью производить преобразование измененных коэффициентов в М измененные отсчеты; при этом Q является коэффициентом разрешения по частоте меньшим, чем порядок транспозиции Т. Кроме того, система может содержать модуль окна синтеза, сконфигурированный с возможностью применять окно синтеза длиной LS отсчетов к М измененным отсчетам, создавая таким образом фрейм выходного сигнала.
М может зависеть от коэффициента разрешения по частоте Q. В частности, разность между М и средней длиной окна анализа и окна синтеза (612) может быть пропорциональна (Q-1). В одном из вариантов воплощения, значение М больше или равно (QLA+LS)/2. В воплощении, значение М может быть меньше, чем (TLA+LS)/2.
Согласно одному из признаков изобретения, предлагается способ создания высокочастотной компоненты сигнала из низкочастотной компоненты этого сигнала. Данный способ может включать в себя этап обеспечения набора сигналов субполос анализа из низкочастотной компоненты сигнала, используя банк фильтров анализа, содержащий модуль преобразования анализа с разрешением по частоте Δf и окно анализа длительностью DA. Кроме того, способ может включать в себя этап определения набора сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т. Наконец, способ может включать в себя этап обеспечения высокочастотной составляющей сигнала по набору сигналов субполос синтеза, используя банк фильтров синтеза, содержащий модуль преобразования синтеза с разрешением по частоте QΔf и окно синтеза длительностью DS. Q является коэффициентом разрешения, причем Q≥1 и меньше, чем порядок транспозиции Т. Значение произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа может выбираться исходя из коэффициента разрешения по частоте Q.
В соответствии со следующим признаком изобретения, предлагается способ создания выходного сигнала, содержащего высокочастотную компоненту, из входного сигнала, содержащего низкочастотную компоненту, при помощи транспозиции порядка Т. Данный способ может включать в себя применение окна анализа длиной LA отсчетов, таким образом получая фрейм входного сигнала, и преобразование этого фрейма из LA отсчетов входного сигнала в М комплексных коэффициентов, используя преобразование анализа порядка М с разрешением по частоте Δf. Кроме того, данный способ может включать в себя этап изменения фазы комплексных коэффициентов, используя порядок транспозиции Т.
Изменение фазы может производиться в соответствии со способами, изложенными в настоящем документе. Кроме того, способ может включать в себя этапы по преобразованию измененных коэффициентов в М измененных отсчетов, используя преобразование синтеза порядка М с разрешением по частоте QΔf, причем Q является коэффициентом разрешения по частоте меньшим, чем порядок транспозиции Т; и применению окна синтеза длиной LS отсчетов к М измененным отсчетам, создавая таким образом фрейм выходного сигнала. М может зависеть от коэффициента разрешения по частоте Q.
Согласно следующему признаку изобретения, предлагается способ проектирования гармонического транспозитора для создания высокочастотной компоненты сигнала из низкочастотной компоненты этого сигнала. Данный способ может включать в себя этап создания банка фильтров анализа, содержащего модуль преобразования анализа с разрешением по частоте Δf; а также окно анализа длительностью DA, причем банк фильтров анализа сконфигурирован с возможностью создавать набор сигналов субполос анализа из низкочастотной компоненты сигнала. Кроме того, способ может включать в себя этап создания модуля нелинейной обработки, сконфигурированного с возможностью определять набор сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т. Кроме того, способ может включать в себя этап создания банка фильтров синтеза, состоящего из модуля преобразования синтеза с разрешением по частоте QΔf; а также окна синтеза длительностью DS; причем банк фильтров синтеза сконфигурирован с возможностью создавать высокочастотную компоненту сигнала по набору сигналов субполос синтеза; при этом Q является коэффициентом разрешения по частоте, таким, что Q≥1 и меньше порядка транспозиции Т. Кроме того, способ может включать в себя этап выбора значения произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа, исходя из коэффициента разрешения по частоте Q.
В соответствии со следующим признаком изобретения, предлагается способ проектирования транспозитора, сконфигурированного с возможностью создания выходного сигнала, содержащего высокочастотную компоненту, из входного сигнала, содержащего низкочастотную компоненту, при помощи транспозиции порядка Т. Способ может включать в себя этапы создания модуля окна анализа, при этом модуль сконфигурирован с возможностью применения окна анализа длиной LA отсчетов, и получения таким образом фрейма входного сигнала; создания модуля преобразования анализа порядка М с разрешением по частоте Δf, при этом модуль сконфигурирован с возможностью преобразования LA отсчетов в М комплексных коэффициентов. Более того, способ может включать в себя этап создания модуля нелинейной обработки, сконфигурированного с возможностью изменения фазы комплексных коэффициентов, используя порядок транспозиции Т. Кроме того, способ может включать в себя этапы по созданию модуля преобразования синтеза порядка М с разрешением по частоте QΔf, при этом модуль сконфигурирован с возможностью преобразования измененных коэффициентов в М измененных отсчетов; причем Q является коэффициентом разрешения по частоте, меньшим, чем порядок транспозиции Т; а также создания модуля окна синтеза, при этом модуль сконфигурирован с возможностью применения окна синтеза длиной LS отсчетов к М измененным отсчетам, создавая таким образом фрейм выходного сигнала. Наконец, способ может включать в себя этап выбора параметра М, исходя из коэффициента разрешения по частоте Q.
Необходимо отметить, что данные способы и системы, включая предпочтительные варианты воплощения, изложенные в настоящей патентной заявке, могут быть использованы отдельно или в сочетании с иными способами и системами, описанными в этом документе. Более того, все разновидности способов и систем, описанных в настоящей патентной заявке, могут быть скомбинированы произвольным образом. В частности, признаки, входящие в формулу изобретения, могут быть скомбинированы друг с другом произвольным образом.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Настоящее изобретение будет описано на иллюстративных примерах, не ограничивающих объема и сущности изобретения, со ссылками на прилагаемые чертежи, где:
Фиг.1 представляет принцип действия образца гармонического транспозитора частотной области первого порядка;
Фиг.2 представляет принцип действия образца гармонического транспозитора, использующего несколько порядков транспозиции;
Фиг.3 представляет принцип действия образца прототипа гармонического транспозитора, использующего несколько порядков транспозиции, с применением общего банка фильтров анализа;
Фиг.4 представляет принцип действия образца прототипа гармонического транспозитора, использующего несколько порядков транспозиции, с применением общего банка фильтров синтеза;
Фиг.5 представляет принцип действия образца гармонического транспозитора, использующего несколько порядков транспозиции, с применением общего банка фильтров синтеза и общего банка фильтров анализа;
Фиг.5b представляет образец отображения сигналов субполос для схемы параллельного транспозитора в соответствии с Фиг.5;
Фиг.6а представляет образец параллельного транспозитора порядка T=2, 3, 4, использующего общий банк фильтров анализа и раздельные банки фильтров синтеза;
Фиг.6b представляет образец параллельного транспозитора порядка T=2, 3, 4, использующего общий банк фильтров анализа и общий банк фильтров синтеза;
Фиг.7 представляет образец отображения сигналов субполос для схемы параллельного транспозитора в соответствии с Фиг.6b;
Фиг.8 представляет частный случай расположения функции Дирака в окнах анализа и синтеза гармонического транспозитора;
Фиг.9 представляет частный случай другого расположения функции Дирака в окнах анализа и синтеза гармонического транспозитора; и
Фиг.10 представляет случай расположения функции Дирака согласно Фиг.9 при использовании передискретизации частотной области.
ОПИСАНИЕ ПРЕДПОЧТИТЕЛЬНЫХ ВАРИАНТОВ ВОПЛОЩЕНИЯ
Описанные ниже варианты воплощения приводятся с единственной целью, а именно для наглядной иллюстрации принципов настоящего изобретения по передискретизации в комбинированном банке фильтров транспозитора. При этом понимается, что модификации и вариации механизмов и подробностей, описанных в настоящем документе, являются очевидными для специалистов данной отрасли техники. Следовательно, ограничительным является только объем прилагаемой формулы изобретения, а не конкретные подробности, упомянутые в настоящем документе в ходе описания и объяснения вариантов воплощения.
Фиг.1 представляет принцип действия гармонического транспозитора частотной области 100. В своей основе, гармонический транспозитор T-го порядка представляет собой модуль, который сдвигает все сигнальные компоненты H(f) входного сигнала, т.е. субполосу сигнала в частотной области, в H(Tf). To есть, частотная компонента H(f) входного сигнала сдвигается в сторону в Т раз более высокой частоты. Для осуществления такой транспозиции в частотной области, банк фильтров анализа 101 преобразует входной сигнал из временной области в частотную область и формирует выходной сигнал, представляющий собой комплексные субполосы или сигналы субполос, также называемые субполосами анализа или сигналами субполос анализа. Банк фильтров анализа, как правило, включает в себя преобразование анализа, например БПФ, ДПФ или вейвлет-преобразование, и скользящее окно анализа. Сигналы субполос анализа подвергаются нелинейной обработке 102, в процессе которой изменяется фаза и/или амплитуда в соответствии с выбранным порядком транспозиции Т. Как правило, результатом нелинейной обработки является некоторое число сигналов субполос, равное числу сигналов входных субполос, т.е. числу сигналов субполос анализа. Модифицированные субполосы или сигналы субполос, также называемые субполосами синтеза или сигналами субполос синтеза, подаются на вход банка фильтров синтеза 103, который преобразует сигналы субполос из частотной области во временную, и его выходной сигнал представляет собой транспонированный сигнал временной области. Банк фильтров синтеза 103, как правило, включает в себя обратное преобразование, например обратное БПФ, обратное ДПФ или обратное вейвлет-преобразование, и скользящее окно синтеза.
Как правило, каждый банк фильтров имеет физическое разрешение по частоте Δf, измеряемое в герцах, и физический параметр временного шага Δt, измеряемый в секундах, причем физическое разрешение по частоте Δf обычно ассоциируется с разрешением по частоте функции преобразования, а физический параметр временного шага Δt обычно ассоциируется с временным интервалом между последовательными оконными функциями. Эти два параметра, т.е. разрешение по частоте и временной шаг, определяют дискретно-временные параметры банка фильтров при заданной частоте дискретизации. Выбирая физические параметры временного шага, т.е. параметр временного шага, измеренный в единицах времени, например секундах, для банков фильтров анализа и синтеза идентичными, на выходе транспозитора 100 можно получить выходной сигнал, имеющий ту же частоту дискретизации, что и входной сигнал. Более того, при исключении этапа нелинейной обработки 102, на выходе может быть получено полное восстановление входного сигнала. Для этого необходим тщательный расчет банков фильтров анализа и синтеза. С другой стороны, в случае, когда на выходе и входе выбраны различные частоты дискретизации, может быть получено преобразование частоты дискретизации. Такой режим работы может быть необходим в случае, когда требуемая ширина полосы частот выходного сигнала y больше, чем половина частоты дискретизации входного сигнала x, т.е. когда требуемая ширина полосы частот выходного сигнала превосходит частоту Найквиста для входного сигнала.
Фиг.2 представляет принцип работы параллельного транспозитора или системы параллельных транспозиторов 200, содержащей некоторое число гармонических транспозиторов 201-1, …, 201-Р различного порядка. Входной транспонируемый сигнал пропускается через банк из Р отдельных транспозиторов 201-1, 201-2, …, 201-Р. Отдельные транспозиторы 201-1, 201-2, …, 201-Р выполняют гармоническую транспозицию входного сигнала в соответствии с пояснением к Фиг.1. Как правило, каждый отдельный транспозитор 201-1, 201-2, …, 201-Р выполняет гармоническую транспозицию различного порядка транспозиции Т. Например, транспозитор 201-1 может выполнять транспозицию порядка Т=1, транспозитор 201-2 может выполнять транспозицию порядка T=2, … и транспозитор 201-Р может выполнять транспозицию порядка T=P. Однако, в общем случае, любой из транспозиторов 201-1, …, 201-Р может выполнять гармоническую транспозицию произвольного порядка транспозиции Т. Вклады, т.е. выходные сигналы отдельных транспозиторов 201-1, 201-2, …, 201-Р, могут суммироваться в сумматоре 202 для формирования суммарного выходного сигнала транспозитора.
Требуется отметить, что для каждого транспозитора 201-1, 201-2, …, 201-Р необходим банк фильтров анализа и синтеза, как представлено на Фиг.1. Более того, в типовой реализации отдельные транспозиторы 201-1, 201-2, …, 201-Р, как правило, изменяют частоту дискретизации обрабатываемого входного сигнала на различную величину. Например, частота дискретизации выходного сигнала транспозитора 201-Р может быть в Т раз выше, чем частота дискретизации входного сигнала транспозитора 201-Р, причем T является порядком транспозиции, осуществляемой транспозитором 201-Р. Это может происходить благодаря коэффициенту расширения полосы пропускания Т, примененного в транспозиторе 201-Р, т.е. благодаря использованию банка фильтров синтеза, имеющего в Т раз больше подканалов, чем банк фильтров анализа. В результате частота дискретизации и частота Найквиста возрастает с коэффициентом T. В результате чего, для отдельных сигналов временной области может потребоваться редискретизация с целью создания возможности для суммирования различных выходных сигналов в сумматоре 202. Редискретизация сигналов частотной области может быть произведена на входе или на выходе каждого независимого транспозитора 201-1, 201-2, …, 201-Р.
Фиг.3 представляет образец конфигурации параллельного гармонического транспозитора или системы параллельных транспозиторов 300, выполняющей несколько порядков транспозиции и использующей общий банк фильтров анализа 301. Отправным пунктом для расчета параллельного транспозитора 300 может быть такой расчет отдельных транспозиторов 201-1, 201-2, …, 201-Р, представленных на Фиг.2, при котором банки фильтров анализа (ссылочная позиция 101 на Фиг.1) всех транспозиторов 201-1, 201-2, …, 201-Р являются идентичными и могут быть заменены единым банком фильтров анализа 301. Вследствие этого, сигналы временной области преобразовываются в один набор сигналов субполос частотной области, т.е. один набор сигналов субполос анализа. Эти сигналы субполос поступают на входы различных модулей нелинейной обработки 302-1, 302-2, …, 302-Р для различных порядков транспозиции. Как изложено выше в контексте Фиг.1, каждый модуль нелинейной обработки производит изменение фазы и/или амплитуды сигналов субполос, и для различных порядков транспозиции эти изменения различны. В дальнейшем, различным образом модифицированные сигналы субполос или субполосы поступают на входы различных банков фильтров синтеза 303-1, 303-2, …, 303-Р, соответствующие различным модулям нелинейной обработки 302-1, 302-2, …, 302-Р. Результатом этого являются Р по-разному транспонированные выходные сигналы временной области, которые суммируются в сумматоре 304, образуя суммарный выходной сигнал транспозитора.
Следует отметить, что если банки фильтров синтеза 303-1, 303-2, …, 303-Р, соответствующие различным порядкам транспозиции, функционируют с различными частотами дискретизации, например, используя различные степени расширения полосы частот, то для выходных сигналов временной области различных банков фильтров синтеза 303-1, 303-2, …, 303-Р потребуется различная редискретизация для привязки Р выходных сигналов к общей временной сетке перед суммированием в сумматоре 304.
Фиг.4 представляет принцип работы параллельного гармонического транспозитора 400, использующего несколько порядков транспозиции, использующего при этом общий банк фильтров синтеза 404. Исходным пунктом расчета такого параллельного транспозитора 400 может быть такой расчет отдельных транспозиторов 201-1, 201-2, …, 201-Р, представленных на Фиг.2, при котором банки фильтров синтеза всех транспозиторов являются идентичными и могут быть заменены единым банком фильтров синтеза 404. Следует отметить, что аналогично ситуации, представленной на Фиг.3, здесь используются различные модули нелинейной обработки 402-1, 402-2, …, 402-Р для каждого порядка транспозиции. Более того, для различных порядков транспозиции используются различные банки фильтров анализа 401-1, 401-2, …, 401-Р. В связи с этим, набор из Р банков фильтров анализа 401-1, 401-2, …, 401-Р определяет Р наборов сигналов субполос анализа. Эти Р наборов сигналов субполос анализа поступают на входы соответствующих модулей нелинейной обработки 402-1, 402-2, …, 402-Р для формирования Р наборов модифицированных сигналов субполос. Эти Р наборов сигналов субполос могут быть просуммированы в частотной области в сумматоре 403 для формирования суммарного набора сигналов субполос, подаваемого на вход единого банка фильтров синтеза 404. Такое суммирование в сумматоре 403 может включать в себя подачу различным образом обработанных сигналов субполос в различные интервалы субполос и/или суперпозицию вкладов сигналов субполос на перекрывающиеся интервалы субполос. Другими словами, различные сигналы субполос анализа, обработанные путем транспозиции различных порядков, могут занимать перекрывающиеся интервалы частот. Например, транспозитор второго порядка может транспонировать субполосу анализа [2А, 2В] в интервал субполос [4А, 4В]. В то же время, транспозитор четвертого порядка может транспонировать субполосу анализа [А, В] в тот же интервал субполос [4А, 4В]. В таких случаях, накладывающиеся вклады могут суммироваться, т.е. складываться и/или усредняться, сумматором 403. Выходной сигнал временной области параллельного транспозитора 400 получают на выходе общего банка фильтров синтеза 404. Аналогично вышеизложенному, в случае, когда банки фильтров анализа 401-1, 401-2, …, 401-Р работают с различными частотами дискретизации, для входных сигналов временной области различных банков фильтров анализа 401-1, 401-2, …, 401-Р может понадобиться редискретизация для привязки выходных сигналов различных модулей нелинейной обработки 402-1, 402-2, …, 402-Р к той же временной сетке.
Фиг.5 представляет принцип действия гармонического параллельного транспозитора 500, использующего несколько порядков транспозиции, в состав которого входит единый общий банк фильтров анализа 501 и единый общий банк фильтров синтеза 504. В этом случае, отдельные транспозиторы 201-1, 201-2, …, 201-Р, представленные на Фиг.2, должны быть рассчитаны таким образом, чтобы как банки фильтров анализа, так и банки фильтров синтеза всех Р гармонических транспозиторов были идентичны. При соблюдении условия идентичности банков фильтров анализа и синтеза для Р различных гармонических транспозиторов, идентичные банки фильтров могут быть заменены единым банком фильтров анализа 501 и единым банком фильтров синтеза 504. Модули усовершенствованной нелинейной обработки 502-1, 502-2, …, 502-Р формируют на выходах различные вклады в частично перекрывающиеся интервалы частот, которые суммируются в сумматоре 503, образуя суммарный входной сигнал для соответствующих субполос банка фильтров синтеза 504. Подобно параллельному гармоническому транспозитору 400, представленному на Фиг.4, суммирование в сумматоре 503 может включать в себя подачу различных выходных сигналов совокупности модулей нелинейной обработки 502-1, 502-2, …, 502-Р в различные интервалы субполос и суперпозицию вкладов параллельных выходных сигналов на перекрывающиеся интервалы субполос.
Как уже показано выше, на выходе процесса нелинейной обработки 102, как правило, формируется некоторое число субполос, соответствующее числу субполос на входе. Нелинейная обработка 102, как правило, изменяет фазу и/или амплитуду субполосы или сигнала субполосы в соответствии с основным порядком транспозиции Т. Например, субполоса на входе преобразовывается в субполосу на выходе, имеющую в Т раз большую частоту, т.е. субполоса на входе нелинейной обработки 102, т.е. субполоса анализа,
может быть транспонирована в субполосу на выходе процесса нелинейной обработки 102, т.е. субполосу синтеза,
, где k является индексом субполосы и Δf является разрешением по частоте банка фильтров анализа. Для того, чтобы использовать общие банки фильтров анализа 501 и общие банки фильтров синтеза 504, один или более модулей усовершенствованной обработки 502-1, 502-2, …, 502-Р могут быть сконфигурированы с возможностью формировать некоторое число субполос выхода, которое может отличаться от числа субполос входа.
Далее будут изложены принципы усовершенствованной нелинейной обработки в модулях нелинейной обработки 502-1, 502-2, …, 502-Р. Для этой цели предполагается, что
- банк фильтров синтеза и банк фильтров анализа имеют одинаковый физический параметр временного шага Δt,
- банк фильтров анализа имеет физическое разрешение по частоте Δf,
- банк фильтров синтеза имеет физическое разрешение по частоте QΔf, где коэффициент разрешения Q≥1 принимает целочисленные значения.
Кроме того, предполагается, что эти банки фильтров набраны четным образом, т.е. центром субполосы с нулевым индексом является нулевая частота, так что центральные частоты банка фильтров анализа имеют вид kΔf, где индекс субполосы анализа k=1, …, KA-1 и KA является числом субполос банка фильтров анализа. Центральные частоты банка фильтров синтеза имеют вид kQΔf, где n=1, …, NS-1 - индекс субполосы синтеза и NS - число субполос банка фильтров синтеза.
При выполнении традиционной транспозиции целого порядка Т≥1, как представлено на Фиг.1, коэффициент разрешения Q выбирается как Q=T, и нелинейно обработанная субполоса анализа k отображается в субполосу синтеза с таким же индексом n=k. Нелинейная обработка 102, как правило, включает в себя умножение фазы субполосы или сигнала субполосы на коэффициент Т. То есть, для каждого отсчета в субполосах банка фильтров можно записать
где θA(k) - фаза (комплексного) отсчета в субполосе анализа k и θS(k) - фаза (комплексного) отсчета в субполосе синтеза k. Величина или амплитуда отсчета субполосы может оставаться неизменной или же может увеличиваться или изменяться на постоянный коэффициент усиления. Вследствие того, что Т - целое число, операция уравнения (1) не зависит от определения фазового угла.
В традиционных параллельных транспозиторах коэффициент разрешения Q банка фильтров анализа/синтеза выбирается равным порядку транспозиции Т соответствующего транспозитора, т.е. Q=Т. В этом случае, разрешение по частоте банка фильтров синтеза равно TΔf и, следовательно, зависит от порядка транспозиции T. Следовательно, для различных порядков транспозиции Т необходимо использовать различные банки фильтров как на этапе анализа, так и на этапе синтеза. Это следует из того факта, что порядок транспозиции T определяет отношение физических разрешений по частоте, т.е. отношение физического разрешения по частоте Δf банка фильтров анализа и разрешения по частоте TΔf банка фильтров синтеза.
Для того, чтобы использовать общий банк фильтров анализа 501 и общий банк фильтров синтеза 504 для совокупности различных порядков транспозиции T, предлагается установить разрешения по частоте банка фильтров синтеза 504 равным QΔf, т.е. предлагается сделать разрешения по частоте банка фильтров синтеза 504 независимым от порядка транспозиции T. Тогда возникает вопрос о том, как выполнить транспозицию порядка Т в случае, когда коэффициент разрешения Q, т.е. коэффициент Q физического разрешения по частоте банков фильтров синтеза и анализа, не всегда удовлетворяет условию Q=Т.
Как изложено выше, принцип гармонической транспозиции состоит в том, что входной сигнал n-ой субполосы банка фильтров синтеза с центральной частотой nQΔf определяется исходя из субполосы анализа, имеющего в Т раз более низкую центральную частоту, т.е. центральную частоту, равную nQΔf/T. Центральные частоты субполос анализа выражаются с помощью индекса субполосы анализа k как kΔf. Оба выражения для центральной частоты индекса субполос анализа, т.е. nQΔf/T и kΔf, могут быть приравнены друг к другу. С учетом того, что индекс n принимает целочисленные значения, выражение
является рациональным числом, которое может быть выражено в виде суммы целочисленного индекса субполосы анализа k и остатка r∈{0,1/Т, 2/Т, …, (Т-1)/Т}, так что
В связи с этим, можно наложить условие, что входной сигнал субполосы синтеза с индексом субполосы синтеза n может быть получен при помощи транспозиции Т-го порядка из субполосы анализа с индексом k, задаваемого уравнением (2). В виду того факта, что
является рациональным числом, остаток r может быть неравен нулю, и значение k+r может быть больше, чем индекс субполосы анализа k и меньше, чем индекс субполосы анализа k+1, т.е. k≤k+r≤k+1. Следовательно, входной сигнал субполосы синтеза с индексом субполос синтеза n следует получать, используя транспозицию T-го порядка из субполос анализа с индексами субполос анализа k и k+1, где k задается уравнением (2). Другими словами, входной сигнал субполосы синтеза может быть получен из двух последовательных субполос анализа.
В результате приведенного выше рассуждения, процесс усовершенствованной нелинейной обработки, производимой модулем нелинейной обработки 502-1, 502-2, …, 502-Р, может включать в себя этап рассмотрения двух соседних субполос анализа с индексами k и k+1 с целью формирования выходного сигнала для субполосы синтеза n. Для транспозиции порядка Т изменение фазы, выполняемое модулем нелинейной обработки 502-1, 502-2, …, 502-Р, может, например, определяться по правилу линейной интерполяции,
где θA(k) - фаза отсчета в субполосе анализа k, QA(k+1) - фаза отсчета в субполосе анализа k+1 и θS(n) - фаза отсчета в субполосе синтеза n. Если остаток r близок к нулю, т.е., если значение k+r близко к k, то основной вклад в фазу отсчета субполосы синтеза определяется фазой отсчета в субполосе анализа k. С другой стороны, если остаток r близок к единице, т.е. если значение k+r близко к k+1, то основной вклад в фазу отсчета в субполосе синтеза определяется фазой отсчета в субполосе анализа k+1. Следует отметить, что фазовые множители Т(1-r) и Tr имеют такие целочисленные значения, что фазовые изменения уравнения (3) строго определены и независимы от определения фазового угла.
Что касается величины отсчетов субполос, то для определения величины отсчетов в субполосах синтеза может быть выбрано следующее среднее геометрическое значение,
где aS(n) обозначает величину отсчета в субполосе синтеза n, aA(k) обозначает величину отсчета в субполосе анализа k и aA(k+1) обозначает величину отсчета в субполосе анализа k+1. Следует отметить, что могут предполагаться и иные правила интерполяции для фазы и/или величины.
В случае, когда банк фильтров набран нечетным образом, и центральные частоты банка фильтров анализа имеют вид
, где k=1, …, KA-1, и центральные частоты банка фильтров синтеза имеют вид
, где n=1, …, NS-1, уравнение, соответствующее уравнению (2), может быть получено путем приравнивания транспонированной центральной частоты банка фильтров синтеза
к центральной частоте банка фильтров анализа
. Полагая индекс k целочисленным и остаток r∈[0,1[, для банков фильтров, набранных нечетным образом, может быть выведено следующее уравнение:
Для специалиста ясно, что если Т-Q, т.е. разность между порядком транспозиции и коэффициентом разрешения, четная, то T(1-r) и Tr выражаются целыми числами, и могут быть применены правила интерполяции, выраженные уравнениями (3) и (4).
Отображение субполос анализа в субполосы синтеза представлено на Фиг.5b. Фиг.5b представляет четыре диаграммы для различных порядков транспозиции от Т=1 до Т=4. Каждая диаграмма иллюстрирует, как исходные элементы дискретизации 510, т.е. субполосы анализа, отображаются в целевые элементы дискретизации 530, т.е. субполосы синтеза. Для упрощения иллюстрации предполагается, что коэффициент разрешения Q равен единице. Другими словами, Фиг.5b иллюстрирует отображение сигналов субполос анализа в сигналы субполос синтеза согласно уравнениям (2) и (3). В приведенном примере банк фильтров анализа/синтеза набран четным образом, с Q=1 и максимальным порядком транспозиции Т=4.
В данном случае уравнение (2) может быть записано как
. Следовательно, для порядка транспозиции T=1, субполоса анализа с индексом k отображается в соответствующую субполосу синтеза n, и остаток r всегда равен нулю. Это представлено на Фиг.5b, где, например, исходный элемент дискретизации 511 взаимно-однозначно отображается в целевой элемент дискретизации 531.
В случае, когда порядок транспозиции Т=2, остаток r принимает значения 0 и ½, и исходный элемент дискретизации отображается в несколько целевых элементов дискретизации. Другими словами, можно утверждать, что каждый целевой элемент дискретизации 532, 535 принимает вклад от не более чем двух исходных элементов дискретизации. Это представлено на Фиг.5b, где целевой элемент дискретизации 535 принимает вклад от исходных элементов дискретизации 512 и 515. Однако, целевой элемент дискретизации 532 принимает вклад только от исходного элемента дискретизации 512. Если предполагается, что целевой элемент дискретизации 532 имеет четный индекс n, например, n=10, то из уравнения (2) следует, что целевой элемент дискретизации 532 принимает вклад от исходного элемента дискретизации 512 с индексом k=n/2, например k=5. Остаток r равен нулю, т.е. при этом отсутствует вклад от исходного элемента дискретизации 515 с индексом k+1, например k+1=6. Это не верно для целевого элемента дискретизации 535 с нечетным индексом n, например, n=11. В этом случае, согласно уравнению (2), целевой элемент дискретизации 535 принимает вклады от исходного элемента дискретизации 512 (индекс k=5) и исходного элемента дискретизации 515 (индекс k+1=6). Аналогичным образом это относится и к более высоким порядкам транспозиции Т, например Т=3 и T=4, как представлено на Фиг.5b.
Дальнейшее объяснение процесса вышеописанной усовершенствованной нелинейной обработки может быть следующим. Под усовершенствованной нелинейной обработкой может подразумеваться комбинация транспозиции заданного порядка Т в сигналы промежуточных субполос сетки промежуточных частот TΔf, и последующее отображение сигналов промежуточных субполос в сетку частот, определяемую общим банком фильтров синтеза, т.е. сетку частот QΔf. Для иллюстрации этого объяснения снова потребуется обращение к Фиг.5b. Однако, для такой иллюстрации исходные элементы дискретизации 510 рассматриваются как промежуточные субполосы, которые получаются из субполос анализа, при использовании порядка транспозиции Т. Сетка частот этих промежуточных субполос имеет вид ТΔf. Для того чтобы создать сигналы субполос синтеза в заданной сетке частот QΔf определяемой целевыми элементами дискретизации 530, исходные элементы дискретизации 510, т.е. промежуточные субполосы с сеткой частот TΔf, должны быть отображены в заданную сетку частот QΔf. Это может быть выполнено путем определения целевого элемента дискретизации 530, т.е. сигнала субполосы синтеза в сетке частот QΔf, путем интерполяции одного или двух элементов дискретизации 510, т.е. сигналов промежуточных субполос в сетке частот TΔf. В предпочтительном варианте воплощения используется линейная интерполяция, причем весовые коэффициенты интерполяции обратно пропорциональны разности между центральной частотой целевого элемента дискретизации 530 и соответствующего исходного элемента дискретизации 510. Например, если эта разность равна нулю, то весовой коэффициент равен 1, и если эта разность равна TΔf, то весовой коэффициент равен 0.
Таким образом, описан способ нелинейной обработки, который позволяет производить определение вкладов в субполосу синтеза путем транспозиции нескольких субполос анализа. Данный способ нелинейной обработки допускает использование единых общих банков фильтров субполос синтеза и анализа для различных порядков транспозиции, благодаря чему значительно снижается вычислительная сложность параллельных гармонических транспозиторов.
Фиг.6а и 6b иллюстрируют образец банков фильтров анализа/синтеза, использующих М=1024 точечное БПФ/ДПФ (быстрое преобразование Фурье или дискретное преобразование Фурье) для параллельной транспозиции порядков Т=2, 3, 4. Фиг.6а иллюстрирует пример традиционного параллельного гармонического транспозитора 600, использующего общий банк фильтров анализа 601 и отдельные банки фильтров синтеза 602, 603, 604 для каждого порядка транспозиции Т=2, 3, 4. Фиг.6а представляет окно анализа νA 611 и окна синтеза νS 612, 613, 614, примененные в банке фильтров анализа 601 и банках фильтров синтеза 602, 603, 604, соответственно. В иллюстрируемом примере, окно анализа νA 611 имеет длину LA=1024, которая равна размеру М БПФ или ДПФ банков фильтров анализа/синтеза 601, 602, 603, 604. Аналогичным образом, окна синтеза νS 612, 613, 614 имеют длину LS=1024, которая равна размеру М БПФ или ДПФ.
Фиг.6а также иллюстрирует величину сдвига ΔsA, используемую банком фильтров анализа 601, и величину сдвига ΔsS, используемую банками фильтров синтеза 602, 603, 604, соответственно. Величина сдвига Δs соответствует числу дискретных значений сигнала, на которое соответствующие окна 611, 612, 613, 614 перемещаются между последовательными этапами преобразования. Величина сдвига Δs связана с физическим временным шагом Δt через частоту дискретизации основного сигнала, т.е. Δs=fSΔt, где fS - частота дискретизации.
Как показано, окно анализа 611 перемещается на величину сдвига 621 из 128 отсчетов. Окно синтеза 612, соответствующее транспозиции порядка T=2, перемещается на величину сдвига 622 из 256 отсчетов, т.е. величину сдвига 622, которая в два раза больше величины сдвига 621 окна анализа 611. Как изложено выше, это приводит к изменению временного масштаба сигнала с коэффициентом T=2. Кроме того, если предполагается в T=2 раза более высокая частота дискретизации, разница между величиной сдвига анализа 621 и величиной сдвига синтеза 622 приводит к гармонической транспозиции порядка Т=2. Т.е. изменение временного масштаба на порядок T может быть преобразовано в гармоническую транспозицию путем выполнения преобразования частоты дискретизации порядка Т.
Аналогичным образом, можно показать, что величина сдвига синтеза 623, соответствующая гармоническому транспозитору порядка T=3, в T=3 раз больше, чем величина сдвига анализа 621, а величина сдвига синтеза 624, соответствующая гармоническому транспозитору порядка T=4, в T=4 раза больше, чем величина сдвига анализа 621. Для привязки частот дискретизации транспозитора 3го порядка и транспозитора 4го порядка к выходной частоте дискретизации транспозитора 2го порядка, в состав транспозитора 3го порядка и транспозитора 4го порядка включены дециматор 633 с коэффициентом 3/2 и дециматор 634 с коэффициентом 2, соответственно. В общем случае, транспозитор Тго порядка должен содержать дециматор с коэффициентом T/2, если требуемая частота дискретизации выходного сигнала в 2 раза выше частоты дискретизации входного сигнала. Т.е. для гармонического транспозитора порядка T=2 децимация не требуется.
Наконец, Фиг.6а иллюстрирует отдельные модули изменения фазы 642, 643, 644 для транспозиции порядка T=2, 3, 4, соответственно. Такие модули изменения фазы 642, 643, 644 производят умножение фазы соответствующих сигналов субполос на порядок транспозиции T=2, 3, 4, соответственно (см. уравнение (1)).
Эффективная структура комбинированного банка фильтров для транспозитора может быть получена путем ограничения параллельного транспозитора, представленного на фиг.6а, единым банком фильтров анализа 601 и единым банком фильтров синтеза 602. Гармоники 3го и 4го порядка в этом случае создаются модулем нелинейной обработки 650 в банке фильтров 2го порядка, как представлено на Фиг.6b. Фиг.6b представляет банк фильтров анализа, содержащий модуль 1024 точечного прямого БПФ 601 и окно анализа 611, которое применяется к входному сигналу x, с величиной сдвига анализа 621. Банк фильтров синтеза содержит модуль 1024 точечного обратного БПФ 602 и окно синтеза 612, которое применяется с величиной сдвига синтеза 622. В приведенном примере, величина сдвига синтеза 622 в два раза больше величины сдвига анализа 621. Более того предполагается, что частота дискретизации выходного сигнала y в два раза больше частоты дискретизации входного сигнала x.
Банк фильтров анализа/синтеза, представленный на Фиг.6b, содержит единый банк фильтров анализа и единый банк фильтров синтеза. При использовании усовершенствованной нелинейной обработки 650 в соответствии со способами, изложенными в контексте Фиг.5 и Фиг.5b, т.е. усовершенствованную нелинейную обработку, выполняемую модулями 502-1, …, 502-Р, такой банк фильтров анализа/синтеза может быть использован для создания параллельного транспозитора, т.е. гармонического транспозитора для совокупности порядков транспозиции Т.
Как было изложено в контексте Фиг.5 и 5b, взаимно однозначное отображение субполос анализа на соответствующие субполосы синтеза, включающее в себя умножение фазы сигналов субполос на соответствующий порядок транспозиции Т, может быть обобщено на правила интерполяции (см. Уравнения (3) и (4)), включающие один или более сигналов субполос. Было показано, что если физический интервал QΔf субполос банка фильтров синтеза в Q раз больше физического интервала Δf банка фильтров анализа, то входной сигнал полосы синтеза с индексом n получается из полос анализа с индексами k и k+1. Взаимосвязь между индексами n и k описывается Уравнением (2) или (5), в зависимости от того, четным или нечетным образом набраны банки фильтров. Геометрическая интерполяция для величин производится с показателями степеней 1-r и r (Уравнение (4)), и линейная комбинация фаз производится с весовыми коэффициентами T(1-r) и Tr (уравнение (3)). В описываемом примере, где Q=2, отображения фазы для каждого коэффициента транспозиции графически представлены на Фиг.7.
В порядке, аналогичном случаю для Q=1, представленному на Фиг.5, целевая субполоса или целевой элемент дискретизации 730 принимает вклады от не более чем двух субполос или исходных элементов дискретизации 710. В случае, когда T=Q=2, для каждого исходного элемента дискретизации с измененной фазой 711 определяется соответствующий целевой элемент дискретизации 731. Для транспозиции более высоких порядков, когда T>Q, целевой элемент дискретизации 735 может быть получен из одного соответствующего исходного элемента дискретизации с измененной фазой 715. Это имеет место в случае, когда остаток r, полученный из Уравнения (2) или (5), равен нулю. В противном случае, целевой элемент дискретизации 732 получается путем интерполяции двух исходных элементов дискретизации с измененной фазой 712 и 715.
Вышеописанная нелинейная обработка выполняется модулем параллельного транспозитора 650, который определяет целевые элементы дискретизации 730 для различных порядков транспозиции T=2, 3, 4, используя модули усовершенствованной нелинейной обработки 502-2, 502-3, 502-4. Впоследствии соответствующие целевые элементы дискретизации 730 суммируются модулем сумматора 503, в результате чего получается единый набор сигналов субполос синтеза, поступающий на вход банка фильтров синтеза. Как изложено выше, модуль сумматора 503 сконфигурирован с возможностью производить суммирование совокупности вкладов в перекрывающиеся частотные интервалы, поступающие от различных модулей нелинейной обработки 502-2, 502-3, 502-4.
Далее описывается гармоническая транспозиция импульсных сигналов с помощью гармонических транспозиторов. В данном контексте требуется отметить, что гармоническая транспозиция порядка Т, использующая банки фильтров анализа/синтеза, может быть представлена как изменение временного масштаба базового сигнала на целочисленный коэффициент транспозиции Т с последующей децимацией и/или преобразованием частоты дискретизации. Изменение временного масштаба выполняется таким образом, что частоты синусоидальных сигналов, составляющих входной сигнал, остаются неизменными. Такое изменение временного масштаба может быть произведено при помощи банка фильтров анализа/синтеза в сочетании с промежуточным изменением фаз сигналов субполос в зависимости от порядка транспозиции Т. Как изложено выше, банк фильтров анализа может представлять собой банк фильтров, использующий оконное ДПФ, с окном анализа νA, а банк фильтров синтеза может представлять собой банк фильтров, использующий оконное ДПФ, с окном синтеза νS. Такое преобразование анализа/синтеза может также называться короткое преобразование Фурье (КПФ).
Входной сигнал временной области x подвергается короткому преобразованию Фурье для получения последовательности перекрывающихся спектральных фреймов. Для минимизации возможного влияния соседних полос необходим правильный выбор соответствующих окон анализа/синтеза, например, окон Гаусса, косинус-окон, окон Хэмминга, окон Ханна, прямоугольных окон, окон Бартлетта, окон Блэкмана и других. Временная задержка, с которой каждый спектральный фрейм выбирается из входного сигнала x, называется величина сдвига Δs или физический временной шаг Δt. Результатом КПФ входного сигнала x, называемого этапом анализа, является представление входного сигнала x в частотной области. Представление в частотной области включает в себя совокупность сигналов субполос, в которой каждый сигнал субполосы представляет собой определенную частотную компоненту входного сигнала.
Для изменения временного масштаба входного сигнала может быть изменен временной масштаб каждого сигнала субполосы, например, путем задержки отсчетов сигнала субполосы. Это может быть достигнуто при использовании большей величины сдвига синтеза, чем величина сдвига анализа. Сигнал временной области может быть восстановлен при помощи обратного (быстрого) преобразования Фурье с дальнейшим последовательным накоплением фреймов. Такая операция стадии синтеза называется операцией перекрытия с накоплением. Результирующий выходной сигнал является разновидностью входного сигнала с измененным масштабом времени, содержащей те же частотные компоненты, что и входной сигнал. Другими словами, результирующий выходной сигнал имеет тот же спектральный состав, что и входной сигнал, но он медленнее, чем входной сигнал, т.е. его изменение растянуто во времени.
Транспозиция к более высоким частотам может быть произведена на последующих этапах, или же интегрированным образом, путем децимации сигналов с измененным масштабом времени или путем преобразования частоты дискретизации выходного сигнала с измененным масштабом времени. В результате транспонированный сигнал имеет временную длительность исходного сигнала, но содержит частотные компоненты, которые сдвинуты вверх с заданным коэффициентом транспозиции.
В виду вышесказанного гармоническая транспозиция импульсных сигналов с помощью гармонических транспозиторов будет описана исходя из изменения временного масштаба эталонного импульсного сигнала, т.е. дискретного импульса Дирака в момент времени t=t0,
Преобразование Фурье для такого импульса Дирака имеет величину единицы и линейную фазу с угловым коэффициентом, пропорциональным t0:
где
- центральная частота сигнала mой субполосы КПФ анализа, и М размер дискретного преобразования Фурье (ДПФ). Такое преобразование Фурье может рассматриваться как этап анализа вышеописанного банка фильтров анализа, в котором использовано плоское окно анализа νA бесконечной длительности. Для создания выходного сигнала y, масштаб времени которого изменен с коэффициентом Т, т.е. импульса Дирака δ(t-Tt0) в момент времени t=Tt0, фаза сигналов субполос анализа должна быть умножена на коэффициент T для получения такого сигнала субполосы синтеза Y(Ωm)=ехр(-jΩmTt0), который обеспечит требуемый импульс Дирака δ(t-Tt0) в качестве результата обратного преобразования Фурье.
Однако следует отметить, что вышеприведенные соображения относятся к этапу анализа/синтеза, в котором используются окна анализа и синтеза бесконечной длины. Действительно, теоретический транспозитор с окном бесконечной длительности в состоянии выполнить точное изменение масштаба для импульса Дирака δ(t-t0). Для анализа, использующего окно конечной длительности, ситуация осложняется тем фактом, что каждый блок анализа должен быть интерпретирован как один периодический интервал периодического сигнала с периодом, равным размеру ДПФ.
Это иллюстрирует Фиг.8, представляющая анализ и синтез 800 импульса Дирака δ(t-t0). Верхняя часть Фиг.8 представляет входной сигнал этапа анализа 810, а нижняя часть Фиг.8 представляет выходной сигнал этапа синтеза 820. Верхняя и нижняя диаграммы представляют временную область. Стилизованное окно анализа 811 и окно синтеза 821 изображены как треугольные окна (Бартлетта). Импульс входного сигнала δ(t-t0) 812 в момент времени t=t0 изображен на верхней диаграмме 810 вертикальной стрелкой. Предполагается, что блок преобразования ДПФ имеет размер М=L=LA=LS, т.е. размер ДПФ выбран равным величине окна. Умножение фазы сигналов субполос на коэффициент T произведет ДПФ анализ импульса Дирака δ(t-Tt0) в момент времени t=Tt0, однако, в виде последовательности импульсов Дирака с периодом L. Это является следствием конечной длины применяемого окна и преобразования Фурье. Эта последовательность периодических импульсов с периодом L изображена пунктирными стрелками 823, 824 на нижней диаграмме.
В реальных системах последовательность импульсов фактически состоит всего лишь из нескольких импульсов (в зависимости от коэффициента транспозиции) - одного основного импульса, т.е. искомого члена, и нескольких предшествующих и последующих импульсов, т.е. нежелательных членов. Предшествующие и последующие импульсы возникают вследствие периодической природы ДПФ (с периодом L). В случае, когда импульс расположен в пределах окна анализа так, что при умножении на T происходит перенос комплексного значения фазы (т.е. импульс сдвигается в сторону конца окна и переносится к его началу), нежелательный импульс возникает внутри окна синтеза. Нежелательные импульсы могут иметь (необязательно) ту же полярность, что и входные импульсы, в зависимости от положения окна анализа и коэффициента транспозиции.
В примере на Фиг.8, в процессе оконного синтеза используется конечное окно νS 821. Конечное окно синтеза 821 выбирает требуемый импульс δ(t-Tt0) в момент времени t=Tt0, обозначенный непрерывной стрелкой 822, и исключает прочие нежелательные вклады, изображенные пунктирными стрелками 823, 824.
По мере продвижения этапов анализа и синтеза по временной оси в соответствии с коэффициентом сдвига Δs или временным шагом Δt, импульс δ(t-t0)812 будет иметь иное положение относительно центра соответствующего окна анализа 811. Как было изложено выше, операция по изменению временного масштаба заключается в перемещении импульса 812 в положение, соответствующее T его положениям относительно центра окна. Если это положение находится в пределах окна 821, то в процессе операции изменения временного масштаба все вклады гарантированно добавятся к единому импульсу с измененным масштабом времени δ(t-Tt0) в момент времени t=Tt0.
Однако, для ситуации, представленной на Фиг.9, проблема возникает в том случае, когда импульс δ(t-t0)912 перемещается за пределы блока ДПФ. Фиг.9 иллюстрирует ту же конфигурацию анализа/синтеза 900, что и Фиг.8. На верхней диаграмме 910 показан входной сигнал этапа анализа и окно анализа 911, а нижняя диаграмма 920 представляет выходной сигнал этапа синтеза и окно синтеза 921. При изменении временного масштаба входного импульса Дирака 912 с коэффициентом Т, импульс Дирака с измененным временным масштабом 922, т.е. δ(t-Tt0), оказывается вне окна синтеза 921. В то же время, другой импульс Дирака 924 из последовательности импульсов, т.е. δ(t-Tt0+L) в момент времени t=Tt0-L, выбирается окном синтеза. Другими словами, импульс Дирака 912 не задерживается до момента времени, в T раз более позднего, а передвигается вперед к моменту времени, предшествующему импульсу Дирака на входе 912. В результате в звуковом сигнале появляется опережающее эхо на временном промежутке, характерном для достаточно длительных окон транспозитора, т.е. в момент времени t=Tt0-L, который на величину L-(T-1)t0 опережает входной импульс Дирака 912.
Принцип решения этой задачи будет описан с помощью Фиг.10. Фиг.10 иллюстрирует сценарий анализа/синтеза 1000, аналогичный Фиг.9. Верхняя диаграмма 1010 представляет входной сигнал этапа анализа с окном анализа 1011, а на нижней диаграмме 1020 показан выходной сигнал этапа синтеза с окном синтеза 1021. Размер ДПФ адаптируется так, чтобы избежать появления опережающего эхо. Этого можно достичь, задавая такой размер М ДПФ, при котором окно синтеза не выбирает нежелательных копий импульса Дирака из получающейся последовательности импульсов. Величина ДПФ 1001 увеличивается до M=FL, где L - длина оконной функции 1002, а коэффициент F - коэффициент передискретизации частотной области. Другими словами, размер ДПФ 1001 выбирается большим, чем величина окна 1002. В частности, размер ДПФ 1001 может быть выбран больше, чем величина окна 1002 синтеза. Вследствие увеличенного размера 1001 ДПФ, период последовательности импульсов, содержащий импульсы Дирака 1022, 1024, равен FL. Выбирая достаточно большое значение F, т.е. выбирая достаточно большой коэффициент передискретизации частотной области, можно исключить нежелательные вклады в изменение временного масштаба импульса. Это представлено на Фиг.10, где импульс Дирака 1024 в момент времени t=Tt0-FL находится за пределами окна синтеза 1021. Следовательно, импульс Дирака 1024 не выбирается окном синтеза 1021, в результате чего можно избежать появления опережающего эхо.
Требуется отметить, что в предпочтительном варианте воплощения окно синтеза и окно анализа имеют равную "номинальную" длину (измеренную числом отсчетов). Однако, при использовании неявной редискретизации выходного сигнала путем удаления или вставки отсчетов в частотные диапазоны преобразования или банка фильтров, величина окна синтеза (измеренная числом отсчетов), как правило, будет отличаться от величины анализа, в зависимости от коэффициента редискретизации и/или транспозиции.
Минимальное значение F, т.е. минимум коэффициента передискретизации частотной области, можно получить, воспользовавшись Фиг.10. Условие для пропуска нежелательных копий импульса Дирака может быть сформулировано следующим образом: для любого входного импульса δ(t-t0) в положении
, т.е. для любого входного импульса, находящегося внутри окна анализа 1011, нежелательная копия δ(t-Tt0+FL) в момент времени t=Tt0-FL должна находиться слева от левой границы окна синтеза
. Эквивалентным образом, должно выполняться условие
, которое приводит к формуле
Как видно из формулы (6), минимальный коэффициент передискретизации частотного домена F является функцией порядка транспозиции Т. Точнее говоря, минимальный коэффициент передискретизации частотного домена F пропорционален порядку транспозиции Т.
Повторив вышеприведенный ход рассуждения для случая, когда окна анализа и синтеза имеют различные длины, получают более общую формулу. Пусть LA и LS - длины окон анализа и синтеза (измеренные числом отсчетов), соответственно, и пусть М - размер используемого ДПФ. В таком случае, формулу (6) можно обобщить в виде
То, что эта формула действительно является обобщением выражения (6), может быть проверено путем подстановки М=FL и LA=LS=L в выражение (7) и делением на L обеих сторон получившегося уравнения.
Рассмотренный выше анализ произведен для довольно специфичной модели переходного импульса, т.е. импульса Дирака. Однако это рассуждение может быть продолжено для того, чтобы продемонстрировать, что при использовании вышеописанной схемы изменения временного масштаба и/или гармонической транспозиции, входные сигналы, имеющие практически плоскую огибающую спектра и стремящиеся к нулю за пределами временного интервала [а, b], будут преобразованы в выходные сигналы, малые по величине за пределами интервала [Ta, Tb]. При изучении спектрограмм реальных звуковых и/или речевых сигналов также может быть подтверждено отсутствие опережающего эхо в сигналах с измененным масштабом времени или в транспонированных сигналах при соблюдении вышеописанного условия для выбора соответствующего коэффициента передискретизации частотной области. Углубленный количественный анализ также показывает, что опережающие эхо по-прежнему подавляются при использовании коэффициентов передискретизации частотного домена с величиной, несколько меньшей, чем значение, определяемое формулой (6) или (7). Это объясняется тем фактом, что типовые оконные функции νS имеют малые значения вблизи границ, тем самым ослабляя нежелательные опережающие эхо, которые находятся вблизи границ оконных функций.
Таким образом, описан способ для улучшения переходной характеристики гармонических транспозиторов частотного домена, или преобразователей временного масштаба, при помощи введения преобразования передискретизации, причем степень передискретизации является функцией выбранного коэффициента транспозиции. Улучшенная переходная характеристика транспозитора получена за счет передискретизации частотного домена.
В параллельном транспозиторе, представленном на Фиг.6, передискретизация частотной области может производиться с помощью ядер ДПФ 601, 602, 603, 604 длиной 1024F и заполнения окон анализа и синтеза нулями симметричным относительно этой длины образом. Следует отметить, что из соображений сложности, предпочтительны малые степени передискретизации. Если для параллельного транспозитора, представленного на Фиг.6, используется формула (6), то коэффициент передискретизации F=2.5 необходимо применять ко всем коэффициентам транспозиции T=2, 3, 4. Однако можно показать, что применение коэффициента F=2.0 уже приводит к значительному повышению качества реальных звуковых сигналов.
Далее будет описано использование передискретизации частотной области в контексте комбинированных банков фильтров анализа/синтеза, например, описанных в контексте Фиг.5 или Фиг.6b.
В общем случае, для комбинированного банка фильтров транспозиции, у которого физический интервал QΔf субполос банка фильтров синтеза в Q раз больше, чем физический интервал Δf банка фильтров анализа, и физическая длительность окна анализа DA (измеренная в единицах времени, например секундах) также в Q раз больше, чем у банка фильтров синтеза, DA=QDS, вышеописанный анализ импульса Дирака применим для всех коэффициентов транспозиции Т=Q, Q+1, Q+2, …, также как для T=Q. Другими словами, формула для степени передискретизации частотного домена, необходимой для комбинированного банка фильтров транспозиции, имеет вид
В частности, требуется отметить, что для Т>Q, коэффициент передискретизации частотной области
является достаточным, поскольку обеспечивает подавление искажений импульсных сигналов, возникающих в процессе гармонической транспозиции порядка Т. Т.е при использовании вышеприведенных формул передискретизации для комбинированного банка фильтров можно убедиться в том, что даже при использовании более высоких порядков транспозиции T>Q, дальнейшее увеличение коэффициента передискретизации F не требуется. Как видно из уравнения (6b), для того, чтобы избежать появления опережающего эхо, в реализации комбинированного банка фильтров, представленного на Фиг.6b, достаточно использовать коэффициент передискретизации F=1.5. Эта величина меньше, чем коэффициент передискретизации F=2.5, который необходим для параллельного транспозитора, представленного на Фиг.6. Следовательно, сложность выполнения передискретизации частотной области для улучшения переходной характеристики параллельных гармонических транспозиторов может быть дополнительно снижена за счет использования комбинированного банка фильтров анализа/синтеза (вместо использования отдельных банков фильтров анализа/синтеза для различных порядков транспозиции).
В более общем варианте, физическая временная длительность окон анализа и синтеза DA и DS, соответственно, могут выбираться произвольно. В таком случае физический интервал Δf субполос банка фильтра анализа должен удовлетворять условию
для того чтобы избежать вышеописанных искажений, возникающих в процессе гармонической транспозиции. Требуется отметить, что длительность окна D как правило, отличается от длины окна L. В то время как длина окна L соответствует числу отсчетов сигнала, накрываемых окном, длительность окна D соответствует временному интервалу сигнала, накрываемого окном. Как представлено на Фиг.6а, окна 611, 612, 613, 614 имеют равную длину из L=1024 отсчетов. Однако, длительность DA окна анализа 611 в T раз больше длительности DS окна синтеза 612, 613, 614, где Т - соответствующий порядок транспозиции и коэффициент разрешения соответствующего банка фильтров синтеза. Аналогично, длительность DA окна анализа 611 на Фиг.6b в Q раз больше длительности DS окна синтеза 612, где Q - коэффициент разрешения банка фильтров синтеза. Длительность окна D связана с длиной окна L через частоту дискретизации fS, а именно
. Аналогичным образом, разрешение по частоте преобразования Δf связано с числом точек или длиной М преобразования через частоту дискретизации fS, a именно
. Более того, физический временной шаг Δt банка фильтров связан с величиной сдвига Δs банка фильтров через частоту дискретизации fS, а именно
.
Используя вышеприведенные соотношения, уравнение (6b) может быть записано в виде
т.е. произведение разрешения по частоте и длины окна для банка фильтров анализа и/или произведение разрешения по частоте и длины окна для банка фильтров синтеза необходимо выбирать меньшим или равным
. Для случая, когда T>Q, произведение ΔfDA и/или QΔfDS может быть выбрано большим, чем
, при этом вычислительная сложность банков фильтров снижается.
В настоящем документе описаны различные методы выполнения гармонической транспозиции сигналов, преимущественно звуковых и/или речевых сигналов. Особое внимание было уделено вычислительной сложности параллельных гармонических транспозиторов. В этом контексте был описан параллельный гармонический транспозитор, сконфигурированный с возможностью производить несколько порядков транспозиции, используя комбинированные банки фильтров анализа и синтеза, т.е. банк фильтров, содержащий единый банк фильтров анализа и единый банк фильтров синтеза. Параллельный транспозитор, использующий комбинированный банк фильтров анализа/синтеза, имеет меньшую вычислительную сложность по сравнению с традиционным параллельным транспозитором. Более того, в контексте комбинированных банков фильтров анализа/синтеза описана передискретизация частотной области. Передискретизация частотной области может быть использована для того, чтобы подавить или устранить искажения, возникающие вследствие гармонической транспозиции импульсных сигналов. Показано, что передискретизация частотной области может быть реализована в комбинированных банках фильтров анализа/синтеза с меньшей вычислительной сложностью по сравнению с традиционными вариантами исполнения параллельного транспозитора.
Хотя в настоящем документе описаны конкретные варианты воплощения и применения настоящего изобретения, для специалиста средней квалификации очевидно, что в варианты воплощения и применения, описанные в настоящем документе, могут быть внесены многочисленные изменения, не выходящие за пределы объема изобретения, описанного и заявленного в настоящем документе. Совершенно ясно, что хотя были показаны и описаны конкретные конфигурации настоящего изобретения, это изобретение не ограничено конкретными описанными и показанными вариантами воплощения или конкретными описанными способами.
Способы и системы, описанные в настоящем документе, могут быть реализованы в виде программного обеспечения, микропрограммы и/или аппаратного обеспечения. Некоторые компоненты могут быть реализованы, например, в виде программного обеспечения, выполняемого цифровым сигнальным процессором или микропроцессором. Другие компоненты могут быть реализованы, например, в виде аппаратного обеспечения или в виде специализированной интегральной микросхемы (ASIC). Сигналы, встречающиеся в описанных способах и системах, могут быть сохранены на носителе информации, например, в запоминающем устройстве с произвольным доступом или оптическом носителе информации. Они могут передаваться по сетям, таким как радиосети, спутниковые сети, беспроводные сети или проводные сети, например Интернет. Типичными устройствами, в которых могут найти применение способы, описанные в настоящем документе, являются, например, мультимедийные проигрыватели или абонентские приставки, декодирующие звуковые сигналы. Для целей кодирования эти системы и способы могут использоваться, например, станциями вещания и производителями мультимедийной продукции.
Claims (25)
1. Система, сконфигурированная для генерирования высокочастотной компоненты сигнала из низкочастотной компоненты этого сигнала, содержащая:
- банк фильтров анализа (501), содержащий модуль преобразования анализа (601) с разрешением по частоте Δf и окно анализа (611) длительностью DA; причем банк фильтров анализа (501) сконфигурирован с возможностью создания набора сигналов субполос анализа из низкочастотной компоненты сигнала;
- модуль нелинейной обработки (502, 650), сконфигурированный с возможностью определения набора сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т, и
- банк фильтров синтеза (504), содержащий модуль преобразования синтеза (602) с разрешением по частоте QΔf и окно синтеза (612) длительностью DS; причем банк фильтров синтеза (504) сконфигурирован с возможностью создания высокочастотной компоненты сигнала из набора сигналов субполос синтеза;
при этом Q является коэффициентом разрешения по частоте, с величиной Q≥1 и меньшей, чем порядок транспозиции Т, а значение произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа выбирается исходя из коэффициента разрешения по частоте Q.
- банк фильтров анализа (501), содержащий модуль преобразования анализа (601) с разрешением по частоте Δf и окно анализа (611) длительностью DA; причем банк фильтров анализа (501) сконфигурирован с возможностью создания набора сигналов субполос анализа из низкочастотной компоненты сигнала;
- модуль нелинейной обработки (502, 650), сконфигурированный с возможностью определения набора сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т, и
- банк фильтров синтеза (504), содержащий модуль преобразования синтеза (602) с разрешением по частоте QΔf и окно синтеза (612) длительностью DS; причем банк фильтров синтеза (504) сконфигурирован с возможностью создания высокочастотной компоненты сигнала из набора сигналов субполос синтеза;
при этом Q является коэффициентом разрешения по частоте, с величиной Q≥1 и меньшей, чем порядок транспозиции Т, а значение произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа выбирается исходя из коэффициента разрешения по частоте Q.
5. Система по по п.4, отличающаяся тем, что значение произведения ΔfDA для банка фильтров анализа (501) равно значению произведения QΔfDS для банка фильтров синтеза (504).
6. Система по по п.5, отличающаяся тем, что
- модуль преобразования анализа (601) сконфигурирован с возможностью выполнения одного из следующих преобразований: преобразование Фурье, быстрое преобразование Фурье, дискретное преобразование Фурье, вейвлет-преобразование; и
- модуль преобразования синтеза (602) сконфигурирован с возможностью выполнения соответствующего обратного преобразования.
- модуль преобразования анализа (601) сконфигурирован с возможностью выполнения одного из следующих преобразований: преобразование Фурье, быстрое преобразование Фурье, дискретное преобразование Фурье, вейвлет-преобразование; и
- модуль преобразования синтеза (602) сконфигурирован с возможностью выполнения соответствующего обратного преобразования.
7. Система по по п.6, отличающаяся тем, что окно анализа (611) и/или окно синтеза (612) является одним из следующих окон:
- окно Гаусса;
- косинус-окно;
- окно Хэмминга;
- окно Ханна;
- прямоугольное окно;
- окно Бартлетта;
- окно Блэкмана.
- окно Гаусса;
- косинус-окно;
- окно Хэмминга;
- окно Ханна;
- прямоугольное окно;
- окно Бартлетта;
- окно Блэкмана.
8. Система по по п.7, отличающаяся тем, что дополнительно содержит:
- второй модуль нелинейной обработки (502), сконфигурированный с возможностью определять второй набор сигналов субполос синтеза из набора сигналов субполос анализа, используя второй порядок транспозиции Т2; причем второй набор сигналов субполос синтеза определяется исходя из части набора сигналов субполос анализа, сдвинутой по фазе на второй порядок транспозиции Т2; при этом порядок транспозиции Т и второй порядок транспозиции Т2 различны; и
- суммирующий модуль (503), сконфигурированный с возможностью производить суммирование набора сигналов субполос синтеза и второго набора сигналов субполос синтеза, получая таким образом суммарный набор сигналов субполос синтеза, поступающий на вход банка фильтров синтеза (602).
- второй модуль нелинейной обработки (502), сконфигурированный с возможностью определять второй набор сигналов субполос синтеза из набора сигналов субполос анализа, используя второй порядок транспозиции Т2; причем второй набор сигналов субполос синтеза определяется исходя из части набора сигналов субполос анализа, сдвинутой по фазе на второй порядок транспозиции Т2; при этом порядок транспозиции Т и второй порядок транспозиции Т2 различны; и
- суммирующий модуль (503), сконфигурированный с возможностью производить суммирование набора сигналов субполос синтеза и второго набора сигналов субполос синтеза, получая таким образом суммарный набор сигналов субполос синтеза, поступающий на вход банка фильтров синтеза (602).
9. Система по п.8, отличающаяся тем, что
- суммирующий модуль (503) сконфигурирован с возможностью совмещать сигналы субполос синтеза набора сигналов субполос синтеза и второго набора сигналов субполос синтеза в соответствии с перекрывающимися интервалами частот.
- суммирующий модуль (503) сконфигурирован с возможностью совмещать сигналы субполос синтеза набора сигналов субполос синтеза и второго набора сигналов субполос синтеза в соответствии с перекрывающимися интервалами частот.
10. Система по по п.9, отличающаяся тем, что
- банк фильтров анализа (501) имеет некоторое число КА субполос анализа, где КА>1, и k является индексом субполосы анализа, принимающим значения k=0,…, КА-1; и
- банк фильтров синтеза (504) имеет некоторое число NS субполос синтеза, где NS>0, и n является индексом субполосы синтеза, принимающим значения n=0,…, NS-1.
- банк фильтров анализа (501) имеет некоторое число КА субполос анализа, где КА>1, и k является индексом субполосы анализа, принимающим значения k=0,…, КА-1; и
- банк фильтров синтеза (504) имеет некоторое число NS субполос синтеза, где NS>0, и n является индексом субполосы синтеза, принимающим значения n=0,…, NS-1.
11. Система по п.10, отличающаяся тем, что модуль нелинейной обработки (502, 650) сконфигурирован с возможностью определения сигнала n-й субполосы синтеза из набора сигналов субполос синтеза, исходя из сигнала k-й субполосы анализа и сигнала (k+1)-й субполосы анализа из набора сигналов субполос анализа.
12. Система по п.11, отличающаяся тем, что модуль нелинейной обработки (502, 650) сконфигурирован с возможностью
- определения фазы сигнала n-й субполосы синтеза в виде суммы фазового сдвига сигнала k-й субполосы анализа и фазового сдвига сигнала (k+1)-й субполосы анализа; и/или
- определения величины сигнала n-й субполосы синтеза в виде произведения возведенной в степень величины сигнала k-й субполосы анализа и возведенной в степень величины сигнала (k+1)-й субполосы анализа.
- определения фазы сигнала n-й субполосы синтеза в виде суммы фазового сдвига сигнала k-й субполосы анализа и фазового сдвига сигнала (k+1)-й субполосы анализа; и/или
- определения величины сигнала n-й субполосы синтеза в виде произведения возведенной в степень величины сигнала k-й субполосы анализа и возведенной в степень величины сигнала (k+1)-й субполосы анализа.
14. Система по п.13, отличающаяся тем, что модуль нелинейной обработки (502, 650) сконфигурирован с возможностью
- определения фазы сигнала n-й субполосы синтеза в виде суммы фазы сигнала k-й субполосы анализа, умноженной на T(1-r) и фазы сигнала (k+1)-й субполосы анализа, умноженной на Т(r); и/или
- определения величины сигнала n-й субполосы синтеза в виде произведения величины сигнала k-й субполосы анализа, возведенной в степень (1-r) и величины сигнала (k+1)-й субполосы анализа, возведенной в степень r.
- определения фазы сигнала n-й субполосы синтеза в виде суммы фазы сигнала k-й субполосы анализа, умноженной на T(1-r) и фазы сигнала (k+1)-й субполосы анализа, умноженной на Т(r); и/или
- определения величины сигнала n-й субполосы синтеза в виде произведения величины сигнала k-й субполосы анализа, возведенной в степень (1-r) и величины сигнала (k+1)-й субполосы анализа, возведенной в степень r.
15. Система по по п.14, отличающаяся тем, что
- банк фильтров анализа (501) и банк фильтров синтеза (504) набраны четным образом так, что центральная частота субполосы анализа имеет вид kΔf, а центральная частота субполосы синтеза имеет вид nQΔf.
- банк фильтров анализа (501) и банк фильтров синтеза (504) набраны четным образом так, что центральная частота субполосы анализа имеет вид kΔf, а центральная частота субполосы синтеза имеет вид nQΔf.
16. Система по любому из пп.1-14, отличающаяся тем, что
- банк фильтров анализа (501) и банк фильтров синтеза (504) набраны нечетным образом так, что центральная частота субполосы анализа имеет вид
, а центральная частота субполосы синтеза имеет вид
, и
- разность между порядком транспозиции Т и коэффициентом разрешения Q выражается четным числом.
- банк фильтров анализа (501) и банк фильтров синтеза (504) набраны нечетным образом так, что центральная частота субполосы анализа имеет вид
- разность между порядком транспозиции Т и коэффициентом разрешения Q выражается четным числом.
17. Система по по п.16, отличающаяся тем, что
- частота дискретизации для низкочастотной компоненты равна fA;
- модуль преобразования анализа (601) является дискретным М-точечным преобразованием; при этом окно анализа (611) имеет длину LA отсчетов; и окно анализа (611) сдвигается по низкочастотной компоненте на величину сдвига анализа, равную ΔsA отсчетам;
- разрешение по частоте имеет вид
; длительность имеет вид
; физический временной шаг банка фильтров анализа (501) имеет вид
;
- частота дискретизации для высокочастотной компоненты равна fS=QfA;
- модуль преобразования синтеза (602) является дискретным М-точечным преобразованием; окно синтеза (612) имеет длину LS отсчетов; и окно синтеза (612) сдвигается по высокочастотной компоненте на величину сдвига синтеза, равную ΔsS отсчетам; и
- разрешение по частоте имеет вид
; длительность имеет вид
; физический временной шаг банка фильтров синтеза (504) имеет вид
.
- частота дискретизации для низкочастотной компоненты равна fA;
- модуль преобразования анализа (601) является дискретным М-точечным преобразованием; при этом окно анализа (611) имеет длину LA отсчетов; и окно анализа (611) сдвигается по низкочастотной компоненте на величину сдвига анализа, равную ΔsA отсчетам;
- разрешение по частоте имеет вид
- частота дискретизации для высокочастотной компоненты равна fS=QfA;
- модуль преобразования синтеза (602) является дискретным М-точечным преобразованием; окно синтеза (612) имеет длину LS отсчетов; и окно синтеза (612) сдвигается по высокочастотной компоненте на величину сдвига синтеза, равную ΔsS отсчетам; и
- разрешение по частоте имеет вид
18. Система для создания выходного сигнала, содержащего
высокочастотную компоненту, из входного сигнала, содержащего низкочастотную компоненту, при помощи порядка транспозиции Т, включающая:
- модуль окна анализа, сконфигурированный с возможностью применять окно анализа (611) длиной LA отсчетов для получения фрейма входного сигнала.
- модуль преобразования анализа (601) порядка М, имеющий разрешение по частоте Δf и сконфигурированный с возможностью преобразования LA отсчетов в М комплексных коэффициентов;
- модуль нелинейной обработки (643, 644, 650), сконфигурированный с возможностью производить изменение фазы комплексных коэффициентов, используя порядок транспозиции Т;
- модуль преобразования синтеза (602) порядка М, имеющий разрешение по частоте QΔf и сконфигурированный с возможностью преобразования измененных коэффициентов в М измененных отсчетов;
причем коэффициент разрешения по частоте Q меньше, чем порядок транспозиции Т; и
- модуль окна синтеза, сконфигурированный с возможностью применять окно синтеза (612) длиной LS отсчетов к М измененным отсчетам, создавая таким образом фрейм выходного сигнала;
при этом М зависит от коэффициента разрешения по частоте Q.
высокочастотную компоненту, из входного сигнала, содержащего низкочастотную компоненту, при помощи порядка транспозиции Т, включающая:
- модуль окна анализа, сконфигурированный с возможностью применять окно анализа (611) длиной LA отсчетов для получения фрейма входного сигнала.
- модуль преобразования анализа (601) порядка М, имеющий разрешение по частоте Δf и сконфигурированный с возможностью преобразования LA отсчетов в М комплексных коэффициентов;
- модуль нелинейной обработки (643, 644, 650), сконфигурированный с возможностью производить изменение фазы комплексных коэффициентов, используя порядок транспозиции Т;
- модуль преобразования синтеза (602) порядка М, имеющий разрешение по частоте QΔf и сконфигурированный с возможностью преобразования измененных коэффициентов в М измененных отсчетов;
причем коэффициент разрешения по частоте Q меньше, чем порядок транспозиции Т; и
- модуль окна синтеза, сконфигурированный с возможностью применять окно синтеза (612) длиной LS отсчетов к М измененным отсчетам, создавая таким образом фрейм выходного сигнала;
при этом М зависит от коэффициента разрешения по частоте Q.
19. Система по п.18, отличающаяся тем, что разность между М и средней длиной окна анализа (611) и окна синтеза (612) пропорциональна (Q-1).
20. Система по п.19, отличающаяся тем, что М больше или равно (QLA+LS)/2.
21. Система по любому из пп.19 и 20, отличающаяся тем, что М больше или равно (TLA+LS)/2.
22. Способ создания высокочастотной компоненты сигнала из низкочастотной компоненты этого сигнала, включающий:
- создание набора сигналов субполос анализа из низкочастотной компоненты сигнала, используя банк фильтров анализа (501), содержащий модуль преобразования анализа (601) с разрешением по частоте Δf и окно анализа (611) длительностью DA;
- определение набора сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т; и
- создание высокочастотной компоненты сигнала по набору сигналов субполос синтеза, используя банк фильтров синтеза (504), содержащий модуль преобразования синтеза (602) с разрешением по частоте QΔf и окно синтеза (612) длительностью DS;
при этом Q является коэффициентом разрешения по частоте с величиной Q≥1 и меньшей, чем порядок транспозиции Т; и значение произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа выбирают исходя из коэффициента разрешения по частоте Q.
- создание набора сигналов субполос анализа из низкочастотной компоненты сигнала, используя банк фильтров анализа (501), содержащий модуль преобразования анализа (601) с разрешением по частоте Δf и окно анализа (611) длительностью DA;
- определение набора сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т; и
- создание высокочастотной компоненты сигнала по набору сигналов субполос синтеза, используя банк фильтров синтеза (504), содержащий модуль преобразования синтеза (602) с разрешением по частоте QΔf и окно синтеза (612) длительностью DS;
при этом Q является коэффициентом разрешения по частоте с величиной Q≥1 и меньшей, чем порядок транспозиции Т; и значение произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа выбирают исходя из коэффициента разрешения по частоте Q.
23. Способ создания выходного сигнала, содержащего высокочастотную компоненту, из входного сигнала, содержащего низкочастотную компоненту, при помощи транспозиции порядка Т, включающий:
- применение окна анализа (611) длиной LA отсчетов для получения фрейма входного сигнала;
- преобразование фрейма, состоящего из LA отсчетов входного сигнала, в М комплексных коэффициентов, используя преобразование анализа порядка М с разрешением по частоте Δf;
- изменение фазы комплексных коэффициентов, используя порядок транспозиции Т;
- преобразование измененных коэффициентов в М измененных отсчетов, используя преобразование синтеза порядка М с разрешением по частоте QΔf; при этом коэффициент разрешения по частоте Q меньше порядка транспозиции Т; и
- применение окна синтеза (612) длиной LS, отсчетов к М измененным отсчетам для создания фрейма выходного сигнала.
при этом М зависит от коэффициента разрешения по частоте Q.
- применение окна анализа (611) длиной LA отсчетов для получения фрейма входного сигнала;
- преобразование фрейма, состоящего из LA отсчетов входного сигнала, в М комплексных коэффициентов, используя преобразование анализа порядка М с разрешением по частоте Δf;
- изменение фазы комплексных коэффициентов, используя порядок транспозиции Т;
- преобразование измененных коэффициентов в М измененных отсчетов, используя преобразование синтеза порядка М с разрешением по частоте QΔf; при этом коэффициент разрешения по частоте Q меньше порядка транспозиции Т; и
- применение окна синтеза (612) длиной LS, отсчетов к М измененным отсчетам для создания фрейма выходного сигнала.
при этом М зависит от коэффициента разрешения по частоте Q.
24. Способ проектирования гармонического транспозитора, сконфигурированного с возможностью создания высокочастотной компоненты сигнала из низкочастотной компоненты этого сигнала, включающий:
- обеспечение банка фильтров анализа (501), содержащего модуль преобразования анализа (601) с разрешением по частоте Δf и окно анализа (611) длительностью DA; причем банк фильтров анализа (501) сконфигурирован с возможностью создавать набор сигналов субполос анализа из низкочастотной компоненты сигнала;
- обеспечение модуля нелинейной обработки (502, 650), сконфигурированного с возможностью определять набор сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т;
- обеспечение банка фильтров синтеза (504), состоящего из модуля преобразования синтеза (602) с разрешением по частоте QΔf и окна синтеза (612) длительностью DS; при этом банк фильтров синтеза (504) сконфигурирован с возможностью создавать высокочастотную компоненту сигнала по набору сигналов субполос синтеза; причем Q является коэффициентом разрешения по частоте таким, что Q≥1 и меньше порядка транспозиции Т; и
- выбор значения произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа исходя из коэффициента разрешения по частоте Q.
- обеспечение банка фильтров анализа (501), содержащего модуль преобразования анализа (601) с разрешением по частоте Δf и окно анализа (611) длительностью DA; причем банк фильтров анализа (501) сконфигурирован с возможностью создавать набор сигналов субполос анализа из низкочастотной компоненты сигнала;
- обеспечение модуля нелинейной обработки (502, 650), сконфигурированного с возможностью определять набор сигналов субполос синтеза по части набора сигналов субполос анализа, причем эта часть набора сигналов субполос анализа сдвинута по фазе на порядок транспозиции Т;
- обеспечение банка фильтров синтеза (504), состоящего из модуля преобразования синтеза (602) с разрешением по частоте QΔf и окна синтеза (612) длительностью DS; при этом банк фильтров синтеза (504) сконфигурирован с возможностью создавать высокочастотную компоненту сигнала по набору сигналов субполос синтеза; причем Q является коэффициентом разрешения по частоте таким, что Q≥1 и меньше порядка транспозиции Т; и
- выбор значения произведения разрешения по частоте Δf и длительности DA для банка фильтров анализа исходя из коэффициента разрешения по частоте Q.
25. Способ проектирования транспозитора, сконфигурированного с возможностью создания выходного сигнала, содержащего высокочастотную компоненту, из входного сигнала, содержащего низкочастотную компоненту, при помощи транспозиции порядка Т, включающий в себя:
- обеспечение модуля окна анализа, при этом модуль сконфигурирован с возможностью применять окно анализа (611) длиной LA отсчетов для получения фрейма входного сигнала.
- обеспечение модуля преобразования анализа М-го порядка (601) с разрешением по частоте Δf, при этом модуль сконфигурирован с возможностью производить преобразование LA отсчетов в М комплексных коэффициентов;
- обеспечение модуля нелинейной обработки (643, 644, 650), сконфигурированного с возможностью производить изменение фазы комплексных коэффициентов, используя порядок транспозиции Т;
- обеспечение модуля преобразования синтеза (602) порядка М, при этом модуль имеет разрешение по частоте QΔf и сконфигурирован с возможностью производить преобразование измененных коэффициентов в М измененных отсчетов; а коэффициент разрешения по частоте Q меньше порядка транспозиции Т;
- обеспечение модуля окна синтеза, при этом модуль сконфигурирован с возможностью применять окно синтеза (612) длиной LS отсчетов к М измененным отсчетам, создавая таким образом фрейм выходного сигнала;
- выбор значения М исходя из коэффициента разрешения по частоте Q.
- обеспечение модуля окна анализа, при этом модуль сконфигурирован с возможностью применять окно анализа (611) длиной LA отсчетов для получения фрейма входного сигнала.
- обеспечение модуля преобразования анализа М-го порядка (601) с разрешением по частоте Δf, при этом модуль сконфигурирован с возможностью производить преобразование LA отсчетов в М комплексных коэффициентов;
- обеспечение модуля нелинейной обработки (643, 644, 650), сконфигурированного с возможностью производить изменение фазы комплексных коэффициентов, используя порядок транспозиции Т;
- обеспечение модуля преобразования синтеза (602) порядка М, при этом модуль имеет разрешение по частоте QΔf и сконфигурирован с возможностью производить преобразование измененных коэффициентов в М измененных отсчетов; а коэффициент разрешения по частоте Q меньше порядка транспозиции Т;
- обеспечение модуля окна синтеза, при этом модуль сконфигурирован с возможностью применять окно синтеза (612) длиной LS отсчетов к М измененным отсчетам, создавая таким образом фрейм выходного сигнала;
- выбор значения М исходя из коэффициента разрешения по частоте Q.
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US25377509P | 2009-10-21 | 2009-10-21 | |
US61/253,775 | 2009-10-21 | ||
US33078610P | 2010-05-03 | 2010-05-03 | |
US61/330,786 | 2010-05-03 | ||
PCT/EP2010/057156 WO2011047887A1 (en) | 2009-10-21 | 2010-05-25 | Oversampling in a combined transposer filter bank |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2494478C1 true RU2494478C1 (ru) | 2013-09-27 |
Family
ID=42542966
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2012116009/08A RU2494478C1 (ru) | 2009-10-21 | 2010-05-25 | Передискретизация в комбинированном банке фильтров транспозитора |
Country Status (11)
Country | Link |
---|---|
US (10) | US8886346B2 (ru) |
EP (7) | EP4276823B1 (ru) |
JP (3) | JP5345737B2 (ru) |
KR (1) | KR101309671B1 (ru) |
CN (3) | CN103559889B (ru) |
BR (2) | BR112012009375B1 (ru) |
ES (6) | ES2936307T3 (ru) |
HK (2) | HK1168189A1 (ru) |
PL (5) | PL3723090T3 (ru) |
RU (1) | RU2494478C1 (ru) |
WO (1) | WO2011047887A1 (ru) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2756042C2 (ru) * | 2013-09-30 | 2021-09-24 | Конинклейке Филипс Н.В. | Передискретизация звукового сигнала для кодирования/декодирования с малой задержкой |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7928310B2 (en) * | 2002-11-12 | 2011-04-19 | MediaLab Solutions Inc. | Systems and methods for portable audio synthesis |
PL4231291T3 (pl) * | 2008-12-15 | 2024-04-15 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Dekoder powiększania szerokości pasma audio, powiązany sposób oraz program komputerowy |
US8971551B2 (en) | 2009-09-18 | 2015-03-03 | Dolby International Ab | Virtual bass synthesis using harmonic transposition |
EP4276823B1 (en) * | 2009-10-21 | 2024-07-17 | Dolby International AB | Oversampling in a combined transposer filter bank |
KR101483157B1 (ko) * | 2010-03-09 | 2015-01-15 | 프라운호퍼 게젤샤프트 쭈르 푀르데룽 데어 안겐반텐 포르슝 에. 베. | 오디오 신호들의 대역폭 연장에 기반한 위상 보코더의 개선된 크기 응답과 시간적 정렬을 위한 방법과 장치 |
PL2532002T3 (pl) * | 2010-03-09 | 2014-06-30 | Fraunhofer Ges Forschung | Urządzenie, sposób i program komputerowy do przetwarzania sygnału audio |
EP2545553B1 (en) | 2010-03-09 | 2014-07-30 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Apparatus and method for processing an audio signal using patch border alignment |
BR112012026502B1 (pt) * | 2010-04-16 | 2022-10-18 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V | Aparelho, método para gerar um sinal de banda larga usando extensão de largura de guiada e extensão de largura de banda cega |
US8958510B1 (en) * | 2010-06-10 | 2015-02-17 | Fredric J. Harris | Selectable bandwidth filter |
SG10201505469SA (en) * | 2010-07-19 | 2015-08-28 | Dolby Int Ab | Processing of audio signals during high frequency reconstruction |
US12002476B2 (en) | 2010-07-19 | 2024-06-04 | Dolby International Ab | Processing of audio signals during high frequency reconstruction |
EP2777042B1 (en) | 2011-11-11 | 2019-08-14 | Dolby International AB | Upsampling using oversampled sbr |
EP2717261A1 (en) | 2012-10-05 | 2014-04-09 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Encoder, decoder and methods for backward compatible multi-resolution spatial-audio-object-coding |
US9494035B2 (en) | 2012-11-06 | 2016-11-15 | Evolution Engineering Inc. | Fluid pressure pulse generator and method of using same |
US10753201B2 (en) | 2012-12-17 | 2020-08-25 | Evolution Engineering Inc. | Mud pulse telemetry apparatus with a pressure transducer and method of operating same |
CA2895346C (en) | 2012-12-17 | 2018-10-23 | Evolution Engineering Inc. | Downhole telemetry signal modulation using pressure pulses of multiple pulse heights |
CA3036490C (en) | 2012-12-17 | 2021-08-03 | Evolution Engineering Inc. | Mud pulse telemetry apparatus with a pressure transducer and method of operating same |
EP2830061A1 (en) | 2013-07-22 | 2015-01-28 | Fraunhofer Gesellschaft zur Förderung der angewandten Forschung e.V. | Apparatus and method for encoding and decoding an encoded audio signal using temporal noise/patch shaping |
US9577798B1 (en) * | 2014-04-30 | 2017-02-21 | Keysight Technologies, Inc. | Real-time separation of signal components in spectrum analyzer |
US9306606B2 (en) * | 2014-06-10 | 2016-04-05 | The Boeing Company | Nonlinear filtering using polyphase filter banks |
CA2895681A1 (en) | 2014-06-27 | 2015-12-27 | Evolution Engineering Inc. | Fluid pressure pulse generator for a downhole telemetry tool |
US9631487B2 (en) | 2014-06-27 | 2017-04-25 | Evolution Engineering Inc. | Fluid pressure pulse generator for a downhole telemetry tool |
US9670774B2 (en) | 2014-06-27 | 2017-06-06 | Evolution Engineering Inc. | Fluid pressure pulse generator for a downhole telemetry tool |
EP2980791A1 (en) | 2014-07-28 | 2016-02-03 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Processor, method and computer program for processing an audio signal using truncated analysis or synthesis window overlap portions |
US10847170B2 (en) | 2015-06-18 | 2020-11-24 | Qualcomm Incorporated | Device and method for generating a high-band signal from non-linearly processed sub-ranges |
US9837089B2 (en) * | 2015-06-18 | 2017-12-05 | Qualcomm Incorporated | High-band signal generation |
US9704497B2 (en) * | 2015-07-06 | 2017-07-11 | Apple Inc. | Method and system of audio power reduction and thermal mitigation using psychoacoustic techniques |
CN107919136B (zh) * | 2017-11-13 | 2021-07-09 | 河海大学 | 一种基于高斯混合模型的数字语音采样频率估计方法 |
DE102018206159B3 (de) * | 2018-04-20 | 2019-09-19 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Paket-Detektor/Decoder für ein Funkübertragungssystem |
KR20210005164A (ko) * | 2018-04-25 | 2021-01-13 | 돌비 인터네셔널 에이비 | 고주파 오디오 재구성 기술의 통합 |
US11322127B2 (en) * | 2019-07-17 | 2022-05-03 | Silencer Devices, LLC. | Noise cancellation with improved frequency resolution |
US11889280B2 (en) * | 2021-10-05 | 2024-01-30 | Cirrus Logic Inc. | Filters and filter chains |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998057436A2 (en) * | 1997-06-10 | 1998-12-17 | Lars Gustaf Liljeryd | Source coding enhancement using spectral-band replication |
US20030187663A1 (en) * | 2002-03-28 | 2003-10-02 | Truman Michael Mead | Broadband frequency translation for high frequency regeneration |
RU2256293C2 (ru) * | 1997-06-10 | 2005-07-10 | Коудинг Технолоджиз Аб | Усовершенствование исходного кодирования с использованием дублирования спектральной полосы |
RU2006118682A (ru) * | 2003-10-30 | 2007-12-10 | Конинклейке Филипс Электроникс Н.В. (Nl) | Кодирование иди декорирование аудиосигнала |
Family Cites Families (140)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4016540A (en) * | 1970-12-28 | 1977-04-05 | Gilbert Peter Hyatt | Apparatus and method for providing interactive audio communication |
US4060848A (en) * | 1970-12-28 | 1977-11-29 | Gilbert Peter Hyatt | Electronic calculator system having audio messages for operator interaction |
JPS5237465B2 (ru) | 1973-03-28 | 1977-09-22 | ||
JPS5237465U (ru) | 1975-09-09 | 1977-03-16 | ||
US4395700A (en) * | 1980-08-15 | 1983-07-26 | Environmental Research Institute Of Michigan | Image analyzer with variable line storage |
FR2628918B1 (fr) * | 1988-03-15 | 1990-08-10 | France Etat | Dispositif annuleur d'echo a filtrage en sous-bandes de frequence |
US5357594A (en) * | 1989-01-27 | 1994-10-18 | Dolby Laboratories Licensing Corporation | Encoding and decoding using specially designed pairs of analysis and synthesis windows |
US5297236A (en) * | 1989-01-27 | 1994-03-22 | Dolby Laboratories Licensing Corporation | Low computational-complexity digital filter bank for encoder, decoder, and encoder/decoder |
US5235623A (en) * | 1989-11-14 | 1993-08-10 | Nec Corporation | Adaptive transform coding by selecting optimum block lengths according to variatons between successive blocks |
US5408580A (en) * | 1992-09-21 | 1995-04-18 | Aware, Inc. | Audio compression system employing multi-rate signal analysis |
US5732389A (en) * | 1995-06-07 | 1998-03-24 | Lucent Technologies Inc. | Voiced/unvoiced classification of speech for excitation codebook selection in celp speech decoding during frame erasures |
US5890106A (en) * | 1996-03-19 | 1999-03-30 | Dolby Laboratories Licensing Corporation | Analysis-/synthesis-filtering system with efficient oddly-stacked singleband filter bank using time-domain aliasing cancellation |
US6073100A (en) * | 1997-03-31 | 2000-06-06 | Goodridge, Jr.; Alan G | Method and apparatus for synthesizing signals using transform-domain match-output extension |
US6026356A (en) * | 1997-07-03 | 2000-02-15 | Nortel Networks Corporation | Methods and devices for noise conditioning signals representative of audio information in compressed and digitized form |
JP3164038B2 (ja) * | 1997-11-05 | 2001-05-08 | 日本電気株式会社 | 音声帯域分割復号装置 |
US6449590B1 (en) * | 1998-08-24 | 2002-09-10 | Conexant Systems, Inc. | Speech encoder using warping in long term preprocessing |
US7072832B1 (en) * | 1998-08-24 | 2006-07-04 | Mindspeed Technologies, Inc. | System for speech encoding having an adaptive encoding arrangement |
US6493665B1 (en) * | 1998-08-24 | 2002-12-10 | Conexant Systems, Inc. | Speech classification and parameter weighting used in codebook search |
US6363338B1 (en) * | 1999-04-12 | 2002-03-26 | Dolby Laboratories Licensing Corporation | Quantization in perceptual audio coders with compensation for synthesis filter noise spreading |
US6246345B1 (en) * | 1999-04-16 | 2001-06-12 | Dolby Laboratories Licensing Corporation | Using gain-adaptive quantization and non-uniform symbol lengths for improved audio coding |
US6298322B1 (en) * | 1999-05-06 | 2001-10-02 | Eric Lindemann | Encoding and synthesis of tonal audio signals using dominant sinusoids and a vector-quantized residual tonal signal |
US6978236B1 (en) * | 1999-10-01 | 2005-12-20 | Coding Technologies Ab | Efficient spectral envelope coding using variable time/frequency resolution and time/frequency switching |
SE0001926D0 (sv) * | 2000-05-23 | 2000-05-23 | Lars Liljeryd | Improved spectral translation/folding in the subband domain |
DE10030583A1 (de) * | 2000-06-21 | 2002-01-10 | Marconi Comm Gmbh | Vorrichtung und Verfahren zum Verarbeiten von Frequenzsignalen |
US20020016698A1 (en) * | 2000-06-26 | 2002-02-07 | Toshimichi Tokuda | Device and method for audio frequency range expansion |
SE0004818D0 (sv) * | 2000-12-22 | 2000-12-22 | Coding Technologies Sweden Ab | Enhancing source coding systems by adaptive transposition |
US7644003B2 (en) * | 2001-05-04 | 2010-01-05 | Agere Systems Inc. | Cue-based audio coding/decoding |
WO2002093546A2 (en) * | 2001-05-16 | 2002-11-21 | Telefonaktiebolaget Lm Ericsson (Publ) | A method for removing aliasing in wave table based synthesisers |
WO2002101726A1 (en) * | 2001-06-08 | 2002-12-19 | Stmicroelectronics Asia Pacific Pte Ltd | Unified filter bank for audio coding |
WO2003019532A1 (en) * | 2001-08-21 | 2003-03-06 | Koninklijke Philips Electronics N.V. | Audio coding with non-uniform filter bank |
US6895375B2 (en) * | 2001-10-04 | 2005-05-17 | At&T Corp. | System for bandwidth extension of Narrow-band speech |
EP1440432B1 (en) * | 2001-11-02 | 2005-05-04 | Matsushita Electric Industrial Co., Ltd. | Audio encoding and decoding device |
FI20012313A (fi) * | 2001-11-26 | 2003-05-27 | Genelec Oy | Menetelmä matalataajuista ääntä muokkaavan modaalisen ekvalisaattorin suunnittelemiseksi |
DE60202881T2 (de) * | 2001-11-29 | 2006-01-19 | Coding Technologies Ab | Wiederherstellung von hochfrequenzkomponenten |
US7146313B2 (en) * | 2001-12-14 | 2006-12-05 | Microsoft Corporation | Techniques for measurement of perceptual audio quality |
US20030138117A1 (en) * | 2002-01-22 | 2003-07-24 | Goff Eugene F. | System and method for the automated detection, identification and reduction of multi-channel acoustical feedback |
FI20020865A (fi) * | 2002-05-07 | 2003-11-08 | Genelec Oy | Menetelmä matalataajuisen modaalisen ekvalisaattorin suunnittelemiseksi erityisesti lähekkäin sijaitsevia moodeja varten |
CA2388352A1 (en) * | 2002-05-31 | 2003-11-30 | Voiceage Corporation | A method and device for frequency-selective pitch enhancement of synthesized speed |
EP1527441B1 (en) * | 2002-07-16 | 2017-09-06 | Koninklijke Philips N.V. | Audio coding |
CN1328707C (zh) * | 2002-07-19 | 2007-07-25 | 日本电气株式会社 | 音频解码设备以及解码方法 |
JP3881943B2 (ja) * | 2002-09-06 | 2007-02-14 | 松下電器産業株式会社 | 音響符号化装置及び音響符号化方法 |
JP3879922B2 (ja) * | 2002-09-12 | 2007-02-14 | ソニー株式会社 | 信号処理システム、信号処理装置および方法、記録媒体、並びにプログラム |
KR100728428B1 (ko) * | 2002-09-19 | 2007-06-13 | 마츠시타 덴끼 산교 가부시키가이샤 | 오디오 디코딩 장치 및 오디오 디코딩 방법 |
US7248711B2 (en) | 2003-03-06 | 2007-07-24 | Phonak Ag | Method for frequency transposition and use of the method in a hearing device and a communication device |
JP4076887B2 (ja) * | 2003-03-24 | 2008-04-16 | ローランド株式会社 | ボコーダ装置 |
DE10328777A1 (de) * | 2003-06-25 | 2005-01-27 | Coding Technologies Ab | Vorrichtung und Verfahren zum Codieren eines Audiosignals und Vorrichtung und Verfahren zum Decodieren eines codierten Audiosignals |
US20050018796A1 (en) * | 2003-07-07 | 2005-01-27 | Sande Ravindra Kumar | Method of combining an analysis filter bank following a synthesis filter bank and structure therefor |
FI118550B (fi) * | 2003-07-14 | 2007-12-14 | Nokia Corp | Parannettu eksitaatio ylemmän kaistan koodaukselle koodekissa, joka käyttää kaistojen jakoon perustuvia koodausmenetelmiä |
EP1507334A1 (fr) * | 2003-08-12 | 2005-02-16 | STMicroelectronics S.A. | Composant électronique permettant notamment le décodage de signaux modulés par une modulation numérique en quadrature sur un grand nombre de porteuses orthogonales |
KR20050027179A (ko) * | 2003-09-13 | 2005-03-18 | 삼성전자주식회사 | 오디오 데이터 복원 방법 및 그 장치 |
KR20050064442A (ko) * | 2003-12-23 | 2005-06-29 | 삼성전자주식회사 | 이동통신 시스템에서 입체음향 신호 생성 장치 및 방법 |
KR20060131798A (ko) * | 2004-01-20 | 2006-12-20 | 돌비 레버러토리즈 라이쎈싱 코오포레이션 | 블록 그룹화에 기반한 오디오 코딩 |
US6980933B2 (en) * | 2004-01-27 | 2005-12-27 | Dolby Laboratories Licensing Corporation | Coding techniques using estimated spectral magnitude and phase derived from MDCT coefficients |
CA2457988A1 (en) * | 2004-02-18 | 2005-08-18 | Voiceage Corporation | Methods and devices for audio compression based on acelp/tcx coding and multi-rate lattice vector quantization |
CN1662100B (zh) * | 2004-02-24 | 2010-12-08 | 三洋电机株式会社 | 低音强调电路以及低音强调处理方法 |
US8983834B2 (en) * | 2004-03-01 | 2015-03-17 | Dolby Laboratories Licensing Corporation | Multichannel audio coding |
US7587254B2 (en) * | 2004-04-23 | 2009-09-08 | Nokia Corporation | Dynamic range control and equalization of digital audio using warped processing |
KR100608062B1 (ko) * | 2004-08-04 | 2006-08-02 | 삼성전자주식회사 | 오디오 데이터의 고주파수 복원 방법 및 그 장치 |
FR2875351A1 (fr) * | 2004-09-16 | 2006-03-17 | France Telecom | Procede de traitement de donnees par passage entre domaines differents de sous-bandes |
US7937271B2 (en) * | 2004-09-17 | 2011-05-03 | Digital Rise Technology Co., Ltd. | Audio decoding using variable-length codebook application ranges |
US7630902B2 (en) * | 2004-09-17 | 2009-12-08 | Digital Rise Technology Co., Ltd. | Apparatus and methods for digital audio coding using codebook application ranges |
JP4580210B2 (ja) * | 2004-10-19 | 2010-11-10 | ソニー株式会社 | 音声信号処理装置および音声信号処理方法 |
US8326606B2 (en) * | 2004-10-26 | 2012-12-04 | Panasonic Corporation | Sound encoding device and sound encoding method |
CN101048649A (zh) * | 2004-11-05 | 2007-10-03 | 松下电器产业株式会社 | 可扩展解码装置及可扩展编码装置 |
KR100657916B1 (ko) * | 2004-12-01 | 2006-12-14 | 삼성전자주식회사 | 주파수 대역간의 유사도를 이용한 오디오 신호 처리 장치및 방법 |
US20070147518A1 (en) * | 2005-02-18 | 2007-06-28 | Bruno Bessette | Methods and devices for low-frequency emphasis during audio compression based on ACELP/TCX |
NZ562182A (en) * | 2005-04-01 | 2010-03-26 | Qualcomm Inc | Method and apparatus for anti-sparseness filtering of a bandwidth extended speech prediction excitation signal |
KR100818268B1 (ko) * | 2005-04-14 | 2008-04-02 | 삼성전자주식회사 | 오디오 데이터 부호화 및 복호화 장치와 방법 |
EP1722360B1 (en) * | 2005-05-13 | 2014-03-19 | Harman Becker Automotive Systems GmbH | Audio enhancement system and method |
DK1742509T3 (da) * | 2005-07-08 | 2013-11-04 | Oticon As | Et system og en fremgangsmåde til eliminering af feedback og støj i et høreapparat |
FR2888699A1 (fr) * | 2005-07-13 | 2007-01-19 | France Telecom | Dispositif de codage/decodage hierachique |
US7197453B2 (en) * | 2005-07-29 | 2007-03-27 | Texas Instruments Incorporated | System and method for optimizing the operation of an oversampled discrete Fourier transform filter bank |
JP4701931B2 (ja) | 2005-09-02 | 2011-06-15 | 日本電気株式会社 | 信号処理の方法及び装置並びにコンピュータプログラム |
US20070083377A1 (en) * | 2005-10-12 | 2007-04-12 | Steven Trautmann | Time scale modification of audio using bark bands |
WO2007043811A1 (en) * | 2005-10-12 | 2007-04-19 | Samsung Electronics Co., Ltd. | Method and apparatus for encoding/decoding audio data and extension data |
CN101297356B (zh) * | 2005-11-04 | 2011-11-09 | 诺基亚公司 | 用于音频压缩的方法和设备 |
US7774396B2 (en) * | 2005-11-18 | 2010-08-10 | Dynamic Hearing Pty Ltd | Method and device for low delay processing |
KR100717058B1 (ko) * | 2005-11-28 | 2007-05-14 | 삼성전자주식회사 | 고주파 성분 복원 방법 및 그 장치 |
WO2007080211A1 (en) * | 2006-01-09 | 2007-07-19 | Nokia Corporation | Decoding of binaural audio signals |
US8081762B2 (en) * | 2006-01-09 | 2011-12-20 | Nokia Corporation | Controlling the decoding of binaural audio signals |
FR2897977A1 (fr) * | 2006-02-28 | 2007-08-31 | France Telecom | Procede de limitation de gain d'excitation adaptative dans un decodeur audio |
US8920343B2 (en) * | 2006-03-23 | 2014-12-30 | Michael Edward Sabatino | Apparatus for acquiring and processing of physiological auditory signals |
US20070299655A1 (en) * | 2006-06-22 | 2007-12-27 | Nokia Corporation | Method, Apparatus and Computer Program Product for Providing Low Frequency Expansion of Speech |
US7876904B2 (en) * | 2006-07-08 | 2011-01-25 | Nokia Corporation | Dynamic decoding of binaural audio signals |
WO2008035949A1 (en) * | 2006-09-22 | 2008-03-27 | Samsung Electronics Co., Ltd. | Method, medium, and system encoding and/or decoding audio signals by using bandwidth extension and stereo coding |
WO2008045950A2 (en) * | 2006-10-11 | 2008-04-17 | Nielsen Media Research, Inc. | Methods and apparatus for embedding codes in compressed audio data streams |
CN101163240A (zh) * | 2006-10-13 | 2008-04-16 | 国际商业机器公司 | 一种滤波装置及其方法 |
US7953595B2 (en) * | 2006-10-18 | 2011-05-31 | Polycom, Inc. | Dual-transform coding of audio signals |
US8041578B2 (en) * | 2006-10-18 | 2011-10-18 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Encoding an information signal |
US8126721B2 (en) * | 2006-10-18 | 2012-02-28 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Encoding an information signal |
US7647229B2 (en) * | 2006-10-18 | 2010-01-12 | Nokia Corporation | Time scaling of multi-channel audio signals |
KR101056253B1 (ko) * | 2006-10-25 | 2011-08-11 | 프라운호퍼 게젤샤프트 쭈르 푀르데룽 데어 안겐반텐 포르슝 에. 베. | 오디오 서브밴드 값을 생성하는 장치 및 방법과 시간 영역 오디오 샘플을 생성하는 장치 및 방법 |
US20100017197A1 (en) * | 2006-11-02 | 2010-01-21 | Panasonic Corporation | Voice coding device, voice decoding device and their methods |
US7756350B2 (en) * | 2006-11-13 | 2010-07-13 | Global Ip Solutions, Inc. | Lossless encoding and decoding of digital data |
US20080243518A1 (en) * | 2006-11-16 | 2008-10-02 | Alexey Oraevsky | System And Method For Compressing And Reconstructing Audio Files |
JP5103880B2 (ja) * | 2006-11-24 | 2012-12-19 | 富士通株式会社 | 復号化装置および復号化方法 |
JP4967618B2 (ja) * | 2006-11-24 | 2012-07-04 | 富士通株式会社 | 復号化装置および復号化方法 |
JPWO2008066071A1 (ja) * | 2006-11-29 | 2010-03-04 | パナソニック株式会社 | 復号化装置および復号化方法 |
KR101377135B1 (ko) * | 2007-01-02 | 2014-03-21 | 삼성전자주식회사 | 오디오 신호의 저주파 및 중주파 성분 보강 방법 및 그장치 |
SG144752A1 (en) * | 2007-01-12 | 2008-08-28 | Sony Corp | Audio enhancement method and system |
JP4708446B2 (ja) * | 2007-03-02 | 2011-06-22 | パナソニック株式会社 | 符号化装置、復号装置およびそれらの方法 |
JP4984983B2 (ja) * | 2007-03-09 | 2012-07-25 | 富士通株式会社 | 符号化装置および符号化方法 |
GB0710211D0 (en) * | 2007-05-29 | 2007-07-11 | Intrasonics Ltd | AMR Spectrography |
US9653088B2 (en) * | 2007-06-13 | 2017-05-16 | Qualcomm Incorporated | Systems, methods, and apparatus for signal encoding using pitch-regularizing and non-pitch-regularizing coding |
EP2015293A1 (en) * | 2007-06-14 | 2009-01-14 | Deutsche Thomson OHG | Method and apparatus for encoding and decoding an audio signal using adaptively switched temporal resolution in the spectral domain |
US7885819B2 (en) * | 2007-06-29 | 2011-02-08 | Microsoft Corporation | Bitstream syntax for multi-process audio decoding |
EP2026597B1 (en) * | 2007-08-13 | 2009-11-11 | Harman Becker Automotive Systems GmbH | Noise reduction by combined beamforming and post-filtering |
US7522074B2 (en) * | 2007-09-17 | 2009-04-21 | Samplify Systems, Inc. | Enhanced control for compression and decompression of sampled signals |
WO2009038056A1 (ja) * | 2007-09-20 | 2009-03-26 | National University Corporation University Of Toyama | 信号解析方法、信号解析装置、及び信号解析プログラム |
US8588427B2 (en) * | 2007-09-26 | 2013-11-19 | Frauhnhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Apparatus and method for extracting an ambient signal in an apparatus and method for obtaining weighting coefficients for extracting an ambient signal and computer program |
JP5098569B2 (ja) * | 2007-10-25 | 2012-12-12 | ヤマハ株式会社 | 帯域拡張再生装置 |
KR101373004B1 (ko) * | 2007-10-30 | 2014-03-26 | 삼성전자주식회사 | 고주파수 신호 부호화 및 복호화 장치 및 방법 |
US8515767B2 (en) * | 2007-11-04 | 2013-08-20 | Qualcomm Incorporated | Technique for encoding/decoding of codebook indices for quantized MDCT spectrum in scalable speech and audio codecs |
ATE456908T1 (de) * | 2007-11-12 | 2010-02-15 | Harman Becker Automotive Sys | Mischung von ersten und zweiten tonsignalen |
US8560307B2 (en) * | 2008-01-28 | 2013-10-15 | Qualcomm Incorporated | Systems, methods, and apparatus for context suppression using receivers |
KR101413967B1 (ko) * | 2008-01-29 | 2014-07-01 | 삼성전자주식회사 | 오디오 신호의 부호화 방법 및 복호화 방법, 및 그에 대한 기록 매체, 오디오 신호의 부호화 장치 및 복호화 장치 |
US8000487B2 (en) * | 2008-03-06 | 2011-08-16 | Starkey Laboratories, Inc. | Frequency translation by high-frequency spectral envelope warping in hearing assistance devices |
US8060042B2 (en) * | 2008-05-23 | 2011-11-15 | Lg Electronics Inc. | Method and an apparatus for processing an audio signal |
US8880410B2 (en) * | 2008-07-11 | 2014-11-04 | Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. | Apparatus and method for generating a bandwidth extended signal |
KR20100008749A (ko) * | 2008-07-16 | 2010-01-26 | 삼성전자주식회사 | 스위칭 파워 증폭 장치 및 그 제어 방법 |
JP4818335B2 (ja) * | 2008-08-29 | 2011-11-16 | 株式会社東芝 | 信号帯域拡張装置 |
US8532983B2 (en) * | 2008-09-06 | 2013-09-10 | Huawei Technologies Co., Ltd. | Adaptive frequency prediction for encoding or decoding an audio signal |
US8831958B2 (en) * | 2008-09-25 | 2014-09-09 | Lg Electronics Inc. | Method and an apparatus for a bandwidth extension using different schemes |
PL4231291T3 (pl) * | 2008-12-15 | 2024-04-15 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Dekoder powiększania szerokości pasma audio, powiązany sposób oraz program komputerowy |
CN101430882B (zh) * | 2008-12-22 | 2012-11-28 | 无锡中星微电子有限公司 | 一种抑制风噪声的方法及装置 |
CA3162807C (en) * | 2009-01-16 | 2024-04-23 | Dolby International Ab | Cross product enhanced harmonic transposition |
CA2750065C (en) * | 2009-01-20 | 2013-12-03 | Med-El Elektromedizinische Geraete Gmbh | High accuracy tonotopic and periodic coding with enhanced harmonic resolution |
CA2749239C (en) * | 2009-01-28 | 2017-06-06 | Dolby International Ab | Improved harmonic transposition |
ES2826324T3 (es) * | 2009-01-28 | 2021-05-18 | Dolby Int Ab | Transposición armónica mejorada |
JP4945586B2 (ja) * | 2009-02-02 | 2012-06-06 | 株式会社東芝 | 信号帯域拡張装置 |
US8359195B2 (en) * | 2009-03-26 | 2013-01-22 | LI Creative Technologies, Inc. | Method and apparatus for processing audio and speech signals |
US8526650B2 (en) * | 2009-05-06 | 2013-09-03 | Starkey Laboratories, Inc. | Frequency translation by high-frequency spectral envelope warping in hearing assistance devices |
US8971551B2 (en) * | 2009-09-18 | 2015-03-03 | Dolby International Ab | Virtual bass synthesis using harmonic transposition |
TWI556227B (zh) * | 2009-05-27 | 2016-11-01 | 杜比國際公司 | 從訊號的低頻成份產生該訊號之高頻成份的系統與方法,及其機上盒、電腦程式產品、軟體程式及儲存媒體 |
ES2673637T3 (es) * | 2009-06-23 | 2018-06-25 | Voiceage Corporation | Cancelación prospectiva de solapamiento en dominio de tiempo con aplicación en dominio de señal ponderada u original |
AU2010309838B2 (en) * | 2009-10-20 | 2014-05-08 | Dolby International Ab | Audio signal encoder, audio signal decoder, method for encoding or decoding an audio signal using an aliasing-cancellation |
EP2491555B1 (en) * | 2009-10-20 | 2014-03-05 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Multi-mode audio codec |
CN102648495B (zh) * | 2009-10-21 | 2014-05-28 | 杜比Ab国际公司 | 用于利用适应性过取样产生高频音频信号的装置及方法 |
EP4276823B1 (en) * | 2009-10-21 | 2024-07-17 | Dolby International AB | Oversampling in a combined transposer filter bank |
RU2526745C2 (ru) * | 2009-12-16 | 2014-08-27 | Долби Интернешнл Аб | Низведение параметров последовательности битов sbr |
-
2010
- 2010-05-25 EP EP23194739.1A patent/EP4276823B1/en active Active
- 2010-05-25 EP EP21210177.8A patent/EP3998606B8/en active Active
- 2010-05-25 EP EP17192422.8A patent/EP3291231B1/en active Active
- 2010-05-25 EP EP22204414.1A patent/EP4152320B1/en active Active
- 2010-05-25 BR BR112012009375-4A patent/BR112012009375B1/pt active IP Right Grant
- 2010-05-25 PL PL20172571T patent/PL3723090T3/pl unknown
- 2010-05-25 PL PL14173151T patent/PL2800094T3/pl unknown
- 2010-05-25 CN CN201310423103.4A patent/CN103559889B/zh active Active
- 2010-05-25 PL PL17192422T patent/PL3291231T3/pl unknown
- 2010-05-25 EP EP20172571.0A patent/EP3723090B1/en active Active
- 2010-05-25 BR BR122020007866-9A patent/BR122020007866B1/pt active IP Right Grant
- 2010-05-25 JP JP2012530183A patent/JP5345737B2/ja active Active
- 2010-05-25 CN CN2010800434072A patent/CN102576541B/zh active Active
- 2010-05-25 CN CN201310423105.3A patent/CN103559890B/zh active Active
- 2010-05-25 US US13/499,893 patent/US8886346B2/en active Active
- 2010-05-25 RU RU2012116009/08A patent/RU2494478C1/ru active
- 2010-05-25 ES ES21210177T patent/ES2936307T3/es active Active
- 2010-05-25 PL PL22204414.1T patent/PL4152320T3/pl unknown
- 2010-05-25 ES ES14173151.3T patent/ES2656668T3/es active Active
- 2010-05-25 EP EP10730734.0A patent/EP2491557B1/en active Active
- 2010-05-25 ES ES17192422T patent/ES2805349T3/es active Active
- 2010-05-25 EP EP14173151.3A patent/EP2800094B1/en active Active
- 2010-05-25 ES ES10730734.0T patent/ES2507165T3/es active Active
- 2010-05-25 PL PL21210177.8T patent/PL3998606T3/pl unknown
- 2010-05-25 ES ES22204414T patent/ES2963061T3/es active Active
- 2010-05-25 ES ES20172571T patent/ES2906085T3/es active Active
- 2010-05-25 KR KR1020127007844A patent/KR101309671B1/ko active IP Right Grant
- 2010-05-25 WO PCT/EP2010/057156 patent/WO2011047887A1/en active Application Filing
-
2012
- 2012-09-06 HK HK12108706.4A patent/HK1168189A1/xx unknown
-
2013
- 2013-08-14 JP JP2013168635A patent/JP5613802B2/ja active Active
-
2014
- 2014-09-08 JP JP2014181997A patent/JP5950974B2/ja active Active
- 2014-10-03 US US14/505,739 patent/US9384750B2/en active Active
- 2014-11-28 HK HK14112055.1A patent/HK1198561A1/xx unknown
-
2016
- 2016-05-26 US US15/165,735 patent/US9830928B2/en active Active
-
2017
- 2017-10-25 US US15/792,956 patent/US10186280B2/en active Active
-
2018
- 2018-12-18 US US16/223,112 patent/US10584386B2/en active Active
-
2020
- 2020-03-05 US US16/810,786 patent/US10947594B2/en active Active
-
2021
- 2021-03-15 US US17/201,107 patent/US11591657B2/en active Active
-
2023
- 2023-01-19 US US18/156,701 patent/US11993817B2/en active Active
-
2024
- 2024-04-17 US US18/637,832 patent/US20240294987A1/en active Pending
- 2024-04-17 US US18/637,814 patent/US20240271217A1/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998057436A2 (en) * | 1997-06-10 | 1998-12-17 | Lars Gustaf Liljeryd | Source coding enhancement using spectral-band replication |
RU2256293C2 (ru) * | 1997-06-10 | 2005-07-10 | Коудинг Технолоджиз Аб | Усовершенствование исходного кодирования с использованием дублирования спектральной полосы |
US20030187663A1 (en) * | 2002-03-28 | 2003-10-02 | Truman Michael Mead | Broadband frequency translation for high frequency regeneration |
RU2006118682A (ru) * | 2003-10-30 | 2007-12-10 | Конинклейке Филипс Электроникс Н.В. (Nl) | Кодирование иди декорирование аудиосигнала |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2756042C2 (ru) * | 2013-09-30 | 2021-09-24 | Конинклейке Филипс Н.В. | Передискретизация звукового сигнала для кодирования/декодирования с малой задержкой |
Also Published As
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2494478C1 (ru) | Передискретизация в комбинированном банке фильтров транспозитора | |
KR101858948B1 (ko) | 고조파 전위에 기초하여 개선된 서브밴드 블록 |