求助：轴承信号仿真

linlyjiang

　　数学模型：
　　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 编辑 ]