matlab求sin60度,matlab数学实验作业

深圳大学考试答题纸

(以论文、报告等形式考核专用)

二○一 ○ ～二○一 一 学年度第 一 学期

课程编号 23160008

课程名称 Matlab 与数学实验

主讲教师

仇高新

评分

学 号 2008160120 姓名 黄日才 专业年级

电子学院微电子二班

题目：

《Matlab 与数学实验》期末作业

1. 用数值积分公式计算 (结果保留小数点后8位)：

θ

d θsin 15.0120

22?

-=π

S

(1) 取积分步长2/π=h , 用梯形公式计算S= 6.247641317417333 。

>> x=0:pi/2:2*pi;

>> trapz(0:pi/2:2*pi, sqrt(1-(0.15.*sin(x)).^2))

(2) 要求相对误差为10-6, 用Simpson 公式S= 6.247691887569109 ，Matlab 命令是__

quad('sqrt(1-(0.15.*sin(x)).^2)',0,2*pi,1e-6)_.

2. 设,0)0(,1)0(,0sin )()(='==-''y y x x y x y 用数值解法算出 y (1)= 1.163536246222507 ,你用的方法是 Runga-Kutta 方法 ，

调用的 Matlab 命令[t,x]=ode45(@verderpol2,[0:0.1:1.5], [1;0]) ， 算法精度为 4阶 。

解：先编写函数文件verderpol2.m ： function xprime=verderpol2(t,x) xprime=[x(2); x(1)*sin(t)];

调用命令：

[t,x]=ode45(@verderpol2,[0:0.1:1.5], [1;0])；

3. 设,0)0(,1)0(,0sin )()(='==+'-''y y ye x x y x y x

用数值解法算出y (1)= 0.2714 (精确到4位小数), 你用的方法是 Runga-Kutta 方法 ，调用的 Matlab 命令是 [t,x]=ode45(@verderpol3, [0:0.1:1.5], [1; 0]) ，算法精度为 4-5阶 。

解：先编写函数文件verderpol3.m ： function xprime=verderpol3(t,x)

xprime=[x(2); x(2)*sin(t)-x(1)*exp(t)];

调用命令：[t,x]=ode45(@verderpol3, [0:0.1:1.5], [1; 0])