声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1193|回复: 0

[其他] 求助:轴承信号仿真

[复制链接]
发表于 2009-6-12 11:47 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 wdhd 于 2016-9-18 14:52 编辑

  数学模型:
  x(t)=∑m(kT)exp[-a(t-kT)]sin[2*pi*f0(t-kT)+r]U(t-kT)+n(t);
  ∑求和,k=-∞到∞;
  其中:m(kT)为kT时刻的冲击幅值,m(t)为幅值调制函数,T为故障特征周期,f0为轴承座-传感器系统的某一固有频率,a为对应系统的阻尼比,U(t)为单位阶跃函数,n(t)为噪声信号。
  设仿真参数为:m(kT)=5.0,T=1/65,f0=6450Hz,a=0.04,r=0,叠加白噪声的方差为1.05;
  我写了一段程序,调试不出来:
  syms k;
  fs=16384;
  t=0:1/fs:0.5-1/fs;
  T=1/65;
  f0=6450;
  a=0.004;
  N=length(t);
  n=randn(1,N);
  n=n-mean(n);
  n=sqrt(1.05)*n./sqrt(var(n));%方差为1.05的白噪声
  sig=5*exp(-a(t-k/65)).*sin(2*pi*f0.*(t-k/65)).*heaviside(t-k/65)+n;
  sig=symsum(sig,k,-inf,inf);
  figure(1);
  plot(t,sig)
  向各为高手求助,请各位高手指正,不胜感激。
  [ 本帖最后由 linlyjiang 于 2009-6-12 12:02 编辑 ]
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 04:56 , Processed in 0.076662 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表