автореферат диссертации по приборостроению, метрологии и информационно-измерительным приборам и системам, 05.11.16, диссертация на тему:Алгоритмы и устройства корреляционного анализа неравномерно-дискретизированных сигналов. Приложения

Из фондов Российской государственной библиотеки

Белолипецкий, В.Н.

Алгоритмы и устройства корреляционного анализа неравномерно-дискретизированных сигналов. Приложения

Москва

Российская государственная библиотека diss.rsl.ru 2007

Белолипецкий, В.Н.

Алгоритмы и устройства корреляционного анализа неравномерно-дискретизированных сигналов. Приложения [Электронный ресурс] : Дис. ... канд. технические науки : 05.11.16. - М.: РГБ, 2007. - (Из фондов Российской государственной библиотеки)

Информационно-измерительные системы

Текст воспроизводится по экземпляру, находящемуся в фонде РГБ:

Белолипецкий, В.Н.

Алгоритмы и устройства корреляционного анализа неравномерно-дискретизированных сигналов. Приложения

О, [б. г.]

Российская государственная библиотека, 2007 (электронный текст)

iV,.: .: -v

IL

ют-ss

И e^j Vt и ^ ^ Й , ft / 7 Л

v^ <J Ъо чу* fa Jj '4- & ОУ с fa?

v n* . A.«.. -. . _ гч А гх , Л

bxJ & i01/ о 6 6Ш a J ^ь kep cii ил Ju&p^^.

