尺度-小波能量谱程序的问题，求解答！

lonsvm

仿照《尺度- 小波能量谱在滚动轴承故障诊断中的应用》程军圣，一文编写程序。
其中用的小波是自定义小波，想做小波能量谱，但现在遇到的问题是，对自定的小波做连续小波变换报错。

顺便问一句，尺度-小波能量谱的程序到底应该怎么编写？

这是自定义的小波函数，后面程序调用时，总是报错，提示omiga未定义

1. function [psi,t] = pulse(LB,UB,N,sigma,omiga)
   
2. %自定义小波函数  
   
3. %sigma改变变换尺度,omiga改变振荡周期
   
4. t = linspace(LB,UB,N);   %参考morlet小波的形式写的
   
5. t1=t(1:1:N/2-1);       %函数分段，t<0
   
6. psi1=exp(-sigma*omiga.*t1).*sin(omiga.*t1);
   
7. t2=t(N/2:1:end);     %t>0

复制代码

psi=[psi2 psi1];
end       %为什么这两行总是自动被取消了。

将其添加到小波工具箱

1. clc
   
2. clear all
   
3. close all
   
4. % 创建一个第4类型的小波
   
5. familyName      = 'Mywave_Pulse';  %小波全名
   
6. familyShortName = 'puls';          %小波简称,缩写名，不大于4个字母的字符串
   
7. familyWaveType  = 4;  %小波类型
   
8. familyNums      = '';                       %只有一个小波，序列号留空
   
9. fileWaveName    = 'pulse';           %创建小波的函数文件名
   
10. familyBounds    = [-4 4];               %定义有效支撑集
    
11. % 将新创建的小波添加到matlab中
    
12. wavemngr('add',familyName,familyShortName,familyWaveType,...
    
13.     familyNums,fileWaveName,familyBounds);
    
14. % 查看自定义小波是否创建成功.
    
15. wavemngr('read')

复制代码

用自定义小波，执行小波连续变换。 用morl小波，和shannon小波可以执行。自定义出错。

1. %调用编写的小波基函数pulse.m，已用add_custom_wave.m文件将其添加到小波族函数中
   
2. clc; clear;close all;
   
3. lb = -4; ub = 4;
   
4. n=100;
   
5. sigma=0.1 ; omiga=7;
   
6. [psi2,x2] = pulse(lb,ub,n,sigma,omiga);   %自定义小波-脉冲小波
   
7. figure(1);
   
8. plot(x2,psi2); title('me');grid on;
   
9. 
   
10. load sumsin; x = sumsin;
    
11. scales = 1:64;
    
12. 
    
13. %x=load('org_1.txt');  %1024
    
14. 
    
15. wname='morl';
    
16. figure(2)
    
17. coefs1 = cwt(x,scales,wname,'plot');
    
18. 
    
19. wname='shan1-1.5';
    
20. figure(3)
    
21. coefs2 = cwt(x,scales,wname,'plot');
    
22. 
    
23. wname='puls';
    
24. figure(4)
    
25. coefs = cwt(x,scales,wname,'plot');

复制代码

以下是出错信息。
??? Input argument "omiga" is undefined.

Error in ==> pulse at 18
psi1=exp(-sigma*omiga.*t1).*sin(omiga.*t1);

??? Input argument "omiga" is undefined.

Error in ==> pulse at 18
psi1=exp(-sigma*omiga.*t1).*sin(omiga.*t1);

寂寞的部落

lonsvm

有大神帮我看看吗

xiaosong02

韵天之色