连续函数最佳平方逼近的算法程序.docVIP

一. 任务： 用MATLAB语言编写连续函数最佳平方逼近的算法程序（函数式M文件）。并用此程序进行数值试验，写出计算实习报告。 二. 程序功能要求： 在书中Page355或Page345的程序leastp.m（见附一）的基础上进行修改，使其更加完善。要求算法程序可以适应不同的具体函数，具有一定的通用性。所编程序具有以下功能： 用Lengendre多项式做基，并适合于构造任意次数的最佳平方逼近多项式。 可利用递推关系 被逼近函数f(x)不用内联函数构造，而改用M文件建立数学函数。这样，此程序可通过修改建立数学函数的M文件以适用不同的被逼近函数（要学会用函数句柄）。 要考虑一般的情况。因此，程序中要有变量代换的功能。 计算组合系数时，计算函数的积分采用变步长复化梯形求积法（见附三）。 程序中应包括帮助文本和必要的注释语句。另外，程序中也要有必要的反馈信息。 程序输入：（1）待求的被逼近函数值的数据点(可以是一个数值或向量) （2）区间端点：a,b。 7. 程序输出：（1）拟合系数： （2）待求的被逼近函数值 三：数值试验要求： 试验函数：；也可自选其它的试验函数。 用所编程序直接进行计算，检测程序的正确性，并理解算法。 分别求二次、三次、。。。最佳平方逼近函数。 分别作出逼近函数和被逼近函数的曲线图进行比较。 （分别用绘图函数plot(,s())和fplot(‘xcosx’,[x1 x2,y1,y2])） 1．本次实验程序如下： Lengendre多项式的自动生成M文件 function P=legendre(N) % evaluate legendre polynomial up to nth order syms t x; for n=1:N PP(n)=diff((t^2-1)^(n-1),n-1); Q(n)=2^(n-1)*prod([1:n-1]); end PP(1)=1; Q=sym(Q); P=PP*(inv(diag(Q))); 二、被逼近函数用M文件建立 function F=creat_fun(x) %creat function use M-file n=length(x); F=x.*cos(x(1:n)); 三、区间变换 function f=convert_interval(a,b,F) %conversion interval [a,b] between [-1,1],and expression of function syms x t; s=2\((b-a)*t+a+b);%区间转换 f=subs(F,x,s);%符号替换 四、变步长复化梯形求积公式的算法 function I=Rctrap(g) %求积分的变步长梯形算法 %f为被积函数； %积分的上下限为【-1,1】； %n为复花区间个数 %T为输出的变步长梯形公式计算值； m=1; h=1-(-1);%所有积分都区间都已转化为【-1,1】； T=zeros(1,100); T(1)=h*(feval(g,-1)+feval(g,1))/2; i=1; %进行变步长计算 while i100 m=2*m; h=h/2; s=0; for k=1:m/2 x=-1+h*(2*k-1); s=s+feval(g,x); end T(i+1)=T(i)/2+h*s; if abs(T(i+1)-T(i))0.00001%满足条件输出停机 I=T(i+1);%把终值返回 break; end i=i+1; end 五、主程序 function [c s]=least_square(a,b,N) %least square fit using Legendre polynomials,s(x)=c1p1(x)+c2p2(x)+c3p3(x)+...; %[a,b] initial interval; %c is coefficient; %N is the degree of polynomial %s is the wanted function syms t x; %sub-function to create function F=creat_fun(x); %sub-function to evaluate legendre polynomial up to nth order P=legendre(N); %sub-function:conversion interval [a,b] between [-1,1],and expression