当前位置：文档之家› 最新重庆大学-数学实验-3微分方程

最新重庆大学-数学实验-3微分方程

重庆大学

学生实验报告

实验课程名称数学实验

开课实验室DS1422

学院

学生姓名

开课时间

数学与统计学院制

开课学院、实验室：数统学院DS1422 实验时间：

-2*x-2+3*exp(x)

作图（y =3*exp(x)-2*x–2）

x=0:0.02:1;

y1=3*exp(x)-2*x-2;

plot(x,y1,'b-')

（2）算法设计

通解y=dsolve('D2y+y*cos(x)=0')

y =

C19*exp(t*(-cos(x))^(1/2)) + C20/exp(t*(-cos(x))^(1/2))

2．用向前欧拉公式和改进的欧拉公式求方程y ’= y - 2x /y , y (0) = 1 (0≤x ≤1,h = 0.1) 的数值解，要求编写程序，并比较两种方法的计算结果，说明了什么问题？程序：

向前欧拉和向后欧拉法:

x1(1)=0;y1(1)=1;y2(1)=1;h=0.1;

for k=1;10

x1(k+1)=x1(k)+h;

y1(k+1)=y1(k)+h*(y1(k)-2*x1(k)/y1(k))

y2(k+1)=y2(k)+h*(y2(k+1)-2*x1(k+1)/y2(k+1)); end x1,y1,y2 ans = 10 y1 =

