bpsk调制解调matlab仿真

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