mulDNewton
本代码为牛顿下山法求解非线性方程组。其调用格式为[r,m]=mulDNewton(F,x0,eps) 其中F:非线性方程组,x0:初始解,eps:解的精度;r:求得的一组解,n:迭代步数。(The code for Newton' s method for solving nonlinear equations. Its call format [r, m] = mulDNewton (F, x0, eps) where F: non-linear equations, x0: initial solution, eps: solution accuracy r: obtain a set of solution, n : the number of iterative steps.)
- 2010-01-01 14:18:33下载
- 积分:1
xiaobosuanfa
采样频率
fs=10000
轴承外环故障信号
fid=fopen( bearingout.dat , r ) 故障
N=1024
xdata=fread(fid,N, int16 )
fclose(fid)
xdata=(xdata-mean(xdata))/std(xdata,1)
时域波形
figure(1)
plot(1:N,xdata)
xlabel( 时间 t/n )
ylabel( 电压 V/v )
db10小波进行4层分解
一维小波分解
[c,l] = wavedec(xdata,4, db10 )
重构第1~4层细节信号
d4 = wrcoef( d ,c,l, db10 ,4)
d3 = wrcoef( d ,c,l, db10 ,3)
d2 = wrcoef( d ,c,l, db10 ,2)
d1 = wrcoef( d ,c,l, db10 ,1)
( Sampling frequency fs = 10000 bearing outer ring fault signal fid = fopen (' bearingout.dat' , ' r' ) failure N = 1024 xdata = fread (fid, N, ' int16' ) fclose (fid ) xdata = (xdata-mean (xdata))/std (xdata, 1) time-domain waveform figure (1) plot (1: N, xdata) xlabel (' Time t/n' ) ylabel ( ' voltage V/v' ) db10 wavelet decomposition 4 layer one-dimensional wavelet decomposition [c, l] = wavedec (xdata, 4, ' db10' ) 1 ~ 4 reconstructed detail signal d4 = wrcoef (' d' , c, l, ' db10' , 4) d3 = wrcoef (' d' , c, l, ' db10' , 3) d2 = wrcoef (' d' , c, l, ' db10' , 2) d1 = wrcoef (' d' , c, l, ' db10' , 1) )
- 2011-05-21 16:48:36下载
- 积分:1