计算机数值方法设计报告.docVIP

《计算机数值方法》 课程设计报告 设计题目：曲线拟合问题的数值试验 学号：XXXXX 姓名：X 电话：1599624XXXX 指导老师：徐元 于2010年11月27日 1、设计准备： 在11月20日，开始查阅有关书籍，学习matlab，在此推荐了一本参考书《MATLAB与计算方法实验》。注册MATLAB中文论坛账号，学习相关的matlab知识。 系统看了一遍插值与拟合的相关知识，掌握最小二乘的确切含义。 于11月24日开始设计。 2、设计题目： 9.2 曲线拟合问题的数值试验 1、对试验数据 （） i=1,2,….,N 利用最小二乘法原理，编制确定拟合多项式 =x+ nN 的数值软件，要求当0时，软件能对拟合多项式中的任意项数=0，能自动调整相应的正则方程组。 现在就选j=2； x=[5 10 15 20 25 30 35 40 45 50]; y=[1.0029 1.0023 1.0000 0.9990 0.9983 0.9979 0.9978 0.9981 0.9987 0.9996]; N=4; B=[10 sum(x) sum(x.^2) sum(x.^3);sum(x) sum(x.^2) sum(x.^3) sum(x.^4); sum(x.^2) sum(x.^3) sum(x.^4) sum(x.^5);sum(x.^3) sum(x.^4) sum(x.^5) sum(x.^6)]; Y=[sum(y);sum(x.^y);sum((x.^2).*y);sum((x.^3).*y)]; j=2; B(j,:)=[]; B(:,j)=[]; Y(j,:)=[]; A=B*B; C=B*Y; X=inv(A)*C; xi=5:5:50; yi=X(1)+X(2)*xi.^2+X(3)*xi.^3; plot(xi,yi,g,x,y,o); xlabel(x); ylabel(y); legend(?-ê?êy?Y,?ao?);（乱码注释：原始数据，拟合） 根据此程序段，已自动调整相应的正则方程组； 这是形如y= 以下是拟合后的图：1.0026 =-1.1563e-05 =2.10479e-07 对比的是：不缺省a1的完全多项式拟合曲线fig2 已知试验数据 x=[5 10 15 20 25 30 35 40 45 50] y=[1.0029 1.0023 1.0000 0.9990 0.9983 0.9979 0.9978 0.9981 0.9987 0.9996] 由于此文件编辑成GUI，M文件代码会比较长，选取其中几段有价值的部分进行叙述： 以下是选取的对上述数据拟合的M文件的部分代码： GUI入口函数 这段是GUI入口的function，将出现的变量x,y,n全部在此赋初值定义。 比如对n赋值2，进行默认进入软件界面时阶数是2。 handles.output = hObject;%选择默认命令行输出 Fig3（x,y,n赋初值图) 2）polyfit按钮下的函数 P=polyfit(x,y,n)是拟合求系数的，ymean是求y的平均值，从而求到偏差平方Rsq，此下均在GUI中显现，故需设定x,y,n。 Strcat()是将参数的结果拷贝至对应的字符下，比如strcat(p=,num2str(p)))就是将求的的p先数组化，后粘贴至’p=’下。 norm() 函数是对其求范数。 Fig4 现在，解释下这题的GUI界面制作：fig5 如上图所示,对x,y,n,系数、偏差平方变量分别用edit，画曲线的位置放置了一axes，polyfit、close全用按钮，x,y,n,系数、偏差平方5个变量后对应的用panel，从而完成软件的制作。 2、拟合形如的函数的一种快速方法是将最小二乘法用于下列问题： (1+cx)f(x)=a+bx 试用者方法拟合表9-2给出的这组数据。 表9-2 世界人口数据 年 人口/亿 1000 3.40 1650 5.45 1800 9.07 1900 16.10 1950 25.09 1970 36.50 1）验证如用此拟合公式，在2010年8月30日这个星期一的下午5时至6时之间，世界人口将变的无穷大。 利用上述程序段，求得的人口变为inf时的时期是： 更精确的值是：2.010737337381618e+03，利用计算机处理折合成2010年8月26日凌晨3—4时，世界人口变得无穷大。 图上X标示在2011年（大约），在此横坐标下 可以算作在误差范围内1）