clc;
close all;
clear;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% 假定:
% 2倍载波频率采样的bpsk信号
% 调制速率为在波频率的 N/2m
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
m=128;
N=512;
n=1:1:N;
N0=0.5*randn(1,N) %噪声
h0=zeros(1,N);
% 30阶低通滤波器 h0
f = [0 0.3 0.3 1]; w0 = [1 1 0 0];
b = fir2(30,f,w0);
[h,w] = freqz(b,1,N/2);
h0(1,1:N/2)=abs(h');
for i=1:N/2
h0(1,N-i+1)=h0(1,i);
end;
%%%%%%%%% 随机序列
a=rand(1,m);
for i=1:m
if(a(1,i)>0.5)
a(1,i)=1;
else
a(1,i)=-1;
end;
end;
a
%%% 生成BPSK信号
bpsk_m=zeros(1,N);
j=1;k=1;
for i=1:N
if(j==(