load?gauss1.txt;
load?gauss2.txt;
load?gauss3.txt;
y1=gauss1;
y2=gauss2;
y3=gauss3;
x=-8:0.1:8;
figure（1）;
subplot（311）;
hist（y1x）;%統計直方圖；
hold?on;
plot（x1000*normpdf（x01）?）;%正態分布概率密度函數曲線
title（‘N（01）‘）;
xlabel（‘x‘）;ylabel（‘PDF‘）;
legend（‘仿真結果‘‘理論結果‘）;
hold?off

subplot（312）;