Приложения

Расчет системы кратных периодов

rogram multidet integer real

real, allocatable real, allocatable integer, allocatable

n=19

i,j,kk,key,n

gg,ll,lhk,tt,r

h(:),lain(:),g(:),glh(:)

hs(:),ls(:),gs(:),left(:),right(:)

k(:),y(:)

allocate (h(n), g(n), gs(0:n), lain(n) , glh(n)) allocate (hs(n), ls(n), k(n), y(n), left(n), right(n))

gs(0)=1.33

open (3, file=multidet.datО do i=l,n

read (3,43el4.5)0 g(i) ,h(i) , lain(i) end do

do i=l,n;, glh(i)=g(i)/(lam(i)*h(i)); end do

call chanlin(n,glh,y) !! Массив glh упорядочен no возрастанию

!! В массиве у --- соответствие номеров

do i=l,n

j=y(i); hs(i)=h(j); ls(i)=lam(j); gs(i)=g(j);

end do

!! Независимая оптимизация 11=0

do i=i,n

ll=ll+sqrt(2*(gs(0)+gs(i))*ls(i)*hs(i)) end do

print *, Независимая оптимизация, 11 = ,11

print *, Общий период gg=gs(0); lhk=0 do i=l,n

gg=gg+gs(i); lhk=lhk+ls(i)*hs(i) end do gg=2*gg

tt=sqrt(gg/lhk); ll=sqrt(gg*lhk)

write (*,4a,el2.5)) Минимальные затраты = , 11 write (*,(a,el2.5)) Базисный период = , tt

!! Начальное приближение tt=sqrt(2*(gs(0)+gs(l))/(ls(l)*hs(l)))

!! Начальная группировка к(1)=1; кк=1

r=2*gs(l)/(kk*(kk+l)*tt**2) do i=2,n

do while (Is(i)*hs(i)<=r)

kk=kk+l; r=2*gs(i)/(kk*(kk+l)*tt**2)

end do

k(i)=kk end do

gg=gs(0); lhk=0 do i=l,n

gg=gg+gs(i)/k(i); lhk=lhk+ls(i)*hs(i)*k(i) end do gg=2*gg

tt=sqrt(gg/lhk); ll=sqrt(gg*lhk)