clear?all
close?all
i=10;%10個碼元
j=5000;
t=linspace（05j）;%0-5之間產生5000個點行矢量，即分成5000份
fc=10;%載波頻率
fm=i/5;%碼元速率
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產生基帶信號
x=（rand（1i））%rand函數產生在0-1之間隨機數，共1-10個
%figure（2）
%plot（x）
a=round（x）;%隨機序列，round取最接近小數的整數
%figure（3）
%stem（a）%火柴梗狀圖
st=t;
for?n=1:10
????if?a（n）<1;
????????for?m=j/i*（n-1）+1:j/i*n
????????????st（m）=0;
????????end
????else
????????for?m=j/i*（n-1）+1:j/i*n
????????????st（m）=1;
????????end
????end
end
figure（1）;
subplot（311）;
plot（tst）;
axis（[05-12]）;
title（‘基帶信號st‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%載波
s1=cos（2*pi*fc*t）;
subplot（312）;
plot（s1）;
title（‘載波信號s1‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調制
e_2ask=st.*s1;
subplot（313）;
plot（te_2ask）;
title（‘已調信號‘）;
noise?=rand（1j）;
e_2ask=e_2ask+noise;%加入噪聲
figure（2）
subplot（411）;
plot（te_2ask）;
title（‘加入噪聲的信號‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解調
at=e_2ask.*cos（2*pi*fc*t）;
at=at-mean（at）;%因為是單極性波形，還有直流分量，應去掉
subplot（412）;
plot（tat）;
title（‘與載波相乘后信號‘）;
[faf]?=?T2F（tat）;%通過低通濾波器
[tat]?=?lpf（faf2*fm）;
subplot（413）;
plot（tat）;
title（‘相干解調后波形‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽樣判決
for?m=0:i-1;
????if?at（1m*500+250）+0.5<0.5;
???????for?j=m*500+1:（m+1）*500;
???????????at（1j）=0;
???????end
????else
???????for?j=m*500+1:（m+1）*500;
???????????at（1j）=1;
???????end
????end
end
subplot（414）;
plot（tat）;
axis（[05-12]）;
title（‘抽樣判決后波形‘）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%頻譜觀察
sf=fftshift（abs（fft（e_2ask）））;
figure（3）;
plot（tsf）
title（‘已調信號頻譜‘）;

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

?????文件???????2954??2016-06-17?16:18??two_DPSK.m

?????文件???????2088??2016-06-17?15:54??two_FSK.m

?????文件???????1934??2016-06-17?15:10??two_PSK.m

?????文件???????1732??2016-06-17?15:00??two_ASK.m

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

?????????????????8708????????????????????4