matlab扩展编程里面的初始化问题
来源:互联网 发布:c语言程序实例代码 编辑:程序博客网 时间:2024/06/07 21:33
function hmm = inithmm(samples, M)
K = length(samples); %语音样本数
N = length(M); %状态数
hmm.N = N;
hmm.M = M;
% 初始概率矩阵
hmm.init = zeros(N,1);
hmm.init(1) = 1;
% 转移概率矩阵
hmm.trans=zeros(N,N);
for i=1:N-1
hmm.trans(i,i) = 0.5;
hmm.trans(i,i+1) = 0.5;
end
hmm.trans(N,N) = 1;
% 概率密度函数的初始聚类
% 平均分段
for k = 1:K
T = size(samples(k).data,1);
samples(k).segment=floor([1:T/N:T T+1]);
end
%对属于每个状态的向量进行K均值聚类,得到连续混合正态分布
for i = 1:N
%把相同聚类和相同状态的向量组合到一个向量中
vector = [];
for k = 1:K
seg1 = samples(k).segment(i);
seg2 = samples(k).segment(i+1)-1;
vector = [vector ; samples(k).data(seg1:seg2,:)];
end
mix(i) = getmix(vector, M(i));
end
hmm.mix = mix;
function mix = getmix(vector, M)
[mean esq nn] = kmeans(vector,M);
% 计算每个聚类的标准差, 对角阵, 只保存对角线上的元素
for j = 1:M
%ind = find(j==nn);
ind = find(j==mean);
tmp = vector(ind,:);
var(j,:) = std(tmp);
end
% 计算每个聚类中的元素数, 归一化为各pdf的权重
weight = zeros(M,1);
for j = 1:M
%weight(j) = sum(find(j==nn));
weight(j) = sum(find(j==mean));
end
weight = weight/sum(weight);
% 保存结果
mix.M = M;
mix.mean = esq; % M*SIZE
mix.var = var.^2; % M*SIZE
mix.weight = weight; % M*1
K = length(samples); %语音样本数
N = length(M); %状态数
hmm.N = N;
hmm.M = M;
% 初始概率矩阵
hmm.init = zeros(N,1);
hmm.init(1) = 1;
% 转移概率矩阵
hmm.trans=zeros(N,N);
for i=1:N-1
hmm.trans(i,i) = 0.5;
hmm.trans(i,i+1) = 0.5;
end
hmm.trans(N,N) = 1;
% 概率密度函数的初始聚类
% 平均分段
for k = 1:K
T = size(samples(k).data,1);
samples(k).segment=floor([1:T/N:T T+1]);
end
%对属于每个状态的向量进行K均值聚类,得到连续混合正态分布
for i = 1:N
%把相同聚类和相同状态的向量组合到一个向量中
vector = [];
for k = 1:K
seg1 = samples(k).segment(i);
seg2 = samples(k).segment(i+1)-1;
vector = [vector ; samples(k).data(seg1:seg2,:)];
end
mix(i) = getmix(vector, M(i));
end
hmm.mix = mix;
function mix = getmix(vector, M)
[mean esq nn] = kmeans(vector,M);
% 计算每个聚类的标准差, 对角阵, 只保存对角线上的元素
for j = 1:M
%ind = find(j==nn);
ind = find(j==mean);
tmp = vector(ind,:);
var(j,:) = std(tmp);
end
% 计算每个聚类中的元素数, 归一化为各pdf的权重
weight = zeros(M,1);
for j = 1:M
%weight(j) = sum(find(j==nn));
weight(j) = sum(find(j==mean));
end
weight = weight/sum(weight);
% 保存结果
mix.M = M;
mix.mean = esq; % M*SIZE
mix.var = var.^2; % M*SIZE
mix.weight = weight; % M*1
阅读全文
0 0
- matlab扩展编程里面的初始化问题
- 关于游戏里面装备有概率掉落、套装收集问题的matlab模拟编程
- matlab cell无法扩展的问题求解:
- matlab 编程遇到的问题
- matlab 编程的细节问题
- C里面的二维数组的初始化问题
- Matlab里面的指针
- Matlab里面的数据类型
- [javabuilder]MATLAB的问题,javabuilder不能初始化类
- matlab C++ 混合编程出现的问题
- MATLAB面向对象编程属性初始化方法
- 语音信号处理 matlab 扩展编程
- 关于编程之美一个房间里面三盏灯的问题
- vc里面字符串变量的初始化
- switch 里面不同分支的不同初始化
- vc里面进行常量数组的初始化
- 扩展修改PHP.INI里面的值
- 我的第一篇编程日志——初始化问题
- 【Scikit-Learn 中文文档】特征选择
- django contentypes 框架
- 关于IO口上拉、下拉、漏极开漏
- numpy中矩阵相乘的用法:dot函数和乘号*
- error: command 'gcc' failed with exit status 1
- matlab扩展编程里面的初始化问题
- 5.10
- Netty源码分析(八)—内存池分析
- [BZOJ3282]LCT模板(单点修改,区间查询)
- HDU 4421 Bit Magic 【2-SAT】
- 数据hive常见问题
- 解决vue2.0路由跳转未匹配相应路由而出现空白页面的问题
- golang中recover和panic用法
- C3P0连接池XML版