clc;
clear;
fc=4800;fs=12000;fb=2400;
%要調制的數字信號
a=randint（1122）;%隨機產生12個“0”，“1”
s=zeros（160）;
for?i=1:12
????for?j=1:5
if（a（i）==0）
????s（j+（i-1）*5）=0;
else
????s（j+（i-1）*5）=1;
end
????end
end


plot（s）;xlabel（‘基帶信號‘）;
figure
pwelch（s）;%功率譜
figure
%波形成形濾波器（平方根升余弦滾降）
h=firrcos（14120012004800‘sqrt‘）;
figure
stem（h）;xlabel（‘成形濾波器的單位沖擊響應‘）;

[HW]=freqz（h1）;
H=abs（H）;
figure
plot（H）;xlabel（‘成形濾波器的頻率響應‘）;
s=fftfilt（hs）;
figure

plot（s）;xlabel（‘通過成形濾波器后的基帶信號‘）;
figure
pwel