求助：轴承信号仿真
本帖最后由 wdhd 于 2016-9-18 14:52 编辑数学模型：
x(t)=∑m(kT)exp[-a(t-kT)]sinU(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 编辑 ]
页:
[1]