如何确定奇异点位置

zsaliwei007

各位老师，我用小波来求奇异点位置，这个例子如下：
load freqbrk; x= freqbrk;
% 对信号进行傅里叶变换
f = fft(x,1024); f = abs(f);
figure(1); subplot(211); plot(x); subplot(212); plot(f);
% 使用DB6小波进行6层分解
[c,l]=wavedec(x,6,'db6');
figure(2);
subplot(8,1,1); plot(x); ylabel('x');
% 对分解的第6层低频系数进行重构
a=wrcoef('a',c,l,'db6',6);
subplot(8,1,2); plot(a); ylabel('a6');
for i=1:6
    % 对分解的第6层到第1层的高频系数分别进行重构
    d=wrcoef('d',c,l,'db6',7-i);
    subplot(8,1,i+2); plot(d); ylabel(['d',num2str(7-i)]);
end
这个程序运行后，从图中可以在第一层高频系数d1和第二层高频系数d2中看出这个信号的不连续点，这个不连续点是第一类间断点。
我现在应该增加什么程序可以让matlab自己计算出间断点的位置呢？
非常感谢各位老师！

[ 本帖最后由 zhangnan3509 于 2010-4-9 16:19 编辑 ]

anfengping

是啊 请各位高人指教！！！

aprilcat

试下模极大值的方法。图中d1，d2分量在某点处出现了局部的极值，通过求得局部极值的方法，可以求出极值点的坐标位置。

tangyuanyuan

既然搂主也在研究信号奇异性问题，我就占个座位，也问问相关问题啊。搂主。
    信号的lipchitz指数怎么求？
    比如信号f(t)=1-abs(0.5-t)^1.1，我们知道在0.5处其lip指数为1.1。但是针对比较复杂的信号(没有具体的表达式)，如何求其lip指数呢？有一种近似方法，通过小波计算的，即通过小波模最大值来计算lip指数。我说的不是这种求法，是直接针对信号来求得lip指数，有大大虾知道么？搂主知道么？