计算Fisher信息之基础矩阵(一)

import numpy as np
import scipy as sp
import matplotlib.pyplot as plt
stepsize = 0.01
'''
这里有个问题 就是特征值的排列方式可能跟下一个的顺序不一样,所以就会导致求导时候的错位。
这里需要在求导的时候,搞一下吧。
另外,矩阵的特征值,会不会突然变异?有跃迁?需要设置一个,检测的一些。
'''
def EigVect(Matrix):
    eigs,vectors = np.linalg.eig(Matrix)
    return eigs,vectors

def test1ForContinueOfEig(eigss=[]):
    #这里的是两个或者多个特征值,我们假设前面已经将其复位好了
    plt.figure()
    eigss1 = np.array(eigss)
    for i in eigss1.shape[1]:
        plt.plot(eigss1[:,i])
    plt.title('tezhengzhi')
    plt.xlabel('x')
    plt.show()

def EigDaoshu(Eig1=1,Eig2=2,stepsize=0.01):
    if abs((Eig1-Eig2)/Eig1)>5:
      #  print('特征值变化超过10倍,出现奇异的点了')
        raise Exception('特征值变化超过原来的10倍,出现奇异的点了')
    return (Eig1-Eig2)/stepsize

def VectDaoshu(Vect1=[2,2],Vect2=[1,1],stepsize=0.01):
    result = (Vect1-Vect2)/Vect1
    for item in result:
        if abs(item)>5:
            raise Exception('向量变化超过原来的十倍啦,应该是出现奇异的点了')
    return (Vect1-Vect2)/stepsize

def isnear(eigs):
    if abs(eigs[0]-eigs[1])<stepsize:
        raise Exception('两个特征值很相近,无法得出哪个特征值是哪个,换换矩阵吧')
    return True
'''
解决的情形是[0.4,0.6] [0.61,0.39] 这里第一位和第二位的就需要换一个位置
'''
def sorteigs(eig1,eig2,eig3,vect1,vect2,vect3):
    if abs(eig1[0]-eig2[0])<abs(eig1[0]-eig2[1]):
        pass
    else:
        print('需要换位2....')
        eig2[0],eig2[1] = eig2[1],eig2[0]
        vect2[0],vect2[1] = vect2[1],vect2[0]
    if abs(eig2[0]-eig3[0])<abs(eig2[0]-eig3[1]):
        pass
    else:
        print('需要换位3....')
        eig3[0],eig3[1] = eig3[1],eig3[0]
        vect3[0],vect3[1] = vect3[1],vect3[0]
    return eig1,eig2,eig3,vect1,vect2,vect3

函数说明

函数函数名称备注
EigVect求特征值和特征向量
test1ForContinueOfEig测试特征值是否连续还没用到
EigDaoshu特征值导数
VectDaoshu特征向量导数
isnear是否相邻为下面函数所用
sorteigs特征值从小到大排序使得特征值变化连续
已标记关键词 清除标记
©️2020 CSDN 皮肤主题: 编程工作室 设计师:CSDN官方博客 返回首页