clear?all
close?all
i=10;
j=5000;
t=linspace（010j）;%???linspace（x0x1n）其中n代表的是點的數目，即分成n-1等分。步長應當是（x1-x0）/（n-1）?
fc=10;
SNR=10;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%產生基帶信號
a=round（rand（1i））;%隨機序列
plot（rand（1i））
st=t;
for?n=1:10
????for?m=j/i*（n-1）+1:j/i*n
??????????if?a（n）<1
??????????????st（m）=0;
??????????else
??????????????st（m）=1;
??????????end
??????end
??end?
figure（1）
subplot（311）;
plot（tst）;
axis（[010-12]）;
title（‘基帶信號‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%載波
s1=cos（2*pi*fc*t）;
subplot（312）;
plot（ts1）;
axis（[0?10?-1.2?1.2]）;
title（‘載波信號‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%調制
e_2ask=st.*s1;
subplot（313）;
plot（te_2ask）;
axis（[?0?10?-2?2]）;
title（‘已調信號‘）;
awgn_2ask=awgn（e_2askSNR‘measured‘‘dB‘）;%%%%%%%%%%%%%%%加入高斯白噪聲
figure（2）
subplot（411）;
plot（tawgn_2ask）;
axis（[?0?10?-2?2]）;
title（‘加入噪聲的信號‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解調
[B1A1]=butter（1[8/（500/2）?12/（500/2）]）;%帶通濾波器
awgn_ask1=filter（B1A1awgn_2ask）;
at=awgn_ask1.*cos（2*pi*fc*t）;
subplot（412）;
plot（tat）;
axis（[?0?10?-2?2]）;
title（‘帶通濾波后與和載波同頻同相的信號相乘‘）;
[BA]=butter（12/（500/2）‘low‘）;
at=filter（BAat）;
subplot（413）;
plot（tat）;
title（‘解調后波形‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%抽樣判決
at1=zeros（1j）;
for?m=1:i
????for?j=（m-1）*500+1:m*500
????????if?at（1（m-1）*500+250）+0.5<0.5
????????????at1（1j）=0;
????????else
????????????at1（1j）=1;
????????end
????end
end
subplot（414）
plot（tat1）;
axis（[010-12]）;
title（‘抽樣判決后波形‘）
%%%%%%%%%%%%%%%%%%%%%%%%%%%?????相干解調后誤碼率統計
for?m=1:10
????for?j=（m-1）*500+1:m*500
????????if?at（1j）+0.5<0.5
????????????at2（1j）=0;
????????else??at2（1j）=1;
????????end
????end
end
error=st-at2;
error=abs（error）;
errorrate=sum（error）/5000
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%??過程中各波形頻譜
figure（3）
subplot（511）
[fstf]=T2F（tst）;
plot（fstf）;
ylabel（‘（dB/Hz）‘）;
axis（[-15?15?-2?6]）;
title（‘輸入信號頻譜‘）
subplot（512）
[fs1f]=T2F（ts1）;
plot（fs1f）;
ylabel（‘（dB/Hz）‘）;
axis（[-15?15?-2?6]）;
title（‘載波信號頻譜‘）
subplot（513）
[fe_2askf]=T2F（te_2ask）;
plot（fe_2askf）;
ylabel（‘（dB/Hz）‘）;
axis（[-15?15?-2?6]）;
title（‘已調信號頻譜‘）
subplot（514）
[fawgn_2askf]=T2F（tawgn_2ask）;
plot（fawgn_2askf）;
ylabel（‘（dB/Hz）‘）;
axis（[-15?15?-2?6]）;
title（‘加入高斯白噪聲后已調信號頻譜‘）
subplot（515）
[fat1f]=T2F（tat1）;
plot（fat1f）;
ylabel（‘（dB/Hz）‘）;
axis（[-15?15?-2?6]）;
title（‘解調之后進行抽樣判決后信號頻譜‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2594??2010-06-10?19:04??ASK.m