clc;
clear;
%%%%%%%%%%%%%%?調(diào)制?%%%%%%%%%%%%%%%%%%%

M?=?4;%調(diào)制元數(shù)
n?=?1000;%符號個數(shù)
Is?=?randi（M1n）;%輸入的碼元
d?=?1;
width?=?1/1000;%碼元寬度
fs?=?10/width?*（1+M）;
ts?=?1/fs;?%采樣時間間隔
t?=?0:?ts?:?n*width-ts;%符號寬度
T?=?0?:?ts?:?width-ts;
Send?=?zeros（1length（t））;
lenT?=?length（T）;
for?i?=?1:1:n
????Send（（i-1）*lenT+1?:?i*lenT）?=?cos（pi/width*（1+Is（i））*T）;
end
figure（1）
subplot（211）
plot（tSend）
title（‘已調(diào)信號時域波形‘）;
xlabel（‘時間t/s‘）;
ylabel（‘已調(diào)信號幅度‘）;
%?已調(diào)信號頻域波形
Nf=1024;
%YSend=abs（fft（SendNf））;
YSend=20*log10（abs（fft（SendNf）））;
f=0:Nf/2;
subplot（212）
plot（f/Nf*fsYSend（f+1））
axis（[0?20000?0?150]）;
title（‘已調(diào)信號頻域波形‘）;
xlabel（‘頻率f/Hz‘）;
ylabel（‘已調(diào)信號幅度‘）;

%%%%%%%%%%%%%%?信道?%%%%%%%%%%%%%%%%%%%
SNR?=?（0:1:20）‘;
Es?=?1/2;
Esdb?=?10*log10（Es）*ones（length（SNR）1）;
Endb?=?Esdb?-?SNR;
En?=?10.^（Endb/10）;
awgn?=repmat（sqrt（26*En/2）1length（t））?.*?repmat（（randn（1length（t）））length（SNR）1）;
%%%%%%%%%%%%%%?收端?%%%%%%%%%%%%%%%%%%%
rt?=?zeros（length（SNR）length（t））;
for?i?=?1?:?1?:?length（SNR）
????rt（i:）?=?Send?+?awgn（i:）;
end
rt1?=?rt（10:）;
figure（2）
subplot（211）
plot（trt1）
title（‘接收信號時域波形‘）;
xlabel（‘時間t/s‘）;
ylabel（‘接收信號幅度‘）;
%接收信號頻域波形
Nf=1024;
Yr=20*log10（abs（fft（rt1Nf）））;
%Yr=abs（fft（rt1Nf））;
f=0:Nf/2;
subplot（212）
plot（f/Nf*fsYr（f+1））
title（‘接收信號頻域波形‘）;
xlabel（‘頻率f/Hz‘）;
ylabel（‘接收信號幅度‘）;
%%%%%%%%%%%%%%?解調(diào)?%%%%%%%%%%%%%%%%%%%
Sm?=?zeros（Mlength（T））;%M個符號
for?m?=?1?:?M?
????Sm（m:）?=?cos（pi/width*（1+m）*T）;%符號波形
end
C?=?zeros（length（SNR）M）;%判決準則
Ir?=?zeros（length（SNR）n）;%解調(diào)符號
for?j?=?1?:?1?:?length（SNR）
for?i?=?1:1:n
????for?m?=?1?:?M
????????C（jm）?=?rt（jwidth/ts*（i-1）+1:1:width/ts*i）?*?Sm（m:）‘?*?ts?-?0.5?*?Sm（m:）?*?Sm（m:）‘?*?ts;?%判決準則
????end
????[CmaxIr（ji）]?=?max（C（j:））;
????C?=?zeros（M1）;
end
end

%誤比特率
ne?=?zeros（1length（SNR））;%錯誤的碼元數(shù)
PM?=?zeros（1length（SNR））;%誤碼率
Pb?=?zeros（1length（SNR））;%誤比特率
for?j?=?1?:?1?:?length（SNR）
??for?i?=?1:1:n
???????if（Ir（ji）?~=?Is（i））
????????????ne（j）?=?ne（j）+1;
????????end
???end
????PM（j）?=?ne（j）/n;
????Pb（j）?=?PM（j）/2;
end


figure（3）
semilogy（SNR‘Pb‘-o‘）
title（‘誤比特率曲線‘）;
xlabel（‘SNR（dB）‘）;
ylabel（‘誤比特率‘）;
grid?on
saveas（gcf‘MFSK.fig‘）

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

?????文件???????2381??2014-03-22?23:27??4FSK.m

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

?????????????????2381????????????????????1