双音多频Goertzel算法Matlab仿真

function Analyser_3D（keytone）
%UNtitleD Summary of this function goes here
%   Detailed explanation goes here
figure（‘Name‘‘Pe-SNR-N关系3D图‘）;
nX_SNR=0;
nY_nGoertzel=0;
for nGoertzel=2:10:800
    nY_nGoertzel=nY_nGoertzel+1;
	for SNR=-10:10
        nX_SNR=nX_SNR+1;
        cnt=0;
        for n=1:10
            noisetone= Noise_Adder（toneSNR）;
            [keyDetected spectrumGeortzel timeGoertzel] = DTMF_Goertzel_Detector（noisetonenGoertzel）;
            if keyDetected==key
                cnt=cnt+1;
            end 
        end
        Pe（nY_nGoertzelnX_SNR）=1-cnt/10;
    end
    nX_SNR=0;
end
[SNRnGoertzel]=meshgrid（-10:102:10:800）;
surf（SNRnGoertzelPe）;
title（‘Pe-SNR-N关系3D图‘）;
xlabel（‘SNR‘）;
ylabel（‘N‘）;
zlabel（‘Pe‘）;
end