close?all
clc
clear

%設置仿真參數
global?dt?t?f?df?N?T
N=2^16+2^9?;
L=64;
M=N/L;
Rb=2;
Ts=1/Rb;
dt=Ts/L;
df=1/（N*dt）;
T=N*dt;
Bs=N*df/2;
Na=4;
alpha=0.8;
t=linspace（-T/2T/2N）;
f=linspace（-BsBsN）+eps;
P_s0=zeros（1N）;
P_s1=zeros（1N）;
P_s2=zeros（1N）;
P_s3=zeros（1N）;
Eb_N0=zeros（110）;
eb_n0=zeros（110）;
Pe=zeros（110）;


%升余弦滾降函數
hr1=sin（pi*t/Ts）./（pi*t/Ts）;
hr2=cos（alpha*pi*t/Ts）./（1-（2*alpha*t/Ts）.^2）;
hr=hr1.*hr2;
HR=abs（t2f（hr））;
GT=sqrt（HR）;
GR=GT;
%余弦載波信號
A=1.7;
fc=20;
ca=A*cos（2*pi*fc*t）;
sa=A*sin（2*pi*fc*t）;


%QAM信號通過頻帶系統
for?lp1=1:10
????Eb_N0（lp1）=（lp1-1）;
????eb_n0（lp1）=10^（Eb_N0（lp1）/10）;
????Eb=1;n0=Eb/eb_n0（lp1）;?
????n0=sqrt（n0*Bs）;
????n_err=0;?????				
????for?lp2=1:10
????????a=round（rand（1M））;%輸入的碼元信號0、1序列
????????as=a（1:2:M）;%串并變換
????????ac=a（2:2:M）;
????????as=Graycode2（as）;
????????ac=Graycode2（ac）;%電平變化
????????imps=zeros（1N）;??
????????impc=zeros（1N）;
????????imps（3*L:6*L:N）=as/dt;
????????impc（3*L:6*L:N）=ac/dt;
????????IMPS=t2f（imps）;
????????IMPC=t2f（impc）;????????
????????bfs=IMPS.*GT;?
????????bfc=IMPC.*GT;????????
????????bs=real（f2t（bfs））;
????????bc=real（f2t（bfc））;
????????sts=bs.*sa;
????????stc=bc.*ca;
????????sc=stc-sts;?%QAM???????
????????
????????Ps0=bfc.*conj（bfc）/T;
????????P_s0=P_s0+Ps0;
????????
????????S1=t2f（stc）;
????????Ps1=S1.*conj（S1）/T;
????????P_s1=P_s1+Ps1;
????????
????????S2=t2f（sc）;
????????Ps2=S2.*conj（S2）/T;
????????P_s2=P_s2+Ps2;
????????
????????nr=n0*randn（1N）;???????????
????????sr=sc+nr;
????????src=sr.*ca;
????????srs=-sr.*sa;
????????Yc=t2f（src）.*GR;????????
????????Ys=t2f（srs）.*GR;
????????yc=real（f2t（Yc））/A^2;
????????ys=real（f2t（Ys））/A^2;
????????
????????Ps3=Yc.*conj（Yc）/T;
????????P_s3=P_s3+Ps3;
????????
????????yc1=yc（3*L:6*L:N）;
????????ys1=ys（3*L:6*L:N）;
????????yoc=Graydeco2（yc1）;
????????yos=Graydeco2（ys1）;%解碼過程
????????yo=reshape（[yos;yoc]1M）;%回復原始序列???????????????????????????????
????????n_err=n_err+length（find（yo~=a））;???????
????end
????Pe（lp1）=n_err/（M*lp2）;
end
%同相支路調制信號波形與平均功率譜
figure（1）;
subplot（211）;
plot（tbc）;
title（‘同相支路調制信號波形‘）;
axis（[-8?8?-5?5]）;
subplot（212）;
P_s0=P_s0/（lp1*lp2）;
plot（f30+10*log10（P_s0））;
title（‘同相支路調制信號功率譜‘）;
axis（[-25?25?-50?30]）;
xlabel（‘f?（MHz）‘）;
ylabel（‘Ps（f）?（dBm/MHz）‘）;


