clear?all;

fs=107;%采樣頻率
Ts=1;%基帶信號(hào)碼元周期
N_sample=fs;%碼元采樣點(diǎn)數(shù)
N=1000;%碼元數(shù)目

%?------?時(shí)域參數(shù)?------?%
dt=1/N_sample;
t=0:dt:（N*N_sample-1）*dt;

%?------?繪制升余弦滾降濾波器的沖激響應(yīng)函數(shù)?------?%
a=1;
t0=-3:dt:3;
ht=sin（pi*t0/Ts）./（pi*t0./Ts）.*cos（a*pi*t0/Ts）./（1-4*a.^2*t0.^2/Ts.^2）;

ht（（length（t0）+1）/2）=1;

figure（1）;
plot（t0ht‘LineWidth‘1.5）;
axis（[-3?3?-0.2?1.2]）;
title（‘升余弦濾波器的沖擊響應(yīng)函數(shù)h（t）‘）;
xlabel（‘時(shí)間（s）‘）;
ylabel（‘電壓值（V）‘）;
grid?on;

%?------?初始化二進(jìn)制數(shù)字信號(hào)碼元?------?%
bc0=[-1zeros（1106）];
bc1=[1zeros（1106）];

%?------?生成二進(jìn)制隨機(jī)基帶信號(hào)?------?%
rand_code=randi（[01]1N）;

%?------?初始化二進(jìn)制基帶脈沖信號(hào)?------?%
mt=zeros（1N*N_sample-1）;

%?------?計(jì)算二進(jìn)制基帶脈沖信號(hào)?------?%
for?i=1:N
????if?rand_code（i）==1
????????mt（（i-1）*N_sample+1:i*N_sample）=bc1;
????else
?????????mt（（i-1）*N_sa