Lipschitz指数matlab程序

Lipschitz指数的matlab程序设计

function sig = GetSignal（）
	% 引用wavelab中的数据即A Wavelet Tour of Signal Processing（2nd edition）中fig6.6的信号
    N = 64;
	i = （1:N）./N;
	x = （1-sqrt（i）） + i/2 -.5;
	M = 8*N;
    load tt.out;
	sig = tt（:2）;
	sig（M-1.5.*N+1:M-.5*N） = x;
	sig（M-2.5*N+2:M-1.5.*N+1） = x（length（x）:-1:1）;
	sig（3*N-30:3*N + N） = .5*ones（1N+31）;
	
	t = （1:M）/M;
	sigma = 0.01;
	g = exp（-.5.*（abs（t-.5）./sigma）.^2）./sigma./sqrt（2*pi）;
	g = fftshift（g）;
	
	sig2 = iconv（g‘sig）‘/M; 
	
	m = max（sig2）;
	j = find（sig2 == m）;
	sig（8*N/2+1:12*N/2） =  sig2（j:j+2*N-1）.*.5./m;
	sig（1:2.5*N） = sig2（M-2.5*N+1:M）;
	sig（M-2.5*N+2:M-1.5.*N+1） = x（length（x）:-1:1）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        660  2009-08-08 14:12  Lipschitz\GetSignal.m

     文件        900  1999-10-03 00:00  Lipschitz\iconv.m

     文件        595  2005-07-04 21:38  Lipschitz\ImageWT.m

     文件       1195  2005-06-30 23:52  Lipschitz\MM_WT.m

     文件       1567  2005-06-22 23:35  Lipschitz\SkelMap.m

     文件       2383  2005-07-04 21:44  Lipschitz\smallExp.m

     文件      24344  2009-08-08 09:32  Lipschitz\tt.out

     目录          0  2009-08-08 14:07  Lipschitz

----------- ---------  ---------- -----  ----

                31644                    8