%%%%功能：BPSk信號的調制和解調
function??bpsk（Nfdfcfc_sample）
%?clc;
%?clear?all;

%*******************??參數設置??*************************
%?N=200;
%?fd=32000;%對輸入信號的符號速率
%?fc=10.24e6;%載波頻率
%?fc_sample=8;%每個載波的采樣點

%%以上參數可以更改來產生不同要求的信號

Tb=1/fd;
fs=fc*fc_sample;%采樣速率
M=fs/fd;%過采樣率
t=0:1/fs:N*Tb-1/fs;
%%%輸入需要產生的碼元數目產生雙極性比特序列
rand_data=rand（1N）;
for?i=1:N
????if?rand_data（i）>=0.5
????????rand_data（i）=1;
????else
????????rand_data（i）=-1;
????end
end

figure（1）;
subplot（411）;
stem（rand_data）;
title（‘原始信號輸出‘）;


%%%根據送入的雙極性比特序列，產生成型信號
%?fd=20.48e3;%對輸入信號的符號速率
%?fs=1024e3;%對輸出信號的采樣頻率
%?M=fs/fd;%過采樣率
R=0.35;%升余弦滾降系數α
delay=3;
A=1;
xt=rcosflt（rand_datafdfs‘fir/sqrt‘Rdelay）;%根升余弦濾波器，長度為Fs/Fd?*?（length（X）?+?2?*?DELAY）
%將xt函數移位，消除delay影響
xt=（xt（（delay*M+1）:（end-delay*M）））‘;
subplot（412）;
plot（t/Tbxt）;
title（‘成型濾波之后的基帶信號‘）;


%%%bpsk信號調制
%?fc=20.48e3;
L=length（xt）;%測量xt序列長度
%?t=0:1/fs:L/fs-1/fs;
St=A*xt.*cos（2*pi*fc*t）;
subplot（413）;
plot（t

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-11-18?21:12??調制方式總結\
?????目錄???????????0??2018-11-18?21:12??調制方式總結\bpsk產生方案\
?????文件????????2490??2018-08-01?20:03??調制方式總結\bpsk產生方案\bpsk.m
?????目錄???????????0??2018-11-18?21:12??調制方式總結\fm產生方案\
?????文件????????1429??2018-08-01?16:07??調制方式總結\fm產生方案\fm1.m
?????目錄???????????0??2018-11-18?21:12??調制方式總結\gmsk產生方式\
?????文件????????6194??2018-08-01?16:49??調制方式總結\gmsk產生方式\gmsk.m
?????文件?????????256??2018-07-31?09:59??調制方式總結\gmsk產生方式\sigexpand.m
?????目錄???????????0??2018-11-18?21:12??調制方式總結\oqpsk產生方案\
?????文件????????4072??2018-08-01?16:31??調制方式總結\oqpsk產生方案\oqpsk3.m
?????目錄???????????0??2018-11-18?21:12??調制方式總結\qpsk產生方案\
?????文件????????3768??2018-08-01?19:39??調制方式總結\qpsk產生方案\qpsk3.m
?????目錄???????????0??2018-11-18?21:12??調制方式總結\upsk產生信號\
?????文件????????4082??2018-08-01?16:38??調制方式總結\upsk產生信號\uqpsk1.m
?????文件?????1047177??2018-08-01?20:07??調制方式總結\數傳信號產生說明.docx