clear?all
close?all
clc

%全局變量
%N=2^12?????????????????????????????????%寬度，或稱為點數(shù)
Ts=0.01?????????????????????????????????%時間分辨率（時間間隔），抽樣時間間隔
t=0:Ts:10???????????????????????????????%時間坐標
N=length（t）
fs=1/Ts
df=fs/（N-1）?????????????????????????????%df為頻率分辨率（頻率間隔）
n=-（N-1）/2:（N-1）/2
f=n*df??????????????????????????????????%頻率坐標

%基帶信號
figure
p=[1?0?1?1?0?1?0?0?1?1]????????????????%待傳送的編碼串可用randint（1n）產(chǎn)生信源；若用rand（1nM），則范圍是1~M-1
d1=[0:1:t（end）-1]??????????????????????%每個編碼的發(fā)送延遲時間
p1=[pzeros（length（d1）-length（p））]?????%在發(fā)送串后補零
d=[d1;p1]‘?????????????????????????????%產(chǎn)生d矩陣
m=pulstran（t-0.5d‘rectpuls‘）
subplot（211）
plot（tm）
axis（[t（1）?t（end）?-0.5?1.5]）
grid?on
title（‘基帶信號p‘）
%基帶信號頻譜
M=fft（m）
M=fftshift（M）
subplot（212）
plot（fabs（M）/N）
title（‘基帶信號頻譜‘）

%抽取IQ信號??BPSK只映射一路I。QPSK兩個基帶符號為一組，映射到一對IQ
figure
pQI=reshape（p25）
pQ=pQI（1:）
pI=pQI（2:）
d2=[0:2:t（end）-2]
dQ=[d2;pQ]‘
dI=[d2;pI]‘
mI=pulstran（t-1dI‘rectpuls‘2）
mQ=pulstran（t-1dQ‘rectpuls‘2）
subplot（211）
plot（tmI）
axis（[t（1）?t（end）?-0.5?1.5]）
grid?on
title（‘I路抽取‘）
subplot（212）
plot（tmQ）
axis（[t（1）?t（end）?-0.5?1.5]）
grid?on
title（‘Q路抽取‘）
%IQ值映射???“0”映射成1，“1”映射成-1
figure
I（pI==1）=-1
I（pI==0）=1
Q（pQ==1）=-1
Q（pQ==0）=1
dImap=[d2;I]‘
dQmap=[d2;Q]‘
mImap=pulstran（t-1dImap‘rectpuls‘2）
mQmap=pulstran（t-1dQmap‘rectpuls‘2）
subplot（211）
plot（tmImap）
axis（[t（1）?t（end）?-1.5?1.5]）
grid?on
title（‘映射后的I路信號‘）
subplot（212）
plot（tmQmap）
axis（[t（1）?t（end）?-1.5?1.5]）
grid?on
title（‘映射后的Q路信號‘）

%載波信號
fc=5?????????????????%載波頻率
c1=cos（2*pi*fc*t）
c2=sin（2*pi*fc*t）
%載波信號頻譜
figure
C1=fft（c1）
C1=fftshift（C1）
subplot（211）
plot（fabs（C1）/N）
title（‘載波cos頻譜‘）
C2=fft（c2）
C2=fftshift（C2）
subplot（212）
plot（fabs（C2）/N）
title（‘載波sin頻譜‘）

%已調(diào)信號
figure
s1=mImap.*c1
subplot（311）
plot（ts1）
axis（[t（1）?t（end）?-1.5?1.5]）
title（‘I路調(diào)制cos‘）
grid?on
s2=mQmap.*c2
subplot（312）
plot（ts2）
axis（[t（1）?t（end）?-1.5?1.5]）
title（‘Q路調(diào)制sin‘）
grid?on
s=s1-s2
subplot（313）
plot（ts）
axis（[t（1）?t（end）?-1.5?1.5]）
title（‘進入信道的信號s=s1-s2‘）
grid?on

%已調(diào)信號頻譜
figure
S1=fft（s1）
S1=fftshift（S1）
subplot（311）
plot（fabs（S1）/N）
title（‘s1頻譜‘）
S2=fft（s2）
S2=fftshift（S2）
subplot（312）
plot（fabs（S2）/N）
title（‘s2頻譜‘）
S=fft（s）
S=fftshift（S）
subplot（313）
plot（fabs（S）/N）
title（‘s頻譜‘）

%通過信道后接收到信號
figure
r=awgn（s10）
subplot（211）
plot（tr）
axis（[t（1）?t（end）?-1.5?1.5]）
title（‘接收信號時域‘）
%接收信號的頻譜
R=fft（r）
R=fftshift（R）
subplot（212）
plot（fabs（R）/N）
title（‘接收信號頻域‘）

%解調(diào)
figure
y1=r.*c1
subplot（211）
plot（ty1）
axis（[t（1）?t（end）?-1.5?1.5]）
title（‘I路解調(diào)‘）
y2=r.*（-c2）
subplot（212）
plot（ty2）
axis（[t（1）?t（end）?-1.5?1.5]）
title（‘Q路解調(diào)‘）


%解調(diào)后從時域積分（求和）判決
y1sum=0
y2sum=0
for?i=1:5
????y1sum（i）=sum（y1（200*（i-1）+1:200*i））
????y2sum（i）=sum（y2（200*（i-1）+1:200*i））
end?
y1d=sign（y1sum）
y2d=sign（y2sum）

%反映射
y1rmap（y1d==1）=0
y1rmap（y1d==-1）=1
y2rmap（y2d==1）=0
y2rmap（y2d==-1）=

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

?????文件???????3514??2011-11-02?16:25??IQ\3971002my_QPSKmod_201107\my_QPSKmod_201107.m

?????文件???????5398??2008-03-07?14:44??IQ\40474840lfm_iq\lfm_iq.m

?????文件???????1203??2009-04-18?21:44??IQ\69154523IQ\IQ.m

?????文件???????5811??2009-06-22?14:43??IQ\786579139IQ\IQ.m

?????文件???????1457??2010-04-22?21:28??IQ\877815358IQ\IQ.m

?????目錄??????????0??2012-03-30?10:59??IQ\3971002my_QPSKmod_201107

?????目錄??????????0??2012-03-30?10:59??IQ\40474840lfm_iq

?????目錄??????????0??2012-03-30?10:59??IQ\69154523IQ

?????目錄??????????0??2012-03-30?10:59??IQ\786579139IQ

?????目錄??????????0??2012-03-30?11:00??IQ\877815358IQ

?????目錄??????????0??2012-03-30?11:00??IQ

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

????????????????17383????????????????????11