Ф (х) — функция Лапласа; &==2 5...%‚; а г— коэффици- ент корреляции между равномерными величинами. При других способах моделирования будут получать- ся векторы © другими совместными „плотностями (так как коэффициент корреляции и одномерные плотности не определяют однозначно распределение). Идеальным было бы при моделировании реального процесса зада- вать совместную плотность, исходя из физических со- ображений. Однако до тех пор пока о совместной плот- ности известно очень мало, естественно предположить из общих соображений, что эта плотность непрерывна отно- сительно своих аргументов. В этом смысле метод, пред- ложенный Г. Г. Сванидзе [Л. 87] с введением в выраже- ние плотности функции Дирака 6, что делает ее раз- рывной, является более далеким от действительно происходящего в реках процесса, чем наш или метод Е. Г. Блохинова и О. В. Сарманова [Л. 10], которые раскладывают совместную плотность по полиномам Ла- герра — (что можно сделать с любой непрерывной функцией). Указанный нами способ получения коррелированного случайного вектора является универсальным, т. е. его. можно применять для моделирования любого распреде- ления, предварительно получив функцию р(г) (аналити- ческую или в виде таблицы достаточного объема). Составление таблицы р(г) требует довольно значи- тельных затрат машинного времени. Например, на ЭЦВМ «Урал-2» с быстродействием 6000 операций в секунду расчет для одной пары Р(х) и Р;(х) при 7г от —0,9 до +0,9 с шагом 0,2 занимает 1 4. Но это затраты единовременные, и при существующих темпах развития вычислительной техники (увеличении быстро- действия ЭЦВМ) предложенный метод моделирования случайных — коррелированных — величин представляется перспективным и пригодным для практики. 5-4. ОБЩАЯ СХЕМА МОДЕЛИРОВАНИЯ ГИДРОЛОГИЧЕСКИХ РЯДОВ В НЕСКОЛЬКИХ СТВОРАХ ПРИ РАЗЛИЧНЫХ ВРЕМЕННЫХ ИНТЕРВАЛАХ ОСРЕДНЕНИЯ В & 5-2 были приведены иллюстративные расчеты многолетнего регулирования при использовании искус- ственных годовых стоков в нескольких створах Ангаро- Енисейского каскада ГЭС. Подобные расчеты можно 205