資源簡介
用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
- 上一篇:DFT的matlab實現的 源程序
- 下一篇:matlab 數字信號處理函數
評論
共有 條評論