formants.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 32 行

%formants.m

function [formc]=formants(area,leng)
% arrange the area function and section length, call bform to calculate
% the lossless formants, then correct the formants with loss

%	clear all;
%	load arle.mat;
	form=zeros(1,5);
	aa=zeros(1,60);
	xx=zeros(1,60);
	Inn = 60;
	Nform = 4;

	for j = 1:Inn,
		aa(1,j) = area(1,Inn+1-j);
		xx(1,j) = leng(1,Inn+1-j);
	end

	for j = 1:Inn,
		if (aa(1,j) <= 0.0001)
			aa(1,j) = 0.0001;
		end
	end

% call bform to calculate the lossless formants */
	form=bform(aa, xx, Inn);
% correction for loss */
	for j = 1:Nform,
		formc(j) = sqrt(form(j)*form(j) + 36100);
	end