第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) % 根据多项式系数构造多项式 【例7-2】已知一组测量数据如下表所示。 给定一组拟合基函数y=1,y=x,y=x2,y=cosx,y=ex,y=sinx，试求其最小二乘拟合函数。 解：编写如下程序： x=0:0.5:3; % x轴数据 y=[0 0.4794 0.8415 0.9815 0.9126 0.5985 0.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) % 最小二乘拟合求解 syms x digits(4) % 设定精度 Phifun=[1;x;x.^2;cos(x);exp(x);sin(x)]; y=vpa(p*Phifun) % 最小二乘拟合解函数 【例7-6】多项式拟合示例。 利用1到8次多项式拟合下表中的数据。 编写如下程序： x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.07 7.34 7.66 9.56 9.48 9.30 11.2]; xi=linspace(0,1); for k=1:8 p=polyfit(x,y,k); % 多项式拟合 yp=polyval(p,xi); % 求节点处的拟合值 plot(x,y,ro,xi,yp,k) % 绘制图形 for i=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) % 添加标题 end pause(1) % 暂停1秒 end 【例7-7】一元非线性曲线拟合。 解：根据上述公式编写如下程序： t=0:5:60; y=[0.2 0.4 0.5 0.9 1.5 2.4 3.1 3.8 4.1 4.2 4.5 4.4 4.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)) % 绘制图形 hold on % 图形保持 plot([39.3,38.5,38.5,37.7,38.5,38.5,39.3],... [2.1,2.05,1.55,1.5,1.45,0.95,0.9],k) % 绘制左括弧 h(1)=text(40,2,[{\ita}=,num2str(x(1))]); % 添加标注 h(2)=text(40,1.5,[{\itb}=,num2str(x(2))]); % 添加标注 h(3)=text(40,1,[{\itc}=,num2str(x(3))]); % 添加标注 set(h,fontname,times,fontsize,12) % 设置标注的字号与字形 【例7-8】多元非线性拟合。 编写如下程序： load reaction beta = nlinfit(reactants,rate,@hougen,beta) 另外上述结果也可以由界面操作得到： nlintool(reactants,rate,@hougen,beta,0.01, {氢气;n-戊烷;异戊烷},反应速率 ) 曲线拟合工具箱操作 Step 1：输入观测点数据 t=0:5:60; y=[0.2 0.4 0.