clear?all
nsymbol=1000;?????????%每種信噪比下的發送符號數

T=1;????????????????????????%符號周期
fs=100;?????????????????????%每個符號的采樣點數
ts=1/fs;????????????????????%采樣時間間隔
t=0:ts:T-ts;????????????????%時間向量
fc=10;??
Fd=1;
Fc=0;
Fs=Fd*8;%載波頻率



c=sqrt（2/T）*exp（j*2*pi*fc*t）;???%載波信號
c1=sqrt（2/T）*cos（2*pi*fc*t）;????%同相載波
c2=-sqrt（2/T）*sin（2*pi*fc*t）;???%正交載波?????

M=4;??%8-PSK
%graycode=[0?1?2?3?6?7?4?5];?????????????%Gray編碼規則??????????????????????
EsN0=0:10;??????????????????????????????%信噪比，Es/N0
snr1=10.^（EsN0/10）;?????????????????????%信噪比轉換為線性值


msg=randint（1nsymbolM）;???????????????%消息數據
%msg1=graycode（msg+1）;???????????????????%Gray映射
msgmod1=pskmod（msgM）.‘;????????????????%基帶Q-PSK調制


%?經過濾波器???成型濾波器
rolloff=0.2;
delay=3;
msgmod=rcosflt（msgmod1FdFs‘fir/sqrt‘rolloffdelay）;

msgmod=msgmod（（delay*Fs/Fd+1:（end?-?delay*Fs/Fd））1）;

scatterplot（msgmod）

?%載波調制
tx=real（msgmod*c）;?????????????????????
tx1=reshape（tx.‘1length（msgmod）*length（c））;???

spow=norm（tx1）.^2/nsymbol;??????????????%求每個符號的平均功率

for?indx=1:length（EsN0）
????sigma=sqrt（spow/（2*snr1（indx）））;????????????????%根據符號功率求噪聲功率
????rx=tx1+sigma*randn（1length（tx1））;??????????????%加入高斯白噪聲
????rx1=reshape（rxlength（c）length（msgmod））;?????
????
????r1=（c1*rx1）/length（c1）;?????????????????????????%相關運算
????r2=（c2*rx1）/length（c2）;
????r=r1+j*r2;
????eyediagram（r2）;
????scatterplot（r）;
????
????
?????%接受端??信號經過低通濾波器??匹配濾波器
?????rolloff=0.2;?????
?????%delay=3;
?????k=rcosflt（rFdFs‘fir/sqrt/Fs‘rolloff）;
?????%?g=?downsample（k8）;
?????k=k（（delay*Fs/Fd+1:（end?-?delay*Fs/Fd））1）;
%?????h?=?resample（k?1000?8048）;

????y=pskdemod（kM）;????????%PSK解調;
????%h?=?resample（y?1000?8000）;
?????g=?downsample（y8）;
????
????u=g‘;
?
????%[errber（indx）]=biterr（msgylog2（M））;?????%誤比特率
????[errser（indx）]=symerr（msgu）;?????????????%誤符號率
end
figure（）;
ser1=2*qfunc（sqrt（2*snr1）*sin（pi/M））;???????????????%理論誤符號率
%ber1=1/log2（M）*ser1;????????????????????????????????%理論誤比特率
semilogy（EsN0ser‘-k*‘EsN0ser1‘-k.‘）;
title（‘Q-PSK載波調制信號在AWGN信道下的性能‘）
xlabel（‘Es/N0‘）;ylabel（‘誤比特率‘）
legend（‘誤比特率‘‘理論誤比特率‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????43475??2014-04-11?00:09??QQ截圖20140411000906.png
?????文件????????2341??2014-04-10?23:45??qpsk過高斯信道.m
?????文件???????14579??2014-04-11?00:07??QQ截圖20140411000753.png