Кривые роста с насыщением: аппроксимация Паде

 Где начало того конца, которым оканчивается начало?

Афоризм №78 из "Плодов раздумья" Козьмы Пруткова

Интерлюдия

 Декан - бесцветная жидкость с характерным запахом.

       Из международной карты химической безопасности (ICSC 0428)

С теми, кто избежал пропуска пары, было обсуждение разных вопросов математического анализа, полезных для расчетов физических свойств веществ и происходящих с ними процессов.

Определение параметров светодиода, часть 3: нелинейная регрессия

И у черта, и у Бога

На одном, видать, счету

Ты, российская дорога,

Семь загибов на версту!

Из песни из т/ф "Гардемарины, вперед!"

В том случае, когда экспериментальные данные надо приблизить функцией, отличной от полинома, или же линейной комбинации функций, используется nlinfit, реализующая нелинейную регрессию, которая ищет параметры, дающие наименьшее среднеквадратичное отклонение от заданных точек. 

Основной принцип ее применения:

задать функцию в виде y=fun(параметры,переменная), например, как

    y=@fun(param,x) выражение для функции 

задать начальное приближение вектора переменных

    param0=[param0(1),param0(2), и так далее];

посчитать оптимизированные параметры

    param=nlinfit(вектор экспериментальных х, вектор экспериментальных y,fun,param0)

Примечание: при расчетах на GNU Octave для использования nlinfit надо сначала подгрузить пакет с функциями оптимизации (pkg load optim).

Пример соответствующей процедуры для фиттинга ВАХ светодиода, выраженной через функцию Ламберта:     

Определение параметров светодиода, часть 2: прояснение специфики реального диода

...а сову эту мы разъясним...

Цитата из повести М.А. Булгакова "Собачье сердце"

Для многих реальных полупроводниковых диодов их вольтамперная характеристика в полулогарифмических координатах не ложится на прямую (как следовало бы из формулы Шокли), так как они обладают еще и реактивным внутренним сопротивлением. 

Разъяснение соответствующей физики и уточненной процедуры фиттинга данных:

Графики, иллюстрирующие описанную выше в рукописном файле двухшаговую процедуру фиттинга и получившуюся аналитическую аппроксимамацию экспериментальных данных W-функцией Ламберта в обычных и полулогарифмических коодинатах (в MATLAB и GNU Octave вычисляется при помощи функции lambertw):

Примечание: при использовании GNU Octave, чтобы функция lambertw стала доступной, необходимо перед запуском программы загрузить пакет symbolic (pkg load symbolic); если он еще не установлен - то установить (pkg install symbolic).

Дополнительно почитать про свойства W-функции Ламберта:

• Corless R.M., Gonnet G.H., Hare D.E., Jeffrey D.J.,  Knuth, D. . (1996). On the Lambert W function. Advances in Computational Mathematics, 5 (1996) 329-359 [ссылка].
• Stewart S. A new elementary function for our curricula?. Australian Senior Mathematics Journal, 19 (2005) 8-26 [ссылка].
• Kesisoglou I., Singh G., Nikolaou M. The Lambert function should be in the engineering mathematical toolbox. Computers & chemical engineering, 148 (2021) 107259 [ссылка]

Специально в приложении к фотовольтаике:

• Nassar-Eddine I., Obbadi A., Errami Y., Agunaou M.  Parameter estimation of photovoltaic modules using iterative method and the Lambert W function: A comparative study. Energy Conversion and Management, 119 (2016) 37-48 [ссылка].

Практическое занятие: определение параметров светодиода, часть 1: выбор и предварительный анализ данных

 светить — 

            и никаких гвоздей!

Из стихотворения В.В. Маяковского

 Задание: 

• выбрать они какой-нибудь светодиод (выбранные должны быть разными у всех студентов в обоих подгруппах) на сайте OSRAM: https://ams-osram.com/products/leds; 
• на странице описания выбранного светодиода скачать файл "Datasheet", там найти график вольт-амперной характеристики прямого тока (IF от VF), оцифровать ее при помощи WebPlotDigitizer и построить соответствующий график (внешний вид графика должен совпадать с тем, который приведен в Datasheet; линии сетки добавляются командой grid on);
• Ориентируясь на уравнение Шокли для диода, построить регрессию оцифрованных данных (подсказка: для достаточно больших напряжений можно пренебречь вычитаемой единицей ).
• Проверить путем построения графика тока в полулогарифмических координатах, ложатся ли данные на прямую, если да - то они соответствуют формуле Шокли и можно непосредственно  найти величины тока насыщения и величины, обратной произведению коэффициента идеальности на тепловое напряжение, найденные значения параметров привести, сопроводив их величиной 95% доверительного интервала.
• Если данные не ложатся на прямую в полулогарифмических координатах - то диод далек от идеального, как работать в этом случае будет рассмотрено на следующем занятии.  

При оформлении графиков - не забывать подписать в заголовке, какой именно диод исследуется.