(j/U '

iO(Uiso > ы i тем - йь^ы/^

os. /i /4

ПРИЛОЖЕНИЯ

У/

ПРИЛОЖЕНИЕ I

SUBROUTINE RAND(N»XiIY)

ГЕНЕРИРОВАНИЕ П"П С РАВНОМЕРНЫМ распределением В ДИАПАЗОНЕ 0,1

IV - ПАРАПЕТ? ГЕНЕРАТОРА N - 05'эЕМ ЗЫВОР^И X{N) - ВЫХОДНОЙ МАССИВ

DIMENSION XtN> On 10 l=t,M

I ^ = IY

Call RANOU (I x, I Y ♦ YELI 10 X { I ) Sy -U

R6TURM END

SUBROUTINE RANOt(R5 ,RF,N,Y)

ГЕНЕРИРОВАНИЕ ПОП С РАВНОМЕРНОМ РАСПРЕДЕЛЕНИЕМ В ДИАПАЗОНЕ *$iRP

RS - наименьшее значение ПСП

RF - наибольшее значение ПСП N - ОВэЕМ зызорки Y ! N' - ВУХЗДНОй МАССИВ

dimension Yin)

lY=22?70il25 DO 10 1 = 1 >>1 I w = I Y

call ramouUX»iy.yfu

RffYFL

10 Y(I)Я(4F-R5)*R + RS R f! T U R м

end

SUBROUTINE P!UEH(Ml*0tT«4tXfMiXK,G)

ГЕНЕРИРОВАНИЕ ПСП С КФ EXP(~A3S(TAY)/Тj

N1 - КОЛИЧЕСТВО ОРДИНАТ ИМПУЛЬСНОЙ ПЕРЕХОДНОЕ ХАРАКТЕРИСТИК/1 D - ИНТЕРВАЛ ДИСКР5ТИЗАЦИИ Т - ПАРАМЕТР Кф

X(N' - ВХОДНОЙ МАССИВ ПСЧ - "аЕЛЫЙ ШУМ»

N - Размер массива

ХК(М) - ЗЫ*0ДН0Й МАССИВ Г)ОЧ

MaN-Nl - РАЗМЕР ВЫХОДНОГО МАССИВА

G < N1i » МАССИВ ОРДИНАТ МПУЛЬСНОй ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ

DIMENSION K(N)» х К(N > * G(N 1 >

JsNUi

A=£XP{-О/Т)

Н* 1 .

do 2г i = 11 ^ 1 H = H * A 0 S ( I )sy

DO 10 '<=! Ss0 . DO 30

И S = S + X(J-I)*G ( !)

Хк(К)з(s+X(J))*(D/T)

10 J*J+l RgTURy

EnD

SUBROUTINE FILE 1M (N1 ,OELT'T ,N .X',M »XK ,C}

ГЕНЕРИРОВАНИЕ ПОП С КФ EXP(-А3${ТAY>/Т)*(1♦A3S(ТАY)/Т)

N1 - КОЛИЧЕСТВО ОРДИНАТ ИМПУЛЬСНОЙ ПЕРЕХОДНОЕ ХАРАКТЕРИСТИКИ DHLT - ИНТЕРВАЛ ДИСКРЕТИЗАЦИИ Т - ПАРАМЕТР Кф

X(N>-3K0flH0p| МАССИВ ПОЧ"ЬЕЛЫЯ ШУМ"

N"PA3M£P массива

ХК(м)-ЗЫХОЛНОЙ массив псч

М г N-N 1 " РАЗМЕР ВЫХОДНОГО МАСЗИВА

С!N1) - МАССИВ ОРДИНАТ ИмПУЛЬЗНОй ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ

DIMENSION X(N),X К(М) • С(N t'

JsNl+i

BsDEUT/T

AsEXp < -В J

H 1 = 1 .

DO 10 I si Hi H 1 = н i * 0 G{I)=Hl*FlPAT(I) DO 20 <=Ы Ss0.

DO 30 I = 1 , Mi 0 S=S*X(I)*G{T)

ХК(K)s5«B**2 0 JsJ+1

RETURN End

SUBROUTINE FILE2H(N1 ,OELT»TiN»X-,M»XKfC) ГЕНЕРи'ОВА^ИЕ ПОП С кф EXP(-A3s(TAY>/Т)*«1-A3$(JAY)/Т)

Nl - КЗЛИЧ5СТВО ОРДИНАТ ИМПУЛЬСНОЙ ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ OPL Т ; ИНТЕРВАЛ ДИСКРЕТИЗАЦИИ Т - ПАЭАНЕТР Кф

X!N' - ВХОДНОЙ МАССИВ псч - "БЕЛЫЙ ШУИ» n - Размер массива

Xh-fM) - ВЫГОДНОЙ МАССИВ ПСЧ MnN-Ni - PА3МНР ВУХОдНОгО МДССиВА

s(n 1) - массив ординат импульсной ПЕРЕХСДнОЙ характеристики

DIMENSION Х(М), X К ( М )» G { N I)

JsNl+i

A=DEUT/T

СгЕХр{- A )

8st.

00 10 Isj ,Vi

а?в + с

Ъ G(I1,-FUOAT{I!*A> 00 20 <=l,* S*0»

DO 30 ! = 1, N 1 Ъ S*S + Xf>I)*G(I)

X К ( К ) E ( s ♦ X I 0 H * A

3 JpJ*i R e T U R ч End

subroutine fiue3m(n 1,0,t,n > x,m,xk , g » генерирование псп с кф

exp(-s3s(tay)/t)*(l+a8s(tay!/r+tay**2/(2*t**2i)

n1 - количество 3 3 д и и а т импульсной переходное хар а.ктерист ик а

d - интервал дискрет иsauии т - параметр кф

x f n > - входной массив псч - "5елый шум» n - размер массива xк < м) - выгодной массив псч MbN-NI - размер выходного массива

g(n1) - массиз ординат ИМПУЛьЗНОй переходной характеристики

01 MENS ION KfMJ 9ХК(М)«C(N1>

JnNl+i

b = d/t

A»EXP{-В)

нз 1 ,

OS 20 1=1,MJ нвн* a

0 С(I)»H*FLOAT(I*«2) DO 10 Sa0 .

DO 30 1=1,Mi 0 SsS+X{J-I 1*G <I)

xa(k)=$*b**3 lo JKO+1

RETURN EMD

SUBROUTINE FICl4(4bOELT,f,OMEGiNiXiM,XK»C) ГЕНЕРИРОВАНИЕ ПСП С КФ EXP(-A3s<ТAY5/т)*COsС 3MEG*TАу)

N! ~ КЗПИЧССТВО ОРДИНАТ ИМПУЛЬСНОЙ ПЕРЕХОДНОЙ ХАРАКТЕРИСТИК А 05:LТ я иНТ;РОАЛ ДИСКРЕТИЗАЦИИ Т - ПАРАМЕТР Кф

ОMEG - ПАРАМЕТР КФ

X f N > « ВХОДНОЙ МАССИВ ПСЧ - "5ЕЛЫЙ ШУМ"

N " РАЗМЕР МАССИВА

ХК(М) - ВЫГОДНОЙ МАССИВ печ

HsN-Nl - РАЗМЕР ВЫХОДНОГО МАССИВА

G(nI) - маЗсиз ординат импульзнои переходное характеристик^

DIMENSION х(N)» X К(М) » G{N I)

J«NJ * !

B-DELt/T

A«EXp{-В )

В 1=T*oiEG

B2=OME3*0EwT

T0=(SsU(l»+Bl**2)-l.)/Bl

S? = SIm(B2 )

CCsCOi(B2 )

H 1= 1 .

C 1 = 1.

S1=0,

do 10 1 = 1,ml H i = H 1 a A

Cscuca-s H*s0

G( I )4Hi*(C + T0*S)

с i=c ъ si=s

do 20 <=1,1 SUM=0,

do 30 1=1,ml ? SUM=SU1+X!J-I'*G(I) XK(K)i(SUM + XIJ) )* 3 Ъ j=j+1

RETURm

end

SUBROUTINE FlCsiMlNt,OELT «Т,01EC»N|X|M|XK,C> ГЕНЕРИРОВАНИЕ ПСП С КФ

EXP'-A3S«T*Y)/T>*'COS(0MEG*TAY)+1/(Т*ОМЕС)*SIN(OMEG«ABS(ТAY)i)

N1 - КОЛИЧЕСТВО ОРДИНАТ ИМПУЛЬСНОЙ ПЕРЕХОДНОЕ- ХАРАКТЕРИСТИКИ DELT . ИНТЕРВАЛ ДИСКРЕТИЗАЦИИ Т - ПАРАМЕТР Кф OHEG - ПАРАМЕТР КФ

XtN> - ВХОДНОЙ МАССИВ ПСЧ - "БЕЛЫЙ ШУМ" N - размер массиза Xк СМ) - ВЫХОДНОЕ МАССИВ ПСЧ MsN-Nl - РАЗМЕР ВЫХОДНОГО МАССИВА

G(Nl) - МАССИВ ОРДИНАТ ИМПУЛЬСНОЙ ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ

01 HEMS ION *tN>,XK(M),G(Nl>

JsNl* 1

BSOELT'T

А*ЕХР<-В>

в1-0мё»*0е^т Вг = Т*£)1ЕС

t0 = sa^ t(1.*b2**21*b/32 sa=sjM(B1) c3=ca*(в i) H i = 1 , С ! = 1 »

sj = 0.

OQ 10 1 = 1,«I Hl=HUft C=C1*03-S l*S0 s=s i*c?+c i*sa g ( 1 > s h 1 * s С J =C 5 S1 = S

Do 20 <=i>M

SUHsjl,

Do 30 1 = 1, Ml ъ SUH=SU1+X(J-I'»C{I)

X К ( к ) S 5 и м * т 0 5 JsJ+1 RETURN END

subroutine FIcs2m(NI ,DELT»T,OIEGINIXIM,-*K,G>

ГЕНЕРИРОВАНИЕ ПСП С КО

EXP<-A3S(TAY)/T)*(COS(OMEG*TAY)-1/(T*OMEG}*SIN(QMEG*ABS(TAY))>

N1 - КОЛИЧЕСТВО ОРДИНАТ ИМПУЛЬСНОЙ ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ DELТ . ИНТЕРВАЛ ДИСКРЕТИЗАЦИИ Т - ПА'АМЕТР Кф OMEG ; ПАРАМЕТР КФ

X(N> - 8Х0ДН0Й МАССИВ ПСЧ - "БЕЛЫЙ ШУМ»

N - размер массива

ХК(М) -> ВЫХОДНОЙ МАССИВ псч

HsN-Nl - РАЗМЕР ВЫХОДНОГО МАССИВА

С f N1) - МАССИВ ОРДИНАТ ИМПУЛЬСНОЙ П^РЕХСДНОЙ ХАРАКТЕРИСТИКИ

DIMENSION X(N>.х К(М ) •G(N 1 )

J я N 1 ♦ 1 в1=DELT/T

АсЕХр ( -В1 ) BrOMEctOElT

тг=т»0чЕс s ® = s i ы ( в) C0=COS(В) Н 1 = 1, С1=1,

s i=0.

On 10 1 = 1 Hi

HlsHUA

сяснсг-s i*s0

SPSl»(i2 + ol*S0 G ( I >sHlMC-S/T0l С 1 =c 0 st=s

Do 20 <=1,fl SUM=0.

00 30 I = 1 , Ml

SUM = SUi + X(V-I)*G ( I)

хк {k ) j- ( sum + x { 0) ) #9 1

JsJM

RfTURM

e^d

SUBROUTINE XGEN(X,N,C1,C2>

ГЕНЕРИРОВАНИЕ ПОП С Кф ЕХР<-AIPHA*ABs!ТAY>)

X(N> - ВХОДНОЙ МАССИВ - "6ЕЛЫ* ШУМ" »ВЫХОДНОЕ МАССИВ N - РАЗМЕР МАССИВА С1 - ПАРАМЕТР АЛГОРИТМА С2 - ПАРАМЕТР АЛГОРИТМА

DIMENSION X{N) DO 10 1=2 Н 10 х (i )=ai*x < I - Р + С2+Х (i) EMD

SUBROUTINE XGf:N2(N,X,Y,A0»Al,3i;B2j ГЕНЕРИРОВАНИЕ ПСП С КФ

ЕХР(-А,РНА»АВЛТАГ))« ( 1 + ALPHA *ABS < Т A Y ) ) » EXPUA„PHA*ABS(TAy) )*COS(ОМЕС*ТАу> t

EXP <-A.PHA*A3S(TAY))«(COs(OHES*TAY)+ALPbA/0M£G*SIN(0MEG*A8S(T*Y)i) EXP(-A.PHA*A3S(TAY) )* (COS <OMES*TAY)-ALPbA/0MEG*SUM0MEC*A3S(TAY)j >

X{N > - ВХОДНОЙ МАССИВ - "6ЕЛЫЯ ШУМ"

N - РАЗМЕР МАССИВА

Y{N' - ВЫХОДНОЙ МАССИВ

А0 - ПАРАМЕТР АЛГОРИТМА

А} - ПАРАМЕТР АЛГОРИТМА

В1 - ПАРАМЕТР АЛГОРИТМА

В2 ~ ПАРАМЕТР АЛГОРИТМА

ЗНАЧЕНИЯ ПАРАМЕТРОВ ЗАВИСЯТ ОТ ВИДА МОДЕЛИ КФ

DIMENSION *(N!,Y(N) Y(1)в* ( 1 j Y f 2 > ev(2 )

00 10 1=3.4

10 Y ( I ) вда*х ( П + U*X ( I - i ) +B1*Y ( 1 - 1) + 32*Y ( 1-2 ) RFTURv ewd

subroutine bpscl ( alpha.Ofiec.of a0, a1 |b1,e2}

ПРдПрО-РдМЧд ВЫЧИСЛЕНИЯ ПАРАМЕТРОВ МОДЕЛИРОВАНИЯ ПСп С КФ EXP(-A.PHA*AB>(TAY))*(C0s«0ME3*TAY)-AuPhA/0HEG*SIN(0MEG*A9S(TftY)!)

ALPHA - ПА?АМ5ТР КФ

О НЕ С ПАРЛНЕ'Р КФ

О - ИНТЕРВАЛ ДИСКРЕТИЗАЦИИ

ая,а1,31,в2 - параметры алгоритма моделирования - subroutine: *сеч2

GI=ALP4A*D C5 = 0hf?:*0 рсе*р«-с1) н s = s i w ( g 2 ) HC=C0S(G2)

his I-P*M + 4*P**2* ALPHA/OMEG*HS*HC

H2=P*<3**2"l)^HC-ALPHA/0MEG+{l+P**2)»P*bS B2=-P**2 В Js2#p*HC

H3 = SQpr( H 1 * * 2 *H2 **2) IMH1.3T.H3) -30 TO 10 A*. = SQR r ( ( Hl + HJ) /2 > GO TO 20 10 AesSfiPT((H1-И 3)/2) a 1sh2/a 0 RETURN emd

SUBROUTINE AVX{N,X,AX)

ПОДПРО-РАМ-ЧА ОЦЕНКИ МАТ, ОЖИДАНИЯ ПСП

X ( N ' - ВХОДНОЕ МАССИВ

N - РАЗМЕР МАССИВА

АХ - 04ЕНКА МАТ, ОЖИДАНИЙ ПСП

DIMENSION X(N> Ss0 .

DO 10 1=1,N S-S+X(I) SsS/Ft.OAT ( M ) AX = S

return

esjd

SUBROUTINE СЕ-ЧТ(N « X * Л X)

подпро-рамна цент°ир0ванич псп

X f N > - ВХОДНОЙ И ВЫХОДНОЙ МАССИВЫ N - РАЗМЕР МАССИВА

ах - оценка мат, ожидания псп

DIMENSION *(N> On 10 1 = 1 Н X ! I ) s к { I ) - A X rgturw E N' D

SUBROUTINE XH0M(N.X|VAR,SIGX,3A$|CEC)

ПОДПРО'РАМ^А ОЦЕНКИ МОМЕНТНЫХ ХАРАКТЕРИСТИК ПСП

X(N) - ВХОДНОЙ МАССИВ (МЕНТРи'ОВАННЫй) N - РАЗМЕР МАССИВА VAR « ДИСП=РСИЯ ПСП SIGX _ СКО ПСП

С AS « КОЭФФИЦИЕНТ АСИММЕТРИИ СЁС - КОЭФФИЦИЕНТ ЭКСЦЕССА

DIMENSION X(N)

V A R = 0 .

С A S = 0 ,

срс = й ,

DO 10 1=1,4

VAR=VA*+X(I)**2

CAS=OA5+X{I)**3

СЕС=ска+х < n**4

ZaFLOAT(N)

VAR = VM/{Z'i)

SjGXssiRTl^ar)

CAS=(CAS/Z?/3IGX**3

CEC=(CEC/Z)/VAr**2"3,

return

END

subroutine 05sp(n.x ,var,sigx)

подпрограмма оценки дисперсии и ско псп

x{n ? - входной массив (центрированный) n " размер массива var П. дисперсия псп sigх _ ско пап

dimension к(я) v?r=0.

s t gxcff. 01 40 1=1,4 40 var=va*+x(i)**2 var=var/fl3at(n-1) sigxesirr (v'ar) rp.turk;

end

subroutine dn0rm(4,x,sigx) п0дпр0-рам1а нормирования псп

x(n) - входное массив

n - размер массива sigx - параметр нормирования

DIMENSION dr- 10 i = 1,m Xfl)/SIGX RETURN

end

subroutine с 0 r f(n » m , x , с 0 r i с 0 r1)

подпрограмма оценки автокорреляционной функции реп

x с n) - входной массив (центрированный) n - размер массива cor (м} - массив нкф

м - размер массива

c0r1 - дисперсия псп

dimension м n),с о s(м) Jk0 s = 0 • len-j

do 20 1=1,w sr.s + x<i ) *x(ю) SsS/tl-t)

jsj+1

Cor(j j =s

i г{j,l t.m ) go to 10 CpR lpfOR ( 1) do j=1,1 c0r(0?=c0r(j)/c0r1 r p t u r ы end

subroutine corr(n » x ,h ,cor1) подпрограмма Оценки нкф релейным нетодог-

x(n) - 8x0дной массив (центрированный)

n - рдзмер массива cor 1 (м) - массив нкф м - Размер массива

41mensi0n 7 ( n ) » с or 1 (м) j = 0

s50 . inn-m

do ui 1 = 1, ц! а гтх < i * j )

IP(a,FJ,0,) go to 10 CONTinje j = j + 1 cor 1 (jj! =s

i г(0» lt,m) go to 5 ci=cosl(1) do 20 j=i ,.1 c0r1(j)=c0.41 ! j ) /с 1

r f: turk: emd

s u b r 0 u гi ne ccf(x,r,n,m,cp)

подпро-рам.ча оценки взаимной <ч>

x f n) - входной массив y(n > - входной массив n - размер массива cf(м) - массив кф

м - размер массива dimension x(n),y(м),df(м)

jя0 . вк0 . с = 0 • 40 Ая0. lsn™j

do 10 1 = 1.1: 10 а = а +х { i ) *y ( i +0) jrj + 1

CF(J)=A/FLOAT(L) if(j.lt.m) go to 40 do 2i) i = l,4 b = b + x(i) 20 сcc +y1i >

ввв/fi oat(n) csc/fl3at(n) do 50 3=1,1 50 cf<j)bsf(j>-b+c r г t u r emd

ПРИЛОЖЕНИЕ 2

SUBROUTINE ERRtM.r ,Y1 ,S.E'

a

q подпро-раша вычисления погрешностей

с '

q y(.m) « входной пассив

1 y 1 (м) - входной массив

о м - размер входных массивов

с S - СР;ДНЕ<ВАДРАТИЧЕСКОЕ ЗНАЧЕНИЕ ПОГРЕШНОСТИ

д е - максимальное значение модуля абсолютной погрешности

с

dimension г (м) ,y мм ) et;abs ( г ( i)-ук 1 ) ) s к е * #

Si^Y j 1)**2 00 20 1=2,1 Е 1 ssABs ( Y С I » -Y 1 < I > ) IF(E«E1) 15,17,17 \Ъ E=Et {7 S'sS + E 1 *»2

2 0 S1 = 51 * Г(I)**2

s rtsqr Т ( S/s 1 )

returh

end

SUBROUTINE SERR(M,C0R,C0RU,SA:R»RSAER)

с

0 ПОДПРОГРАММА вычисления погрелностей f*

d COR(M) - ВХОДНОЙ МАССИВ

С CORUI'M - ЗХОДНОЙ МАССИВ

с М - РАЗМЕР ВХОДНЫХ МАССИВОВ

a SAER - КВАДРАТИЧЕСКОЕ ЗНАЧЕНИЕ ПОГРЕШНОСТИ

о RSAER - СРЕДНЕКВАДРАТИЧЕСКОЕ ЗНАЧЕНИЕ ПОГРЕШНОСТИ

С

DIMENSION COR(Ml,CORU(М) SUM=0. SUD=0. 00 10 1=1,4 SUM=SU1*C0^(I>**2 10 SUDSSU3*C0*0(I)**2 SAERsS JM/SJD RSAERiSQRT <$ AER) RETURN END

subroutine apexp(m,tay»rx»ersialpha,err)

Й

с АППР0КЗИМА4ИЯ КФ ФУНКЦИЕЙ EXP ( ~А|_РНА*АВ$ ( TAY > )

С

с rx(m> 5 массив значений к'

о т A y(и} - МАССИВ аргументов кф

с м - р&змер массива

с EPS « ПОГРЕШНОСТЬ ОЦЕНКИ ALPHA

с alpha - паэаметр аппроксимируощей функции

с err * <вадлдтическдя погрешность аппроксимации кф

3

DIMENSION Т A Y (М) » R X (М) S 3 = 0, S4 = 0, 5 АМ=ALPHA S 1 = 0, 52 = 0» 00101=1,1 А 1 =ЕХР(~AN*TAY(I))

• RsRX(I)-Л1

АсТАу (I)*А 1 b=r * а

С s В * Т а у ( i ) D«А**? - ■ S1=SU3 10S2=S2+D-C

ALPHA = AN-S 1/S2

I г(ABS(ALP iA~ AN).CT.EPS) GO T3 5 Da 20 1=1,1

S3=S3+(RX(I)-EXP(-TAY ( i ) *ALPHA) )**2 20 S4=S4**X(I)**2 EpR=S3/S4

RETURm

end

SUBROUTINE APEXP2(M,TAY .RX,EPS,ALPHA,ERR)

a

а АППРОКСИМАЦИИ КФ ФУНКЦИЕЙ EXP(-ALPHA*abs(TAY)) *(1 + AlPHA*a3S f T a Y ) '

с

с TAY < M) - МАССИВ АРГУМЕНТОВ КФ

a rx(M) » МАССИВ ЗНАЧЕНИЙ КФ

с м - РАЗМЕР МАССИВА

с EPS - ПОГРЕШНОСТЬ ОЦЕНКИ ALPHA

a alpha „ параметр аппрокси"ируощей функции

с ERR ™ <ВАД3АТИЧЕСКДЯ ПОГРЕШНОСТЬ АППРОКСИМАЦИИ КФ

с

DIMENSION UY(M),RX(M>

53 = 0,

54 = 0,

5 АИ = ALP НА

51 = 0,

52 = 0,

DO 10 1=1,1

АЯ=ехр(-А1)

A1=AN*TAV(U

R"RX(f)-A2*(l.+Al)

A3=TAv(I)**2

А в A 3 # a 2

BsR#AN*A

CRR*(1,-A1)

D = A2*(A 1#*2)

S1=S1+B

10 S 2 = S2 * A*(0*0) ALPHASAN-S1/S2

IF(AB.(;(ALP",AIAN)..ST,EPS) GO TO 5 DO 20 1=1,1

S3*S3* IRX(I)-EXP(-TAY(I)nALPHA)*(It + TAY(I J *AlPHA> > #*2'

20 S4 = S4 + 3 X(I)**2 ERR = SVS4 RFTURm END

SUBROUTINE APEXP3(M,TAY »RX,EPS,A^PHA.ERP)

О

a АППРОКЗИМА4ИЯ КФ ФУНКЦИЕЙ

с EXP(-A.PHA*ABS(TAY))*( 1-ALPHA*А8$ (TAY))

a

О TAY(M) - МАССИВ АРГУМЕНТОВ КФ

О RX(Mj в МАССИВ ЗНАЧЕНИЙ К*

О М - РАЗМЕР МАССИВА

О EPS - ТОРРьШИОСТЬ ОЦЕНКИ ALPHA С ALPHA - ПАРАМЕТР АППРОКСИМИРУОЩЕй ФУНКЦИИ

3 ERR * <ВАДРАТИЧЕСКДЯ ПОГРЕШНОСТЬ АППРОКС ИМAUИИ «Ф б

OIMENSION Гду (М) (М )

53 = 0,

54 = 0,

5 AK'sALP4A

51 = 0.

52 = 0,

DO 10 1=1,1 A1 = AN*TAY ( I ) А2=ЕХр{-А1)

A3=A2*(i.-Al) R!5RX{|)~A3

Si=51*A24TAY(I)»(2,-Al)*R 10 S2 = S2*A2*(Г А у(I J **2)*(A2# ((2,»Al)**2)-R* (3,-Aj } ) ALPHAKAN-S 1/S2

IF(ABs(ALPHA-AN).GT.EPS) GO TO 5 DO 20 1=1,1

S3=S3MRX(I>-EXP(-TAY(I)*ALPHA)*{1,-TAY(I)*ALPHA))**2

20 S4 = S4*U ( I ) **? EPR=S3/S4 RTTURkj

end

subroutine aexcos(cor ,m,alf,01eg.eps,t)

С

с подпрограмм аппроксимации нкф функцией

о вида FXP(-Alf*ABS(TAy))*COS(01EG*TAY) с

с cor си j - массив нкф

с м, - размер массива

a alf . омеg - начальное и конечное значения параметра

- 248 - .......

С АППРОКСИМИРУЮЩЕЙ МОДЕЛИ

С EPS " АБСОЛЮТНАЯ ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЯ ПАРАМЕТРОВ МОДЕЛИ

С ТСМ) « МАССИВ АРГ/МЕНТОВ НКФ

DIMENSION -PR (М) ,Т(М)

9 АН = 0« А = 0 . Н = 0 •

В = 0 . Сз0 .

00 IB >1 T1=T(J)

ys£xp ( -т l*al,f ) Pr.Y*OOS (f I * О ИЕс ) 0BY*SIM(T1*0MEG) XsCQRfJJ-P

AM=AM+X*P*TI

HsH+X*3*T1 A =; A * T J **2*<J* (P-X )

B»B+T 1 (p-x)

10 C=C+T 1 **2*<Q**2+X*P) Q с A * С « 3 * * 2

OsALF PcOMEC

ALF=AL"-{C*AM-B*H)/Q ' OMEGsOIEG-<sB*A4+A*H)/Q G И 0 - A L -ЕяР-OMcG GsABS f 2) ЕяABS(;)

IF((G-5ps).GT.0) GO TO 9 IF ( IE-EPS)•GT.0) GO TO 9 RFTURN EmQ

SUBROUTINE COR&0(N,NK,X.OT,I TAU , RX,DISP)

a

а ВЫЧИСЛЕНИЕ НКФ HPД ПРОЦЕССА с КУООЧНО - ПОСТОЯННЫМ

с ВОССТдНОВЛЁНИЕМ - ПОСЛЕДОВАТЕЛЬНЫЙ АЛГОРИТМ

с

с n - число существенных отсчетов

з N к - РХОДНЗй ПАРАМЕТР, ОПРЕДЕЛЯЕТСЯ ИЗ УСЛОВИЯ jt<N к)< = jt<n)~*

С X(N> - МАССИВ ОТСЧЕТОВ НРД ПРОЦЕССА

С jТ(N) - МАССИВ МЕТОК ВРЕМЕНИ

б . i т а и * номер ОРДИНАТЫ К* о RX - ЗНАЧЕНИЕ КФ В ТОЧКЕ IТ A(J

С D!Sp * ДИСПЕРСИЯ НРД ПРОЦЕССА

ft

dimension *(N).JT(N) j i = I т a j - i R X s 0 , Ks 1

ni=nk'«i

dc 20 1=1,41 I0 = JT(I)*J 1 n = JT{I + l)+Jl

8 KrK + 1

I F (JT(<)«I 3) 8,10,9

9 КзК-1

10 LaK

11 UL*1

I г(JT i ) - l 1 ) 11.12,12

12 1=1-1 Y п 0 .

DO 15 Ч1яК»1 15 YsY + XMlj*<jT(Ml + l)-JT(Ml> )

Y»Y-X{<)# <I0-JT (К) )-X(L> * < JT()1> 20 RY=RX*K(!)*Y

RK=RX/-LQAT(JT(N)-J1)

RX=RX/DISP

RETURN

END

SUBROUTINE CORQl (N ,N1 ,x . JT ,H, rx'idisp)

ВЫЧИСЛЕНИЕ НКФ НРД ПРОЦЕССА С КУСОЧНО-ЛУНЕйНЫМ ВОССТАНОВЛЕНИЕ4 -ПОСЛЕДОВАТЕЛЬНЫЙ АЛГОРИТМ

N- ЧИСЛО СОБСТВЕННЫХ ОТСЧЕТОВ

NK - ВХОДНОЙ ПАРАМЕТР ОПРЕДЕЛЯЕТСЯ ИЗ УСЛОВИЯ JT ( N К ) < s J Т ( М ) - М

X (N) - МАССИВ ОТСЧЕТОВ НРД ПРОЦЕССА

J Т < N) - МАЗвИВ МЕТОК ВРЕМЕНИ

J - НОМЕР ОРДИНАТЫ КФ

RX - ЗНАЧЕНИЕ КФ 3 ТОЧКЕ J

DISP г. ДИСПЕРСИЯ НРД ПРОЦЕССА

DIMENSION * ( N) , JT(N 1) tRX(M) ОП 2.3 Jci.l j 1 = j-i

Y я 0 •

1 s 0 . ' Kb 1

7 IbI + 1 IffsJT(I)*Jl ii=OTfi + u+Ji

8 KsK+l

IF(JT(<)-IЪ) 8,10 , 9

9 КгК - 1

10 LkK

11 L=L+1

IF (L,qE,N) CO TO 13 IF(JT(.)-11)11,12,12

12 1=1-1

GO TO 14

13 jt<u1) = i1

14 Y 1 = 0,

Y 2 = 0 »

DO 15 41-aKil

DJ = <Xn + t>"XU)>/(OT(l + l)-OT(I)) Dk=<X{<+1>*X(K)>/(JT(K+l)-JT(<) ) DL=<X!.+1)"X(U))/(JT(L+1)"JT(L)) Yl=(X(11)+XfHUl))*(JT(Ml + l)-JT(Ml))/2.

Y2=X{Hl)*(JT(Ml*l)#*2-5*jT(MU*JT(Ml*t)-4*(jT(Ml)**2))/6.

15 y? = Y2*X(M1 + 1)'42*(JT(M1M>**2) 'JT (Ml) *jT(Ml+l?~jT(Mi)**2)/6, Y1 = Y1«(X(K>-DK*JT(K))*<JT(I) + J-JT <K>)

Y i = Y J<OK* ( ( JT (I ) + J ) **?-JT f К > **2)/2 ,

Y 1=Y Ю1* < JT (L+l ) **2-(JT ( 1 + 1 ) +J) **2 >/2,

Y ? = Y 2 (X (К' - QK#JT (К ) >*<(0ni)+J)**2-JT{H>**2>/2, Y2 = Y2«DK*MjT(i)+J)**3-JT<K)**3)/3,

Y5 = Y2-tXtL?-DL*JT(L)J*(JT(L*i)**2-{JT(I-»l)4-0)**2)/2. Y? = Y2*OL*(JT(L + 1)* * 3- ( JT{ 1 + 1)+ J)**3 '/3,

Y s Y 1 * { X (i)*DI *(JT ( I )-J))+Y2*D I IF ( I UJT (N> ) 7,20, 20

20 Rx(0)er/F13IT(JT(N)-Jl) DISPBRX(I) DO J=t,1 30 Ry(J}~RX{JJ/DISP RPTURfi EMD

SUBROUTINE RCORF (X ,JT ,N 1 ,N2 ,S «II ,M|M1)

С

а ВЫЧИСЛЕНИЕ НКФ ПО РЕЛЕЙНОМУ АЛГОРИТМУ С ИСПОЛЬЗОВАНИЕМ

(. ИИТЕРВАЛЬНЗЙ КОРРЕЛЯЦИИ (ПОСЛЕДОВАТЕЛЬНЫ?!)

о

3 XINM - МАССИВ ЗНАЧЕНИЙ Н^Д ПэОЦЕССА

g JT(N1 j - МАССИВ МЕТОК ВРЕМЕНИ

о N1 - РАЗМЕР МАССИВА

О N2 - 1МФРО30Й ЭКВИВАЛЕНТ ДЛИТЕЛЬНОСТИ РЕАЛИЗАЦИИ

О S'Ml} - МАССИВ ЗНАЧЕНИЙ НКФ

О II (МП - МАССИВ НЕУСРЕДН'ЕННЫХ ЗНАЧЕНИЙ

С Mi - РАЗМЕ? МАССИВА

ё М - ЦИФРОВОЙ ЭКВИВАЛЕНТ TAYKMAX

a

DIMENSION X(Nl),jT(Nl).S(M1),11(МИ do 5 l = 1 , m j. S (L ) й 5 . 5 I I (L|

DO 10 JsJ.I N3sN2«J+l" 00 10 1 = 1, 43

IK»jt(i+j"i>-JT сi)*i

IF(X(f+J-l),EEl.0.) GO TO 7 S tIK)=5{iK)+ X(I)*X(I + J-l)/ABS(X(I + J-l)) 7 I ! ( IК ) s I i ( I К ) + 1 10 С n N T J w J E

DO 20 .=1,11

IF(II ( .).E3.0) GO TO 20

siu =s (l) /float ! i i (d ) 20 CQNTinjE s i = s 11) DO 15 .=1,11 15 S(l)=5(l) /Si RPTUR'-; E')D

SUBROUTINE pb0RF(X,JT,K|l,N2,H,II,M,Ml)'

ОПРЕДЕЛЕНИЕ НКф НРД ПСП ПОЛЯРНЫМ МЕТОДОВ ПРИ АДАПТИВНОЙ ДИСКРРТИЗА4ИИ БЕЗ ВОССТАНОВЛЕНИЯ ПРОМЕЖУТОЧНЫХ ОТСЧЕТОВ (ПОСЛЕЦОВАТЕЛЬНЫЙ АЛГОРИТМ)

X ( N 1 ) - ЗНАЧЕНИЯ НРД ПСП JT(N1) - МЕТКИ ВРЕМЕНИ N1 ~ РАЗМЕР ЗХОДНЫХ МАССИВОВ N2 - ДЛИТЕЛЬНОСТЬ РЕАЛИЗАЦИИ Н ! М 1) - МАССИВ НКФ

II (МП - МАССИВ НЕУСРЕДНЕННЫХ ЗНАЧЕНИЙ г^КФ Ml - Р А 3 М Е J МАССИВОВ М - Ч^ЗЛО ОРДИНАТ НА ТКМАХ

DIMENSION <(Nl),jT(Nl),H(Ml),II(Mn DO 5 i. = l ,П1-H I L ) = ? . 5 I I

H( 1 )яр.ОАТ(N2) II U > = v| 2 DO 13 0=2,1 N3 = N2 и 3 + 1 DO 10 1 = 1 , M3 IK = JT(I+ J -1)-ОТ(I) + 1 IF (X ( П*X<I+0-1) .EQ.0.) GO TO 7 н{IKj =4(IK? *|X(I)*X( I + O-l >)/ABS(X(I)#x(i + o-l)) 7 Ц(IK! = 1I(IK > + 1

10 continue do 20 .=i,ii

IP ( I I ! .) .E3.0) GO TO 20 H(L)=-C0S(5,Hl59*H(L)/Fl,OAT(II (U ) > 20 CONTlNJE RFTURM EMD

SUBROUTINE PCORF 1 (L ,JT ,N,Z,M,H1)

ВЫЧИСЛЕНИЕ НКФ ПРИ СПОРАДИЧЕСКОЙ ДИСКРЕТИЗАЦИИ п - ПОЛЯэНЫЙ МЕТОД РЕАЛЬНОГО ВРЕМЕНИ

JT(L) - МАССИВ ИНТЕРВАЛОВ ВРЕМЕНИ L - РАЗМЕР МАССИВА

N - ЧИСЛО ОТСЧЕТОВ , ОБРАБАТЫВАЕМЫХ В ОГНОМ 4ИМЕ ОПРОСА Z(N> - РАБОЧИЙ МАССИВ ДЛЯ ХРАНЕНИЯ ИНТЕРВАЛЬНЫХ СУММ Н(М) « МАССИВ НКФ

НИМ) - МАССИВ ИНТЕРВАЛЬНЫХ К Ф • (НЕ у СРЕ ДН Е ННЫ X )

М - РАЗМЕР МАССИВОВ

DIMENSION JT <U .Z с N >|H(M> |H1 (1) DO 10 1=1»^ 10 Z ( I )=c .

DO 2 0 1=1,1 H! I >=p, 20 HI(I)ad, Kp0

DO 55 1 = 1 , V . KoK+ 1

I F ( К,G T,N) K = 1 Z ! К ) я ?■, S»-l » jsi К 1 =K •25 К1 = КЫ

I F ( К 1 t £ (J , О I K1 = N IF (Kl.EQ.K) GO TO 55 Ss-S

z <K1}-I(K1 I+JT( I )

35 INJ,(:T,Z«4! » SO TO 25 HIJ)eH(Jj + s HI(J)sH1(J) + 1. JrJM GO TO 35 5 5 CONTIMJE

DO 60 J=l,1

IF (HI { J) ) GO TO 60

H(J)sslNi1.5708*H(J)/HI(JM 6 0 CSNTINJE RFTURM END

SUBROUTINE MCORD (NiNK ,X .JT . ITAU*,RX)

вычисление одного значения НКФ для нрд 5 Г С П П»И КУСОЧНЗ Г- постоянном восстановлении ПРОиЕССА

Мк - и/|ФРОЗОЙ ЭКВИВАЛЕНТ ДЛИТЕЛЬНОСТИ НАБЛЮДЕНИЯ ПСП

К;<= (Т^ТКМАХ )/DELTA ( ОЕ LT А-ИН ТЕРВ А Л ДИСКРЕТИЗАЦИИ) X(N) - МАССИВ ЗНАЧЕНИЙ НРД ПСП J Т ( N ) - МА-СИВ МЕТОК ВРЕМЕНИ НРД ПСП N - РАЗМЕР МАССИВА

ITAU - АРГУМЕНТ НКФ (НОМЕР ОРДИНАТЫ) R X " ЗНАЧЕНИЕ НКФ.ДЛЯ АРГУМЕНТА IТ AU

dimension x(n),jt(n) J 1=IT 4J-I RX = 0»

тх=0.

Кя I

N 1 = NK » 1

DO 20 1 = 1 ,.41

I?=JT?I)*J1

Il = JTU + l)*Jl

8 KcKM

IF(ОТ{<)-I a > 8,10,9

9 K=K-1

10 LoK

11 L«LM IF(JTf.Nli) 11.12,12

12 LssL-l YP0 .

DO 15 11 gК IL 15 YsY+X Ml ) * ( JT (Ml + 1 )-JT (Ml > )

YsY-X{<)#(I 0-J T(К) >-X(L)*(JT(l*1)-I1) RXsRXAX(j ) *Y 20 TX = TX*X(I) * * 2'*F L 0 A T ( JT ( I + 1) ~ J T ( I ) } R X = RX/Т X RFTURk;

end

SUBROUT I NE PCF(N,NK,X,jT,ITAUiY'»NN)

d

а ВЫЧИСЛЕНИЕ НКФ ПРИ СПОРАдИЧЕС<ОЙ ДИСКРЕТИЗАЦИИ

CJ ПОЛЯРНОМ МЕТОДОМ (ПОСЛЕДОВАТЕЛЬНЫЙ АЛГОРИТМ)

d

d N* - ОПРЕДЕЛЯЕТСЯ ИЗ УСЛОВИЯ Т(ыЮ<-Т(к)-ТКМАХ

й X(N J » МАССИВ ЗНАКОВ НРД ПРОЦЕССА ПРИ СГОРдДИЧЕСКОй ДИСКРЕТИЗАЦИИ

a JT(N) - МАССИВ МЕТОК ВРЕМЕНИ

a N - ЧИСЛО 'СУЩЕСТВЕННЫХ ОТСЧЕТОВ

о ITAU - НОМЕР ОРДИНАТЫ нКФ

о Y - ЗНАЧЕНИЕ ОРДИНАТЫ НКФ В ТОЧКЕ |ТАи

С NM " иИФРОЗОй ЭКВИВАЛЕНТ ДЛИТЕЛЬНОСТИ РЕАЛИЗАЦИИ

Я

dimension x<n)•jt с n > I TAU lei T'AU-1

Y s 0 « Kr 1

DO 22 I s 1 , ЧК

I e = JT fI>*IГди t I 5 = JT<I + i)*lTAul

8 КвК + 1

IF (JT{О-I 3)8,10,9

9 KsK-l

10 LBK

11 LbL*1

IF(JT(.)- I 1 ) 11,14,12

12 L=L-1

14 Y 1 =0 »

00 15 HsKiL

15 Y 1=Y UX (M 1> * < JT (M l + l )-JT (Ml ) ) Y1=Y|«X(K)*{I0-JT(K))-X(l,)*{JT(L*l>-H)

22 Y = Y+X(I)*Y*

YbSIN t1.5708*Y/FLOAT(NN))

RETURm

end

SUBROUTINE PC0R0(NI,NK,X.jr, I T А и V R X )

ri

а ВЫЧИСЛЕНИЕ ОДНОГО. ЗНАЧЕНИЯ HK * для HP д ГСП ПОЛЯРНЫМ А ЛГОР ИТ 10 *

о ДЛЯ АЛ ЭКСТрАПОЛЯТОРОМ НуЛЕВ0г0 ПОРЯДКА

с

в • NK= (Т-ТКМА*)/0Е1.ТА - ЦИФРОВОЙ ЭКВИВАЛЕНТ О ДЛИТЕЛЬНОСТИ НАБЛЮДЕНИЯ ПСП

О DELTA - ИНТЕРВАЛ ДИСКРЕТИЗАЦИИ

о X ! N ' - МАССИВ ЗНАЧЕНИЙ НРД ПСП

О J Т(N) . МАССИВ МЕТОК ВРЕМЕНИ НРД ПСП

б N - РАЗМЕР МАССИВА

с ITAU - АРГУМЕНТ НКф (НОМЕР ОРДИНАТЫ)

3 RX - ЗНАЧЕНИЕ НКФ ДЛЯ АРГУМЕНТА IТ A U

8

DIMENSION * (N) ,-JT (N) J! = IT/1J-1 RX = 0, К а 1

N!=NK«J DO 20 I = i ,Mj I 0sJf(I)*J1

iI=jtti+i)

8 КЙКМ

IF(JT|<)„IЪ) 8,10,9

9 KsK-1 10 L=K П LsL+i

I F ( J T ( •) -SI I ) 11,12,12 12 LeL-1 Y и 0 t

DO 15 11 sI' , L 15 YsYMJT(M14)-jT(f11))*X(M1)/A3S(X(M1))

Y=Y-n3-JT(K))*X(K)/A8S(X?K))'(JT<L+l)-Il)*X(L)/ABS(X{L)) 20 RX=RX*r*X(I)/ABS(X(I)> PI=4,*ATAN<1,>

RX = SI»j(RX*Pl/(2t*FL0AT{JT«N)-Jl)))

RrTURv END

SUBROUTINE RCORD (N,NK,X,JT,ITAU',RX)

fl

о ВЫЧИСЛЕНИЕ ОДНОГО ЗНАЧЕНИЯ НКФ ДЛЯ НРД >г СП РЕЛЕЙНЫМ АЛГОРИТМОВ

а ДЛЯ А о ЭКСТРАПОЛЯТОРОН ну ЛЕВО г0 ПОРЯДКА

3

a NK=(Т-ТКМАХ)/DELTA - ЦИФРОВОЙ ЭКВИВАЛЕНТ

а ЛПИТЕПЬНОСТИ НАБЛЮДЕНИЯ ПСП

С DELTA - ИНТЕРВАЛ ДИСКРЕТИЗАЦИИ в X(N) - МАССИВ ЗНАЧЕНИЙ НРД ПСП

О ОТ(N} - МАССИВ МЕТОК ВРЕМЕНИ НРД ПСП

С N - РАЗМЕР МАССИВА

С ITAU « АРГУМЕНТ НК«? (НОМЕР ОРДИНАТЫ)

О RX - ЗНАЧЕНИЕ НКФ ДЛЯ АРГУМЕНТА jTAU

С

DI-ENSION X(N),JT(N) Ji = iT £ J-i

rx = 0.

ТХ = 0, К в 1

n ! = n к r> 1 00 2 0 I в 1,М J iii)*j1 ij =