samp10_5.m

MATLAB7.x数字信号处理 光盘内容

%Samp10_5
x=zeros(1,500);   %设置空矢量
%给出检测数据
for t=1:1:500  
    if (t<200)
     x(t)=50.*exp(t/300).*sin(2*pi*0.01*t);
 elseif (t>=200)&(t<300)
   x(t)=50.*exp(t/300).*sin(2*pi*0.01*200)+30;
elseif ((t>=300)&(t<400))
    x(t)=50.*exp(t/300).*sin(2*pi*0.01*200)+10;
else
    x(t)=50.*exp(t/300).*sin(2*pi*0.01*200)+10+250*sin(2*pi*0.003*t);
end
end
%结束检测数据
x=x';   %将检测数据转置以备后面应用
x=x-mean(x);   %去掉平均值
    a=1;   %给出小波尺度为1
for n=(1:1:500)
    t0=n*1;
    [g,fai,v]=gdttest(a,t0,500);   %调用小波函数的计算
    xx(n)=sum(x.*g');    %一阶导数的计算（墨西哥草帽波）
    yyy(n)=sum(x.*(fai)');%高斯函数的计算
    zzz(n)=sum(x.*v');   %二阶导数的计算
end
subplot(3,1,1),plot(x);   %绘出原始数据波形
ylabel('检测数据')
subplot(3,1,2),plot(yyy);   %绘出高斯函数为基函数的小波变换
ylabel('高斯函数结果')
subplot(3,1,3),plot(zzz); %绘出高斯函数的一阶导数为基函数的小波变换
ylabel('一阶导数结果')
xlabel('时间/s')