求助：关于小波信号奇异点的检测问题！

猴王

首先我做的工作是，模拟一个单层框架的加速度，问题是，无论我怎样改变我的刚度值，和分解的层数，都无法识别出信号的奇异点，请高手指点，感激不尽，时间紧迫，请高手自己看看给出自己见解，小弟谢谢了，以下是我的m文件
function xdot=danziyoudukuangjia(t,x)
F=20*sin(25*t);
  M=31.2;
  C=0.01;
  K=67.51;
xdot=[x(2);(1/M)*(F-C*x(2)-K*x(1))];%对单层框架的速度和位移进行求解的方程式，主要根据结构动力微分方程
执行文件
clc;clear all;
close all;
t0=0;
   tf=60;
   x0=[0,0]';

   [t,x]=ode23('danziyoudukuangjia',[0:0.005:60],x0);
   d=x(:,1);
   v=x(:,2);
   
t=0:0.005:60;

F=20*sin(25*t);
M=31.2;
C=0.01;
K=67.51;
s1=(1/M)*(F-C*v.'-K*d.');%求解加速度的方程
s=s1+randn(1,12001);%加入噪声

%subplot(311);plot(t,s);
%AXIS([0 60 -10 10]);%单层框架加速度信号输出


%下面进行离散的单尺度小波变换并生成，各尺度上的信号


[c,l]=wavedec(s,5,'bior6.8');%对第一信号进行3尺度一维离散小波分解，采用墨西哥小帽函数

%提取结构的低频和高频信号
ca3=appcoef(c,l,'bior6.8',5);%提取第三尺度系数的低频
[cd1,cd2,cd3,cd4,cd5]=detcoef(c,l,[1,2,3,4,5]);%提取第一、二、三尺度系数的高频

%重构信号的低频和高频部分
a3=wrcoef('a',c,l,'bior6.8',3);
d1=wrcoef('d',c,l,'bior6.8',1);
d2=wrcoef('d',c,l,'bior6.8',2);
d3=wrcoef('d',c,l,'bior6.8',3);
d4=wrcoef('d',c,l,'bior6.8',4);
d5=wrcoef('d',c,l,'bior6.8',5);

%显示多尺度一维信号的分解结果


subplot(611);plot(t,a3);title('信号的低频细节部分');%AXIS([0 60 -5 5]);

subplot(612);plot(t,d1);title('Detailsd1');%AXIS([0 60 -5 5]);

subplot(613);plot(t,d2);title('Detailsd2');%AXIS([0 60 -5 5]);

subplot(614);plot(t,d3);title('Detailsd3');%AXIS([0 60 -5 5]);

subplot(615);plot(t,d4);title('Detailsd4');%AXIS([0 60 -5 5]);

subplot(616);plot(t,d5);title('Detailsd5');%AXIS([0 60 -5 5]);

[ 本帖最后由 xinyuxf 于 2007-6-18 13:24 编辑 ]

VibrationMaster

这种问题根本就没有所谓奇异点问题，要说有奇异点，那就是开始的那一时刻和最后那一时刻，这样的东西检测出来有什么意义。
告诉你的老师（不要痛苦），这种题目根本就是错误的，也别指望发文章。反正到我手里是四拉丝拉的干活。

猴王

我明白您的意思，但现实情况好像不允许我像您说的那样做吧，我也是没办法，希望你能帮忙，:handshake

[ 本帖最后由 zhlong 于 2007-6-20 11:43 编辑 ]

VibrationMaster

我记得国外有人用小波求振动响应的问题。你要想混过去，只能把问题搞的复杂一些。
1）激励不要用连续sin cos，想办法用带脉冲或方波之类的东西，这样响应的奇异实际就是激励的奇异点，还可以模拟一点噪声。
2）弄点非线性进来。
3）这种近伪科学的东西（和伪导师给的伪命题），只有把问题搞复杂了，才能把学术混混整晕，从而蒙混过关。

wy558558558

我觉得还是要结合一些实际问题来用小波．如果出这种文章也就能投个垃圾学校的学报了．

zhaopeng80

我也做小波信号奇异点的检测，也碰到这种情况！就是检测不出信号的奇异点！

猴王

我现在用的是在原始信号中加入冲击信号的方法，对于非线性的问题，我可以采用多自由度来解决，但我还有个问题，请各位执教，我用的matlab编译的程序，现在出现了一个问题，如果各位可以帮忙看看，我感激不尽
m文件：
  function xdot=danziyoudukuangjia(t,x)
F=20*sin(25*t);
  M=31.2;
  C=0.01;
  K=67.51;
xdot=[x(2);(1/M)*(F-C*x(2)-K*x(1))];



执行文件：
clc;clear all;
close all;
t0=0;
   tf=60;
   x0=[0,0]';
   [t,x]=ode23('danziyoudukuangjia',[0:0.005:60],x0);
   d=x(:,1);
   v=x(:,2);
   
t=0:0.005:60;
F=20*sin(25*t);                                
M=31.2;
C=0.01;
K=67.51;
s1=(1/M)*(F-C*v.'-K*d.');
T=12001;
t=0:0.002:60;
x2=zeros(1,T);
x2(6001)=1;
s=s1+randn(1,12001)+x2;

%subplot(311);plot(t,s);
%AXIS([0 60 -10 10]);%单层框架加速度信号输出

%下面进行离散的单尺度小波变换并生成，各尺度上的信号

[c,l]=wavedec(s,5,'bior6.8');%对第一信号进行3尺度一维离散小波分解，采用墨西哥小帽函数
%提取结构的低频和高频信号
ca3=appcoef(c,l,'bior6.8',5);%提取第三尺度系数的低频
[cd1,cd2,cd3,cd4,cd5]=detcoef(c,l,[1,2,3,4,5]);%提取第一、二、三尺度系数的高频
%重构信号的低频和高频部分
a3=wrcoef('a',c,l,'bior6.8',3);
d1=wrcoef('d',c,l,'bior6.8',1);
d2=wrcoef('d',c,l,'bior6.8',2);
d3=wrcoef('d',c,l,'bior6.8',3);
d4=wrcoef('d',c,l,'bior6.8',4);
d5=wrcoef('d',c,l,'bior6.8',5);
%显示多尺度一维信号的分解结果

subplot(611);plot(t,a3);title('信号的低频细节部分');%AXIS([0 60 -5 5]);
subplot(612);plot(t,d1);title('Detailsd1');%AXIS([0 60 -5 5]);
subplot(613);plot(t,d2);title('Detailsd2');%AXIS([0 60 -5 5]);
subplot(614);plot(t,d3);title('Detailsd3');%AXIS([0 60 -5 5]);
subplot(615);plot(t,d4);title('Detailsd4');%AXIS([0 60 -5 5]);
subplot(616);plot(t,d5);title('Detailsd5');%AXIS([0 60 -5 5]);