СПОСОБ И УСТРОЙСТВО ДЛЯ СЖАТИЯ ДАННЫХ, ПРЕДСТАВЛЯЮЩИХ ЗАВИСЯЩИЙ ОТ ВРЕМЕНИ СИГНАЛ
METHOD AND APPARATUS FOR COMPRESSING DATA DEPENDING ON TIME SIGNAL

Настоящее изобретение относится к способу и устройству обработки сигналов для сжатия данных, представляющих зависящий от времени сигнал.
Сжатие данных играет центральную роль в управлении колоссальным объемом данных. Например, многие современные инженерные, научные, и относящиеся к статистике виды деятельности требуют сбора больших объемов электронных данных для обработки и хранения информации. Например, электронные данные могут быть представлением изменяющегося со временем сигнала. Управление данными является в высшей степени востребованной задачей, при этом управление также может содержать хранение и передачу данных. Эффективное управление данными достижимо с помощью интеллектуальных способов сжатия данных, вследствие чего пространство хранения, требуемое для хранения информации, присутствующей в данных, и/или ширина полосы передачи, требуемая для передачи данных, в значительной степени уменьшаются.
В текущее время, когда собираются колоссальные объемы данных, например, в случае непрерывного потока данных, представляющих изменяющийся со временем сигнал, собираемый в течение выполнения непрерывного наблюдения, выполняемого над машиной, т.е. наблюдения вибраций машины в течение ее работы, вычисление и хранение преобразования Фурье, т.е. спектральных данных изменяющегося со временем сигнала, особенно в их полноте, делает необходимым огромный объем пространства хранения.
Нормальные способы сжатия спектральных данных включают в себя либо использование универсальных инструментов сжатия без потерь (ZIP), либо использование оптимизированного в отношении памяти кодирования. Однако отношения сжатия (отношение пространства хранения/ресурсов, требуемых, чтобы хранить сжатые данные, к пространству хранения/ресурсам, требуемым, чтобы хранить фактические данные) для вышеупомянутых способов являются заметно низкими и являются неэффективными способами сжатия данных. Таким образом, пространство хранения и ширина полосы передачи данных, по существу, не уменьшаются.
Целью настоящего изобретения является обеспечить эффективные способ и устройство обработки сигналов для сжатия данных, представляющих зависящий от времени сигнал.
Вышеупомянутая цель достигается посредством способа для сжатия данных согласно пункту 1 формулы изобретения и устройства обработки сигналов согласно пункту 13 формулы изобретения.
i
i
i
i
i
j
j,i
j
j,i
j
j
i
i
j,i
j,i
j
j,i
j,i
v
j
j,i
j
Лежащей в основе целью изобретения является сжатие данных D(nT), представляющих зависящий от времени сигнал A(t). Зависящий от времени сигнал A(t) содержит множество зависящих от времени частичных сигналов A(t). Способ сжатия начинается с этапа, на котором принимается множество спектров S(f). Каждый спектр S(f) однозначно соответствует одному из зависящих от времени частичных сигналов A(t) и каждый спектр S(f) содержит множество частот f и множество амплитуд a множества частот f. Амплитуды a однозначно назначаются f и являются мерами вкладов f в A(t). На последующем этапе, для каждого спектра S(f) множество амплитуд a нормализуются для получения множества нормализованных амплитуд b. На последующем этапе, для каждой частоты f, нормализованные амплитуды b обрабатываются на основе распределения нормализованных амплитуд b по отношению к пороговому значению (T), и один или более соответствующих параметров модели (MP) определяются для представления распределения нормализованных амплитуд b. На последующем этапе, генерируется набор сжатых данных (CDS). CDS содержит, по меньшей мере, упомянутые один или более параметров модели и частоту f. Поэтому CDS обеспечивает возможность сжимать огромные объемы данных.
i
i
i
i
i
i
j,i
j,i
j,i
Согласно одному варианту осуществления настоящего изобретения, для каждого частичного сигнала A(t), генерируется временная метка t. Временная метка t представляет момент времени, в который частичный сигнал A(t) был получен. Временные метки t являются предпочтительными для идентификации последовательности появления A(t), для индексирования по спектрам множества амплитуд a и множества нормализованных амплитуд b, и для обработки множества нормализованных амплитуд b для получения CDS.
j
j,i
v
j
j,i
j
1
j
j
j
j
Согласно другому варианту осуществления настоящего изобретения, для каждой частоты f, если каждая нормализованная амплитуда b меньше, чем пороговое значение (T), то MP, определенный из этапа обработки, будет содержать среднее значение (AV) нормализованных амплитуд b для f. Дополнительно, CDS дополнительно содержит первую временную метку t. Хотя вклад f в D(nT) меньше, чем требуемое значение (т.е. пороговое значение), т.е. f не является значимой частотой. Однако вклад f может моделироваться и сохраняться посредством сохранения среднего значения f в D(nT) для получения достоверного представления D(nT) посредством CDS.
j
j,i
v
j,i
v
j
j,i
j
1
j
j
j
j
Согласно еще другому варианту осуществления настоящего изобретения, для каждой частоты f, если каждая нормализованная амплитуда b больше, чем пороговое значение (T), и если разность  между экстремальными значениями  нормализованных амплитуд b меньше, чем пороговое значение (T), то MP, определенный из этапа обработки, будет содержать среднее значение (AV) нормализованных амплитуд b для f. Дополнительно, CDS дополнительно содержит первую временную метку t. Вклад f в D(nT) является более большим, чем требуемое значение (т.е. пороговое значение), тем самым делая f значимой частотой. Вклад f может моделироваться и поддерживаться посредством сохранения среднего значения f в D(nT) для получения достоверного представления D(nT) посредством CDS.
j
j,i
v
j,i
v
j,i
j
j
k
Согласно еще другому варианту осуществления настоящего изобретения, для каждой частоты f, если, по меньшей мере, одна из нормализованных амплитуд b не меньше, чем пороговое значение (T), и если разность между экстремальными значениями  нормализованных амплитуд b не меньше, чем пороговое значение (T), то на этапе обработки выполняется алгоритм аппроксимации для моделирования распределения нормализованных амплитуд b частоты f. Таким образом, MP, определенный из этапа обработки, будет содержать один или более параметров модели алгоритма аппроксимации. Посредством этого, алгоритм аппроксимации упрощает моделирование усложненных распределений. Согласно варианту вышеупомянутого варианта осуществления, CDS также содержит соответствующие временные метки t.
Согласно варианту вышеупомянутого варианта осуществления, алгоритм аппроксимации является полиномиальной подгонкой, которая является простым способом моделирования распределения.
Согласно другому варианту вышеупомянутого варианта осуществления, алгоритм аппроксимации является ступенчатой функцией, которая увеличивает скорость моделирования нелинейного распределения.
Согласно предпочтительному варианту вышеупомянутого варианта осуществления, алгоритм аппроксимации является итеративным алгоритмом подгонки на основе конечных точек, например, алгоритмом Рамера - Дугласа - Пекера. Таким образом, получается меньший набор точек для определения распределения, содержащего большой набор точек. Таким образом, достигаются очень высокие отношения сжатия и точность сжатия.
i
i
j,i
j,i
Согласно еще другому варианту осуществления настоящего изобретения, CDS содержит для каждого спектра S(f), сумму SF множества амплитуд a. Это является предпочтительным, так как это помогает в восстановлении a, так как обработка выполняется над нормализованными амплитудами.
j
f
j,i
Согласно еще другому варианту осуществления настоящего изобретения, пороговое значение является пропорциональным обратному количественного числа (NoP) множества частот f и также коэффициенту допустимого отклонения (T). Это обеспечивает переменную порога, и помогает в оказании влияния на распределение нормализованных амплитуд b.
Согласно еще другому варианту осуществления настоящего изобретения, способ дополнительно содержит этап, следующий за этапом обработки, на котором строятся восстановленные данные D'(nT) посредством обработки CDS. На последующем этапе, определяется коэффициент корреляции (CC) между данными D(nT) и восстановленными данными D'(nT). Порог изменяется в ответ на CC, и вышеупомянутые этапы повторяются до тех пор, когда получается удовлетворительный CC. Посредством этого, точность CDS улучшается.
i
j,i
j
i
j,i
j,i
v
j,i
j
j,i
v
j
j,i
Здесь раскрывается устройство обработки сигналов, сконфигурированное с возможностью сжимать данные D(nT) согласно любому из вышеупомянутых вариантов осуществления. Устройство обработки сигналов содержит модуль приема спектральных данных, модуль нормализации амплитуд, модуль сравнения, модуль параметров и блок памяти. Модуль приема спектральных данных принимает множество спектров S(f). Модуль нормализации амплитуд нормализует множество амплитуд a для каждой из множества частот f каждого из множества спектров S(f). Таким образом, получается множество нормализованных амплитуд b. Модуль сравнения сравнивает каждую из множества нормализованных амплитуд b с пороговым значением (T). Модуль параметров обрабатывает нормализованные амплитуды b частоты f на основе распределения нормализованных амплитуд b по отношению к пороговому значению (T). Таким образом, определяются упомянутые один или более соответствующих параметров (MP) для представления распределения нормализованных амплитуд b. Блок памяти сохраняет набор сжатых данных (CDS).
Согласно одному варианту осуществления настоящего изобретения, устройство обработки сигналов дополнительно содержит модуль корреляции. Модуль корреляции определяет коэффициент корреляции (CC) между данными D(nT) и восстановленными данными D'(nT). Посредством этого, точность CDS для представления данных D(nT) улучшается.
i
i
i
Согласно одному варианту осуществления настоящего изобретения, устройство обработки сигналов дополнительно содержит модуль спектра. Модуль спектра принимает каждый из зависящих от времени частичных сигналов A(t) и вычисляет соответствующий спектр S(f). Модуль спектра может быть интегрированным с вышеупомянутыми модулями в каскаде. Посредством этого, гибкость в применении устройства обработки сигналов улучшается, так как спектры S(f) легко вычисляются и являются доступными для сжатия данных D(nT).
Вышеупомянутые и другие варианты осуществления изобретения, относящиеся к способу и устройству обработки сигналов для сжатия данных, представляющих зависящий от времени сигнал, будут теперь описываться со ссылкой на сопровождающие чертежи настоящего изобретения. Проиллюстрированные чертежи и варианты осуществления предназначены, чтобы иллюстрировать, но не ограничивать изобретение. Сопровождающие чертежи содержат следующие фигуры, на которых сходные ссылочные позиции указывают на сходные части, всюду в описании и чертежах.
Фигуры схематическим образом иллюстрируют дополнительные примеры вариантов осуществления изобретения, на которых:
Фиг. 1 изображает систему наблюдения состояния, содержащую модуль получения данных и устройство обработки сигналов, сжимающее данные, представляющие зависящий от времени сигнал,
Фиг. 2 изображает различные модули устройства обработки сигналов, показанного на фиг. 1, при этом покомпонентные изображения множества спектров, матрица амплитуд, и матрица нормализованных амплитуд схематически иллюстрируют операции модуля спектра, модуля приема спектральных данных, модуля нормализации амплитуд соответственно,
Фиг. 3 изображает иллюстративную матрицу амплитуд, и иллюстративную матрицу нормализованных амплитуд, упомянутые на фиг. 2, и
Фиг. 4 изображает иллюстративный способ построения набора сжатых данных для восстановления данных для сжатия данных, представляющих зависящий от времени сигнал, упомянутый на Фиг. 1, и
Фиг. 5 изображает блок-схему последовательности операций способа для сжатия данных, представляющего сжатие данных, представляющих зависящий от времени сигнал.
Система 10 наблюдения состояния, соединенная с мотором 20, для наблюдения степени исправности мотора 20 изображена на Фиг. 1.
Система 10 наблюдения состояния содержит модуль 30 получения данных (в дальнейшем упоминаемый как "DAQ 30") и устройство 40 обработки сигналов в соответствии с одним вариантом осуществления изобретения. DAQ 30 получает зависящий от времени сигнал A(t) 45, который является непрерывным сигналом (характеризующимся независимой переменной времени "t"), который имеет отношение к вибрациям мотора 20. DAQ 30 обеспечивает данные D(nT) 65 в устройство 40 обработки сигналов. D(nT) 65 являются представлением дискретного времени (характеризующимся независимой переменной "nT") сигнала A(t) 45, которое впоследствии обрабатывается для сжатия D(nT) 65 для достижения цели настоящего изобретения.
DAQ 30 содержит датчик 31, модуль 32 обработки сигналов, и аналогово-цифровой преобразователь ADC 33 (ADC). Датчик 31 получает A(t) 45 и обеспечивает его в модуль 32 обработки сигналов. Модуль 32 обработки сигналов обрабатывает A(t) 45 и обеспечивает A'(t) в ADC 33. ADC 33 оцифровывает A'(t), тем самым создавая D(nT) 65, которые обеспечиваются в устройство 40 обработки сигналов, которое эффективно представляет A(t).
i
i
i
i
i
i
A(t) 45 может непрерывно получаться посредством DAQ 30 и может толковаться как содержащий множество зависящих от времени частичных сигналов A(t) (с i=1, 2, …, NoS) 46-50. A(t) 46-50 может толковаться как A(t) 45, полученный в течение множества индивидуальных временных интервалов Δt (с i=1, 2, …, NoS) 51-55, соответственно, или как A(t) 45, подвергнутый оконной обработке временной области, при этом интервалы окон соответствуют Δt 51-55. Альтернативно, DAQ 30 может получать A(t) 45 в течение Δt 51-55 посредством оконной обработки временной области, при этом A(t) 45 может толковаться как последовательный сбор A(t) 46-50.
i
i
i
Продолжительности Δt 51-55 могут определяться пользователем или могут быть фиксированными или переменными в зависимости от типа A(t) 45, требований обработки сигналов, характеристик DAQ 30 и/или устройства 40 обработки сигналов, и так далее. Продолжительности могут быть микросекундами, миллисекундами, или в секундах. Δt 51-55 являются, в общем, смежными. Однако Δt 51-55 могут быть перекрывающимися или разделенными некоторыми промежутками времени, и так далее.
i
Здесь, признак "частичные сигналы" определяется как часть A(t) 45, полученная посредством DAQ в течение соответствующего Δt 51-55.
i
i
i
Здесь "NoS" является безразмерной сущностью и указывает на количественное число сигналов A(t) 46-50, полученных в течение Δt 51-55, которые в полноте обеспечивают формирование A(t) 45. Т.е. NoS представляет количество частичных сигналов A(t).
i
i
i
i
1
1
Множество временных меток t (с i=1, 2, …, NoS) 56-60 определяются из A(t), при этом каждая временная метка t 56-60 представляет время начала соответствующего частичного сигнала A(t) 46-50, полученного в течение временного интервала Δt 51-55. Например, временная метка t 56 обозначает время начала частичного сигнала A(t) 46.
i
i
i
i
i
i
i
Данные D(nT) 65 являются представлением дискретного времени сигнала A(t) 45. Таким образом, A(t) 46-50 дает результатом множество частичных сигналов дискретного времени D(nT) (с i=1, 2, …, NoS; n=1, 2, …, NoP) 66-70 соответственно, при этом D(nT) 66-70 соответствуют D(nT) 65, присутствующему в Δt 51-55. D(nT) 66-70 также может получаться посредством оконной обработки области дискретного времени сигнала A(t) 45. Поэтому D(nT) 65 содержит представления дискретного времени множества частичных сигналов A(t) 51-55, т.е. D(nT) 66-70. Дополнительно, D(nT) 65 также может быть цифровым представлением эквивалента дискретного времени сигнала A(t).
i
i
Здесь, "NoP" является безразмерной сущностью и указывает на количественное число выборок, содержащихся в каждом частичном сигнале дискретного времени D(nT) 66-70. "NoP" также равно количественному числу частот, определимых из NoP-точечного дискретного преобразования Фурье (DFT), выполненного над D(nT) 66-70. Значимость признака "NoP" будет объяснена в последующих разделах.
i
i
"NoP" для D(nT) 66-70 может изменяться посредством изменения частоты дискретизации в ADC 33, т.е. посредством субдискретизации или передискретизации. "NoP" также может модифицироваться посредством дополнения нулями, или в другом альтернативном аспекте, может изменяться продолжительность Δt 51-55, чтобы оказывать влияние на изменение "NoP". Изменение "NoP" соответственно изменяет количественное число частот, определимых из NoP-точечного DFT.
i
i
i
i
i
i
i
i
i
i
i
Здесь, временные метки t 56-60 также представляют время начала для D(nT) 66-70 соответственно. Поэтому временные метки t 56-60 содержат информацию относительно времени получения D(nT) 66-70, т.е. времени начала Δt 51-55, и являются предпочтительными для сжатия D(nT) 65, также как для восстановления D(nT) 65. В альтернативном аспекте, t 56-60 могут быть любыми моментами времени соответствующих Δt 51-55, из которых вычисляется информация относительно времени получения сигнала A(nT) 66-70. Другими словами, временные метки t 56-60 представляют момент времени, в который соответствующий частичный сигнал A(t) и D(nT), соответственно, были получены.
Устройство 40 обработки сигналов, содержащее процессор 90 и блок 100 памяти для сжатия D(nT) 65, изображено на Фиг. 2. На Фиг. 2 иллюстративно изображено NoP=6 и NoS=5.
Процессор 90 содержит модуль 105 спектра, модуль 130 приема спектральных данных, модуль 140 нормализации амплитуд, модуль 150 сравнения, модуль 160 параметров, и модуль 170 корреляции. Вышеупомянутые модули 105, 130, 140, 150, 160, 170 сконфигурированы с возможностью сжимать данные D(nT).
i
i
Модуль 105 спектра принимает D(nT) 66-70 и вычисляет множество спектров S(f) (с i=1, 2, …, NoS) 106-110.
i
i
j
j,i
j
j,i
j
j
i
j,i
j
Модуль 105 вычисляет S(f) 106-110, которые могут быть множеством NoP-точечных (дискретных преобразований Фурье) DFT. S(f) 106-110 содержит множество частот f (с j=1, 2 …, NoP) 111-116 и множество амплитуд a (с j=1, 2 …, NoP; с i=1, 2 …, NoS) 121-126 частот f 111-116. Здесь, вычисленные амплитуды a 121-126 однозначно назначаются f 111-116 и являются мерами вкладов f 111-116 в A(t) 46-50. Это назначение амплитуд a 121-126 и частот f 111-116 находится в соответствии с хорошо известным определением спектра зависящего от времени сигнала, т.е. амплитуды являются мерами вклада назначенной частоты в лежащий в основе зависящий от времени сигнал.
i
j
j,i
j
Схематическая работа модуля 105 изображена как покомпонентное изображение "105". Здесь, для каждого из S(f) 106-110, горизонтальная ось "f" представляет "частоту" (f 111-116), и вертикальная ось "A" представляет "амплитуду", то есть амплитуды (a 121-126) частоты f 111-116 лежащего в основе зависящего от времени сигнала.
i
i
i
i
i
i
i
Здесь объясняется иллюстративный способ для реализации работы модуля 105. Например, D(nT) 66-70 принимаются устройством 40 обработки сигналов и они буферизуются как блоки данных конкретной длины в блоке 100 памяти. Каждый блок данных соответствует каждому из D(nT) 66-70, полученному в течение каждого из Δt 51-55. Каждый блок данных обеспечивается временной меткой с помощью соответствующего t 56-60 для идентификации времени начала D(nT) 66-70. Блоки, соответствующие D(nT) 66-70, затем извлекаются и обрабатываются для вычисления S(f) 106-110.
j,i
j
i
j
i
Здесь, a 121-126 и f 111-116 спектра S(f) 106-110 могут быть отображены как двухмерная матрица 128 амплитуд, расположенная в строках и столбцах, например, с NoP строками (представляющими NoP частот f 111-116) и NoS столбцами (представляющими NoS сигналов D(nT) 66-70).
i
i
Каждый из S(f) 106-110 вычисляется из соответствующего D(nT) 66-70. Поэтому "NoS" также обозначает "количество спектров".
i
j,i
j
i
j,i
i
Модуль 130 приема спектральных данных принимает спектры S(f) 106-110 от модуля 105 спектра. Модуль 130 затем размещает вектор из a 121-126 и вектор из f 111-116, присутствующих в S(f) 106-110. Таким образом, вектор a 121-126 назначается соответствующей временной метке t 56-60.
i
i
j
1
j,i
j
1
1
j,i
1
i
В матрице 128 амплитуд, каждому столбцу "i" из NoS столбцов назначена соответствующая временная метка t 56-60, которая соответствует соответствующему D(nT) 66-70. Аналогично, каждой строке "j" из NoP строк назначена соответствующая частота f 111-116. Например, первому столбцу назначена временная метка t 56 и он содержит все амплитуды a 121-126 частот f 111-116, вычисленные из S(f) 106. Аналогично, первой строке назначена f 111 и она содержит амплитуды a 121-126 частоты f 111, вычисленные из количества NoS спектров S(f) 106-110.
i
j
1,1
6,1
1
6
1
1
1
3,1
3,5
3
1
5
1
5
1
5
Покомпонентное изображение "130" иллюстрирует пример матрицы 128 амплитуд для пяти спектров S(f) 106-110 (NoS=5) и шести частот f 111-116 (NoP=6). Например, значения a 121 по a 126 в столбце представляют значения шести частот f 111 по f 116, определенные из S(f) 106, вычисленного для D(nT) 66, несущего временную метку t 56. Значения a 123 по a 123 в строке представляют значения третьей частоты "f 113", определенные из S(f) 106 по S(f) 110, вычисленных соответственно для D(nT) 66 по D(nT) 70, несущих соответствующие временные метки t 56 по t 60.
i
i
j,i
j
Дополнительно, для каждого из S(f) 106-110 модуль 130 вычисляет алгебраическую сумму SF (с i=1 по NoS) 135 амплитуд a 121-126 частот f 111-116.
Например,
, , и так далее.
i
i
i
1
1
3
3
Таким образом, SF 135 также является вектором длины "NoS". Дополнительно, каждая из SF 135 может однозначно соответствовать каждой из NoS временных меток t 56-60. Например, SF соответствует t, SF соответствует t, и так далее.
i
j,i
j,i
j,i
Модуль 140 нормализации амплитуд принимает матрицу 128 амплитуд от модуля 130 приема спектральных данных. Для каждого из S(f) 106-110, амплитуды a 121-126 нормализуются посредством модуля 140. Таким образом, получается множество нормализованных амплитуд b (с j=1, 2 …, NoP; с i=1, 2 …, NoS) 141-146. Каждая из b 141-146 вычисляется в соответствии со следующей формулой:
Т.е.
Например,
, , и так далее.
j,i
i
i
Здесь, b 141-146 также располагаются в матрице 148 нормализованных амплитуд по отношению к S(f) 106-110, несущим временные метки t 56-60, соответственно. Матрица 148 нормализованных амплитуд снова содержит NoP строк и NoS столбцов. Схематическая работа модуля 140 изображена как покомпонентное изображение "140".
i
j,i
j
j,i
j
i
j
j
В вышеупомянутом контексте, признак "нормализация" определяется как, для каждого S(f) 106-110, обработка деления индивидуальных амплитуд a 121-126 частоты f 111-116 на алгебраическую сумму амплитуд a 121-126. Поэтому посредством нормализации, могут измеряться относительные вклады или процентные вклады каждой из f 111-116, присутствующей в каждом из A(t) 46-50. Таким образом, вклад некоторой частоты f 111-116 может оцениваться и может измеряться для определения того, как эта конкретная частота f 111-116 может представляться для сжатия D(nT) 65, что будет объяснено далее.
j,i
j
i
j,i
NoP строк и NoS столбцов матрицы 128 амплитуд и матрицы 148 нормализованных амплитуд всего лишь иллюстрируют и представляют способ, в котором принятые a 121-126 частот f 111-116 каждого из S(f) 106-110 и соответствующие b 141-146 могут размещаться для обработки информации гибким способом. Эти NoP строк и NoS столбцов могут транспонироваться и информация там может представляться и обрабатываться соответственно без потери общности.
Вышеупомянутый пример (NoP=6 и NoS=5) является всего лишь иллюстративным, и не должен толковаться как ограничивающий изобретение. NoP и NoS могут быть любыми двумя положительными целыми числами и величины NoP и NoS могут быть огромными. Например, для 256-точечного FFT, вычисленного для 300 спектров, "NoP" равняется 256, и "NoS" равняется 300.
j
i
i
i
i
j
i
i
i
Здесь, количественное число частот f 111-116, т.е. NoP, определимое из S(f) 106-110, равняется количественному числу выборок, т.е. снова NoP, присутствующих в каждом из D(nT) 66-70. Например, количественное число равняется "16" для 16-точечного DFT, "32" для 32-точечного DFT, и так далее. Однако для вычисления S(f) 106-110 из D(nT) 66-70 эффективным способом, т.е. для настройки количества частот f 111-116, определимых из каждого из S(f) 106-110, D(nT) 66-70 могут дискретизироваться с повышением, дискретизироваться с понижением, дополняться нулями, и так далее, как упомянуто выше без потери общности для D(nT) 66-70.
i
j
i
j,i
i
j
i
Здесь, каждый S(f) 106-110 содержит один и тот же вектор частот f 111-116, т.е. другой S(f) 106-110 может отличаться только в векторе амплитуд a 121-126, вычисленных для другого D(nT) 66-70. Здесь, каждая частота f 111-116 содержит один и тот же вектор временных меток t 56-60.
i
j,i
v
Модуль 150 сравнения принимает матрицу 148 нормализованных амплитуд. Для каждого спектра S(f) 106-110, модуль 150 сравнивает нормализованные амплитуды b 141-146 с пороговым значением T.
v
j
v
T пропорционален обратному "NoP", т.е. количественному числу (NoP) частот f 111-116. Здесь, отношение между T и NoP может быть показано следующим образом:
j,i
i
v
j,i
j
v
Посредством сравнения каждой нормализованной амплитуды b 141-146 каждого спектра S(f) 106-110 с пороговым значением T, модуль 150 сравнения генерирует распределение нормализованных амплитуд b 141-146 частоты f по отношению к T.
j,i
j,i
v
j,i
i
После этого, посредством обработки матрицы 148 нормализованных амплитуд, содержащей b 141-146, и распределения b 141-146 по отношению к T в модуле 160 параметров, является возможным моделировать распределение b 141-146 по S(f) 106-110.
j,i
j,i
j
v
Модуль 160 параметров принимает матрицу 148 нормализованных амплитуд, содержащую b 141-146, и распределение b 141-146 для каждой частоты f 111-116 по отношению к T.
j
j,i
j,i
v
j
j
j,i
1,1
1,6
1
1
4,1
4,6
4
4
Для каждой частоты f 111-116, нормализованная амплитуда b 141-146 и распределение b 141-146 по отношению к T, соответственно, обрабатываются, чтобы получать один или более параметров модели (MP) (с j=1 по NoP). Параметры модели MP являются подходящими для моделирования распределения нормализованных амплитуд b 141-146. Например, b по b 141 обрабатываются для получения MP для моделирования f 111, и b по b 144 обрабатываются для получения MP для моделирования f 114, и так далее.
j
j
j,i
j
v
j
j
j,i
Для каждой из f 111-116, MP зависит от тенденции распределения b 141-146 частоты f 111-116 по отношению к T. MP вычисляется посредством моделирования частот f 111-116 на основе разных тенденций распределения соответствующих b 141-146.
Для этого, выделяются три разных случая:
j,i
j
j,i
v
j
j,i
j
j,i
j
i
j
j,i
j
i
Случай 1) В случае, когда тенденция распределения соответствующих b 141-146 конкретной частоты f 111-116 является такой, что каждая из b 141-146 меньше, чем T, MP для представления нормализованных амплитуд b этой конкретной частоты f 111-116 является средним значением (AV) нормализованных амплитуд b 141-146 конкретной частоты f 111-116 спектров S(f). Среднее значение (AV) для этой конкретной частоты f 111-116 вычисляется посредством обработки нормализованных амплитуд b этой конкретной частоты f 111-116 по S(f) 106-110, как показано ниже:
1
1,1
1,5
v
1
Например, если для f 111, если все из b по b 141 меньше, чем T, то AV для f 111 вычисляется, как показано ниже:
j,i
j
j,i
j,i
j
v
j
j,i
j
j,i
j
Случай 2) В случае, когда тенденция распределения нормализованных амплитуд b 141-146 частоты f 111-116 является такой, что разность  между экстремальными значениями (т.е. разность Δ между максимальным значением  и минимальным значением  для b) нормализованных амплитуд b этой конкретной частоты f 111-116 меньше, чем T, то MP для представления нормализованных амплитуд b 141-146 этой конкретной частоты f 111-116 является опять таки средним значением (AV) нормализованных амплитуд b 141-146 конкретной частоты f 111-116, как вычислено выше.
5
5,4
5,5
5,4
5,5
v
5
Например, если для f 115, если b является максимальным значением  и b является минимальным значением , и если разность Δ между b и b меньше, чем T, то AV для f 115 вычисляется, как показано ниже:
j,i
j
j,i
j
j
Случай 3) В оставшихся случаях, если тенденция распределения b 141-146 конкретной частоты f 111-116 отличается от тенденций вышеупомянутых распределений в случаях 1) и 2), распределение нормализованных амплитуд b 141-146 моделируется и представляется посредством алгоритма аппроксимации. Здесь, MP этой конкретной частоты f 111-116 являются результирующими параметрами модели алгоритма аппроксимации. Алгоритм аппроксимации может быть одним из следующих:
j
j,i
линейная подгонка, в которой результирующий MP, определяющий линейное уравнение для моделирования b 141-146, содержит значение наклона, и значение пересечения оси;
j
j,i
квадратичная подгонка, в которой результирующие MP, определяющие квадратичное уравнение для моделирования b 141-146, являются одним или более коэффициентами, представляющими переменную квадратичного уравнения;
j
j,i
кубическая подгонка, в которой результирующие MP, определяющие кубическое уравнение для моделирования b 141-146, являются одним или более коэффициентами, представляющими переменную кубического уравнения;
j
j,i
полиномиальная подгонка n-ого порядка, в которой результирующие MP, определяющие полиномиальное уравнение n-ого порядка для моделирования b 141-146, являются одним или более коэффициентами, представляющими переменную полиномиального уравнения n-ого порядка;
j
j,i
ступенчатая функция, где результирующий MP, определяющий ступенчатую функцию для моделирования b 141-146, является кумулятивной функцией, определенной как сумма элементарных ступенчатых функций; или
j
j,i
итеративный алгоритм подгонки на основе конечных точек, где результирующие MP, определяющие итеративный алгоритм подгонки на основе конечных точек для моделирования b 141-146, могут получаться из алгоритма Рамера - Дугласа - Пекера, и так далее.
j,i
j
j
Таким образом, в зависимости от соответствующих распределений нормализованных амплитуд b 141-146 для разных частот f 111-116, параметры модели MP вычисляются, в конце концов давая результатом данные, которые являются подходящими для использования в качестве основы для генерирования набора сжатых данных CDS для представления исходных данных D(nT) 65.
CDS, представляющий данные D(nT) 65, содержит следующее:
j
вектор частот f 111-116,
i
вектор временных меток t 56-60,
i
j,i
i
j
j
j
j,i
алгебраическую сумму SF 135 амплитуд a 121-126 каждого спектра S(f), что является предпочтительным для масштабирования MP для каждой из f 111-116, так как MP был создан с использованием только b 141-146, и
j
j,i
j
параметры модели MP для представления тенденций распределения нормализованных амплитуд b 141-148 каждой частоты f 111-116.
j
j
j,i
1
Более того, для тех частот f 111-116, для которых соответствующий параметр модели MP выбирается, чтобы быть средним значением (AV) нормализованных амплитуд b 141-146 (т.е. случаи 1) и 2), как описано выше), дополнительно первая временная метка t 56 назначается AV.
j
j
j
В общем, для всех других случаев (т.е. случаев, нежели случаи 1) и 2)) и для всех частот, соответственно, CDS должен содержать информацию о том, является ли MP средним значением (случаи 1) и 2)), или был смоделирован посредством алгоритма аппроксимации (случай 3)). Если MP для конкретной частоты f 111-116 основывается на алгоритме аппроксимации, то CDS дополнительно содержит следующее:
j,i
j
информацию о типе алгоритма аппроксимации (например, ступенчатая функция, итеративный алгоритм подгонки на основе конечных точек, и так далее), использованного для приближения нормализованных амплитуд b 141-148 конкретной частоты f 111-116, и
j
j
i
j
j,i
j,i
j
MP алгоритма аппроксимации для этой конкретной частоты f 111-116 также будет содержать пары из одного или более значений амплитуд и соответствующих одной или более временных меток t 56-60 этой конкретной частоты f 111-116. Здесь, значения амплитуд могут быть самими одной или более нормализованными амплитудами b 141-146 или средним некоторого количества нормализованных амплитуд b 141-146 этой конкретной частоты f 111-116.
Посредством обработки CDS, является возможным, по существу, восстанавливать данные D(nT).
i
Чтобы оценить качество сжатия перед окончательным сохранением CDS, сигнал данных D'(nT) восстанавливается из набора сжатых данных CDS и сравнивается с исходным сигналом D(nT). Для оценки, модуль 170 корреляции принимает набор сжатых данных CDS и обрабатывает его для построения сигнала данных D'(nT) для соответствующего Δt 51-55.
Восстановление D'(nT) посредством обработки CDS объясняется иллюстративным образом со ссылкой на Фиг. 4.
Модуль 170 корреляции определяет коэффициент корреляции CC посредством сопоставления D'(nT) и D(nT). CC является безразмерной сущностью и является просто качественным индексом для представления сходства между D'(nT) и D(nT).
v
j,i
j
j
На основе определенного CC, можно принимать решение повторять сжатие с другими предварительными условиями. Например, если коэффициент корреляции показывает, что сходство D'(nT) и D(nT) не является достаточным, пороговое значение T должно изменяться. Это может давать результатом измененную тенденцию распределения нормализованных амплитуд b 141-146 частоты f 111-116 и, таким образом, новые параметры модели MP' для этой частоты.
v
v
v
f
v
f
Пороговое значение T может изменяться, чтобы получать измененное пороговое значение T' посредством умножения T на коэффициент допустимого отклонения T. Здесь, отношение между T', NoP и T может быть показано следующим образом:
f
v
j,i
v
j
j
v
f
j
T может изменяться, чтобы изменять T. Фактически, это дает результатом изменение распределения b 141-146 по отношению к T. Тем самым, это дает результатом изменение MP для каждой f 111-116. Например, если CC между D'(nT) и D(nT) находится между "0" и "0,8", можно предполагать, что качество восстановленного сигнала D'(nT) является плохим, и качество может быть улучшено посредством понижения T, т.е. посредством уменьшения T, тем самым выбирая другой MP.
j
Аналогично, если CC равняется "0,8" или больше, можно предполагать, что MP для восстановления D'(nT) являются достаточными. Затем текущий CDS рассматривается как конечный CDS и он сохраняется в памяти.
v
v
f
f
v
f
v
f
v
f
T, T' и T являются безразмерными сущностями. T просто служит, чтобы увеличивать или уменьшать T. Если T выбирается между "0" и "1", то T уменьшается. Альтернативно, если T больше, чем "1", то T увеличивается. Однако предпочтительно T должен быть неравным "1".
v
v
j,i
v
v
v
После изменения порогового значения, заново определенное пороговое значение T' посылается в модуль 150 сравнения. Там, пороговое значение T, которое используется для сравнения с нормализованными амплитудами b 141-146, устанавливается на новое пороговое значение T', т.е. T=T'.
j,i
i
v
j,i
v
v
v
После этого, выполняются такие же этапы, как описано выше, для определения распределения нормализованных амплитуд b 141-146 каждого спектра S(f) 106-110 по отношению к пороговому значению T. Т.е. Модуль 150 сравнения снова генерирует распределение b 141-146 по отношению к T, при этом T соответствует измененному пороговому значению T'.
j,i
i
После определения распределения, снова является возможным моделировать распределение b 141-146 по S(f) 106-110.
j,i
v
Модуль 160 параметров принимает информацию относительно сравнения между каждой из b 141-146 с T, т.е. информацию о вышеупомянутом распределении.
j
j,i
j,i
v
j
j
j
j,i
Снова, для каждой частоты f 111-116, нормализованная амплитуда b 141-146 и распределение b 141-146 по отношению к T, соответственно, обрабатываются, чтобы получать один или более параметров модели MP, как описано выше. MP снова вычисляется посредством моделирования частот f 111-116 на основе разных тенденций распределения соответствующих b 141-146, снова рассматривая вышеупомянутые три разных случая 1), 2) и 3).
j
j
После вычисления параметров модели MP для каждой частоты f 111-116, генерируется новый набор сжатых данных CDS, который содержит:
j
вектор частот f 111-116,
i
вектор временных меток t 56-60,
i
j,i
i
j
j
j
j,i
алгебраическую сумму SF 135 амплитуд a 121-126 каждого спектра S(f), что является предпочтительным для масштабирования MP для каждой из f 111-116, так как MP был создан с использованием только b 141-146, и
j
j,i
j
параметры модели MP для представления тенденций распределения нормализованных амплитуд b 141-148 каждой частоты f 111-116.
j
j
j,i
1
Более того, для тех частот f 111-116, для которых соответствующий параметр модели MP выбирается, чтобы быть средним значением (AV) нормализованных амплитуд b 141-146 (т.е. случаи 1) и 2), как описано выше), дополнительно первая временная метка t 56 назначается AV.
j
j
j
В общем, для всех других случаев (т.е. случаев, нежели случаи 1) и 2)) и для всех частот, соответственно, CDS должен содержать информацию о том, является ли MP средним значением (случаи 1) и 2)), или был смоделирован посредством алгоритма аппроксимации (случай 3)). Если MP для конкретной частоты f 111-116 основывается на алгоритме аппроксимации, то CDS дополнительно содержит следующее:
j,i
j
информацию о типе алгоритма аппроксимации (например, ступенчатая функция, итеративный алгоритм подгонки на основе конечных точек, и так далее), использованном для приближения нормализованных амплитуд b 141-148 конкретной частоты f 111-116, и
j
j
i
j
j,i
j,i
j
MP алгоритма аппроксимации для этой конкретной частоты f 111-116 также будет содержать пары из одного или более значений амплитуд и соответствующих одной или более временных меток t 56-60 этой конкретной частоты f 111-116. Здесь, значения амплитуд могут быть самими одной или более нормализованными амплитудами b 141-146 или средним некоторого количества нормализованных амплитуд b 141-146 этой конкретной частоты f 111-116.
k
k
k
Частоты, для которых параметры модели были вычислены с использованием алгоритмов аппроксимации, являются f и соответствующие временные метки, сохраненные в параметрах модели для f, являются t.
Чтобы обеспечивать достаточное качество сжатия, как описано выше, сигнал данных D'(nT) снова восстанавливается из набора сжатых данных CDS и сравнивается с исходным сигналом D(nT). Модуль 170 корреляции определяет новый коэффициент корреляции CC посредством сопоставления D'(nT) и D(nT).
v
f
j
На основе определенного CC, можно принимать решение повторять сжатие снова с другими предварительными условиями. В этом случае, вышеупомянутая процедура, начинающаяся с изменения порогового значения T с другим коэффициентом допустимого отклонения T, повторяется, давая результатом новое распределение и новые параметры модели MP и т.д.
j,i
j
В заключение, когда CC показывает достаточное сходство между D'(nT) и D(nT), соответствующий набор сжатых данных CDS сохраняется. Так как вместо всех амплитуд a сохраняются только средние значения и/или параметры модели MP, может достигаться эффективное отношение сжатия.
Модуль 105 спектра и модуль 130 приема спектральных данных могут быть объединены, чтобы формировать одиночный модуль с возможностью выполнения функций обоих модулей 105, 130.
Модуль 150 сравнения, модуль 160 параметров, и модуль 180 корреляции могут быть объединены, чтобы формировать одиночный модуль с возможностью выполнения функций модулей 150, 160, 170.
i
i
j,i
j
i
Блок 100 памяти выполнен с возможностью хранения D(nT) 66-70, принятых от DAQ 30, промежуточных данных, полученных в течение разных этапов обработки D(nT) 66-70, посредством вышеупомянутых модулей 105, 130, 140, 150, 160, 170, a 121-126 и f 111-116, вычисленных из S(f) 106-110, и так далее.
Упомянутые один или более вышеупомянутых модулей 105, 130, 140, 150, 160, 170 работоспособным образом соединены с процессором 90 и могут реализовываться как независимые модули или как частично объединенные модули или полностью объединенный модуль, при этом процессор 90 сконфигурирован с возможностью соответственно для выполнения соответствующих функций упомянутых одного или более вышеупомянутых модулей 105, 130, 140, 150, 160, 170. Дополнительно, упомянутые один или более модули 105, 130, 140, 150, 160, 170, могут быть соединенными и могут располагаться внутри или вне процессора 90.
Процессор 90 может быть процессором общего назначения, микроконтроллером, цифровым сигнальным процессором, программируемой пользователем вентильной матрицей (FPGA), FPGA с частично динамически перестраиваемой конфигурацией, специализированной интегральной схемой, и комбинацией перечисленного.
Достаточными модулями для достижения цели изобретения являются модуль 130, модуль 140, модуль 150 и модуль 160. Однако устройство обработки сигналов из Фиг. 2 обеспечено дополнительными модулями, т.е. блоком 100 памяти, модулем 105, и модулем 170.
j,i
j,i
v
i
Пример матрицы 128 амплитуд и матрицы 148 нормализованных амплитуд, где NoP=4 и NoS=10, содержащих иллюстративные значения a 121-126 и b 141-146 соответственно, проиллюстрированы на Фиг. 3. Так как NoP равняется "4", соответствующий T равняется "0,25". Столбцы матрицы 128 амплитуд и матрицы 148 нормализованных амплитуд несут соответствующие временные метки t 56-60.
1
2
3
4
1
10
1
2
3
4
5
6
7
8
9
10
1
10
1
10
Здесь, для f, матрица 128 амплитуд содержит иллюстративные значения амплитуд [5, 6, 6, 7, 5, 6, 1, 2, 5, 3], для f, матрица 128 амплитуд содержит иллюстративные значения амплитуд [54, 56, 52, 54, 52, 58, 59, 55, 53, 51], для f, матрица 128 амплитуд содержит иллюстративные значения амплитуд [10, 20, 30, 10, 20, 30, 0, 20, 30, 30], и для f, матрица 128 амплитуд содержит иллюстративные значения амплитуд [20, 5, 25, 30, 23, 5, 24, 6, 23, 5], полученные из иллюстративных десяти спектров S(f) по S(f). Соответствующим образом, SF равняется 89, SF равняется 87, SF равняется 113, SF равняется 101, SF равняется 100, SF равняется 99, SF равняется 84, SF равняется 83, SF равняется 111, и SF равняется 89. Каждый столбец матрицы 128 амплитуд назначен соответствующей временной метке t по t, сгенерированной из соответствующих D(nT) по D(nT).
1
2
3
4
Здесь, для f, матрица 148 нормализованных амплитуд содержит иллюстративные значения нормализованных амплитуд [0,06, 0,07, 0,05, 0,07, 0,05, 0,06, 0,01, 0,02, 0,05, 0,03], для f, матрица 148 нормализованных амплитуд содержит иллюстративные значения нормализованных амплитуд [0,61, 0,64, 0,46, 0,53, 0,52, 0,59, 0,70, 0,66, 0,48, 0,57], для f, матрица 148 нормализованных амплитуд содержит иллюстративные значения нормализованных амплитуд [0,11, 0,23, 0,27, 0,10, 0,20, 0,30, 0,00, 0,24, 0,27, 0,34], и для f, матрица 148 нормализованных амплитуд содержит иллюстративные значения нормализованных амплитуд [0,22, 0,06, 0,22, 0,30, 0,23, 0,05, 0,29, 0,07, 0,21, 0,06].
1
v
1
1
1
Для f, можно наблюдать отсюда, что каждое из значений нормализованных амплитуд [0,06, 0,07, 0,05, 0,07, 0,05, 0,06, 0,01, 0,02, 0,05, 0,03] ниже T. Поэтому MP, соответствующий f, является средним значением (AV) для [0,06, 0,07, 0,05, 0,07, 0,05, 0,06, 0,01, 0,02, 0,05, 0,03], которое равняется "0,047", и первой временной меткой t.
2
v
v
2
2
1
Для f, можно наблюдать отсюда, что каждое из значений нормализованных амплитуд [0,61, 0,64, 0,46, 0,53, 0,52, 0,59, 0,70, 0,66, 0,48, 0,57] выше T, также как разность между экстремальными значениями ("0,70" и "0,46") равняется "0,24", что ниже T. Поэтому MP, соответствующий f, является средним значением (AV) для [0,61, 0,64, 0,46, 0,53, 0,52, 0,59, 0,70, 0,66, 0,48, 0,57], которое равняется "0,576", и первой временной меткой t.
3
v
v
v
3
3
3
1
6
7
8
10
Для f, можно наблюдать отсюда, что некоторые из значений нормализованных амплитуд [0,11, 0,23, 0,27, 0,10, 0,20, 0,30, 0,00, 0,24, 0,27, 0,34] выше T, тогда как другие ниже T. Также, разность между экстремальными значениями ("0,30" и "0,00") равняется "0,30", что выше T. Поэтому MP, соответствующий f, получаются из алгоритма Рамера - Дугласа - Пекера, при этом MP составляет [(0,11, t), (0,3, t), (0, t), (0,24, t), (0,34, t)].
4
v
v
v
3
4
3
1
2
4
6
7
8
9
10
Для f, можно наблюдать отсюда, что некоторые из значений нормализованных амплитуд [0,22, 0,06, 0,22, 0,30, 0,23, 0,05, 0,29, 0,07, 0,21, 0,06] выше T, тогда как другие ниже T. Также, разность между экстремальными значениями ("0,30" и "0,05") равняется "0,25", что равно T. Поэтому MP, соответствующий f, получаются из алгоритма Рамера - Дугласа - Пекера, при этом MP составляет [(0,22, t), (0,06, t), (0,3, t), (0,05, t), (0,29, t), (0,07, t), (0,21, t), (0,06, t)].
Поэтому CDS содержит следующее:
1
10
вектор частот f по f,
1
10
вектор временных меток t по t,
1
10
1
10
алгебраическую сумму SF по SF амплитуд каждого спектра S(f) по S(f), т.е. [89, 87, 113, 101, 100, 99, 84, 83, 111, 89], и
параметры модели (MP), как показано ниже:
1
1
MP: [(0,047, t)]
2
1
MP: [(0,576, t)]
3
1
6
7
8
10
MP: [(0,11, t), (0,3, t), (0, t), (0,24, t), (0,34, t)], и
4
1
2
4
6
7
8
9
10
MP: [(0,22, t), (0,06, t), (0,3, t), (0,05, t), (0,29, t), (0,07, t), (0,21, t), (0,06, t)].
3
4
k
1
6
7
8
10
k
3
1
2
4
7
8
9
10
k
4
Здесь, f и f составляют f, и временные метки t, t, t, t и t составляют t для MP и временные метки t, t, t, t, t, t и t составляют t для MP.
j
j
j
Фиг. 4 изображает иллюстративную матрицу восстановленных амплитуд для восстановления D'(nT) на основе вышеупомянутого CDS. Для каждой частоты f восстановленная амплитуда частоты f зависит от типа параметра модели MP.
j
j
j
1
j,i
j
i
j
i
j,i
j
1
1
NoS-1
NoS
В вышеупомянутых случае 1) и случае 2), если частота f представлена посредством соответствующего MP, который содержит только среднее значение нормализованных амплитуд частоты f и первую временную метку t, то каждая восстановленная амплитуда a' частоты f' для временной метки t в матрице восстановленных амплитуд будет генерироваться посредством умножения среднего значения нормализованных амплитуд частоты f на соответствующую SF. То есть, восстановленные амплитуды a' частоты f' будут следующими: [AV*SF, AV*SF, ----, AV*SF, AV*SF].
j
j
j,i
j
j
j,i
j
i
В вышеупомянутом случае 3), если частота f представлена посредством MP, полученного из алгоритма аппроксимации, например, алгоритма Рамера - Дугласа - Пекера, то каждая восстановленная амплитуда a' частоты f может генерироваться посредством передачи MP в модуль обратного алгоритма аппроксимации, в этом случае модуль обратного алгоритма Рамера - Дугласа - Пекера для получения каждой из восстановленной амплитуды a' частоты f' для каждой из временной метки t.
j
i
i
Как упомянуто выше, CDS обрабатывается для восстановления матрицы 180 восстановленных нормализованных амплитуд и впоследствии матрицы 190 восстановленных амплитуд частот f. Посредством вычисления обратного DFT каждого из столбцов матрицы 190 восстановленных амплитуд является возможным построить D'(nT), т.е. восстановленные частичные сигналы. Посредством непрерывного и последовательного размещения D'(nT) по отношению к временным меткам, строится D'(nT).
Фиг. 5 изображает резюмирующую блок-схему последовательности операций способа для сжатия D(nT) 65, представляющих A(t) 45.
i
i
i
На этапе 200, принимаются сигналы D(nT) 66-70 с t 56-60 и вычисляются соответствующие S(f) 106-110.
i
j,i
j
i
i
i
i
i
На последующем этапе 210, для каждого из количества NoS спектров S(f) 106-110, количество NoP амплитуд a 121-126, каждая из которых соответствует каждой из f 111-116, принимаются и располагаются, как упомянуто выше. Каждый из S(f) 106-110 обеспечивается временной меткой с помощью соответствующего t 56-60. Таким образом, создается матрица 128 амплитуд с NoP строк и NoS столбцов, при этом каждый столбец назначен временной метке t 56-60. Матрица 128 амплитуд может создаваться при обеспечении выводов модуля 105 в модуль 130. Дополнительно, для каждого из S(f) 106-110, вычисляется соответствующая SF 135, как упомянуто выше.
j,i
i
j,i
j,i
На последующем этапе 220, каждый столбец матрицы 128 амплитуд нормализуется, чтобы получать матрицу 148 нормализованных амплитуд, как раскрыто в предшествующих разделах. Матрица 128 амплитуд может обеспечиваться в модуль 140 для получения b 141-146. Для каждого из S(f) 106-110, соответствующие a 121-126 нормализуются, чтобы получить b 141-146.
j,i
i
v
j
i
j,i
i
v
j,i
v
На последующем этапе 230, каждая из b 141-146 каждого из S(f) 106-110 сравнивается с T, который пропорционален обратному "NoP", т.е. количественного числа частот f 111-116 каждого из S(f) 106-110. Матрица 148 нормализованных амплитуд может обеспечиваться в модуль 150 для сравнения каждой из b 141-146 каждого из S(f) 106-110 с T. Таким образом, распределение b 141-146 по отношению к T может получаться и иллюстрироваться.
j
j,i
j
i
j
j
На последующем этапе 240, для каждой из f 111-116, распределение b 141-146 частоты f 111-116 по S(f) 106-110 принимается для вычисления соответствующего MP для соответствующей частоты f 111-116, как описано выше.
На этапе 250, генерируется набор сжатых данных (CDS), как определено выше.
f
v
v
v
j,i
i
v
v
На этапе 260, D'(nT) восстанавливается из CDS. На последующем этапе 270, D'(nT) сопоставляется с D(nT) 65 для получения коэффициента корреляции CC. На основе определенного CC, как объяснено в предшествующих разделах, T может изменяться, чтобы изменять пороговое значение T. В случае, когда CC показывает, что сходство между D'(nT) и D(nT) не является достаточным, генерируется новое пороговое значение T' посредством изменения последнего порогового значения T. Затем, процедура возвращается на этап 230, при этом каждая из нормализованных амплитуд b 141-146 каждого из спектра S(f) 106-110 сравнивается с пороговым значением T=T'. В случае, когда CC показывает, что сходство между D'(nT) и D(nT) является достаточным, процедура продолжается на этапе 280.
f
Коэффициент корреляции CC является обратно пропорциональным коэффициенту допустимого отклонения T.
На этапе 280, текущий набор сжатых данных CDS сохраняется или, в зависимости от приложения, дополнительно обрабатывается.
В заключение, D(nT) 65 сжимается, т.е. пространство хранения для хранения информации (амплитуды, частоты и времени), содержащейся в D(nT) 65, в значительной степени уменьшается в силу хранения только CDS в блоке 100 памяти. Впоследствии, является возможным использовать его либо для ссылки, либо также для передачи.
В вышеупомянутом способе для сжатия D(nT) 65, этап 200 и этапы 260, 270 являются необязательными этапами для достижения цели изобретения.
Хотя изобретение было описано со ссылкой на конкретные варианты осуществления, это описание не предназначено, чтобы толковаться в ограничивающем смысле. Различные примеры раскрытых вариантов осуществления, также как альтернативные варианты осуществления изобретения, должны быть ясными для специалистов в данной области техники при обращении к описанию изобретения. Поэтому предполагается, что такие модификации могут делаться без отхода от вариантов осуществления настоящего изобретения, как определены.
СПИСОК ССЫЛОЧНЫХ ПОЗИЦИЙ
10 Система наблюдения состояния
20 Мотор
30 Модуль получения данных
31 Датчик
32 Модуль обработки сигналов
33 Аналогово-цифровой преобразователь
40 Устройство обработки сигналов
45 Зависящий от времени сигнал A(t)
i
46-50 Множество зависящих от времени частичных сигналов A(t)
i
51-55 Множество временных интервалов Δt
i
56-60 Множество моментов времени t (временных меток)
65 Данные D(nT)
i
66-70 множество частичных сигналов дискретного времени D(nT)
90 Процессор
100 Блок памяти
105 Модуль спектра
i
106-110 Множество спектров S(f)
j
111-116 Множество частот f
j,i
121-126 Множество амплитуд a
128 Матрица амплитуд
130 Модуль приема спектральных данных
i
131-135 Сумма частот SF
140 Модуль нормализации амплитуд
j,i
141-146 Множество нормализованных амплитуд b
148 Матрица нормализованных амплитуд
150 Модуль сравнения
160 Модуль параметров
170 Модуль корреляции
180 Матрица восстановленных нормализованных амплитуд
190 Матрица восстановленных амплитуд
i
i
200 Этап приема D(nT) и вычисления S(f)
210 Этап создания матрицы амплитуд
220 Этап создания матрицы нормализованных амплитуд
230 Этап сравнения матрицы нормализованных амплитуд с порогом
j,i
j
240 Этап обработки нормализованных амплитуд b для вычисления параметров модели MP
250 Этап генерирования набора сжатых данных
260 Этап построения D'(nT)
270 Этап сопоставления D'(nT) и D(nT)
280 Этап сохранения набора сжатых данных