Python Model.as1DDistr示例

示例#1

文件： diffeq.py 项目： dkasak/pacal

    M.varschange(Y[i], Y[i+1])
print M
M.varschange(K, A)
M.plot()
#M.condition(Y[n], 2)
print M
M.eliminate(K)
print M
for i in xrange(n-1,-1,-1):
    M.eliminate(Y[i])
print M
M.eliminate(A)
print M
figure()
M.plot()
print M.nddistr.pdf(linspace(0,2.5,100))
X0 = BetaDistr(1, 1)
y = X0 * exp(A)
y.summary()
y.plot(label="Y0*exp(A)")
figure()
err = y.get_piecewise_pdf() - M.as1DDistr().get_piecewise_pdf()
err.plot()

show()
#M.varschange(X3, S3)
#print M
#M.condition(S3, 0.8)
#M.varschange(S2, X3)
#M.varschange(X3, X1)

示例#2

文件： diffeq.py 项目： ntankovic/pacal

    M.varschange(Y[i], Y[i + 1])
print M
M.varschange(K, A)
M.plot()
#M.condition(Y[n], 2)
print M
M.eliminate(K)
print M
for i in xrange(n - 1, -1, -1):
    M.eliminate(Y[i])
print M
M.eliminate(A)
print M
figure()
M.plot()
print M.nddistr.pdf(linspace(0, 2.5, 100))
X0 = BetaDistr(1, 1)
y = X0 * exp(A)
y.summary()
y.plot(label="Y0*exp(A)")
figure()
err = y.get_piecewise_pdf() - M.as1DDistr().get_piecewise_pdf()
err.plot()

show()
#M.varschange(X3, S3)
#print M
#M.condition(S3, 0.8)
#M.varschange(S2, X3)
#M.varschange(X3, X1)