??-+=+=--=)('''c x z b z ay

x y z y x

-0.5

0.5

-1

-0.5

0.10.20.30.40.50.6

当a=0.2时

重庆大学数学模型数学实验作业四讲解

开课学院、实验室：数统学院实验时间：2015年10月28日课程名称数学实验实验项目名称种群数量的状态转移—— 微分方程实验项目类型验证演示综合设计其他指导教师肖剑成绩实验目的 [1] 归纳和学习求解常微分方程(组)的基本原理和方法； [2] 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； [3] 熟悉MATLAB软件关于微分方程求解的各种命令； [4] 通过范例学习建立微分方程方面的数学模型以及求解全过程；通过该实验的学习，使学生掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是十分必要的。实验内容 1．微分方程及方程组的解析求解法； 2．微分方程及方程组的数值求解法——欧拉、欧拉改进算法； 3．直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)； 4．利用图形对解的特征作定性分析； 5．建立微分方程方面的数学模型，并了解建立数学模型的全过程。基础实验一、问题重述 1．求微分方程的解析解, 并画出它们的图形， y’= y + 2x, y(0) = 1, 0

重庆大学数值分析试卷

重庆大学数值分析课程试卷 2012 ~2013 学年第 1学期开课学院：数统学院课程号：考试日期：考试方式：考试时间 120 分钟题号一二三四五六七八九十总分得分注：1.大标题用四号宋体、小标题及正文推荐用小四号宋体；2.按A4纸缩小打印一、选择题（3分/每小题，共15分） 1、以下误差公式不正确的是（ A ） A. ()()()1212x x x x εε ε- =- B. ()()()1212x x x x εεε+=+ C ．()()()122112x x x x x x εε ε = + D. ()()2 2 x x x εε = 2、通过点()0 0,x y ，()11,x y 的拉格朗日插值基函数()0l x ，()1l x 满足（C ） A. ()000l x =,()110l x = B. ()000l x =,()111l x = C. ()001l x =,()111l x = D. ()001l x =,()110l x = 3、已知等距节点的插值型求积公式 ()()3 52 k k k f x d x A f x =≈ ∑ ? ，则3 k k A == ∑ （ C ） A. 1 B. 2 C. 3 D. 4 4、解线性方程组A x b =的简单迭代格式() () 1k k x B x f +=+收敛的充要条件是（ B ） A. ()1A ρ< B. ()1B ρ< C. ()1A ρ> D. ()1B ρ> 5、已知差商021[,,]5 f x x x =，402[,,]9f x x x =，234[,,]14f x x x =，032[,,]8f x x x =，则 420[,,]f x x x = （ B ） A. 5 B. 9 C. 14 D. 8 二、填空题（3分/每小题，共15分） 1取 3.141592x =作为数 3.14159265 4...的近似值，则x 有____6____位有效数字 2、Cotes 求积公式的代数精度为 5 学院专业、班年级学号姓名公平竞争、诚实守信、严肃考纪、拒绝作弊封线密

数学软件MATLAB实验作业

数学软件与数学实验作业一.《数学软件》练习题（任选12题，其中19－24题至少选2题）： 3.对下列各式进行因式分解. (1). syms x y >> factor(x^5-x^3) (2). syms x y >> factor(x^4-y^4) (3). syms x >> factor(16-x^4) (4). syms x >> factor(x^3-6*x^2+11*x-6) (5). syms x y >> factor((x+y)^2-10*(x+y)+25) (6). syms x y >> factor(x^2/4+x*y+y^2) (7). syms x y a b >> factor(3*a*x+4*b*y+4*a*y+3*b*x) (8). syms x >> factor(x^4+4*x^3-19*x^2-46*x+120) 5.解下列方程或方程组. (1).solve('(y-3)^2-(y+3)^3=9*y*(1-2*y)') (2). solve('3*x^2+5*(2*x+1)') (3). solve('a*b*x^2+(a^4+b^4)*x+a^3*b^3','x') (4). solve('x^2-(2*m+1)*x+m^2+m','x') (5). [x,y]=solve('4*x^2-9*y^2=15','2*x-3*y=15') 6.计算极限. (1). syms x f=(exp(x)-exp(-x))/sin(x); limit(f,x,0) (2) syms x >> f=(x/(x-1)-1/log(x)); >> limit(f,x,1) (3). syms x >> f=(1-cos(x))/x^2; >> limit(f,x,0)

重庆大学数学实验方程模型及其求解算法参考答案

实验2 方程模型及其求解算法一、实验目的及意义 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程；通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。二、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。四、实验要求与任务基础实验 1．用图形放大法求解方程x sin(x) = 1. 并观察该方程有多少个根。画出图形程序： x=-10:0.01:10; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB运行结果：

-10-8-6-4-20246810 -8-6 -4 -2 2 4 6 8 扩大区间画图程序： x=-50:0.01:50; y=x.*sin(x)-1; y1=zeros(size(x)); plot(x,y,x,y1) MATLAB 运行结果： -50-40-30-20-1001020304050 由上图可知，该方程有偶数个无数的根。

数学实验作业

练习2﹒1 画出下列常见曲线的图形（其中a=1，b=2,c=3）。 1. 立方抛物线y = 解： x=-4:0.1:4; y=x.^(1/3); plot(x,y) -4 -3-2-101234 0.20.40.60.811.21.4 1.6 2.高斯曲线2 x y e -= 解： fplot('exp(-x^2)',[-4,4])

-4 -3 -2 -1 1 2 3 4 00.10.20.30.40.50.60.70.80.9 1 3、笛卡儿曲线23 3 2 2 33,(3)11at at x y x y axy t t = = +=++ 解：ezplot('x^3+y^3-3*x*y',[-4,4])

-4 -3-2-1 01234 -4-3-2-10123 4x y x 3+y 3-3 x y = 0 或：t=-4:0.1:4; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)

-1.5 -1-0.500.51 1.5 00.5 1 1.5 2 2.5 3 4、蔓叶线233 2 2 2 ,()11at at x x y y t t a x = = = ++- 解：t=-4:0.1:4; x=t.^2./(1+t.^2); y=t.^3,/(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)

00.10.20.30.40.50.60.70.80.91 -4 -3-2-10123 4 或： ezplot('y .^2-x.^3/(1-x)',[-4,4])

数学实验作业题目(赛车跑道)

数学实验报告实验题目：赛车车道路况分析问题小组成员：填写日期2012 年 4 月20 日

一．问题概述赛车道路况分析问题现要举行一场山地自行车赛，为了了解环行赛道的路况，现对一选手比赛情况进行监测，该选手从A地出发向东到B，再经C、D回到A地（如下图）。现从选手出发开始计时，每隔15min观测其位置，所得相应各点坐标如下表（假设其体力是均衡分配的）：由D→C→B各点的位置坐标（单位：km）假设：1. 车道几乎是在平原上，但有三种路况（根据平均速度v(km/h)大致区分）：平整沙土路（v>30）、坑洼碎石路（10

2．估计车道的长度和所围区域的面积； 3．分析车道上相关路段的路面状况（用不同颜色或不同线型标记出来）； 4．对参加比赛选手提出合理建议. 二．问题分析 1.模拟比赛车道的曲线：因为赛道散点分布不规则，我们需要用光滑曲线来近似模拟赛道。由于数据点较多，为了避免龙格现象，应采用三次样条插值法来对曲线进行模拟（spline命令）。全程曲线为环路，我们需要对上下两部分分别模拟，设模拟出的曲线为P：。 2.把A到B点的曲线分成若干小段：赛道的路程L：取dL=，对模拟出的整条曲线求线积分，即所围区域的面积：用上下部分曲线的差值对求定积分，即 3.用样条插值法模拟出比赛车道曲线后，根据曲线分别计算出原数据中每两点（）间的路程，即求线积分由于每两点间时间间隔相同且已知（15min），故可求出每段路程的平均速度易知即为的积分中值将此速度近似作为两点间中点时刻的速度，然后再次采用样条插值法，模拟出全过程的图像。而根据求出的与之间的关系，再次采用样条插值法，即可模拟出全过程的图像 4. 由赛道曲线可求出赛道上任一点到点的路程同时图像也可以求出赛道上任一点到点的路程

重庆大学《数值分析》期末考试真题及答案

一.填空题: 1．若求积公式对任意不超过 m 次的多项式精确成立,而对 m+1 次多项式不成立,则称此公式的代数精度为m 次. 2．高斯消元法求解线性方程组的的过程中若主元素为零会发生计算中断 ;. 主元素的绝对值太小会发生误差增大 . 3． ) 4．当A 具有对角线优势且不可约时，线性方程组Ax=b 用简单迭代法和塞德尔迭代法均收敛. 5．求解常微分方程初值问题的欧拉方法是 1 阶格式; 标准龙格库塔法是 4 阶格式. 6．一个n 阶牛顿-柯特斯公式至少有 n 次代数精度,当n 偶数时,此公式可以有 n+1 次代数精度. 、 7．相近数相减会扩大相对误差,有效数字越多,相对误差越大 . 二计算题: 1. 线性方程组: ??? ??-=++-=+-=++5 .1526235.333321 321321x x x x x x x x x 1) ￥ 2) 对系数阵作LU 分解,写出L 阵和U 阵； ???? ? ? ?-=????? ? ?--=79/123/54 1 33 14 /33/113 /11U L 3) 求出此方程组的解. )5.0,1,2('-=x 2. 线性方程组: — ??? ??=++-=++=++3 32212325223321 321321x x x x x x x x x 1)对系数阵作LU 分解,写出L 阵和U 阵；

?? ??? ? ?=?? ?? ? ??=573235223 152321321//////U L 2)求出此方程组的解. ),,(' -=133x 4) # 5) 此方程组能否用用简单迭代法和高斯塞德尔迭代法求解. 073 2 2 232223053 2 2 3 03>=>=>,, A 对称正定,用高斯-塞德尔迭代法收敛; . .,., //////)(,6667033331027 16 3432323232323232131 =-==+-=-?? ?? ? ?? -=+-=-λλλλλJ J B I U L D B 用简单迭代法不收敛 > 3. 设f (x )= x 4, 以-1,0,1,2为插值节点, 1) 试写出f (x )的三次拉格朗日插值多项式P 3(x )及其插值余项R 3(x ); 6 ) 2)(1())()(())()(()(3020103210--- =------= x x x x x x x x x x x x x x x x l

北理工数学实验作业

一． 1. 1/e 2. 3 3.1 4.e3 5. ∞ 6. 0 7.∞ 8.0 9.1/2 10.0 11.e2c12.不存在13. 1/12 Matlab实验过程： 1.1/exp(1) syms n; f=(1-1/n)^n; limit(f,n,inf) ans = 1/exp(1) 2.3 syms n; f=(n^3+3^n)^(1/n); limit(f,n,inf) ans = 3 3. 1 syms n; f=(1+sin(2*n))/(1-cos(4*n)); limit(f,n,pi/4) ans = 1 4.e^3 syms x; f=(1+cos(x))^(3*sec(x)); limit(f,x,pi/2) ans = exp(3) 5.inf syms x; f=(x^2)*exp(1/(x^2));

limit(f,x,0) ans = Inf 6.0 syms x; f=(x^2-2*x+1)/(x^3-x); limit(f,x,1) ans = 7.inf syms x; f=((2/pi)*atan(x))^x; limit(f,x,+inf) ans = Inf 8.0 syms x y; f=(1-cos(x^2+y^2))/((x^2+y^2)*exp(x^2+y^2)); limit(limit(f,x,0),y,0) ans = 9.1/2 syms x; f=(1-cos(x))/(x*sin(x)); limit(f,x,0) ans = 1/2 10.0 syms x;

f=atan(x)/(2*x); limit(f,x,inf) ans = 11.exp(2*c) syms c; f=sym('((x+c)/(x-c))^x'); limit(f,'x',inf) ans = exp(2*c) 12.极限不存在 syms x; f=cos(1/x); limit(f,x,0) ans = limit(cos(1/x), x = 0) 13.1/12 syms x; f=1/(x*log(x)^2)-1/(x-1)^2; limit(f,x,1) ans = 1/12 二．观察函数logbx,当b=1/2,1/3,1/4和b=2,3,4时函数的变化特点，总结logbx的图形特点。

2015年重庆大学数学分析研考题(精)

重庆大学2015年硕士研究生入学考试试题科目代码:621 科目名称:数学分析总分:150 分特别提醒:所有答案一律写在答题纸上,直接写在试题上的不给分。一、计算(6分/每小题,共24分 (1(( (1 2 2lim 111n n x x x -→∞ +++ (1x < (2 (2 1x xe dx x +? (3 2 sin 1cos x x dx x π

+? (4((21 1lim 1n n k nx k nx k n →∞=+++∑ 二、(10分设(f x 在(0,+∞上满足函数方程((2f x f x =,且(0 lim x f x C →=(常数,证明:(f x C ≡,(0,x ∈+∞. 三、(13分若(f x 在(,-∞+∞上可微,且(lim x f x →∞ =-∞,证明:存在(,ξ∈-∞+∞使得(0f ξ'=. 四、(15分设(,α∈-∞+∞,讨论级数????? +∑∞ =n n n n ln 1sin 12 πα 的绝对收敛性与条件收敛性. 五、(13分计算(32sin 2x y z dxdydz Ω ++???,Ω由旋转双曲面2221x y z +-=、平面z H =、z H =-所围成. 六、(15分计算(2 222 axdydz z a dxdy

I x y z ∑ ++=++?? ,其中∑为下半球222 z a x y =---的上侧,0a >. 七、(15分令2 1 sin( (1xt f t dx x +∞ =+?,证明: (1反常积分关于t 在(,-∞+∞上一致收敛; (2函数(f t 在(,-∞+∞上连续,且lim (0t f t →+∞ =. 八、(15分函数(f x 为(,-∞+∞上的单调增加有界函数, (1证明:对于任意(0,x ∈-∞+∞,(0 lim x x f x →+存在; (2讨论(lim x f x →-∞ 的存在性,并说明理由. 九、(15分讨论(肯定,给出证明;否定,举出反例: (1对无穷限反常积分,平方可积与绝对可积之间的关系; (2对无界函数反常积分,平方可积与绝对可积之间的关系. 十、(15分设11a =,21a =,2123n n n a a a ++=+,1n ≥, (1证明{} n a 的通项公式为113(12 n n n a --+-=; (2求

数学实验作业一

数学实验作业一对以下问题,编写M文件: (1)用起泡法对10个数由小到大排序. 即将相邻两个数比较,将小的调到前头. 解：代码如下： zuoye1 clear all;clc; a=[7 2 1 0 9 4 5 -3 8 6]; n=length(a); for ii=1:n-1 if a(ii+1)>=a(ii) t1=a(ii); a(ii)=a(ii+1); a(ii+1)=t1; end for jj=1:n-1 if a(jj+1)>=a(jj) t2=a(jj); a(jj)=a(jj+1); a(jj+1)=t2; end end end a 运行结果显示如下： a = 9 8 7 6 5 4 2 1 0 -3

(2)有一个矩阵,编程求出其最大值及其所处的位置. 解：代码如下：zuoye2.m clear; clc; a=[1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6] max=-1; flage1=0; flage2=0 for i=1:4 for j=1:5 if (a(i,j)>max) t=max; max=a(i ,j); a(i,j)=t; flage1=i; flage2=j ； end end end max flage1 flage2 运行结果显示如下： a = 1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6 flage2 = max = 45′

9 flage1 = 2 flage2 = 5 结果：（3）编程求∑=20 1 !n n 。解：代码如下：zuoye3.m clear; clc; sum=0; for i=2:11 sum=sum+gamma(i); end sum

“2011年重庆大学数学建模竞赛”报名通知

“2011年重庆大学数学建模竞赛”报名通知 2011/3/15 数学建模竞赛能够提高大学生建立数学模型与运用计算机技术解决实际问题的综合能力，拓展知识面，开阔眼界, 培养创新精神及合作意识，磨炼意志。为了推动创新教育的深入开展，让尽可能多的学生参与这项有益的活动，同时也为参加国内外建模竞赛选拔参赛队员，我校每年都要举办“重庆大学数学建模竞赛暨全国大学生数学建模竞赛选拔赛”。2011年重庆大学数学建模竞赛将于5月21日至5月30日举行,我们将从竞赛优胜者中选拔队员参加暑期的培训和9月份的全国竞赛。在3月26日下午起将在D1137举办数学建模竞赛周末讲座（其他具体时间地点请关注重庆大学“数学实验”国家级精品课程网“https://www.doczj.com/doc/5516703763.html,/cmewebhome/”上的“公告”栏目中的通知）。现将有关事宜通知如下：一、报名时间 3月16日——4月20日。二、报名地点各学院学生工作办公室。请各学院将名单汇总后于4月25日前将名单电子稿发送至gongqu@https://www.doczj.com/doc/5516703763.html,。三、对报名学生的要求一至三年级理工科学生。学生自愿组队参加竞赛，每队三人，学校鼓励跨院系组队。四、参加培训学生的待遇 1. 免费听课和上机培训； 2. 参加培训并代表学校参加全国竞赛者给予记2个创新实践学分； 3. 全国数学建模竞赛获奖者有资格代表学校参加次年2月的美国大学生数学建模竞赛； 4. 全国大学生数学建模竞赛获奖者，在推荐免试攻读硕士生时将按《重庆大学推荐优秀应届本科毕业生免试攻读硕士学位研究生和直接攻读博士学位研究实施办法（试行）》（重大校[2009]146号）相关规定获加分奖励。欢迎符合条件的同学踊跃报名，力争在竞赛中取得优秀成绩，为重庆大学争光。重庆大学数学与统计学院 2011.3.15

数值分析实验报告

学生实验报告实验课程名称开课实验室学院年级专业班学生姓名学号开课时间至学年学期

if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); %换行 end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); %消去end end x=zeros(length(b),1); %回代求解 x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end y=x; format short;%设置为默认格式显示,显示5位（2）建立MATLAB界面利用MA TLAB的GUI建立如下界面求解线性方程组：详见程序。五、计算实例、数据、结果、分析下面我们对以上的结果进行测试，求解：

? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - - - 7 2 5 10 13 9 14 4 4 3 2 1 13 12 4 3 3 10 2 4 3 2 1 x x x x 输入数据后点击和，得到如下结果：更改以上数据进行测试，求解如下方程组： 1 2 3 4 43211 34321 23431 12341 x x x x ?? ???? ?? ???? ?? ???? = ?? ???? - ?? ???? - ???? ?? 得到如下结果：

重庆大学数学实验报告七

开课学院、实验室：数统学院DS1421实验时间：2013年03月17日

由于matlab中小数只能是四位，所以我在编程的过程中将距离扩大了1000倍，但是并不会影响我们所求得的结果。运行程序之后我们得到的结果为：我们可以得到当金星与地球的距离（米）的对数值为9.9351799时，只一天恰好是25号。 8.编写的matlab程序如下： x=0:400:2800; y=0:400:2400; z=[1180 1320 1450 1420 1400 1300 700 900 1230 1390 1500 1500 1400 900 1100 1060 1270 1500 1200 1100 1350 1450 1200 1150 1370 1500 1200 1100 1550 1600 1550 1380 1460 1500 1550 1600 1550 1600 1600 1600 1450 1480 1500 1550 1510 1430 1300 1200 1430 1450 1470 1320 1280 1200 1080 940]; [xi,yi]=meshgrid(0:5:2800,0:5:2400); zi=interp2(x,y,z,xi,yi,'cubic'); mesh(xi,yi,zi); xlabel('x'),ylabel('y'),zlabel('高程'); title('某山区地貌图'); figure(2); contour(xi,yi,zi,30); 运行程序我们得到的结果如下所示：山区的地貌图如下所示：

等高线图如下所示：三、附录（程序等） 6. y=18:2:30;

数学实验8月13日作业

1.取不同的初值计算下列平方和形式的非线性规划，尽可能求出所有局部极小点，进而找出全局极小点，并对不同算法（搜索方向、搜索步长、数值梯度与分析梯度等）的结果进行分析、比较。 (2). ( )( ) 2 2 2 22 121212min 12114949812324681x x x x x x +-++++-, (4).()()212222 23 12123min10010,1x x x x x x θ??????-++-+?????????????? ，其中 ()()()21112211 1 arc ,02,11arc ,0 22tg x x x x x tg x x x π θπ ?>??=??+

重庆大学数学实验一 matlab的基本应用参考答案

《数学实验》第一次上机实验 1．设有分块矩阵?? ? ???= ????22322333S O R E A ，其中E,R,O,S 分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证?? ????+= 22 S 0RS R E A 。程序及结果： E=eye(3); %创建单位矩阵E% R=rand(3,2); %创建随机矩阵R% O=zeros(2,3); %创建0矩阵% S=diag(1:2); %创建对角矩阵% A=[E,R;O,S]; %创建A 矩阵% B=[E,(R+R*S);zeros(2,3),S^2] %计算等号右边的值% A^2 %计算等号左边的值% 运行结果： B = 1.00 0 0 1.63 2.74 0 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.00 ans = 1.00 0 0 1.63 2.74 0 1.00 0 1.81 1.90 0 0 1.00 0.25 0.29 0 0 0 1.00 0 0 0 0 0 4.00 2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。表1.1 1)程序： a=[7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30]; b=[11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50]; c=[568 1205 753 580 395 2104 1538 810 694];

数学实验第二次作业任务常微分方程数值求解

实验4常微分方程数值解实验目的： 1.练习数值积分的计算； 2.掌握用MATLAB软件求微分方程初值问题数值解的方法； 3.通过实例学习用微分方程模型解决简化的实际问题； 4.了解欧拉方法和龙格——库塔方法的基本思想和计算公式，及稳定性等概念。实验内容： 3.小型火箭初始质量为1400kg，其中包括1080kg燃料，火箭竖直向上发射是燃料燃烧率为18kg/s，由此产生32000N的推力，火箭引擎在燃料用尽时关闭。设火箭上升是空气阻力正比于速度的平方，比例系数为0.4kg/m，求引擎关闭瞬间火箭的高度，速度，加速度，及火箭到达最高点是的高度，速度和加速度，并画出高度，速度，加速度随时间变化的图形。解答如下：这是一个典型的牛顿第二定律问题，分析火箭受力情况；先规定向上受力为正数建立数学模型： A燃料未燃尽前，在任意时刻（t<60s）火箭受到向上的-F=32000N，向下的重力G=mg，g=9.8，向下的阻力f=kv^2, k=0.4, v表示此时火箭速度；此时火箭收到的合力为F1=（F-mg-f）；火箭的初始质量为1400kg，燃料燃烧率为-18kg/s; 此刻火箭质量为m=1400-18*t

根据牛顿第二定律知，加速度a=F1/m=(F-mg-f)/(m-r*t)= (32000-(0.4.*v.^2)-9.8.*(1400-18.*t)) 由此可利用龙格-库塔方法来实现，程序实现如下 Function [dx]=rocket[t,x] %建立名为rocket的方程 m=1400;k=0.4;r=-18;g=9.8; %给出题目提供的常数值dx=[x(2);(32000-(k*x(2)^2)-g*(m+r*t))/(m+r*t)]; %以向量的形式建立方程[a]=(32000-(k*x(2)^2)-g*(m+r*t))/(m+r*t); %给出a的表达式 End; ts=0:60; %根据题目给定燃烧率计算出燃料燃尽的时间，确定终点 x0=[0,0]; %输入x的初始值[t,x]=ode15s(@rocket,ts,x0); %调用ode15s计算 [t,x]; h=x(:,1); v=x(:,2); plot(t,x(:,1)),grid; %绘出火箭高度与时间的关系曲线 title('h-t'); xlabel('t/s');ylabel('h/m'),pause; plot(t,x(:,2)),grid ; %绘出火箭速度与时间的曲线关系 title('v-t'); xlabel('t/s');ylabel('v/m/s'),pause; a=(32000-(0.4.*v.^2)-9.8.*(1400-18.*t))/(1400-18.*t);

重庆大学--数学模型--数学实验作业七

开课学院、实验室：数统学院实验时间：2015年11月25日课程名称数学实验实验项目名称医用薄膜渗透率的确定 ——数据拟合实验项目类型验证演示综合设计其他指导教师肖剑成绩实验目的 [1] 了解最小二乘拟合的基本原理和方法； [2] 掌握用MATLAB作最小二乘多项式拟合和曲线拟合的方法； [3] 通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区别。 [4] 了解各种参数辨识的原理和方法； [5] 通过范例展现由机理分析确定模型结构，拟合方法辨识参数，误差分析等求解实际问题的过程；通过该实验的学习，掌握几种基本的参数辨识方法，了解拟合的几种典型应用，观察不同方法得出的模型的准确程度，学习参数的误差分析，进一步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。实验内容 1．用MATLAB中的函数作一元函数的多项式拟合与曲线拟合，作出误差图； 2．用MATLAB中的函数作二元函数的最小二乘拟合，

作出误差图； 3．针对预测和确定参数的实际问题，建立数学模型，并求解。应用实验（或综合实验） 1.旧车价格预测一、问题重述某年美国旧车价格的调查资料如下表，其中xi表示轿车的使用年数，yi表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？表1 x i 1 2 3 4 5 6 7 8 9 10 y i 26 15 19 43 14 94 10 87 76 5 53 8 48 4 29 22 6 20 4 二、数学模型的建立与求解先作出散点图分析其应该是一个二次函数，可以采用polyfit线性拟合。编辑程序Untitled1.m: clc x=[1 2 3 4 5 6 7 8 9 10]; y=[2615 1943 1494 1087 765 538 484 290 226 204]; plot(x,y,'+') hold on a=polyfit(x,y,2) y1=polyval(a,x); plot(x,y1,'r') t=4.5; cost=polyval(a,t) 三、实验结果及分析 a =1.0e+03* 0.0361 -0.6508 3.1523 t =4.5000

时数分配 - 重庆大学-数学与统计学院-数学实验

2011—2012学年度第二学期教学日历课程名称：数学实验任课教师姓名：龚劬课程类别：（）必修课（）选修课教材名称：数学实验主编姓名刘琼荪出版时间2004.7 授课对象：计算机学院计算机科学1—5班、网络工程1-3班、信息安全1班140 人填表时间：2012 年 3 月

教学日历

数学软件自己动手做实验。第7次教学内容： 1. 应用实例：放射性废物的处理问题问题重述、分析、假设，建立数学模型，模型求解 2．方程和方程组求解的MATLAB命令及其应用。教学方式：多媒体教学2 14 第8次实验内容： 1．使用MATLAB软件求解方程与方程组的练习； 2．应用问题：炮弹发射角的确定。教学方式：学生在教师指导下，借助于计算机和数学软件自己动手做实验。4 3 18 第9次教学内容： 1. 引例：倒葫芦形状容器壁上的刻度问题微分方程模型及其求解方法解析法，数值解法：欧拉方法，梯形法，改进欧拉方法教学方式：多媒体教学2 20 第10次实验内容： 1．使用MATLAB软件求解微分方程（组）的练习； 2．编用向前欧拉公式和改进的欧拉公式求微分方程数值解的MATLAB程序，并观察其迭代过程；教学方式：学生在教师指导下，借助于计算机和数学软件自己动手做实验。4 3 24 第11次教学内容： 1. 求解微分方程（组）的MATLAB命令追击路线问题教学方式：多媒体教学2 26 第12次实验内容： 1．用MATLAB命令求解Rossler微分方程组，并讨论解随参数的变化情况； 2.考虑两相互竞争种群的数量变化模型；4 3 30

重庆大学数学与统计学院

重庆大学数学与统计学院推荐免试攻读硕士研究生实施办法及操作细则根据教育部办公厅《关于进一步完善推荐优秀应届本科毕业生免试攻读研究生工作的通知》（教学厅〔2013〕8号）和学校《重庆大学推荐优秀应届本科毕业生免试攻读硕士学位研究生工作管理办法（试行）》（重大校〔2014〕268号）及相关文件﹑通知的精神，结合我院的实际情况，特制定本实施细则。一﹑推荐免试研究生条件 1. 申请者应符合重庆大学对应届本科毕业生申请免试攻读硕士学位研究生基本条件的规定。 2. 补充业务条件（1）学习态度端正，成绩优异，前三年的平均成绩在同专业学生中排名应处于前列，方可进入推免资格的候选人名单，并予以公布。（2）本科学习阶段内必修课和专选课补考科目不得超过1门（无不及格成绩），特殊情况需经学院推免研究生工作小组讨论研究决定。 3. 可适当突破第2条限制的情况： (1) 在全国性的大学生数学竞赛，数学建模竞赛活动中获国家二等奖以上的学生，直接具有推免资格，但须满足基本条件且复试合格。 (2) 基础课和专业课成绩优异，并且具有浓厚数学兴趣和具有培养潜质者优先推荐。但需要2位专家推荐。二﹑综合成绩计算办法综合成绩：60% A +40% B + C (附加分数) 1、A——平均成绩平均成绩按三年计算，课程只包括必修课和专业选修课（五级记分折算标准：优=95分、良=85、中=75分、及格=65分）。 2、B——按百分制给出的面试成绩，其中1）笔试科目80%，2）专业面试10%，3）英语口语面试10%。 3、C:见加分细则

三、免试研究生的推荐程序 1. 学院组成推免研究生工作小组，由院长为组长，学院党政班子、学院学位委员会、研究生教学工作委员会和教学管理人员为成员。 2. 由学生本人提出申请，报学院推免研究生工作小组。 3. 学院推免研究生工作小组从符合推荐免试研究生基本条件的申请学生中，根据学生平均成绩及优先情况进行排序，并按推荐免试研究生名额的1.5倍比例，确定具有推荐免试研究生面试人选名单，并予以公布。 4. 学院推免研究生工作小组组织专家对初选合格的学生进行面试，根据面试专家个人评分，计算每个学生的平均分。 5．面试包括： 1）笔试科目：数学分析、高等代数； 2）专业面试：面试老师有统一评分标准，对所有专业大方向相同的考生使用相同的面试题目，已面试的考生在所有面试结束前不能离开面试考场； 3）英语口语面试。 6. 学院推免研究生工作小组根据平均成绩、面试成绩和获奖得分，计算综合成绩，进行排序，并向学生公布。 7. 学院推免研究生工作小组根据综合成绩排名，确定获得推免资格的初选学生名单，并张榜公布三天。 8. 公布无异后，初选学生名单经学院推免研究生工作小组审核盖章后上报教务处。四、本实施细则自颁布之日起实行，并由数学与统计学院推免研究生工作小组负责解释。重庆大学数学与统计学院 2016年9月5日

数学实验作业汇总

（1）产生一个5阶魔方矩阵M：M=magic(5) （2）将矩阵M的第3行4列元素赋值给变量t：t=M(3,4) （3）将由矩阵M第2，3，4行第2,5列构成的子矩阵赋给变N：N=M(2:4,2:3:5) （4）将由矩阵M的前3行赋给变量N：N=M(1:3,:) （5）将由矩阵M的后3列赋给变量N：N=M(:,end:-1:end-2) （6）提取M的主对角线元素，并以这些对角线元素构成对角矩阵N：N=diag(diag(M))或N=tril(triu(M)) (7)随机产生1000个100以内的整数赋值给变量t：t=round(rand(1,1000)*100) （8）随机产生100*5个100以内的实数赋值给变量M：M=rand(100,5)*100 （1）删除矩阵M的第7个元素M(7)=[] （2）将含有12个元素的向量t转换成3*4的矩阵：reshape(t,3,4) （3）产生和M同样大小的单位矩阵：eye(size(M)) （4）寻找向量t中非零元素的下标：find(t) （5）逆序显示向量t中的元素：t(end:-1:1) （6）显示向量t偶数位置上的元素：t(2:2:end) （7）利用find函数，将向量t中小于10的整数置为0：t(find(t<10&rem(t,1)==0))=0 （8）不用find函数，将向量t中小于10的整数置为0：t(t<10&rem(t,1)==0)=0 （9）将向量t中的0元素用机器0（realmin）来代替：t(find(t=0))=realmin （10）将矩阵M中小于10的整数置为0：M(find(M<10)&rem(M,1)==0)=0 2、写出完成下列操作的命令及结果。（1）将1~50这50个整数按行优先存放到5*10的矩阵中，求该矩阵四周元素的和； >> t=[1:10]; >> M=[t;t+10;t+20;t+30;t+40] M = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 >> N=M(2:4,2:9) N = 12 13 14 15 16 17 18 19 22 23 24 25 26 27 28 29 32 33 34 35 36 37 38 39 >> sum(sum(M))-sum(sum(n)) ans = 663 2）n取100、1000、10000，求序列1、1/2、1/3……1/n的和。

文档之家