《工程数值分析报告实验龙格库塔最小二乘法》由会员分享，可在线阅读，更多相关《工程数值分析报告实验龙格库塔最小二乘法（10页珍藏版）》请在金锄头文库上搜索。
1、工程数值分析实验报告指导老师班级学号姓名实验一：最小二乘法拟合曲线实验一、实验名称： 最小二乘法拟合曲线实验实验时间：2015-5-14实验地点：主楼机房实验器材： 计算机matlab二、实验目的：学会用最小二乘法求拟合数据的多项式，并应用算法于实际问题。三、实验要求：（1）根据最小二乘法和加权最小二乘法的基本理论，编写程序构造拟合曲线的法方程，要求 可以方便的调整拟合多项式的次数；（2）采用列主元法解（1）中构造的法方程，给出所拟合的多项式表达式；（3）编写程序计算所拟合多项式的均方误差，并作出离散函数 和拟合函数的图形；（4）用MATLAB的内部函数polyfit求解上面最小二乘法曲线拟合
2、多项式的系数及平方误差， 并用MATLAB 的内部函数plot作出其图形，并与（1）的结果进行比较。四、算法描述（实验原理与基础理论）基本原理：从整体上考虑近似函数 同所给数据点（=，yi）（i=0,1,，m）误差r = P（二.）乙（i=0,1,，m） 的大小，常用的方法有以下三种：一是误差 二p（x,） y（. （i=0,1,，m）绝 对值的最大值m瓷匕，即误差向量r=（r0,人rm）t的-范数；二是误差绝对值的和 才 lr 丨、Er 2ci，即误差向量r的1 -范数；三是误差平方和i 的算术平方根，即误差向量r的2-范数；前两种方法简单、自然，但不便于微分运算，后一种方法相当于考虑 2-
3、范数的平乙 r 2方，因此在曲线拟合中常采用误差平方和i=0 i 来 度量误差 ri （i=0, 1, m）的整体大 小。五、实验内容：共有两组给定数据，把给定的数据拟合成多项式。第一组给定数据点如表1所示如下：表1数据表x .100. 50. 60. 70. 80. 91. 0yi11. 751. 962. 192. 442. 713. 00表2数据表x . i00. 50. 60. 70. 80. 91. 0yi11. 751. 962. 192. 442. 713. 00i1236421六、程序流程图开始输入拟合次数N,X,Y的坐标1矩阵A计算X,Y平均值输出获得系数a0,a1获得A的逆
4、矩阵B获得矩阵C解方程组计算相关系数r七、实验结果 zuixiaoerchenfa ans =27-May-2015 ans =7.3611e+05 ans =1.0e+03 *2.0150 0.0050 0.0270 0.0140 0.0010 0.0213xFigure 1Figure 2xFigure 3xFigure 4八、实验结果分析实验程序quxiannihe.m clear all date,now,clockx0=0.0 0.5 0.6 0.7 0.8 0.9 1.0;y0=1 1.75 1.96 2.19 2.44 2.71 3.00;w=ones(size(x0);x=0:
5、0.01:1;%进行五次曲线拟合N=5;for i=1:Na1=LSF(x0,y0,w,i) ;y=polyval(a1,x);figure(i)plot(x0,y0,ok,x,y,r)title（最小二乘法）;legend（，yO,y）;xlabel（x）;ylabel（，y，）;end实验二：4阶经典龙格库塔法解常微分方程一、实验名称：4阶经典龙格库塔法解常微分方程实验时间：2015-5-14实验地点：主楼机房实验器材： 计算机matlab二、实验目的： 学习掌握4阶经典R-K方法，体会参数和步长对问题的影响。三、实验要求：（1）用4阶经典R-K法编写计算程序，要求用法与ode45 一致。
6、并将计算结果画图比较， 并分析步长变化对解的影响。（2）当激励力幅值F 分别按 0.3, 0.33, 0.4， 0.43, 0.54, 0.58, 0.75， 0.84，11.21， 13.34进行计算。每一个数据画出三幅图，分别为时间位移曲线，时间速度曲线和相图。考 察激励力幅值 F 变化引起的系统响应的变化。（3）请采用MATLAB中的内部库函数ode45求解此常微分方程初值问题的解，并与（1） 中的结果进行比较。四、算法描述（实验原理与基础理论）系统方程和表述如下：2他认谢d二则系统的输出按如下求解：Vn+1 = Vn+ （fcl + 2魅 + 2ft3 + fci）其中：h = f 仏
7、、4）这样，下一个值(y )由现在的值(y )加上时间间隔(h)和一个估算的斜率的乘积决定。该斜率 n+1n是以下斜率的加权平均：k是时间段开始时的斜率；k2是时间段中点的斜率，通过欧拉法采用斜率k1来决定y在点tn +h/2的值；k3也是中点的斜率，但是这次采用斜率k2决定y值；k4是时间段终点的斜率，其y值用k3决定。五、实验内容：求解常微分方程初值问题，考虑著名的Duffing方程。GDuffing 在1918 年引入了一个带有立方项的非线性振子来描述出现在许多力学问题中的质量、弹簧、阻尼系统。从那时起，Duffing方程在非线性动力学系统的研究中占有重要 的地位。Duffing方程的标
8、准形式是d2 xdt 2+ c+ f (x) = g (t)其中：f(x)是一个含有三次项的非线性函数，g(t)是一个周期函数。把 f (x)二 x - x3， g (t) = F cos( t)代入上式，可得1)d2xdx+ c+ x 一 x3 = F cos(w t)dt2dt式中：c = 0.301,= 1.201, 步长h = 0.01 ;初值向量为：x0=(0, 0.1)。要求考察激励力幅 值 F 变化引起的系统响应的变化。积分时间区间为：0, 50。六、程序流程图开始输入 a,b,n,x Q a,y QyOk=n结束x U x + h2 了+2 KK U f (x, y)1 y h
9、K U f (x + -, y + K )2 k U , K U f (x + y + ：K )3 21K u f( x+屯yy u y + - (K + 2K + 2K + K )61234k u k +1输出k,x,y七、实验结果i x(i) y(i)1 0.0000 1.00002 0.1000 0.99013 0.2000 0.96154 0.3000 0.91745 0.4000 0.86216 0.5000 0.80007 0.6000 0.73538 0.7000 0.67119 0.8000 0.609810 0.9000 0.552511 L 0000P500012 L 10
10、00 0- 452513 L 2000 0- 4098 flmcr+ion varargour+Hsaxplaxliu (varargin)clrclearX0H0 - xnnL 2 - yonl - hnp1 -y. xvlgs4j (xo. xp yo“ h)- nnlengr+h (X)- fprinr+f (. i x(i) y(i)n.)-for inlm fprinr+f2d 手4 4f 手 4 4f n “ - x (i) “ y (i)-endflmcr+ion znfFy) zh2*x42funcr+iony. xvlgs4j (xo. xp yo. h)XHXOhxF nnlengr+h (X)-yTX-ylaTyo-forKlf(x(i)yl(i)-K2nf (X (i) +h/% yl (i) +h紫KI)-K3nf (X (i) +h/% yl (i) +h、2奠2)-K4nf (x(i)+h. yl (i)+h*K3)- yl(i+1) Hyl(i) +h、6*(Kl+2奠2+2其3+K4)-endynyr
