|
|
  |
Выделение постоянной из смеси пост + экспонента |
|
|
|
Jul 10 2014, 20:59
|
Участник

Группа: Участник
Сообщений: 21
Регистрация: 2-12-04
Пользователь №: 1 296

|
Подскажите пожалуйста в какую сторону ЦОС смотреть. Есть сигнал вида s(t) = A*(1+B*e^(-t/T)) , отсчеты с АЦП Необходимо найти A , не дожидаясь хвоста экспоненты. Сейчас нахожу Т по трем точкам , вычитаю предсказанный сигнал, в идеале все работает, с реальными сигналами точность никакая, все скачет , и Т в том числе.
|
|
|
|
|
Jul 11 2014, 07:16
|
Гуру
     
Группа: Модераторы
Сообщений: 8 752
Регистрация: 6-01-06
Пользователь №: 12 883

|
Цитата(IVANS @ Jul 11 2014, 00:59)  Подскажите пожалуйста в какую сторону ЦОС смотреть. Есть сигнал вида s(t) = A*(1+B*e^(-t/T)) , отсчеты с АЦП Необходимо найти A , не дожидаясь хвоста экспоненты. Сейчас нахожу Т по трем точкам , вычитаю предсказанный сигнал, в идеале все работает, с реальными сигналами точность никакая, все скачет , и Т в том числе. А T и B Вы заранее знаете? Сначала Вы пишете, что нужно найти А, а потом про скачки Т. Как Вас следует правильно понимать?
|
|
|
|
|
Jul 11 2014, 13:08
|
Участник

Группа: Участник
Сообщений: 21
Регистрация: 2-12-04
Пользователь №: 1 296

|
Спасибо за подсказки по МНК, похоже , что это то, что нужно, но пока разбираюсь как это дело применить. Может есть рекоммендации? Цитата(Tanya @ Jul 11 2014, 11:16)  А T и B Вы заранее знаете? Сначала Вы пишете, что нужно найти А, а потом про скачки Т. Как Вас следует правильно понимать? Есть выборка с ацп , достаточно большая, скажем 100+ , период дискретизации и знание о виде сигнала. А,В и Т неизвестны, внутри выборки фиксированы, могут поменяться в другой выборке. При одиночных вычислениях T по 3 точкам, скачет, как правильно подсказали из за шумов. Найти нужно A, но не найдя T у меня не выходило найти А. Отношение сигнал-шум я думаю достаточно, т.к. если подождать 13T , A мерится с удовлетворительной точностью. Но ждать 13T низя, можно 2T.
Сообщение отредактировал IVANS - Jul 11 2014, 13:15
|
|
|
|
|
Jul 11 2014, 13:21
|
Гуру
     
Группа: Модераторы
Сообщений: 8 752
Регистрация: 6-01-06
Пользователь №: 12 883

|
Цитата(IVANS @ Jul 11 2014, 17:08)  Спасибо за подсказки по МНК, похоже , что это то, что нужно, но пока разбираюсь как это дело применить. Может есть рекоммендации?
Отношение сигнал-шум я думаю достаточно, т.к. если подождать 13T , A мерится с удовлетворительной точностью. Но ждать 13T низя, можно 2T. А где Вы собираетесь все это вычислять? В контроллере или машине, которая собирает эти данные? А чему равно 2Т? В секундах и точках.
|
|
|
|
|
Jul 11 2014, 13:41
|
Участник

Группа: Участник
Сообщений: 21
Регистрация: 2-12-04
Пользователь №: 1 296

|
Цитата(Tanya @ Jul 11 2014, 17:21)  А где Вы собираетесь все это вычислять? В контроллере или машине, которая собирает эти данные? А чему равно 2Т? В секундах и точках. Вычислять в контроллере, T [0;10c], самое плохое T 10c, для малых T можно и подождать хвост. 2Tmax = 20c. В точках 300 , можно еще увеличить.
|
|
|
|
|
Jul 11 2014, 14:13
|
Гуру
     
Группа: Модераторы
Сообщений: 8 752
Регистрация: 6-01-06
Пользователь №: 12 883

|
Цитата(IVANS @ Jul 11 2014, 17:41)  Вычислять в контроллере, T [0;10c], самое плохое T 10c, для малых T можно и подождать хвост. 2Tmax = 20c. В точках 300 , можно еще увеличить. Не знаю, потянет ли Ваш контроллер такое... Это стандартная задача из области нелинейного программирования. Например, подгонка под такую функцию - стандартный VI в LabView. Самому писать алгоритм нет смысла. Их много описанных. Влезет ли он в Ваш контроллер по размеру и скорости - вот основной вопрос. Хотя... могу подсказать, как уменьшить размерность и увеличить скорость, но, к сожалению, и объем. Решаем нелинейную задачу (по Т), а минимум по А и В находим на каждом шаге линейным способом. Формулу лучше переписать немного. И экспоненту по ряду Тейлора по опорным точкам из таблицы. Что за процесс у Вас?
|
|
|
|
|
Jul 11 2014, 14:37
|
Участник

Группа: Участник
Сообщений: 21
Регистрация: 2-12-04
Пользователь №: 1 296

|
Цитата(Tanya @ Jul 11 2014, 18:13)  ... Что за процесс у Вас? Спасибо , буду думать. Измеряем малый ток при большой емкостной составляющей.
|
|
|
|
|
  |
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|