%clear?all;clc;清屏
%產生服從正態分布（10，2）的隨機數
m=input（‘請輸入平均值：‘）;
n=input（‘請輸入標準差：‘）;
t=input（‘請輸入數據長度：‘）;
%產生正態分布的隨機數
for?i=1:t
????a=rand;
????b=rand;
????X1（i）=sqrt（（-2）*log（a））*cos（2*pi*b）;
????X2（i）=sqrt（（-2）*log（a））*cos（2*pi*b）;
????Y1=X1*n+m;
????Y2=X2*n+m;
end
disp（Y1）;
%求平均值和標準差
M1=mean（Y1）;
N1=std（Y1）;
disp（M1）;
disp（N1）;
disp（Y1）;
%求平均值和標準差
M2=mean（Y2）;
N2=std（Y2）;
disp（M2）;
disp（N2）;
%?%將數據寫入文本文件
%?fid=fopen（‘xiefei1.dat‘‘w‘