boxmuller方法生成正态分布随机数

%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‘