Определение напряжённости поля антенны считывателя БСК

[podrivnik]

Всегда было интересно узнать напржённость поля своих считывателей до изготовления опытных образцов.
Воспользовался AN3394 от STMicroelectronics и Scilab.
Теперь буду определять параметры антенны до изготовления.
И могу теперь заранее определить необходимый ток через антенну.
А исходя из тока могу понять на какое сопротивлениее антенну согласовывать и нужен ли будет усилитель.

//Антенна БСК
d=0:0.001:0.1;//Расстояние от антенны, в метрах
a=0.05;//Одна сторона прямоугольника антенны, в метрах
b=0.05;//Вторая сторона прямоугольника антенны, в метрах
N=4;//Число витков антенны
i=0.0625;//Ток антенны. Зависит от сопростивления согласования антенны и выхода усилителя, в Амперах

//Дальше вычисление по AN3394 от STMicroelectronics
//"Antenna design and impedance matching guidelines for CR95HF multiprotocol contactless transceiver IC"
Hnumerator1=2*N*i*a*b;
Hdenominator1=%pi*sqrt(a^2 + b^2 + 4*d^2);
Hones=ones(1, size(d,2));
Hdenominator2=a^2 + 4*d^2;
Hdenominator3=b^2 + 4*d^2;
Hnumerator2=Hones./Hdenominator2 + Hones./Hdenominator3;
H=(Hnumerator1./Hdenominator1).*Hnumerator2;//Напряжённость поля А/м
clf();//Сотрём старый график
xtitle("График напряжённости поля от расстояния");
plot(d,H);//Строим график
xgrid();