数学实验1导论.docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1402 学 院 联合学院 年级 2014 专业班电气工程 学 生 姓 名 赖一多 学 号 开 课 时 间 2014 至 2015 学年第 二 学期 总 成 绩 教师签名 数 学 与 统 计 学 院 制 开课学院、实验室： 实验时间 ： 年 月 日 课程 名称 实验项目 名 称 实验项目类型 验证 演示 综合 设计 其他 指导 教师 成 绩 实验目的 5号宋体 实验内容 5号宋体 基础实验 一、问题重述 5号宋体 二、实验过程（一般应包括实验原理或问题分析，变量说明、程序、调试情况记录、图表等， 实验结果及分析） 5号宋体 应用实验（或综合实验） 一、问题重述 5号宋体 二、问题分析 5号宋体 三、数学模型的建立与求解(一般应包括模型、求解步骤或思路，程序放在后面的附录中) 5号宋体 四、实验结果及分析 5号宋体 五、附录（程序等） 5号宋体 总结与体会 设计记录表格，包括碰到的问题汇总及解决情况 注 行距：选最小值16磅，每一图应有简短确切的题名，连同图号置于图下。每一表应有简短确切的题名，连同表号置于表上。图表的题名及其中的文字采用小5号宋体。公式应该有编号，编号靠右端。 1．求微分方程的解析解, 并画出它们的图形， y’= y + 2x, y(0) = 1, 0x1; y’’+ycos(x) = 0, y(0)=1, y’(0)=0; Code: 1.1 clear y=dsolve(Dy=y+2*x,y(0)=1,x); ezplot(y,[0 1]) Result: 1.2 Code: clear [t y]=ode23(wf,[0 10],[0 1]); plot(t,y(:,1)) M文件 function yd = wf( t,y ) yd=[y(2);-y(1)*cos(t)] Result: 以下为部分数值解和图像 3．Rossler微分方程组： 当固定参数b=2, c=4时，试讨论随参数a由小到大变化（如a∈(0,0.65))而方程解的变化情况，并且画出空间曲线图形，观察空间曲线是否形成混沌状？ Code： clear global a for a=0:0.05:0.65; [t x]=ode45(P3,[0 10],[1 1 1]); plot(t,x) hold on end Mfunction equ=fun(t,x) global a b=2;c=4; equ=[0 -1 -1;1 a 0;x(3) 0 -c]*x+[0;0;b]; end Result： ，图像呈混沌状。 4.Apollo卫星的运动轨迹的绘制 Code：clear global u u1 u=1/82.45;u1=1-u; [t x]=ode45(P4,[0 10],[1.2 0 0 -1); plot(t,x) M文件： function eqt=P4(t,x) global u u1 r1=sqrt(((x(1)+u))^2+x(3)^2); r2=sqrt(((x(1)-u1)^2)+x(3)^2); eqt(1,1)=x(2); eqt(2,1)=2*x(4)+x(1)-u1*(x(1)+u)/r1^3-u*(x(1)-u1)/r2^3; eqt(3,1)=x(4); eqt(4,1)=-2*x(2)+x(3)-u1*x(3)/r1^3-u*x(3)/r2^3; Result： 5．盐水的混合问题 一个圆柱形的容器，内装350升的均匀混合的盐水溶液。如果纯水以每秒14升的速度从容器顶部流入，同时，容器内的混合的盐水以每秒10.5升的速度从容器底部流出。开始时，容器内盐的含量为7千克。求经过时间t后容器内盐的含量。 Code： clear [t x]=ode23(P5,0:0.1:10,[7 350]); count=1; for t=0:0.1:10 plot(t,x(count,1)/x(count,2),--);xlabel(ê±