UPF.docVIP

UPF clear; %load data_lpcc_lpc3_b22;%载入LPCC系数 %load data_lpcc_lpc3_dress6; %load data_lpcc_lpc3_dress; %load data_lpcc_lpc3_v4bie; %load data_lpcc_lpc3_v4ben; %load data_lpcc_lpc3_b22; %load data_lpcc_lpc3_b2; load data_lpcc_lpc3_szdx; %load data_lpcc_lpc3_sjj_com_4; %load data_lpcc_lpc3_sjj_vol_com; %load data_lpcc_lpc3_dress_r3; %load data_lpcc_lpc3_dress_r3_1; %load data_lpcc_lpc3_speaker6; %load data_lpcc_lpc3_speaker6_1; %load data_lpcc_lpc3_wyb_com_2_1; %load data_lpcc_lpc3_pt; [sig]=wavread(D:\语音库\实验语音\szdx.wav); %load data_fmt_lsf_b22_071219;%载入用LSF提取的共振峰参数 figure(10);set(gcf,Color,White); specgram(sig(:,1),1024,fs,40,16); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% dimension = 4; %the dimension of the formant vector.(向量共振峰的维数) d2 = dimension * 2; N = 100; % Number of particles.(颗粒的数目) sigma = 1e-4; % Variance of the Gaussian measurement noise.（高斯测量噪声） g1 = 3; % Paramater of Gamma transition prior.（G变换前的参数） g2 = 2; % Parameter of Gamman transition prior. % Thus mean = 3/2 and var = 3/4. T = franum; % Number of time steps.（时间步的数目） R = 1e-5; % EKFs measurement noise variance.（扩展卡尔曼滤波的测量噪声方差） Q = 3/4 * eye(d2); % EKFs process noise variance.（扩展卡尔曼滤波的进程噪声方差） P0 = 3/4*400/6; % EKFs initial variance of the states.（扩展卡尔曼滤波的初始状态方差） resamplingScheme = 1; % The possible choices are % systematic sampling (2), （可能的选择是系统抽样（2），剩余（1）和 % residual (1) 多项式(3)。他们都是O（N）算法） % and multinomial (3). % Theyre all O(N) algorithms. y = lpcc_array_1(:,2:franum+1); %Q_pfukf = 8*3/4 * eye(6); R_pfukf = 1e-1 * eye(16) ; miu(:,1:T)=zeros(16,1,T); Qv = 1e-5 * eye(16) ; alpha = 1; % UKF : point scaling parameter（点缩放参数） beta = 0; % UKF : scaling para