第六章 小波变换与应用.docVIP

第六章 小波变换与应用 张飞 X6-1.某一直升机齿轮箱上的早期损伤的特征是，在伴随振动信号上产生一个可变周期的非平稳扰动信号。要求按此特征构造一组有物理意义的数据仿真早期故障信号，应用小波变换技术分析故障的特征和发生时刻。假定选取高斯小波函数： 其频谱函数为: (1) 证明小波变换为： (2) 证明当尺度因子为时，小波频谱的半功率带宽为1.776; (3) 将表示为尺度因子的线性函数，即，其中和为常数。假定信号的截止频率为，取；试确定常数和，使不同尺 度的小波族完全覆盖被分析信号的频带，并画出小波族的频谱； (4) 按小波变换的离散形式 其中 进行小波变换，其结果用像素的两维等高线来表示 解：(1) 证明： (2) 小波的频谱半功率带宽计算式为： 解之可得： (3)依题意，为使不同尺度的小波族完全覆盖被分析信号的频带，需满足如下方程组： 对，可先取定，再解以上方程组即可得到合适的常数和。 根据上述分析，取，用Matlab绘出小波族的频谱如下图： 可见小波族的频谱很好的覆盖了被分析信号的频带，其实现程序如下： w=1:1000;w0=5;wc=1000;M=10; s=1; sigma=w0/1.1776; d=(wc/M-2*w0)/(M-s); c=w0-d*s; for s=1:M w0=c+d*s for k=1:1000 FY(k)=sqrt(pi/s)/sigma*exp(-(w(k)-s*w0).^2/(4*s^2*sigma^2)); end plot(w,FY);hold on; end （4） 构造一个持续时间为1秒正弦信号，在0.5秒钟之后，加入可变周期的非平稳扰动信号fs=1024; N=1024; T=1;%信号持续时间1s signal=2*sin(2*pi*f1*[0:1/fs:T-1/fs]);%信号长度1s noise=zeros(1,fs*T); t1=[T/2:1/fs:T-1/fs]; no=cos(2*pi*f2*t1.^2);%可变周期的非平稳扰动信号 for i=N*T/2+1:N*T; noise(i)=no(i-N*T/2);%噪声从第0.5s时开始存在 end t=[0:1/fs:T-1/fs]; d=signal; x=signal+noise; figure(1); subplot(2,1,1);plot(t,d); subplot(2,1,2);plot(t,x); X=fft(x,1024); S=fft(d,1024); w0=10;wc=1024;M=10; s=1; sigma=w0/1.1776; d=(wc/M-2*w0)/(M-s); c=w0-d*s; WTs=zeros(M,1024); WTx=zeros(M,1024); for s=1:M w0=c+d*s; phai=exp(-(sigma*s*t).^2).*exp(-j*w0*t); PHAI=fft(phai,1024); WTs(s,:)=ifft(PHAI.*S,1024); WTx(s,:)=ifft(PHAI.*X,1024); end figure(2); subplot(2,1,1);contour(WTs); subplot(2,1,2);contour(WTx); 运行上述程序可得到原始信号和实际观测信号如下图所示，对它们分别进行小波变换后的结果如图-3。由图-3的小波变换结果可以看出，实际信号的后半段出现了很多高频分量，也就是引进的噪声，由此可以知道在0.5秒时系统出现了故障。 图2 图-3