过零法测频率matlab程序

clear all;
clc;

k = 54000;
n=0.001;
dbzy=15;
timestep=0.0000001;
ncy= 30;
tcy=ncy*timestep
fcy=1/tcy;
alltime = 0;
ntime=-1;
stepdbzy=0.1;

% t = 0:timestep:n*pi;
% x = sin（k*2*pi*t）; % Create sawtooth signal.
% y = awgn（xdbzy‘measured‘）; % Add white Gaussian noise.
% plot（txty） % Plot both signals.
% legend（‘Original signal‘‘Signal with AWGN‘）;
% title（‘带有高斯白噪声的回波信号‘）;
% xlabel（‘时间t （seconds）‘）;
% ylabel（‘幅度‘）;

t = 0:timestep:n*pi;
x = sin（k*2*pi*t）;
xx = 0;

    
    
  
    y = awgn（xdbzy‘measured‘）; % Add white Gaussian noise.
    for i=1:ncy:l