function?y=alphacx（alphabetadeltamuN）??%alpha?特征指數；beta?對稱參數；delta?分散系數；mu位置參數?;N輸入參數數量
y=[];
for?k=1:N
v=pi*rand-pi/2;%（-pi/2，pi/2）的均勻分布
w=exprnd（1）;???%平均參數為1的指數分布
if?alpha~=1????
if?alpha>1
????beta_=2*atan（beta*tan（pi*alpha/2））/pi/（alpha-2）;????
????v0=-（pi/2）*beta_*（alpha-2）/alpha;
????x=sin（alpha*（v-v0））/（cos（v））.^（1/alpha）*（cos（v-alpha*（v-v0））/w）.^（（1-