my_alphanoise.rar

用matlab实现alpha稳定分布序列产生，并进行仿真，函数参照论文-Alpha稳定分布的模型仿真及参数估计

%% 标准参数系数下随机变量的生成方法
%% 此函数参考：毕业论文-华中科技大学-Alpha稳定分布的模型仿真及参数估计
%X～S（αβσμ）
%标准参数下服从S（αβ10）分布的随机变量生成
function x=my_alphanoise（alphabetasigmamumsns）
%α=alphaβ=betaσ=sigmaμ=mu
%pi=3.14;
%% alpha!=1：
M1=（1+beta^2*tan（pi*alpha/2）^2）^（alpha/2）;
N1=-atan（beta*tan（pi*alpha/2））/alpha;
p1=（1-alpha）/alpha;
rand（‘state‘sum（100*clock））;
V=pi*（rand（msns））-pi/2;% 产生（-pi/2pi/2）的均匀分布随机变量
rand（‘state‘7*sum（100*clock）+3）;
W=-log（1-rand（msns））;% 产生均值为1的指数分布随机变量
X1=M1*sin（alpha*（V-N1））./（cos（V））.^（1/alpha）.*（cos（V-alpha*（V-N1））./W）.^p1;% 产生nonnegative innovations η（k）
%% alpha==1：
M2=pi/2;
X2=M2*（（pi/2+beta.*V）.*tan（V）-beta*log（W.*cos（V）./（pi/2+beta.*V）））;
%% 产生标准参数系下服从S（alphabetasigmamu）
if（alpha==1）
    x=sigma*X2+2/pi*beta*sigma*log（sigma）+mu;
else
    x=sigma*X1+mu;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        919  2019-05-24 19:44  my_alphanoise\my_alphanoise.m

     文件        206  2019-05-24 19:57  my_alphanoise\my_noise.m

     目录          0  2019-05-24 20:08  my_alphanoise

----------- ---------  ---------- -----  ----

                 1125                    3