BPSK带宽频率估计

%-------------------------------------------------------%
%BPSK信号产生，500个点，100个bit
%对BPSK信号进行带宽估计，重心法频率估计
%-------------------------------------------------------%
clear;

m=100;                      %m个bit
r=5;                            %一个bit长度为5
s=randsrc(m,1);       %随机产生m个bit,+1/-1
s=rectpulse(s,r);       %基带信号产生，每个bit含有5个点
subplot(221);
plot(s);
axis([1 length(s) -1.5 1.5]);title('基带信号');

fc=0.300;                       %载波为0.3，归一化频率
SNR=10;
x=s.*exp(2*1i*pi*fc*(1:length(s))');%x=s.*{cos(2*pi*fc*t)+j*sin(2*pi*fc*t)},t=1:length(s)
noise=sqrt(1/10^(SNR/10)/2)*(randn(size(x))+1i*randn(size(x)));%通过过采样方式加入噪声
x=x+noise;%带噪声的
subplot(222);
plot(linspace(0,1,length(x)),abs(fft(x)));title('FFT');

[Pxx1, f1]=pwelch(x,hamming(100),100/2,length(x),1);%功率谱密度计算
Pxx1=