x=linspace(0,2*pi);
y=sin(x);
plot(x,y);
w=randn(1,250);
w=w/std(w);
w=w-mean(w);
a=0.0128;
b=sqrt(0.9596);
w=a+b*w;
plot(w,'o');
N=0:1000;
fs=1024;
t=N/fs;
y=10*sin(2*pi*t);
x=randn(1,1001);
w=y+x;
subplot(3,1,1),plot(x),xlabel('时间(ms)'),ylabel('幅度'),title('高斯白噪声');
subplot(3,1,2),plot(y),xlabel('时间(ms)'),ylabel('幅度'),title('正弦信号');
subplot(3,1,3),plot(w),xlabel('时间(ms)'),ylabel('幅度'),title('正弦信号叠加高斯白噪声');
f1=500;%滤波器阶数
fs=1000;%滤波器显示频率最大值的两倍
ff=[0 .51 .52 .68 .69 1];%归一化的截止频率
fa=[0 0 1 1 0 0];
h=remez(f1,ff,fa);%生成带通滤波器
figure(4),
freqz(h,1,512,fs);%绘出滤波器的频率响应
r=filter(h,1,e2psk);%信号通过带通滤波器
k=5;
plotspec(r,Ts,k)
num=10; %码元个数
tnum=200;%码元长度
N=num*tnum;%10个码元整体长度
a=randi(1,num,2); %产生1行num列的矩阵，矩阵内0和1随机出现
fc=0.25; %载波频率为0.5
t=0:0.05:9.99;%t从0到9.99，间隔为0.05
s=[];c=[];
for i=1:num %i从1到10循环
if(a(i)==0)
A=zeros(1,tnum); %i=0时，产生一个码元长度为tnum（200）的0码元
else
A=ones(1,tnum); %i=1时，产生一个码元长度为tnum（200）的1码元
end
s=[s A]; %s为随机基带信号
cs=sin(2*pi*fc*t);
c=[c cs]; %c为载波信号
end
xlabel('BPSK载波信号');
ylabel('BPSK载波信号幅值');
plot(e)
clear
fs=8e5;%采样频率
fm=20e3;
n=2*(6*fs/fm);
final=(1/fs)*(n-1);
fc=2e5; % 载波频率
t=0:1/fs:(final);
Wc=2*pi*fm*t;
A=1;
phi=0;%初始相位
x = A * cos(Wc + phi);
x(x>0)=1;
x(x<0)=-1;
y=cos(2*pi*fc*t);%余弦载波波形
bpsk=x.*y;%调制bpsk载波
plot(t,bpsk);
axis([0 100e-6 -2 2]);
xlabel('时间(s)');
ylabel('幅值');
title('bpsk调制信号');
