covdesign.py－20170728

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Thu Jul 27 16:06:30 2017

@author: vicky
"""

import numpy as np
(nx,px)=np.shape(xdesign)

xdesmax=np.zeros((nx,14))
d=np.zeros((nx,14))
for i in range(0,nx):
    xdesmax[i,0]=np.max(xdesign[i,2:4])
    d=np.append(d,np.where(xdesign[i,2:4]==np.max(xdesign[i,2:4])))
    xdesmax[i,1]=np.max(xdesign[i,4:11])
    d=np.append(d,np.where(xdesign[i,4:11]==np.max(xdesign[i,4:11])),axis=0)
    xdesmax[i,2]=np.max(xdesign[i,11:13])
    d[i,3]=np.array(np.where(xdesign[i,11:13]==np.max(xdesign[i,11:13])))[0,0]
#    xdesmax(i,4)=max(xdesign(i,14:15));
#    d(i,4)=find(xdesign(i,14:15)==max(xdesign(i,14:15)));
#    xdesmax(i,5)=max(xdesign(i,16:24));
#    d(i,5)=find(xdesign(i,16:24)==max(xdesign(i,16:24)));
#    xdesmax(i,6)=max(xdesign(i,25:26));
#    d(i,6)=find(xdesign(i,25:26)==max(xdesign(i,25:26)));
#    xdesmax(i,7)=max(xdesign(i,27:28));
#    d(i,7)=find(xdesign(i,27:28)==max(xdesign(i,27:28)));
#    xdesmax(i,8)=max(xdesign(i,29:30));
#    d(i,8)=find(xdesign(i,29:30)==max(xdesign(i,29:30)));
#    xdesmax(i,9)=max(xdesign(i,31:32));
#    d(i,9)=find(xdesign(i,31:32)==max(xdesign(i,31:32)));
#    xdesmax(i,10)=max(xdesign(i,33:34));
#    d(i,10)=find(xdesign(i,33:34)==max(xdesign(i,33:34)));
#    xdesmax(i,11)=max(xdesign(i,35:37));
#    d(i,11)=find(xdesign(i,35:37)==max(xdesign(i,35:37)));
#    xdesmax(i,12)=max(xdesign(i,38:41));
#    d(i,12)=find(xdesign(i,38:41)==max(xdesign(i,38:41)));
#    xdesmax(i,13)=max(xdesign(i,42:43));
#    d(i,13)=find(xdesign(i,42:43)==max(xdesign(i,42:43)));
#    xdesmax(i,14)=max(xdesign(i,44:45))
#    d(i,14)=find(xdesign(i,44:45)==max(xdesign(i,44:45)))