BP-HMM 源代码

function [hmm] = obsinit (X,hmm,options) % function [hmm] = obsinit (X,hmm,options) % % Initialise observation model in HMM % % …

function [hmm]=hmmtrain(X,T,hmm) % function [hmm]=hmmtrain(X,T,hmm) % % Train Hidden Markov Model using Baum Welch/MAP EM algorit…

function [hmm] = obsupdate (X,T,Gamma,Gammasum,hmm,update) % function [hmm] = obsupdate (X,T,Gamma,Gammasum,hmm,update) % % Upda…

% HMM's % HMM_Backward - HMM backward algorithm % HMM_Boltzmann - Find the transition matrices of an HMM using …

load samples; % tic, for i=1:length(samples) sample=[]; for k=1:length(samples{i}) sample(k).wave=samples{i}{k}; …

function [prob,q] = viterbi(hmm, O) %Viterbi算法 %输入: % hmm -- hmm模型 % O -- 输入观察序列, N*D, N为帧数,D为向量维数 %输出: % prob -- 输出概率 …

clc disp('加载初始模型HMM...'); load hmm; fs=11025; disp('准备语音输入...') pause; y=wavrecord(3*fs,fs); disp('正在播放刚录入的声音...'); pause;…

function hmm = inithmm(samples, M) % （ 改动 03－19) %samples=wavread('101.wav'); %M=[3,3,3,3]; K = length(samples); %语音样本数 N = l…

function [hmm, pout] = train(samples, M) %输入: % samples -- 样本结构 % M -- 为每个状态指定pdf个数,如:[3 3 3 3] %输出: % hmm -- 训…

load samples; % tic, for i=1:length(samples) sample=[]; for k=1:length(samples{i}) sample(k).wave=samples{i}{k}; …