Размещено на http://www.allbest.ru/

Министерство образования и науки Российской Федерации

Федеральное государственное автономное образовательное учреждение высшего образования

«Новосибирский национальный исследовательский государственный университет»

(Новосибирский государственный университет)

Структурное подразделение Новосибирского государственного университета

Высший колледж информатики Университета

КАФЕДРА ЕСТЕСТВЕННОНАУЧНЫХ ДИСЦИПЛИН

Расчет схемы амплитудно-фазовой частотной характеристики ЭХС

Введение

Цель: убедится в перспективности способа измерения импеданса ЭХС с предварительной компенсацией сопротивления электролита и емкости двойного электрического слоя. Используемые сокращения для обозначения параметров:

Re - сопротивление электролита;

Rp - сопротивление электродной границы

Cd- емкость двойного электрического слоя и константы стадии диффузии, константы А Варбурга.;

Zdl - импеданса Варбурга

A=1000 - Константа Варбурга, Ом/сек^-1/2;

Методы расчета

Zdl = A./sqrt(w*(1-j)

R2=100; Значение объемного сопротивления электролита.

Cd=10^-5; Значение параметра Cd - емкость двойню эл. слоя.

R3=1000; Значениесопротивления перехода R3=Rp.

y3=j.*w*Cd+1. /(R3+Zdl)

Адмиттанс электродной границы без Ro.

Zio=R2+1. /y3

Импеданс исследуемого ЭХО.R2 - сопротивление электролита

ImZiomin=min(imag(Zio)) поиск экстремального значения

R3e=abs(ImZiomin*2)

dR3=abs(R3-R3e); абсолютная погрешность определения сопротивления Re=R3e.

oR3=(dR3/R3) *100; относительная погрешность определения Re=R3e.

Используя уравнение w*Cd*R3=1, зная частоту w в точке минимума и зная R3,можно определить значение Cd.

Cde=1/(75*R3e)

dCd=abs(Cd-Cde)

oCd=(dCd/Cd) *100

w=10^4

Rk=min(real(Zio)); Rk=100

R2e=abs(Rk);

dR2=abs(R2-R2e);

oR2=(dR2/R2) *100;

y2=1. /Zio;

C0=10^-10;

R0=10^6;

z1=1. /(j.*w*C0+(1/R0))

входной импеданс измерительной цепи. Учтемшунтирующее влияние измерителя на исследуемый импеданс Zio ЭХО.

1. Непосредственно по частотной характеристике можно определить параметры схемы замещения ИО с импедансом Zio=R2+1. /y3. В этом случае сопротивление электролита определяется с погрешностью, порядка 0,02.

2.Полученная оценка значения сопротивления электрохимической реакции R3=Rp имеют погрешностью порядка 12,5, а погрешность в оценке значения емкости двойного электрического слоя Cd - составляет 18. Это достаточно грубые оценки. Поэтому нужно усовершенствовать импедансный метод.

Применение компенсации для повышения точности измерений параметров Rp и Cd

Методика компенсации емкости cd

1. Введем отрицательное значение проводимости - j*w*Ck и определим:

1. Используем на первом шаге численного эксперимента компенсацию сопротивления электролита, тогда Zio1=Zio -Rk. Определим Rk.

В натурном эксперименте это легко осуществить компенсацию Re с помощью мостовой измерительной цепи. В нашем "численном эксперименте" мы это произведем вычитанием из общего импеданса значения сопротивления электролита Re. Численно значение

Re=min(real(Zio)), т.е. Re=|-Rk|= min(real(Zio))

Построим годограф импеданса Zio1 с компенсацией сопротивления электролита

Zio1=Zio-Rk;

ReZio1=real(Zio1);

ImZio1=imag(Zio1);

Попытаемся определить параметры Cd, сопротивления Rp более точно.

Re=R2; сопротивление электролита Rp=R3;

Преобразуем импеданс Zio1 в комплексную проводимость

Yio1=1. /Zio1;

Попытаемcя скомпенсировать емкость двойного слоя. В качестве первогоприближения выбираем

Ck=9.95*10^-6

Yio2=Yio1-j*w*Ck;

Zio2=1. /Yio2;

MIZ2=min(imag(Zio2))

find(imag(Zio2) ==MIZ2) ind=156.

i=1: lw;

wi=w (i (156)) wi= 9.7452e+003.

Определение параметров электродной границы и их погрешностей.

dCk=abs(Cd-Ck), dCk=1.0000e-007

абсолютная погрешность определения Ск.

oCk=100*dCk/10^-5, oCk=1.0000

относительная погрешность определения Ск.

dCk = 5.0000e-008, oCk = 0.5000.

Rp=2*abs(MIZ2)

Drp=abs(R3-Rp)

абсолютная погрешность определения Rp.

Orp=Drp/10- относительная погрешность определения Rp.

Rp= 1.0082e+003, Drp = 8.1738, Orp = 0.8174.

Определим значение константы Варбурга. Для этой цели воспользуемся значениями мнимой компоненты imag(Zio2(1) и imag(Zio2(2) на первой (1) и второй (2) низких частотах.

Im1=abs(imag(Zio2(1)));

Im2=abs(imag(Zio2(2)));

w1=2*pi*f (1);

a1=abs(1/sqrt(w1));

w2=2*pi*(f (2));

a2=abs(1/sqrt(w2));

A=(Im1+Im2)/abs(a1+a2)

dA= A-1000

odA=dA/10

A = 1.0093e+003

odA= 0.9334.

1.Используя импеданс-частотный метод, т.е. метод определения зависимости импеданса ЭХО от частоты, можно с достаточной высокой точностью определитьсопротивление электролита и его скомпенсировать.

2. Подключив дополнительно операцию по компенсации емкости двойного электрического слоя, повысить точность измерения компонентов Cd и Rp, не менее, чем на порядок.

Таким образом, методом математического моделирования на данном примерепоказано, что предлагаемое к реализации устройство с предварительной компенсацией параметров: сопротивления электролита Re и емкости двойного электрического слоя Cd потенциально может обеспечить относительные погрешности (o)измерения:оRe=0,01, oCd=0.50;

Rp= 1.0082e+003, оRp=0.8174.

Относительная погрешность определения константы Варбурга odA= 0.9334

Следовательно, используя сочетание методов:

1. Определения частотной зависимости импеданса исследуемого объекта и предварительной оценки значений искомых параметров Re, Rp и Cd;

2. Компенсации сопротивления электролита Re и емкости двойного электрического слоя Cd;

3. можно практически более чем на порядок повысить точность определения искомых компонентов Rp, Cd и константы Варбурга, а полученные результаты подтверждают перспективность предложенного принципа предварительной компенсации сопротивления электролита Re и емкости двойного электрического слоя Cd для построения "цифрового измерителя параметров двух- и трехполюсных электрических объектов".

Графики

Годограф импеданса Варбурга

X=real(Zdl)

Y=imag(Zdl)

Годограф адмиттанса электродной границы ЭХО без Ro

X=real(y3)

Y=imag(y3)

Годограф импеданса ЭХО с замедленными стадиями диффузии и перехода

X=real(Zio)

Y=imag(Zio)

Годограф мнимого компонента импеданса исследуемого ЭХО

X=w

Y=imag(Zio)

Годограф входного импеданса z1 измерительной цепи

X=real(z1)

Y=imag(z1)

Годограф адмиттанса исследуемого ЭХО

X=real(y2),

Y=imag(y2)

Годограф Zio1 с компенсацией сопротивления электролита

X=ReZio1

Y=ImZio1

Годограф проводимости электродной границы с компенсацией Re

X=real(Yio1)

Y=imag(Yio1)

Годограф Zio2 ЭХО с предварительной компенсацией Re и Cd

X=real(Zio2)

Y=imag(Zio2)

Зависимость imag(Zio2) от частоты при компенсации Re и Cd

X=i

Y=imag(Zio2)

Размещено на Allbest.ru