matlab第7章实例解析.pptVIP

实例解析【例7-1】求函数在区间上的三次最佳平方逼近多项式。解：编写如下程序：wfun=x.^0;%权函数phifun=char(1,x,x.^2,x.^3);%基函数fun=exp(x);%被逼近函数a=-1;b=1;%区间端点[A,B,p]=square_approximation(wfun,phifun,fun,a,b)%求解三次最佳平方逼近多项式系数p2=poly2str(fliplr(p),x)%根据多项式系数构造多项式A=2.00000.00000.6667-0.00000.00000.666700.40000.666700.40000-0.00000.400000.2857B=2.35040.73580.87890.4495p2=0.17614x^3+0.53672x^2+0.99795x+0.99629【例7-2】已知一组测量数据如下表所示。给定一组拟合基函数y=1,y=x,y=x2,y=cosx,y=ex,y=sinx，试求其最小二乘拟合函数。解：编写如下程序：x=0:0.5:3;%x轴数据y=[00.47940.84150.98150.91260.59850.1645];%y轴数据wfun=ones(1,6);%权系数phifun=@(x)[ones(size(x));x;x.^2;cos(x);exp(x);sin(x)];%拟合基函数[A,b,p]=Least_square(wfun,phifun,x,y)%最小二乘拟合求解symsxdigits(4)%设定精度Phifun=[1;x;x.^2;cos(x);exp(x);sin(x)];y=vpa(p*Phifun)%最小二乘拟合解函数xi00.511.522.53yi00.47940.84150.98150.91260.59850.1645运行结果：

y=0.407*x-0.4598*cos(x)+0.0765*exp(x)+0.5653*sin(x)-0.3901*x^2+0.3828A=[ones(size(x))xx.^2cos(x)exp(x)sin(x)];p=A\y%最小二乘法求解【例7-6】多项式拟合示例。利用1到8次多项式拟合下表中的数据。编写如下程序：x=0:0.1:1;y=[-0.4471.9783.286.167.077.347.669.569.489.3011.2];xi=linspace(0,1);fork=1:8p=polyfit(x,y,k);%多项式拟合yp=polyval(p,xi);%求节点处的拟合值plot(x,y,ro,xi,yp,k)%绘制图形fori=1:length(p)digits(5)%设置精度h1=text(0.08*i-0.03,i-2,[p(,num2str(i),)=,char(vpa(p(i)))]);%添加标注h2=text(0.95,i-2,[{\itx}^,num2str(length(p)-i)]);%添加标注set([h1,h2],fontname,times,fontsize,12)%设置标注的字形与字号title([num2str(k),次多项式],fontname,隶书,fontsize,16)%添加标题endpause(1)%暂停1秒endexample_7_6.m【例7-7】一元非线性曲线拟合。解：根据上述公式编写如下程序：t=0:5:60;y=[0.20.40.50.91.52.43.13.84.14.24.54.44.5];f=@(b,x)1./(b(1)+b(2)*exp(b(3)*x));%定义拟合函数模型x=lsqcurvefit(f,[0.5,10,0],t,y)%最小二乘拟合求解plot(t,y,ko,t,f(x,t))