%同相支路已調信號波形與平均功率譜
figure（2）;
subplot（211）;
plot（tstc）;
title（‘同相支路已調信號波形‘）;
axis（[-8?8?-5?5]）;
subplot（212）;
P_s1=P_s1/（lp1*lp2）;
plot（f30+10*log10（P_s1））;
title（‘同相支路已調信號功率譜‘）;
axis（[-40?40?-50?30]）;
xlabel（‘f?（MHz）‘）;
ylabel（‘Ps（f）?（dBm/MHz）‘）;

%16QAM信號波形與平均功率譜
figure（3）;
subplot（211）;
plot（tsc）;
title（‘64QAM信號波形‘）;
axis（[-8?8?-5?5]）;
subplot（212）;
P_s2=P_s2/（lp1*lp2）;
plot（f30+10*log10（P_s2））;
title（‘64QAM信號功率譜‘）;
axis（[-40?40?-50?30]）;
xlabel（‘f?（MHz）‘）;
ylabel（‘Ps（f）?（dBm/MHz）‘）;

%同相支路匹配濾波器解調信號波形與平均功率譜
figure（4）;
subplot（211）;
plot（tyc）;
title（‘同相支路解調信號波形‘）;
axis（[-8?8?-2.5?2.5]）;
subplo

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

?????文件???????3929??2012-05-04?22:50??運行正確的代碼bpsk，qpsk?qam2?-?自己做\64qam.m

?????文件???????2735??2009-05-25?21:02??運行正確的代碼bpsk，qpsk?qam2?-?自己做\basic.m

?????文件???????2941??2012-05-03?13:56??運行正確的代碼bpsk，qpsk?qam2?-?自己做\BPSK.asv

?????文件???????2945??2012-05-03?13:59??運行正確的代碼bpsk，qpsk?qam2?-?自己做\BPSK.m

?????文件?????????76??2009-05-04?20:20??運行正確的代碼bpsk，qpsk?qam2?-?自己做\f2t.m

?????文件????????219??2009-05-17?17:18??運行正確的代碼bpsk，qpsk?qam2?-?自己做\Graycode.m

?????文件????????270??2012-05-04?22:23??運行正確的代碼bpsk，qpsk?qam2?-?自己做\Graycode2.asv

?????文件????????241??2012-05-04?22:25??運行正確的代碼bpsk，qpsk?qam2?-?自己做\Graycode2.m

?????文件????????382??2012-05-04?22:23??運行正確的代碼bpsk，qpsk?qam2?-?自己做\Graydeco.asv

?????文件????????381??2009-05-17?16:32??運行正確的代碼bpsk，qpsk?qam2?-?自己做\Graydeco.m

?????文件????????342??2012-05-04?22:18??運行正確的代碼bpsk，qpsk?qam2?-?自己做\Graydeco2.asv

?????文件????????243??2012-05-04?22:24??運行正確的代碼bpsk，qpsk?qam2?-?自己做\Graydeco2.m

?????文件????????121??2009-05-18?16:42??運行正確的代碼bpsk，qpsk?qam2?-?自己做\lpf.m

?????文件???????3914??2012-05-03?14:44??運行正確的代碼bpsk，qpsk?qam2?-?自己做\QAM.asv

?????文件???????3924??2012-04-20?16:51??運行正確的代碼bpsk，qpsk?qam2?-?自己做\QAM.m

?????文件???????3601??2012-05-03?14:06??運行正確的代碼bpsk，qpsk?qam2?-?自己做\QPSK.asv

?????文件???????3630??2012-05-03?14:10??運行正確的代碼bpsk，qpsk?qam2?-?自己做\QPSK.m

?????文件???????1288??2009-05-25?23:10??運行正確的代碼bpsk，qpsk?qam2?-?自己做\SSB.asv

?????文件???????1288??2009-05-25?22:28??運行正確的代碼bpsk，qpsk?qam2?-?自己做\SSB.m

?????文件?????????77??2012-04-20?16:34??運行正確的代碼bpsk，qpsk?qam2?-?自己做\t2f.asv

?????文件?????????75??2012-04-20?17:33??運行正確的代碼bpsk，qpsk?qam2?-?自己做\t2f.m

?????目錄??????????0??2012-05-04?22:44??運行正確的代碼bpsk，qpsk?qam2?-?自己做

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

????????????????32622????????????????????22