随机相位正弦波仿真实验MATLAB程序

clearclc
%初始化
fs=100;                       %采样频率100
t=0:1/fs:pi;                  %设置时间范围0~pi
B=0;                          %设置初始相位为0
A=1;                          %信号的幅度为1
f=10;                         %信号频率为10Hz

%画随机相位正弦波
B=2*pi*rand（1length（t））;     %求变化范围是0~2pi的随机相位
x=A.*sin（2*pi*f.*t+B）;
e1=mean（x）                    %计算均值
s1=var（x）                     %计算方差 
[rlags]=xcorr（x）;            %计算自相关函数   
subplot（411）
plot（tx）
title（‘随机相位正弦波‘）;xlabel（‘t‘）;ylabel（‘幅值x‘）;
subplot（412）
plot（te1）
title（‘均值‘）;xlabel（‘t‘）;ylabel（‘e‘）;
subplot（413）
plot（ts1）axis（[03.14-0.20.2]）
title（‘方差‘）;xlabel（‘t‘）;ylabel（‘s‘）;
subplot（414）
plot（lags/fsr）axis（[03.1401]）
title（‘自相关函数‘）;xlabel（‘时间间隔‘）;ylabel（‘r‘）;
%幅度，相位，频率还原初始化情况
A=1;                          %信号的幅度为1
B=0;                          %设置初始相位为0
f=10;                         %信号频率为10Hz

%画幅度随机正弦波
A=raylrnd（1:length（t））;       %产生满足瑞利分布的随机幅度
x=A.*sin（2*pi*f*t+B）; 
e2=mean（x）                    %计算均值
s2=var（x）                     %计算方差 
[rlags]=xcorr（x）;            %计算自相关函数   
figure
subplot（411）
plot（tx）
title（‘随机幅度正弦波‘）;xlabel（‘t‘）;ylabel（‘幅值x‘）;
subplot（412）
plot（te2）
title（‘均值‘）;xlabel（‘t‘）;ylabel（‘e‘）;
subplot（413）
plot（ts2）
title（‘方差‘）;xlabel（‘t‘）;ylabel（‘s‘）;
subplot（414）
plot（lags/fsr）
title（‘自相关函数‘）;xlabel（‘时间间隔‘）;ylabel（‘r‘）;
%幅度，相位，频率还原初始化情况
A=1;                          %信号的幅度为1
B=0;                          %设置初始相位为0
f=10;                         %信号频率为10Hz

%画随机频率正弦波
f=20*rand（1length（t））;       %产生均匀分布的随机频率，以10Hz为中心，上下宽度10Hz
x=A.*sin（2*pi*f.*t+B）;
e3=mean（x）                    %计算均值
s3=var（x）                     %计算方差 
[rlags]=xcorr（x）;            %计算自相关函数   
figure
subplot（411）
plot（tx）
title（‘随机频率正弦波‘）;xlabel（‘t‘）;ylabel（‘幅值x‘）;
subplot（412）
plot（te3）
title（‘均值‘）;xlabel（‘t‘）;ylabel（‘e‘）;
subplot（413）
plot（ts3）
title（‘方差‘）;xlabel（‘t‘）;ylabel（‘s‘）;
subplot（414）
plot（lags/fsr）
title（‘自相关函数‘）;xlabel（‘时间间隔‘）;ylabel（‘r‘）;
%幅度，相位，频率还原初始化情况
A=1;                          %信号的幅度为1
B=0;                          %设置初始相位为0
f=10;                         %信号频率为10Hz

%画随机相位频率正弦波
B=2*pi*rand（1length（t））;     %产生0~2pi均匀分布的随机相位
f=20*rand（1length（t））;       %产生均匀分布的随机频率，以10Hz为中心，上下宽度10Hz
x=A*sin（2*pi*f.*t+B）;
e4=mean（x）                    %计算均值
s4=var（x）                     %计算方差 
[rlags]=xcorr（x）;            %计算自相