%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)<max1
%         max2=absu(i);
%     end
% end
% cw(k)=max1/max2;

% end  %對應18行的for
% 
% i=1:k;

% plot(i,a1,'r+')
% plot(i,a2,'r*')
% plot(i,r)
% plot(i,cw,'go')
% end
hold off;
xlabel('信噪比(dB)')
title('BPSK SFFT)')