datafit3.asv

语音信号子带编码源程序

function  datafit3
%《语音信号子带编码》
% 对不同的子带数，选择相似的字长分派方式进行传输处理
% 产生的误差进行统计分析。并且可以选用不同的参数拟
% 和，有多项式拟和和非线性拟和。
close all;
w0=[2 4 8 16];
p0=[0.0149 0.0427 0.2273 0.2337];
plot(w0,p0,'r*');hold on;
%%% 非线性拟合 %%%
x=[2 4 8 16];y=p0;
a0=[-1 0.5];x0=0.0014;
options=optimset('fminsearch');
options.TolX=0.001;
options.Display='off';
a=fminsearch(@logg2,a0,options,x,y);
i=1:0.1:16;
p1=a(1)./(1+(a(1)/x0-1)*exp(-a(2).*i));
% plot(i,p1)

plot(i,p1,'LineWidth',1.8);axis([1 16 0 0.25])
ylabel('\fontname{隶书}\fontsize{12}相对失真'),xlabel('\fontname{隶书}\fontsize{12}子带数/个'),
title('\fontname{隶书}\fontsize{14}非线性数据拟和曲线');
text(8,.17,'\fontname{隶书}\fontsize{12}拟合函数的表达式：');
ch=['y=',num2str(a(1)),'/[1+[',num2str(a(1)),'/',num2str(x0),'-1]*exp(',num2str(a(2)),'*x)]']
text(8,130,ch);
text(8,.13,'\fontname{隶书}\fontsize{12}产生的相对失真：');
tmin=num2str(min(p1));syms l;
tmax=num2str(double(limit(a(1)./(1+(a(1)/x0-1)*exp(-a(2)*l)),inf)));
text(9,.115,'\fontname{隶书}\fontsize{12}最大值为Dmax:');text(13,.115,tmax);
text(9,.10,'\fontname{隶书}\fontsize{12}最小值为Dmin:');text(13,.10,tmin);