www.pudn.com > bpsk_sfft.rar > bpsk_sfft.m

%BPSK的SFFT参数分析 
figure(74) 
hold on; 
% for ii=1:100; 
clear 
df=0.15; 
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; 
 
% for k=1:20   %对应程序结尾的end 
 
%产生白噪声 
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(u0)); 
u1=u0+noise;%加白噪声 
u=u1; 
subplot(211) 
plot(u) 
subplot(212) 
specgram(u,40,1,hanning(40),0); 
 
% %求瞬时平方包络 
% v=blpf(u); 
% %均值 
% miu=mean(v); 
% %方差 
% delt=var(v); 
% %求R参数 
% r(k)=delt/(miu^2); 
 
%  
% % %信号平方 
% % sqrtu=u.*u; 
%  
% %求功率谱 
% nfft=1024; 
% u2=Rx_set(u,nfft);  %这个100是取了100个自相关系数a 
% u3=fftshift(abs(fft(u2)));%求功率谱 
% ff=-1*fs/2:fs/nfft:fs/2; 
%  
% %求AWGN因子 
% e1=mean(u3); 
% j=0; 
%   for i=1:nfft+1 
%       if u3(i)>e1 
%           j=j+1; 
%           ee(j)=u3(i); 
%       end 
%   end 
%  e2=mean(ee);   
%  j=0; 
%   for i=1:nfft+1 
%       if u3(i)>e2 
%           j=j+1; 
%           eee(j)=u3(i); 
%       end 
%   end 
% e3=mean(eee);  
% a1(k)=(e2-e1)/e1; 
% a2(k)=(e3-e2)/e1; 
%  
% %求CW因子 
% %频谱分析 
% [U,u1,df1]=fftseq(u,ts,df); 
% U=U*ts; 
% f=[0:df1:df1*(length(u1)-1)]-fs/2; 
% absu=abs(U); 
% max1=max(absu); 
% max2=0; 
% for i=1:length(U) 
%     if absu(i)>max2 & absu(i)