phase315.m

来自「是信号处理方面欠采样的源程序 ,是很不错的程序,可以这方面的学者一些参考价值」· M 代码 · 共 23 行

clear all;
clc;
close all;
unit=1E6;
f0=150*unit;
fs=100*unit;
del=1E-9;
T=200*1E-6;
t=0:1/fs:T;
 SNR=-10;
x=exp(j*2*pi*f0*t);
x=awgn(x,SNR);
x1=exp(j*2*pi*f0*t+j*2*pi*f0*del);
% 
% y=abs(fft(x));
% y1=abs(fft(x1));
% [peak,ipeak]=max(y);
% [peak1,ipeak1]=max(y1);
r=x.*conj(x1);
r1=sum(r)/length(r);
% r1=x(ipeak)*conj(x1(ipeak1));

f=abs(1/2/pi/del*angle(r1));