clear?all
close?all
clc

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

%基帶信號(hào)
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]?????????;?????????????%每個(gè)編碼的發(fā)送延遲時(shí)間
p1=[pzeros（length（d1）-length（p））]?;????%在發(fā)送串后補(bǔ)零
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（‘基帶信號(hào)p‘）
%基帶信號(hào)頻譜
M=fft（m）;
M=fftshift（M）;
subplot（212）
plot（fabs（M）/N）
title（‘基帶信號(hào)頻譜‘）

%抽取IQ信號(hào)??BPSK只映射一路I。QPSK兩個(gè)基帶符號(hào)為一組，映射到一對(duì)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路信號(hào)‘）
subplot（212）
plot（tmQmap）
axis（[t（1）?t（end）?-1.5?1.5]）
grid?on
title（‘映射后的Q路信號(hào)‘）

%載波信號(hào)
fc=5?????????????????%載波頻率;
c1=cos（2*pi*fc*t）;
c2=sin（2*pi*fc*t）;
%載波信號(hào)頻譜
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)信號(hà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（‘進(jìn)入信道的信號(hào)s=s1-s2‘）
grid?on

%已調(diào)信號(hà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頻譜‘）

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

%解調(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)后從時(shí)域積分（求和）判決
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;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????3557??2014-11-10?10:47??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

?????文件????????355??2014-11-10?10:55??IQ\786579139IQ\Untitled2.m

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

?????目錄??????????0??2014-11-10?10:45??IQ\3971002my_QPSKmod_201107

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

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

?????目錄??????????0??2014-11-10?10:55??IQ\786579139IQ

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

?????目錄??????????0??2014-11-18?16:43??IQ

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

????????????????17781????????????????????12