北科大 matlab大作业,北科大_matlab第六次作业_常微分方程.doc

北科大_matlab第六次作业_常微分方程

《数学实验》报告

实验名称? 常微分方程

学 院 计算机与通信工程学院

专业班级??? ??

姓?? 名?????

学?? 号??? ? ??

??

2015年 6月

【实验目的】

1、学习常微分方程符号解法和数值解法；

2、学习。

【实验任务】

1.P168 T24

求解微分方程。

2.P168 T27

用数值方法求解下列微分方程，用不同颜色和线型将y和y’画在同一个图形窗口里：

初始时间：；终止时间：；初始条件： ，。

3.P190 T15

描绘以下数组的频数直方图：

6.8,29.6,33.6,35.7,36.9,45.2,54.8,65.8,43.4,53.8,63.7,69.9,70.7,79.5,97.9,139.4,157.0

【实验程序】

1.P168 T24

y=dsolve('Dy=x*sin(x)/cos(y)','x')

2.P168 T27

先将高阶微分方程转化为一阶微分方程，令，即原微分方程化为：

写成矩阵形式为：

放入函数exf.m中，命令如下：

[t,x]=ode23('exf',[t0,tf],x0t)

将导数表达式的右端写成exf.m函数文件：

function [xdot ] = exf( t,x )

u=1-2*t;

xdot=[0 1;1 -t]*x+[0 1]'*u;

end

主程序如下：

clf

t0=0,tf=pi,x0t=[0.1,0.2];

[t,x]=ode23('exf',[t0,tf],x0t);

y=x(:,1);y1=x(:,2);

plot(t,y,'r-',t,y1,'b--');

legend('y','y的一阶导数');

3.P190 T15

A=[6.8,29.6,33.6,35.7,36.9,45.2,54.8,65.8,43.4,53.8,63.7,69.9,70.7,79.5,97.9,139.4,157.0]

hist(A,12)

【实验结果】

1.P168 T24

y =

asin(C3 + sin(x) - x*cos(x))

2.P168 T27

3.P190 T15

【实验总结】

1、常微分方程符号解法和数值解法；

2、学。

4