《数学实验》报告
实验名称 matlab作业 6 学 院 专业班级 姓 名 学 号
2011年 11月
一、 【实验目的】
学会运用MATELAB进行微积分求解和进行概率论计算的有关运算.
二、
【实验任务】
（一） P168 习题24
求解微分方程y??xsinx。 cosy（二） P168 习题27
用数值方法求解下列微分方程，用不同颜色和线形将y与y?画在同一个图形窗口里： y???ty??y?1?2t
初始时间：t0?0；终止时间：tf??；初始条件：y|t?0?0.1 y?|t?0?0.2
（三） P190 习题15
描绘以下数组的频率直方图：
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
（四） P190 习题17
若样本为85,86,78,90,96,82,80,74，求样本均值、标准差、中位数、极差和方差。 三、【实验程序】 (一) (二)
>> y=dsolve('Dy=(x*sinx)/cosy')
2
dsolve('D2y+t*Dy-y=1-2*t','y(0)=0.1,Dy(0)=0.2','t')
function xdot=exf(t,x) u=1-2*t;
xdot=[0,1;1,-t]*x+[0,1]'*u; clf; t0=0; tf=2*pi; x0t=[0;0];
[t,x]=ode23('exf',[t0,tf],x0t) y=x(:,1); dy=x(:,1); u=3*cos(t);
plot(t,y,'-',t,u,'bp',t,dy,'o') legend('êy?μ?y·??a','ê?è?á?','?a???a')
(三)
load p190_15.txt hist(p190_15,10)
(四)
load P190_16.txt
M=[mean(P190_16),median(P190_16),range(P190_16),var(P190_16),std(P190_16)]
四、【实验结果】 （一）
y =
x*sinx/cosy*t+C1
（二） t =
3
0 0.0001 0.0005 0.0025 0.0125 0.0355 0.0586 0.0857 0.1203 0.1633 0.2153 0.2768 0.3483 0.4302 0.5226 0.6258 0.7400 0.8651 1.0012 1.1135 1.2223
4
1.3615 1.5313 1.7182 1.9546 2.2277 2.5328 2.8562 3.1894 3.5308 3.8801 4.2380 4.6047 4.9806 5.3666 5.7643 6.1757 6.2832 x =
0
0
5