数学建模实验报告-实验三

数学建模实验报告

              题目：微分方程与插值方法

院系：计算机科学与工程学院

班级：

姓名：

学号：

【实验目的】

1.掌握matlab软件求解常微分方程解析解的方法；

2.掌握matlab软件求解一维插值问题的方法；

3.掌握matlab软件求解二维插值方法。

【实验内容与实现】

1. 解下面的二阶常微分方程：

实验代码：

y=dsolve('D2y-4*Dy+4*y=0','y(0)=4,Dy(0)=1')

实验结果：

2. （二阶变系数）解下面的一阶常微分方程：

实验代码：

y=dsolve('D2y-sin(x)*Dy+y=cos(x)','y(0)=0,Dy(0)=0')

实验结果：

 y =

cos(x) + (exp(t*(sin(x)/2 + ((sin(x) - 2)*(sin(x) + 2))^(1/2)/2))*cos(x)*(sin(x) - ((sin(x) - 2)*(sin(x) + 2))^(1/2)))/(2*((sin(x) - 2)*(sin(x) + 2))^(1/2)) - (exp(t*(sin(x)/2 - ((sin(x) - 2)*(sin(x) + 2))^(1/2)/2))*cos(x)*(sin(x) + ((sin(x) - 2)*(sin(x) + 2))^(1/2)))/(2*((sin(x) - 2)*(sin(x) + 2))^(1/2))

3. 已知数表：

x	10°	11°	12°	13°	14°
sinx	0.1736482	0.190809	0.2079117	0.2249511	0.2419219

试用最邻近插值、分段线性插值、分段三次样条插值和分段三次插值计算sin(12°54')的近似值，并在区间[10,14]上绘制sinx的图像以及插值法得到的图像。要求绘制4张图片，每张图片上分别绘制原始图像和一种插值图像。

实验代码：

x0=linspace(10,14,5);