- 博客(5)
- 资源 (1)
- 收藏
- 关注
原创 Levinson-Durbin算法的实现
在实现维纳滤波器和预测器的时候,需要计算数据的自相关矩阵的逆。但是当数据量比较大的时候,计算矩阵的逆花费的代价比较大,所以需要使用Levinson-Durbin算法来实现系数的求解。一、数据模型k阶前向维纳预测器: 对上述模型进行变换: 所以a(0) ...
2018-04-01 21:28:34
15185
2
原创 混合高斯模型GMMS的EM算法实现
最近学习机器学习课程,需要实现混合高斯模型的求解。因为混合高斯模型的优化函数,似然函数难以求导优化,所以使用EM算法优化。下面是EM算法的实现。主要参考Andrew Ng的讲义上的公式实现的。其中,j代表第几个类别,i代表第几个样本数,phi(j)代表类别出现的概率,即p(z(i)=j; phi)。数据满足的模型是混合高斯模型,即p(x|z)是高斯分布,由多个高斯分布以一定的权重相加。下面...
2018-03-30 22:50:50
2189
原创 matlab关于fft的应用(二)
y = fft(x,n); 其中n是关于x做n点fft,若x小于n点,则对x做补零处理;若x大于n点,则对x做截断处理。下面讨论补零和截断对fft的影响。一、补零 matlab中fft的计算分辨率为fs/n,其中L是计算fft的点数。所以,若对x做补零处理,则会增大fft的计算分辨率。这里需要注意的是,补零并不会提高对信号的实际分辨率。举个例子,原本信号中100hz和110hz的信号...
2018-03-20 22:56:39
2052
原创 matlab关于fft的应用(一)
在做功率谱分析的时候,周期图法会用到fft()。 fft()是matlab对信号做频谱分析的一个函数。 语法:Y=fft(x); 在实际的处理当中,首先是对数据进行采样,此时就需要有采样率fs,那么采样时间T=1/fs。因此数据就可以根据y=f(nT)进行采样。所以做fft时,就会用到如下采样代码:f0 = 100;fs = 500;Ts = 1/fs;n=1:1:...
2018-03-18 23:00:13
19544
5
原创 matlab关于xcorr函数的用法
在做信号处理的时候,经常会对信号做自相关处理,比如对信号做功率估计,或者是参数拟合。在机器学习领域,如wule-walker方程也会遇到自相互函数的处理。1 自相关矩阵的基本概念 首先给出自相关函数的定义: R(s,t)=E[(x(s)x(t))]其中s,t是x(n)的不同时刻的信号值。 在实际的应用当中,我...
2018-03-18 18:36:42
50641
2
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人