Курсовая работа: Разработка програмного обеспечения для расчёта дисперсионной характеристики планарного волновода

aa:= sqrt(1-sqr(sin(tetta)));

aa:= 1/(2*pi*n2*aa);

bb:= sqrt(sqr(n2*sin(tetta))-sqr(n3));

cc:= sqrt(sqr(n2*sin(tetta))-sqr(n1));

dd:= n2*sqrt(1-sqr(sin(tetta)));

dwl:=aa*(arctan(bb/dd)+arctan(cc/dd)+pi*m);

end;

function q(wl,tetta,n2,n1:real):real;

//функция, вычисляющая волновое число q

begin

q:=(2*pi*sqrt(sqr(n2*sin(tetta))-n1*n1))/wl;

end;

function h(wl,tetta,n2:real):real;

//функция, вычисляющая волновое число h

begin

h:=(2*pi*cos(tetta)*n2)/wl;

end;

function p(wl,tetta,n2,n3:real):real;

//функция, вычисляющая волновое число p

begin

p:=(2*pi*sqrt(sqr(n2*sin(tetta))-n3*n3))/wl;

end;

function e1(wl1,d1,n11,n21,n31,tetta1,x1:real):real;

//функция, задающая профиль Еу ТЕ моды при условии х<=0

begin

e1:=exp(q(wl1,tetta1,n21,n11)*x1);

end;

function e2(wl2,d2,n12,n22,n32,tetta2,x2:real):real;

//функция, задающая профиль Еу ТЕ моды при условии 0<х<d