clear?all;

fs=107;%采樣頻率
Ts=1;%基帶信號碼元周期
N_sample=fs;%碼元采樣點數
N=1000;%碼元數目

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

%?------?繪制升余弦滾降濾波器的沖激響應函數?------?%
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（‘升余弦濾波器的沖擊響應函數h（t）‘）;
xlabel（‘時間（s）‘）;
ylabel（‘電壓值（V）‘）;
grid?on;

%?------?初始化二進制數字信號碼元?------?%
bc0=[-1zeros（1106）];
bc1=[1zeros（1106）];

%?------?生成二進制隨機基帶信號?------?%
rand_code=randi（[01]1N）;

%?------?初始化二進制基帶脈沖信號?------?%
mt=zeros（1N*N_sample-1）;

%?------?計算二進制基帶脈沖信號?------?%
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