clear?all;
close?all;
n=8;%隨機(jī)碼元個(gè)數(shù)
N=100;%模擬一個(gè)碼元的點(diǎn)數(shù)
K=3;%figure1畫三個(gè)小圖
a=round（rand（1n））;%碼元生成
bita=[];%定義空數(shù)組，存放基帶信號(hào)
for?i=1:length（a）
????if?a（i）==0
????????bit1=zeros（1N）;
????else?
????????bit1=ones（1N）;
????end
end
figure（1）;
bitRate=1000;%每一個(gè)碼元中采樣點(diǎn)的間隔寬度0.001s
fc=1000;%載頻1KHZ
t=linspace（01/bitRateN）;
sign=[];
????c1=sin（2*pi*t*fc）;%載波
????c2=sin（2*pi*t*fc*2）;%載波
????for?i=1:length（a）
???????if?a（i）==1
??????????sign=[signc1];
???????else
?????????sign=[signc2];
???????end
????end
subplot（K11）;plot（sign‘LineWidth‘1.5）;title（‘2FSK已調(diào)信號(hào)‘）;grid?on;
signal=awgn（sign20‘measured‘）;%加噪
subplot（K12）;plot（signal‘LineWidth‘1.5）grid?on;title（‘信號(hào)混疊噪聲‘）;
pinpu=spectrum（sign）;
subplot（K13）;plot（pinpu‘LineWidth‘1.5）grid?on;title（‘頻譜圖‘）;
h=waitbar（0‘請(qǐng)等待‘）;%創(chuàng)建進(jìn)度條
rs=1000;%時(shí)間軸頻率步進(jìn)
fc=100;%載波頻率100HZ
tzd=100;
j=0;%進(jìn)度條計(jì)數(shù)
N=10;%設(shè)置循環(huán)次數(shù)
t=0:1/rs:（tzd-1/rs）;
for?snrb=0:1:10?%不同信噪比
????ratio=0;%初始誤碼數(shù)設(shè)為0，累計(jì)十次得到總誤碼數(shù)
????for?k=1:N?%N次循環(huán)每次循環(huán)產(chǎn)生1000個(gè)碼元
????????n=1000;%一次產(chǎn)生碼元數(shù)
????????g=round（rand（1n））;%產(chǎn)生1000個(gè)碼元
????????tz1=g（ceil（10*t+（1/rs）））.*cos（2*pi*fc*t）;%載]波頻率為100HZ
????????tz2=~g（ceil（10*t+（1/rs）））.*cos（2*pi*2*fc*t）;%載波頻率為200HZ
????????tz=tz1+tz2;%得到調(diào)制信號(hào)tz，100個(gè)點(diǎn)表示1個(gè)碼元
????????waitbar（j/（10*N））;
????????signal=awgn（tzsnrb）;
????????%======解調(diào)過程==========%
??