《重庆大学 数学实验报告 机械二零一六年机自 实验三》.docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1401 学 院 机械工程学院 年级 09级 专业班 机自6 学 生 姓 名 学 号 2009 开 课 时 间 2010 至 2011学年第 2 学期 总 成 绩 教师签名 数 理 学 院 制 开课学院、实验室： 实验时间 ： 年 月 日 课程 名称 数学实验 实验项目 名 称 微分方程模型及其求解 实验项目类型 验证 演示 综合 设计 其他 指导 教师 何仁斌 成 绩 实验目的 5号宋体 基础实验 一、实验内容 5号宋体 二、实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等， 实验结果及分析） 5号宋体 应用实验（或综合实验） 实验内容 第一题： 1）通解 y=dsolve(Dy=y+2*x ,x) y = -2*x-2+exp(x)*C1 特解 y1=dsolve(Dy=y+2*x,y(0)=1,x) y1 = -2*x-2+3*exp(x) 画图 x=0:0.03:1; y1=-2*x-2+3*exp(x); plot(x,y1,r-*); 2） 通解 y=dsolve(D2y+y*cos(x) ,x) y=dsolve(D2y+y*cos(x) ,x) y = DESol({diff(_Y(x),$(x,2))+cos(x)*_Y(x)},{_Y(x)}) 特解 y1=dsolve(D2y+y*cos(x)=0 ,Dy(0)=0,y(0)=1,x) x=0:0.1:1; y =DESol({diff(_Y(x),$(x,2))+cos(x)*_Y(x)},{_Y(x)}); plot(x,y,’r-*’) 第二题 %随即生成一个11行，3列的矩阵，用于存放数据 a=rands(11,3) a = -0.0871 0.8338 -0.4556 -0.9630 -0.1795 -0.6024 0.6428 0.7873 -0.9695 -0.1106 -0.8842 0.4936 0.2309 -0.2943 -0.1098 0.5839 0.6263 0.8636 0.8436 -0.9803 -0.0680 0.4764 -0.7222 -0.1627 -0.6475 -0.5945 0.6924 -0.1886 -0.6026 0.0503 0.8709 0.2076 -0.5947 %矩阵的第一列存放X的值，先对其赋值 for i=1:11 a(i,1)=i-1 end a = 0 0.8338 -0.4556 1.0000 -0.1795 -0.6024 2.0000 0.7873 -0.9695 3.0000 -0.8842 0.4936 4.0000 -0.2943 -0.1098 5.0000 0.6263 0.8636 6.0000 -0.9803 -0.0680 7.0000 -0.7222 -0.1627 8.0000 -0.5945 0.6924 9.0000 -0.6026 0.0503 10.0000 0.2076 -0.5947 %将a(1,2)赋值为1 a(1,2)=1 a = 0 1.0000 -0.4556 1.0000 -0.1795 -0.6024 2.0000 0.7873 -0.9695 3.0000 -0.8842 0.4936 4.0000 -0.2943 -0.1098 5.0000 0.6263 0.8636 6.0000 -0.9803 -0.0680 7.0000 -0.7222 -0.1627 8.0000 -0.5945 0.6924 9.