怎样计算平面图形的面积 .docxVIP

《数学实验》报告学 院：**************专业班级：************学 号：*************姓 名：**********实验名称：怎样计算平面图形的面积实验日期：********实验四怎样计算平面图形的面积目的：探索曲线拟合的不同方式，使学生了解泰勒公式的意义，并且对运用定积分计算任意平面图形的面积有更深入的认识。能初步运用所学数学知识及数学软件工具matlab解决实际问题。任务1．在同一坐标系内作出区间上对数函数及多项式函数的图像，观察这些多项式函数逼近指数函数的情况。2．在同一坐标系内作出区间上正弦函数及多项式函数3．假定某天的天气变化，试找出这一天的气温变化规律（求气温关于时间的函数）。时刻t(x)0123456789101112温度oC(y)15141414141516182022232528时刻t(x)131415161718192021222324温度oC(y)3132312927252422201817164．比较这两种方法各自的优点和缺点。对于给定的数据，运用Matlab及以上两种不同的方法找出相应的多项式，并作出其图像。你认为哪种方法要好些。你能不能找到其他的办法解决这个问题。5．在平面上任意画一条简单光滑闭曲线，求此闭曲线围成的平面图形的面积。并观察当曲线上取的点逐渐增加时，所求多项式及所求面积值的变化情况。实验过程在同一坐标系内作出区间上对数函数及多项式函数的图像，观察这些多项式函数逼近指数函数的情况。编制程序绘制出它的图形；【思路】使用plot函数在同一坐标中绘出全部曲线，使用legend函数对全部曲线做出注释【程序】x=(-1:0.001:1);y1=log(1+x) ;y2=x ;y3=x-x.^2/2 ;y4=x-x.^2/2+x.^3/3;y5=x-x.^2/2+x.^3/3-x.^4/4;plot (x , y1 , k , x , y2 , m , x , y3 , g , x , y4 , c,x,y5,r);legend(ln(1+x),x,x-x^2/2,x-x^2/2+x^3/3,x-x^2/2+x^3/3,x-x^2/2+x^3/3-x^4/4);【运行结果】【分析】随着多项式级数的增加，多项式逼近拟合的效果越佳2．在同一坐标系内作出区间上正弦函数及多项式函数【思路】使用plot函数在同一坐标中绘出全部曲线，使用legend函数对全部曲线做出注释【程序】x=(-pi:0.001:pi);y1=sin(x) ;y2=x-x.^3/factorial(3) ;y3=x-x.^3/factorial(3)+x.^5/factorial(5) ;y4=x-x.^3/factorial(3)+x.^5/factorial(5)-x.^7/factorial(7);y5=x-x.^3/factorial(3)+x.^5/factorial(5)-x.^7/factorial(7)+x.^9/factorial(9);plot (x , y1 , k , x , y2 , m , x , y3 , g , x , y4 , c,x,y5,r);legend(sin(x),x-x^3/3!,x-x^3/3!+x^5/5! ,x-x^3/3!+x^5/5!-x^7/7!,x-x^3/3!+x^5/5!-x^7/7!+x^9/9!);【运行结果】　　　【分析】随着多项式级数的增加，多项式逼近拟合的效果越佳3．假定某天的天气变化，试找出这一天的气温变化规律（求气温关于时间的函数）。时刻t(x)0123456789101112温度oC(y)15141414141516182022232528时刻t(x)131415161718192021222324温度oC(y)313231292725242220181716（方法一）拉格朗日插值法【程序】%原始数据-下边界曲线X1=[ 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24 ; 15 ,14 ,14 ,14,14 ,15,16,18,20,22,23,25,28,31,32,31,29,27,25,24,22,20,18,17,16 ] ;t1=X1(1, : ) ; %横坐标y1=X1(2 , : ) ; %竖坐标for i=1:length(X1) %计算范德蒙矩阵b( i , 1 ) = 1 ;for j=2:length(X1)b( i , j ) = X1 ( 1 , i ) .^ ( j - 1 ) ;end ;enda1=inv( b ) * X1( 2 ,:); %解方程组for i=1 : 25%调整h