bpsk_ph.m

BPSK调制程序

%BPSK并频谱分析
%产生频率调制信号
%被调信号为m(t)=sinc(2*pi*20*t),Fm=20,Fc=200,Fs=2000,Kf=100
clear;
df=0.25;
fs=20;%20MHz
fc=2;
rb=0.2;
ts=1/fs;
m=10;%m是码元数，自己随便定义
n=40;%n是每个码元取点数，即对应的正弦波4个周期取多少个点
a=sj(m,n);
%sj(m,n)的解释参见sj.m,sj(m,20)的作用相当于上面的for循环
t=0:1:m*n-1;
Fc1=sin(8*pi*t/n);%码元0发射的载波
Fc2=sin(8*pi*t/n+pi);%码元1发射的载波
b=1-a;
u0=a.*Fc2+b.*Fc1;

%产生白噪声
snr=10;%S/N
snr_lin=10^(snr/10);%线性信噪比
signal_power=(norm(u0)^2)/length(u0);
noise_power=signal_power/snr_lin;
noise_std=sqrt(noise_power);
noise=noise_std*randn(1,length(t));
u1=u0+noise;%加白噪声
u=u1;

phs=xiangweicha(u,ts,fc);
subplot(211)
plot(a);
grid
subplot(212)
plot(phs);
grid