想要生成一个dataset做一个聚类的simulation,X是n行p列的样本数据,p为X的特征个数,设定前50个feature是relevant features,其他的特征都是noise features,总共有三类,每一类的样本数是20,其中\[x_{ij}\]服从均值为\[u_{ij}\],方差为1的正态分布,

取u=0.6和0.7,p=200, 500和1000,求问怎么生成这样的数据集?
我用了MASS包的mvrnorm函数,但是需要有协方差矩阵,如果P=200,协方差矩阵岂不是有40000个值,感觉不对,求教大神指点