bpsk_r.m

BPSK的R参数分析归一化,求瞬时平方包络均值

%BPSK的R参数分析归一化
figure(24)
hold on;
for ii=1:100;
    ii
clear
df=0.15;
fs=20;%20MHz
fc=2;
rb=0.2;
ts=1/fs;
m=256;
n=40
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);

b=1-a;
u0=a.*Fc2+b.*Fc1;

for k=1:20   %对应程序结尾的end
k
%产生白噪声
snr=k;%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(u0));
u1=u0+noise;%加白噪声
uls=length(u1);
maxu1=max(u1);
u=u1/maxu1;


%求瞬时平方包络
v=blpf(u);
%均值
miu=mean(v);
%方差
delt=var(v);
%求R参数
r(k)=delt/(miu^2);


% plot(i,a1,'r+')
% plot(i,a2,'r*')
plot(i,r)
% plot(i,cw,'go')
end
hold off;
xlabel('信噪比(dB)')
title('归一化BPSK R参数(仿真100次)')