????????????????????????%IN?THE?NAME?OF?GOD%
????????%?Computing?of?Moments?and?Cumulants?for?QAM16%

clear?all;
close?all;
M=16;

[a?b]=qaskenco（M）;%qaskdeco?從方形的QASK星座圖反映射到數字信號；qaskenco?把數字信號映射到方形的QASK星座圖
??????????????????%This?is?an?obsolete?function?and?may?be?removed?in?the?future.Please?use?QAMMOD?instead.
aa=a‘;
bb=b‘;
stem（ab‘filled‘）
grid;
axis（[-4?4?-4?4]）;
xlabel（‘Inphase‘）;
ylabel（‘Quadrature‘）;


a2=a.^2;b2=b.^2;?a4=a.^4;b4=b.^4;
a6=a.^6;?b6=b.^6;a8=a.^8;?b8=b.^8;

P=（sum（a2+b2））/M;

m20=（sum（a2-b2））/M;
m21=（sum（a2+b2））/M;



m40=（（sum（a4+b4-6*（a2.*b2））））/M;
m41=（sum（a4-b4））/M;
m42=sum（a4+b4+2*（a2.*b2））/M;



m60=sum（a6-b6+15*（a2.*b4）-15*（a4.*b2））/M;
m61=sum（a6+b6-5*（a2.*b4）-5*（a4.*b2））/M;
m62=sum（a6-b6-（a2.*b4）+（a4.*b2））/M;
m63=sum（a6+b6+3*（a2.*b4）+3*（a4.*b2））/M;



m80=（sum（a8+b8-28*（a6.*b2）+70*（a4.*b4）-28*（a2.*b6）））/M;
m81=（sum（a8-b8-14*（a6.*b2）+14*（a2.*b6）））/M;
m82=（sum（a8+b8-4*（a6.*b2）-10*（a4.*b4）-4*（a2.*b6）））/M;
m83=（sum（a8-b8+2*（a6.*b2）-2*（a2.*b6）））/M;
m84=（sum（a8+b8+4*（a6.*b2）+6*（a4.*b4）+4*（a2.*b6）））/M;

M20=m20/P;
M21=m21/P;

M40=m40/P^2;
M41=m41/P^2;
M42=m42/P^2;

M60=m60/P^3;
M61=m61/P^3;
M62=m62/P^3;
M63=m63/P^3;

M80=m80/P^4;
M81=m81/P^4;
M82=m82/P^4;
M83=m83/P^4;
M84=m84/P^4;


Cum20=M20;
Cum21=M21;

Cum40=（m40-3*m20^2）/P^2;
Cum41=（m41-3*m20*m21）/P^2;
Cum42=（m42-m20^2-2*m21^2）/P^2;

Cum60=（m60-15*m20*m40+30*（m20）^2）/P^3;
Cum61=（m61-10*m20*m41-5*m21*m40+30*（m20）^2*m21）/P^3;
Cum62=（m62-m20*m40-8*m21*m41-6*（m20）*m42+6*m20^3+24*（m21^2）*m20）/P^3;
Cum63=（m63-6*m20*m41-9*m21*m42+18*（m20）^2*m21+12*（m21^3））/P^3;

Cum80=（m80-35*（m40^2）-630*（m20^4）+420*（m20^2）*m40）/P^4;

Cum81=（m81-35*m40*m41-630*（m20^3）*m21+210*m40*m20*m21+210*m20*m41）/P^4;

Cum82=（m82-15*m40*m42-20*（m41^2）+30*m40*（m20^2）+60*m40*（m21^2）+240*m41*m21*m20+90*m42*（m20^2）-90*（m20^4）-540*（m20^2）*（m21^2））/P^4;

Cum83=（m83-5*m40*m41-30*m41*m42+90*m41*（m20^2）+120*m41*（m21^2）+180*m42*m21*m20+30*m40*m20*m21-270*（m20^3）*m21-360*（m21^3）*m20）/P^4;?

Cum84=（m84-m40^2-18*（m42^2）-16*（m41^2）-54*（m20^4）-144*（m21^4）-432*（m20^2）*（m21^2）+12*m40*（m20^2）+96*m41*m21*m20+144*m42*（m21^2）+72*m42*（m20^2）+96*m41*m20*m21）/P^4;?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1552??2012-06-14?09:59??a_mask.m

?????文件???????1140??2012-06-14?15:40??a_MFSK.m

?????文件???????1211??2012-06-14?10:56??a_mpsk.m

?????文件???????1189??2012-06-15?20:47??a_MSK.m

?????文件????????863??2012-06-15?21:10??a_2ask.m

?????文件????????851??2012-06-15?21:46??a_2FSK.m

?????文件???????2238??2012-06-14?16:13??a_16QAM.m

-----------?---------??----------?-----??----

?????????????????9044????????????????????7