exit2.py－20170818

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Fri Aug 18 15:15:43 2017

@author: vicky
"""

[n,p]=np.shape(y) 
   
f=np.zeros((n,p), dtype=np.int)
for i in range(n-3): 
    for j in range (0,p):
        f[i+3,j]=y[i,j]+y[i+1,j]+y[i+2,j]+y[i+3,j]

df1=np.zeros((n,p), dtype=np.int)
df2=np.zeros((n,p), dtype=np.int)
for j in range (p):
    for i in range (n-1):
        df1[i+1,j]=f[i+1,j]-f[i,j]
    for i in range (n-2):
        df2[i+2,j]=f[i+2,j]-2*f[i+1,j]+f[i,j]
        
td2=range(0,n)
for j in range(p):
    fig=plt.plot(td2,f4[:,j])
    plt.title("fig"+np.str(j))
    plt.show()        

t1=np.array([259,259,240])
t2=np.array([336,330,302])
t3=np.array([287,287,287])
#tt=np.hstack((t,t))

for j in range(p):
    fig=plt.plot(td2,f[:,j],t1[j],f[t1[j],j],'x',t2[j],f[t2[j],j],'o',t3[j],f[t3[j],j],'v',td2,ff[:,j])
    plt.show()

fmean=zeros((1,p))
for j in range(p):
    fmean[0,j]=np.mean(f[36:97,j])

ff=f
for j in range(p):
    for i in range(t1[j],t3[j]):
        if ff[i,j]>fmean[0,j]:
            ff[i,j]=ff[i-1,j]*(1+m)

    #for i in range(t3[j]):