cei.m

本代码是关于用matlab算出两路信号的相位差

%clear all;
clc
%------make a signal--------------
fs=100000;  
f0=4900;
n=0:1023;
% xx=exp(j*(2*pi*f0*n/fs+pi/3));
 xx=sin(2*pi*f0*n/fs+0);
 yy=sin(2*pi*f0*n/fs-pi*100/180);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
y=fft(xx);    %--fft
m=abs(y);     %--make the line (a to f)
p=angle(y);   %--make the line (p to f)
pp0=180*p/pi ;  %--把弧度换成角度
leng_fft=length(y)   %--FFT的点数就是信号样本点的个数 
f=(0:length(y)-1)*fs/length(y);  %--把X轴换成频率为单位,最大值是采样频率100K  
s=length(y)*f0/fs                %--信号频率
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
y1=fft(yy);    %--fft
mm=abs(y1);     %--make the line (a to f)
p1=angle(y1);   %--make the line (p to f)
pp1=180*p1/pi ;  %--把弧度换成角度
leng_fft1=length(y1)   %--FFT的点数就是信号样本点的个数 
f1=(0:length(y1)-1)*fs/length(y1);  %--把X轴换成频率为单位,最大值是采样频率100K  
s1=length(y1)*f0/fs                %--信号频率
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 yy1=m(1);
 ttt=leng_fft/2
 for t=1:ttt       %求最大波束位置                   %
     if yy1<m(t)
        yy1=m(t);
        zz1=t;        %波束最大时所对应的波束数
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  yy2=mm(1);                                         %
%  for t=1:leng_fft        %求最大波束位置                   %
%      if yy2<mm(t)
%         yy2=mm(t);
%         zz2=t;        %波束最大时所对应的波束数
%     end
% end
%  %   yy1=B(zz1-1);        %纵座标
%   %  xx1=fai(zz1-1) ;     %横座标
    
     x1=pp0(zz1)    %横座标
     x2=pp1(zz1)    %横座标
     
     aaaa=x2-x1
     
 
     
%     %yy3=B(zz1+1) ;     %纵座标 
%     %xx3=fai(zz1+1) ;   %横座标
%    % 此方法一点好处都没有
%  %  XX=[xx1.^2,xx1,1;xx2.^2,xx2,1;xx3.^2,xx3,1];         %%%%%%%%%%%
%  %  YY=[yy1;yy2;yy3];                                    %%%%%%%%%
%  %  CC=XX\YY;                                             %%%%%%%%
%  %  aa=-CC(2)/CC(1)/2       %以抛物线为模型求解方位%%%%%%%%%%
%  %  bb=abs(fai(zz1)-jiaodu) %最大波束与真实角度的差别   %%%%%%%
%  %  cc=abs(aa-jiaodu)       %与真实角度的差别           %%%%%%%
%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %     figure          
    subplot(3,1,1);
    plot(f,pp0)        
    title('');
    xlabel('');
    ylabel('v');
    %%%%%%%%%%%%%%%%%%
    subplot(3,1,2);
    plot(f,pp1)         
    title('');
    xlabel('');
    ylabel('');
    %%%%%%%%%%%%%%%%%%
%     subplot(3,1,3);
%     plot(f,pp)        
%     title('');
%     xlabel('');
%     ylabel('');
    subplot(3,1,3);
   plot(n,xx,'r*-', n,yy,'bo--')        
    title('');
    xlabel('');
    ylabel('');