数值逼近上机作业3.doc

一、sin(pi*x)的2次最佳平方逼近程序代码： %% 二次逼近主函数 function S=squar_approx(a,b,n) % 全局变量 global i;global j; % 判断输入 if nargin3 n=1; end % 求方程系数矩阵（基函数内积矩阵） Phi2=zeros(n+1); for i=0:n for j=0:n; Phi2(n-i+1,n-j+1)=quad(@rho_phi,a,b); end end % 求方程非齐次项列向量 PhiF=zeros(n+1,1); for i=0:n PhiF(n-i+1)=quad(@fun_phi,a,b); end % 解方程获得多项式系数 S=Phi2\PhiF; % 给出绘图节点 xi=a:(b-a)/100:b; % 作图 plot(xi,polyval(S,xi),r:); hold on fplot(@obj,[a,b],b); legend(逼近多项式,原函数); title(sin\pix的二次平方逼近); hold off %% 基函数内积 function y=rho_phi(x) global i;global j; y=(rho(x).*phi_k(x,i)).*phi_k(x,j); %% 基函数与待逼近函数内积 function y=fun_phi(x) global i; y=(rho(x).*phi_k(x,i)).*obj(x); %% 外部函数 function y=rho(x) y=1; %% 直交函数系 function y=phi_k(x,k) if k==0 y=ones(size(x)); else y=x.^k; end %% 待逼近函数 function y=obj(x) y=sin(pi*x); 运行结果： 二、1/(1+25x2)的2次最佳平方逼近，及其10次拉格朗日插值 程序代码： %% 二次逼近主函数 function S=squar_approx(a,b,n) % 全局变量 global i;global j; % 判断输入 if nargin3 n=1; end % 求方程系数矩阵（基函数内积矩阵） Phi2=zeros(n+1); for i=0:n for j=0:n; Phi2(n-i+1,n-j+1)=quad(@rho_phi,a,b); end end % 求方程非齐次项列向量 PhiF=zeros(n+1,1); for i=0:n PhiF(n-i+1)=quad(@fun_phi,a,b); end % 解方程获得多项式系数 S=Phi2\PhiF; % 给出绘图节点 xi=a:(b-a)/100:b; % 作图 plot(xi,polyval(S,xi),r:); hold on fplot(@obj,[a,b],b); x=a:(b-a)/10:b; y=obj(x); [C,~,~,~]=lagran1(x,y); plot(xi,polyval(C,xi),m:) legend(±??ü?à??ê?,?-oˉêy,10′?à-??àêè?2??μ); title(1/(1+25x^2)μ?10à-??àêè?2??μó?2′?×?????·?±??ü±è??); hold off %% 基函数内积 function y=rho_phi(x) global i;global j; y=(rho(x).*phi_k(x,i)).*phi_k(x,j); %% 基函数与待逼近函数内积 function y=fun_phi(x) global i; y=(rho(x).*phi_k(x,i)).*obj(x); %% 外部函数 function y=rho(x) y=1; %% 直交函数系 function y=phi_k(x,k) if k==0 y=ones(size(x)); else y=x.^k; end %% 待逼近函数 function y=obj(x) y=1./(1+25*x.^2); 运行结果：