f(alpha) !функция
/vab/omega,x1(8), PARAMETER :: I=(0.0,1.0) !мнимая единица(8)
:: f,alpha,x1
complex(8)::eps, eps1, eps2, eps3, ka2, del11, del21, k11, k12, k22, psi, gam1, sig1, sig2, elem !они комплексные
real(8)::mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1 , omega !1 набор
! real::mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1 , omega !2 набор
! real::mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1 , omega !3 набор
! real::mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1 , omega !4 набор
! real::mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=1, ro1=1 , omega !5 набор
! real::mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1,
ro2=1, ro1=1 , omega !6 набор
eps1=(mu2*(h**3))/(6*(1-nu2))
eps2=((h**3)/12)*ro2*(omega**2)
eps3=h*ro2*(omega**2)
ka2=cdsqrt(dcmplx((ro1*(omega**2))/mu1))
eps=(1-2*nu1)/(2-2*nu1)
psi=eps1*(alpha**4)-eps2*(alpha**2)-eps3
gam1=(alpha**2)-0.5*(ka2**2)
sig2=cdsqrt(dcmplx((alpha**2)-(ka2**2)))
sig1=cdsqrt(dcmplx((alpha**2)-eps*(ka2**2)))
k11=(sig2*(ka2**2)*(sig1*sig2*((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)-cdexp(-(sig2*h1)))/2)-(alpha**2)*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2)))/(-4)
k12=(sig1*sig2*(gam1+(alpha**2))*(((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2)-1)-(gam1*(alpha**2)+(sig1**2)*(sig2**2))*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)-cdexp(-(sig2*h1)))/2))/2
k22=(sig1*(ka2**2)*(sig2*sig1*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2)*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)-(alpha**2)*((cdexp(sig2*h1)-cdexp(-(sig2*h1)))/2)*((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)))/(-4)
del21=(sig1*sig2*((1/4)*(ka2**4)-((gam1+(alpha**2))**2)+((1/4)*(ka2**4)+((gam1+(alpha**2))**2))*((cdexp(sig1*h1)+cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2))-(alpha**2)*((gam1**2)+(sig1**2)*(sig2**2))*((cdexp(sig1*h1)-cdexp(-(sig1*h1)))/2)*((cdexp(sig2*h1)-cdexp(-(sig2*h1)))/2))/2
del11=((cdexp(sig2*h1)+cdexp(-(sig2*h1)))/2)
! f=real(((k22/del21)/((k22/del21)*psi+1))*cdexp(-I*alpha*x1)) !для перемещений
f=real(((1)/((k22/del21)*psi+1))*cdexp(-I*alpha*x1)) !для напряжений
end function f
ПРИЛОЖЕНИЕ B
Листинг программы в среде Maple15 (построение графиков нулей и полюсов).
(stats);
:= readdata(`c:\\MOmegaNuli.txt`, 1);:=
readdata(`c:\\MNuli.txt`, 1);:= readdata(`c:\\MOmegaPolusa.txt`, 1);:=
readdata(`c:\\MPolusa.txt`, 1);[scatterplot](x, y);[scatterplot](z, q);
ПРИЛОЖЕНИЕ Г
Листинг программы в среде Maple15 (построение графиков перемещений и напряжений).
(stats);
[listplot](readdata(`c:\\xf19n.txt`, float,
2), color = red, labels = ["x", "u"], labelfont = [TIMES,
BOLDITALIC, 11]);[listplot](readdata(`c:\\yf14.txt`, float, 2), color = red,
labels = ["x", "u"], labelfont = [TIMES, BOLDITALIC, 11]);