Gramian角场生成图像的python代码

以下是一个使用Gramian角场生成图像的Python代码示例：

import numpy as np
import matplotlib.pyplot as plt

def gramian_angular_field(X):
    # 归一化输入数据
    X = (X - np.min(X)) / (np.max(X) - np.min(X))
    
    # 计算sin和cos矩阵
    sin_X = np.sin(X)
    cos_X = np.cos(X)
    
    # 计算Gramian角场矩阵
    GAF = sin_X @ cos_X.T - cos_X @ sin_X.T
    
    return GAF

# 生成示例数据
t = np.linspace(0, 2*np.pi, 100)
X = np.outer(t, np.ones_like(t))

# 计算Gramian角场矩阵
GAF = gramian_angular_field(X)

# 绘制Gramian角场图像
plt.imshow(GAF, cmap='hot', origin='lower')
plt.colorbar()
plt.title('Gramian Angular Field')
plt.show()

这段代码首先定义了一个`gramian_angular_field`函数，用于计算Gramian角场矩阵。然后，生成了一个示例数据`X`，并调用`gramian_angular_field`函数计算Gramian角场矩阵`GAF`。最后，使用`imshow`函数绘制了Gramian角场图像。

请注意，这只是一个简单的示例代码，实际应用中可能需要根据具体需求进行调整和优化。