替换Z2 = mlab.bivariate_normal(X,Y,1.5,0.5,1,1)

mlab.bivariate_normal 函数可以生成一个二元正态分布的概率密度函数。其中的参数分别表示：

- X, Y：生成的网格点，即横纵坐标。
- 1.5, 0.5：正态分布的两个坐标轴的标准差。
- 1, 1：正态分布的两个坐标轴的相关系数。

如果要替换 Z2，可以使用其他的函数来生成概率密度函数。例如，可以使用 numpy 库的 random 函数生成随机概率密度函数。代码如下：

import numpy as np

Z2 = np.random.rand(*X.shape)  # 生成形状与 X 相同的随机矩阵

这样就可以将 Z2 替换为一个随机概率密度函数了。