1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
#! /usr/bin/env python
import openturns as ot
ot.TESTPREAMBLE()
Tmin = 0.0
deltaT = 0.1
N = 11
# Initialization of the RegularGrid timeGrid1
timeGrid = ot.RegularGrid(Tmin, deltaT, N)
# Functional basis
basisDimension = 10
basis = []
coefficients = ot.DistributionCollection(basisDimension)
for i in range(basisDimension):
basis.append(ot.SymbolicFunction("x", "sin(" + str(i) + "*x)"))
coefficients[i] = ot.Normal(0.0, (1.0 + i))
process = ot.FunctionalBasisProcess(ot.JointDistribution(coefficients), ot.Basis(basis))
# Setting the timeGrid
process.setTimeGrid(timeGrid)
# Print of the process
print("process = ", process)
# Initialization of the TimeSeries
timeSerie = process.getRealization()
# Full prints
print("timeSerie = ", timeSerie)
# Some steps further
stepNumber = 4
print("One future=", process.getFuture(stepNumber))
size = 3
print("Some futures=", process.getFuture(stepNumber, size))
|
