多尺度小波模极大值法检测频谱边缘matlab

%求信号x的功率谱
%时域信号--（FFT变换）-->幅度谱--（平方）-->功率谱  直接法
%时域信号--->相关函数--（FFT变换）-->功率谱--（除以频率分辨率）-->功率谱密度  间接法
clear all
%a函数及其频谱
m=-20:0.2:20;
index1=find（m<0）;
a = 10*sin（m（index1））./（m（index1））;
a=[a 10];
index2 = find（m > 0）;
a = [a 10*sin（m（index2））./（m（index2））];
figure（1）;subplot（411）;stem（ma）;title（‘原函数a‘）;
M=length（m）;
p=0:M-1;q=0:M-1;
A=a*exp（-1i*2*pi/M）.^（p‘*q）;
subplot（412）;plot（fftshift（abs（A）））;title（‘原函数a的幅度谱‘）;
%调制函数b
t=0:1/100:2;
b=100*sin（2*pi*20*t）;
f=linspace（0100128）;
subplot（413）;plot（tb）;title（‘其中的一个调制函数b=50sin（2pi20t）‘）;
B=fft（b128）;
subplot（414）;plot（fabs（B））;title（‘其中的一个调制函数b的幅度谱‘）;
%调制函数c
c=80*sin（2*pi*40*t）;
ba=a.*b;ca=a.*c;
x=ba+ca+randn（size（ba））;
X=abs（fft（x128））; 
nfft=512;
index=0:nfft-1;
k=index*100/nfft;%100采样频率
[r tao]=xcorr（x‘unbiased‘）;
subplot（412）;plot（taor）;title（‘x的自相关函数‘）;
Pxx=abs（fft（rnfft））;
Sf=10*log10（Pxx（index+1））;  %功率谱密度函数
n0=length（X）;
n1=length（Sf）;
%  多尺度
m=-4:-1;
for i=1:4
    delta（i）=2.^m（i）;
end
N=16; % 滤波器长度（需要调整，必须是偶数）
A=-1/sqrt（2*pi）;  %  幅度
%  构造高斯函数
for index_x=1:N;
    x=index_x-（N