Fc=10;???%載頻??????
Fs=40;???%系統采樣頻率??????
Fd=1;????%碼速率?????
N=Fs/Fd;????????
df=10;??????
numSymb=25;%進行仿真的信息代碼個數?????
M=2;???????%進制數??????
SNRpBit=60;%信噪比?????
SNR=SNRpBit/log2（M）;%60???
seed=[12345?54321];?????
numPlot=15;????
x=randsrc（numSymb1[0:M-1]）;%產生25個二進制隨機碼?????
figure（1）?????
stem（[0:numPlot-1]x（1:numPlot）‘bx‘）;%顯示15個碼元，桿圖，從x的前十五個隨機數中選取?????
title（‘二進制隨機序列?桿圖‘）??????
xlabel（‘T‘）;?????
ylabel（‘幅度‘）;?

%調制
?y=dmod（xFcFdFs‘fsk‘Mdf）;%數字帶通調制?
?%x是一個隨機序列，如果M=2則x=[1?1?0?0]之類的，Fd為碼元速率，Fc為載頻，Fs為采樣頻率，fsk表示為頻移鍵控調制。?????
?numModPlot=numPlot*Fs;?%15*40?????
?t=[0:numModPlot-1]./Fs;%數組除法（仿真時間）?????
?figure（2）??????
?plot（ty（1:length（t））‘b-‘）;?????
?axis（[min（t）?max（t）?-1.5?1.5]）;??????
?title（‘調制后的信號‘）????
?xlabel（‘Time‘）;?????
?yl