Форум РадиоКот https://radiokot.ru/forum/ |
|
LCR на STM32G431 https://radiokot.ru/forum/viewtopic.php?f=10&t=197692 |
Страница 1 из 3 |
Автор: | linkov1959 [ Пн апр 21, 2025 20:39:01 ] |
Заголовок сообщения: | LCR на STM32G431 |
Плата STM32G431CBU6 (WeAct Studio). Среда разработки CubeIDE. Первый этап - измерение LCR на 10кГц. На выходе двух каналов DAC1 синус в противофазе с буфером на ОУ mc33202 и фильтром на входе 300 Ом 1000рФ. Таблица синуса на 720 точек. На первом этапе нужно выйти на максимальную точность на одном диапазоне. Синус очень хороший. Пробовал много процессоров. Лучший синус на GD32F303, но с CubeIDE он не дружит. С выходов синуса через резисторы 100-200 Ом и щупы Кельвина ток идет на измеряемый LCR, а напряжение c потенциальных контактов щупа на входа ADC через буфера ОУmc33202 . За период измерения заполняем данными ADC массивы и вычисляем Z и RLC. Я еще не сталкивался с такой математикой, надеюсь на помощь. Пытался спрашивать у ИИ, но только потерял время и нервы. ![]() |
Автор: | linkov1959 [ Ср апр 23, 2025 19:55:16 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Получил 3 массива АЦП за период синуса. 1- массив напряжения после буфера 1 канала ЦАП до шунта 150 Ом, 2- массив напряжения после резистора шунта 150 Ом 1 канала ЦАП. 3- массив напряжения 2 канала после резистора 150 Ом. Из 1 и 2 массива получаем массив тока через LCR. Из массивов 2 и 3 по разнице получаем массив напряжения на RLC. Вот что предлагает ИИ Deep Seek. 3. Алгоритм расчёта импеданса Теперь у нас есть: Массив напряжения на нагрузке: V_load[k] = V1[k] - V2[k] Массив тока: I[k] = (2V3[k] - V1[k] - V2[k]) / 300 Шаги расчёта (аналогично предыдущему методу): Применить DFT к V_load и I для получения Re_V, Im_V, Re_I, Im_I. Вычислить импеданс: Z=(ReV+j?ImV)(ReI+j?ImI) Z=(ReI?+j?ImI?)(ReV?+j?ImV?)? Разделить на активную (R) и реактивную (X) части: R=ReVReI+ImVImIReI2+ImI2 R=ReI2?+ImI2?ReV?ReI?+ImV?ImI?? X=ImVReI?ReVImIReI2+ImI2 X=ReI2?+ImI2?ImV?ReI??ReV?ImI?? Определить L или C (как в предыдущем ответе). 4. Плюсы такого подхода ? Уменьшение влияния просадки буферов (используется разностный сигнал). ? Симметричная обработка обоих плеч моста. ? Требуется только один DFT для напряжения и один для тока (вместо трёх). 5. Пример кода (STM32G431) c // Напряжение на нагрузке: V1 - V2 for (int k = 0; k < 720; k++) { V_load[k] = V1[k] - V2[k]; I[k] = (2 * V3[k] - V1[k] - V2[k]) / 300.0f; // Ток через 150 Ом } // DFT для V_load и I dft(V_load, &Re_V, &Im_V); dft(I, &Re_I, &Im_I); // Расчёт Z float denom = Re_I * Re_I + Im_I * Im_I; float R = (Re_V * Re_I + Im_V * Im_I) / denom; float X = (Im_V * Re_I - Re_V * Im_I) / denom; // Определение L или C if (X > 0) { L = X / (2 * M_PI * f0); } else { C = -1 / (2 * M_PI * f0 * X); } 6. Важные замечания Калибровка разности АЦП: Убедитесь, что каналы АЦП синхронизированы (или учитывайте задержку). Если есть смещение нуля, вычитайте его до расчётов. Повышение точности: Используйте усреднение нескольких периодов. Примените цифровую фильтрацию (например, медианный фильтр). Оптимизация для STM32: Используйте DMA для сбора данных с АЦП. CORDIC и FMAC ускорят расчёты DFT. |
Автор: | bob1 [ Чт апр 24, 2025 06:36:27 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
не сталкивался с такой математикой, надеюсь на помощь https://radiokot.ru/circuit/digital/measure/108/ |
Автор: | linkov1959 [ Чт апр 24, 2025 08:53:12 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
bob1, Я не смог его программу перенести в CubeIDE. Там без помощи автора не просто разобраться. Пытаюсь разобраться сам поэтапно, упростив задачу до примитива. Вот пробую включить в STM32G4 функции DSP. Вот уже включил, значит расчеты пойдут быстрее. |
Автор: | bob1 [ Чт апр 24, 2025 10:14:36 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Вы попробуйте во float считать. Расчетов не так много. Это же не весь спектр считать при помощи FFT, а всего 1 точка (частота)!! |
Автор: | linkov1959 [ Чт апр 24, 2025 22:54:16 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
bob1, Сначала нужно данные для расчета подготовить. Это и определяет точность. Добавлено after 4 hours 34 minutes 44 seconds: Стартую синхронно синус ЦАП и АЦП циркулярно от одного таймера. Первый период смазан вначале, брать в расчет можно данные только с второго цикла. После второго цикла останавливаю таймер и можно обрабатывать данные из 3-х буферов АЦП. |
Автор: | Andrey_B [ Пт апр 25, 2025 10:35:47 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
linkov1959, стартовать ЦАП и АЦП можно асинхронно, начальная фаза не важна. Важна синхронная работа. Длина буферов должна быть такой, чтоб в них укладывалось строго целое количество периодов измерительной частоты. Продумывайте алгоритм так, чтобы генерация синуса ЦАП и весь механизм оцифровки не прерывались, без этого высокой точности не достичь. Если не сможете обеспечить обработку всех циклов измерения, можно брать только нечётные, а вместо пропущенных чётных производить вычисления. |
Автор: | bob1 [ Пт апр 25, 2025 11:35:43 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
советует ИИ DeepSeek Куда мир катится....роботы будут думать за вас и вы им верить.// Нормировка (делим на количество точек) Это лишняя операция. Потом все равно числа будите делить (напряжение на ток).Почему не FFT? FFT избыточен для одной частоты. Требует больше ресурсов ЦПУ. Квадратурный метод даёт ту же точность при меньших вычислениях. В примере у вас вычисление FFT в одной точке (для одной частоты). Другие частоты считать не нужно ведь. Они не нужны. Цифровые фильтры тоже не нужны.
|
Автор: | linkov1959 [ Пт апр 25, 2025 17:39:21 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
bob1, ИИ не думает, он сортирует для вас информацию от думающих людей. |
Автор: | linkov1959 [ Сб апр 26, 2025 19:40:55 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Вот уже выходим на ключевые точки. Главная точка - это запуск и остановка таймера ТИМ3, который запускает ADC DMA(3 канала) и DAC DMA(2канала) в циркулярном режиме. По прерыванию ADC DMA два раза за цикл разбираем буфер(массив) ADC DMA(3 канала) 1080 точек по каналам 360 точек на канал и суммируем в 3 отдельных массива 16 раз с накоплением и останавливаем таймер ТИМ3. Накопленные данные делим на 16 и получаем 3 усредненных массивов. 1 и 2 пересчитываем по разнице значений и получаем массив напряжения на RLC. Из массивов 2 и 3 по разнице получаем массив напряжений на опорном резисторе 150 Ом и пересчитываем на массив тока через RLC. Рассчитываем в Экселе и заполняем опорные массивы синуса и косинуса на 360 точек. Вроде уже есть всё для расчета Z, если пренебречь пока калибровкой и задержками. Добавлено after 4 hours 22 minutes 15 seconds: Массивы напряжения и тока получаем в прерывании ADC DMA и поднимаем флаг готовности данных. В основном цикле по этому флагу рассчитываем комплексное Z и снова запускаем ТИМ3 для нового цикла. Усредняем Z за 10 циклов и рассчитываем R,C,L,Q и др. Обнуляем счетчики циклов и запускаем карусель снова. |
Автор: | КРАМ [ Сб апр 26, 2025 21:47:24 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
В примере у вас вычисление FFT в одной точке (для одной частоты) FFT для одной точки не существует. Либо всё, либо ничего. Для одной точки вычисляют DFT. Дискретное преобразование Фурье - это узкополосная фильтрация. И она дает на выходе сразу комплЕксный результат. Обычная квадратурная выборка имеет никакую помехоустойчивость. Те, кто идет таким путем, дальше начинают медленно и печально накапливать скалярные значения вместе с перенесенными в ноль помехами и шумом: Усредняем Z за 10 циклов и рассчитываем R,C,L,Q и др. Обнуляем счетчики циклов и запускаем карусель снова. ![]() |
Автор: | bob1 [ Вс апр 27, 2025 07:53:44 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
FFT Перепутал перевод. Конечно ДПФ.linkov1959 писал(а): По прерыванию ADC DMA два раза за цикл разбираем буфер(массив) ADC DMA(3 канала) 1080 точек по каналам 360 точек ADC в камне 2 штуки. Или выборка 3 каналов идет не одновременно?linkov1959 писал(а): Накопленные данные делим на 16 и получаем 3 усредненных массивов Получили числа с фиксированной точкой 28бит целая часть и 4 бит дробь. Берете и выкидываете дробную часть. Уменьшаете точность результата. Не нужно делить на 16! linkov1959 писал(а): Рассчитываем в Экселе Камень сам не может при старте заполнить таблицы?У Вас какая задача стоит? Сделать замер на 1 частоте или вы хотите анализировать деталь на нескольких частотах сразу и строить графики? |
Автор: | Yuriy_King [ Вс апр 27, 2025 11:25:21 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
У Вас какая задача стоит? Сделать замер на 1 частоте или вы хотите анализировать деталь на нескольких частотах сразу и строить графики? ТС не озвучил свои намерения. Вот для примера измерения требующие более детальной проработки... Спойлер![]() ![]() |
Автор: | linkov1959 [ Вс апр 27, 2025 16:04:11 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Догадливые!:) 1 частота это разминка, 3 частоты это первый подход. А дальше будет интересно... |
Автор: | Andrey_B [ Вс апр 27, 2025 20:57:03 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Создаем 2 массива (на той же частоте 1 кГц): float ref_sin[n] = sin(2 * π * 1000 * n / (720 * fs)); // fs = частота дискретизации float ref_cos[n] = cos(2 * π * 1000 * n / (720 * fs)); Можно обойтись одним массивом, значения второго получить сдвигом индекса первого на 90 градусов. Калибровка нуля: Измерьте сигналы при коротком замыкании (вычтите offset). Калибровка шунта: Учёт точного значения 150 Ом и паразитного RL. Фазовые сдвиги ADC: Если каналы ADC имеют разную задержку, потребуется коррекция Калибровка по трём точкам (open-short-load) скомпенсирует все статические погрешности. linkov1959 писал(а): Догадливые!:) 1 частота это разминка, 3 частоты это первый подход. А дальше будет интересно... balmer это "интересно" уже реализовал. |
Автор: | linkov1959 [ Пн апр 28, 2025 10:12:31 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Andrey_B, Ещё интересней... Беру случайное колечко с платы. Е7-8 - 50mН 11 Ом. RLC-2 - 50mН 10 Ом neekeetos - 47mН 17 Ом RLC 71 - 52mH 25 Ом Частота везде 1кГц. Мост постоянного тока - 0.169 Ом. Ошибка в 100 раз. Единожды совравшему кто поверит?, а мы верим... Не можешь - пиши прочерк или текст. Все LCR врут. Пока это не обсуждаем. Обсуждать будет методы калибровки на одной частоте 1кГц. |
Автор: | КРАМ [ Пн апр 28, 2025 13:58:00 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Мост постоянного тока - 0.169 Ом. А что удивляет? ![]() |
Автор: | linkov1959 [ Пн апр 28, 2025 14:17:31 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
Пора переходить к 2-й точке отсчета. Напомню 1-ю. Запускаем ТИМ3 и получаем готовые данные для расчетов. Даже останавливать не надо. Таймер останавливается сам после обработки нужного количества циклов сбора данных и их усреднений. 2-я точка отсчета это обратная связь - возможность контролировать результат на каждом этапе. В данном случае использую возможности отладки в CubeIDE. Просто останавливаю после заполнения массивов и смотрю их. Первое значение АЦП должно соответствовать началу контрольной таблице синусов для метода квадратурного детектирования. У нас два противофазных канала АЦП и в начале цикла они равны. Подогнал дискретной фазой ЦАП почти равенство в первой точке, остальное выбирается корректировкой синусной таблицы. |
Автор: | bob1 [ Пн апр 28, 2025 15:01:28 ] |
Заголовок сообщения: | Re: LCR на STM32G431 |
ИИ не думает, он сортирует для вас информацию от думающих людей. Исключите меня из списков, чтоб он мне подобное не сортировал ![]() linkov1959 писал(а): Е7-8 - 50mН 11 Ом. RLC-2 - 50mН 10 Ом neekeetos - 47mН 17 Ом RLC 71 - 52mH 25 Ом Ничего не забыли написать? linkov1959 писал(а): Все LCR врут Если вы чего-то не понимаете и не можете объяснить, то при чем ту приборы!
|
Страница 1 из 3 | Часовой пояс: UTC + 3 часа |
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group http://www.phpbb.com/ |