|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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 编辑 ]
|
|