第一篇:数学建模与数学实验教学大纲(工科)
数学建模与数学实验教学大纲(工科)总学分:3 总上课时数:48 或32
一、课程的性质与目的
本课程是面向理工科学生开设的一门选修课。本课程的教学目的是让学生增加一些用数学的感性认识,初步掌握一些基本的建模方法、建模原理和数学软件的应用。学生通过这门课的学习,在数学知识的综合运用,将实际问题转化为数学问题的能力方面、创新能力、自学能力方面、发散性思维能力方面都能得到一定培养。
二、适用专业
数学大类、工科各专业
三、课程内容的教学要求
(1)数学建模与数学实验概述:介绍数学建模与数学实验的基本概念,熟悉建模步骤。
(2)初等模型:掌握用初等函数对实际问题的变化关系作简单的定量分析;熟悉用图示法对实际问题作定性分析。
(3)量纲分析建模:掌握量纲分析原理,学会用量纲分析原理对一些物理问题作一些分析;了解数学中的无量纲化方法;掌握非线性方程求根的常用方法。
(4)代数学模型:介绍矩阵在解决实际问题中的应用,熟悉层次分析法的建模步骤,学会用矩阵思想分析实际问题;掌握线性方程组的数值揭解法和矩阵特征值与特征向量的近似求法。
(5)静态优化模型:了解微积分在解决实际问题中应用,掌握静态优化建模的基本步骤;熟悉微分、积分的数值方法。
(6)数值分析法建模:掌握曲线拟合、插值的基本方法,学会用插值、拟合作数据处理,了解插值、拟合建模的大致过程。
(7)常微分方程模型:熟悉微分方程建模的基本步骤,掌握线性微分方程建模基本方法,了解非线性微分方程模型的一些特殊性质;熟悉微分方程的数值解法。
(8)差分方程模型:了解差分法的基本思想,学会建立实际问题的离散模型,掌握递推、迭代法的求解过程。
(9)统计模型与实验 学习简单的随机模型的建模方法,熟悉Matlab工具箱的应用;
(10)优化模型:了解最优化思想,熟悉优化建模思路,能建立和求解一些简单的优化模型;会在适当的数学软件上实现优化模型。
四、上机要求
学会Matlab的基本操作、学会非线性方程求根,能在该软件平台上进行较大规模的数据处理及求解微分方程及优化问题。能更具体实际问题在软件上实现小规模编程运算。
五、能力培养
1.实际问题分析能力的培养:通过对实际问题的分析,抓住问题本质,才能建立满意的数学模型。
2.实际问题转化为数学问题能力的培养:要求学生通过本课程的学习,初步掌握将实际问题转化为数学问题的方法,能够建立简单的实际问题的数学模型。
3.自学能力、语言表达能力的培养:课程安排了大量自学内容,要求学生通过查阅文献,写论文等形式完成课后作业,使学生自学能力等得到培养。
4.创新能力的培养:课程里许多范例都是来源于实际问题,属于开放型的问题,学生可以充分展开自己的思维,开放式的学习,促使学生独立思考、深入钻研。
六、教材与参考书
1.陈恩水.《数学建模与实验》,自编讲义,2004.2.姜启源编.数学模型.北京,高等教育出版社,1992,第二版.3.郑家茂编.数学建模基础.南京,东南大学出版社,1997.4.朱道元编.数学建模精品案例.南京,东南大学出版社,1999.5.萧树铁主编.数学实验.北京, 高等教育出版社,1998.6.乐经良主编.数学实验.北京, 高等教育出版社,1999.
第二篇:数学建模与数学实验
通过多年来的教学改革与教学实践,教学效果显著,模块化分层次教学、换位式教学和启发式教学的方法得到了学生们的认可。这种方式大大提高了学生们的动手能力,并贯穿于平时的教学实践中,同时也反映出学生撰写科技论文的写作水平,为学生进一步参加数学建模竞赛奠定了良好的基础。该课程的成功经验在我校、市内以及西部地区起到很好的示范辐射作用,得到专家和学生的好评。
校外专家
(一)评价:
刘琼荪(全国数学建模竞赛重庆赛区组委会秘书长,重庆大学教授)
重庆邮电大学是我国最早开设数学建模系列课程的学校之一, 经过十多年的努力,该课程已经建设成为培养学生的创新和竞争能力的优秀课程。该课程在教学环节上充分体现出了教学研究型大学的特色,坚持培养学生“以竞赛为契机,以能力提升为宗旨”的指导思想,在教学内容和教学方式方面进行了大胆、慎重的改革, 把课堂教学、课后实践、在数学建模基地做数学实验、参加讨论班研讨、参加国内外数学建模竞赛结合起来,既激发了学生进一步学习数学的兴趣,又提高了学生的科学素质和能力,收到了很好的效果。该类课程自开设以来,已有逾万名学生学习本课程。全校每年有1000余名学生参加全国或校内竞赛,近三年参加全国大学生数模竞赛中, 获全国奖27项(规定每年一个学校最多10项), 成绩在重庆赛区参赛学校中名列前茅。另外,陈理荣教授等编著的教材《数学建模导论》(北京邮电大学出版社出版)也已为全国20余所大学用作数学建模课程的教材被广泛使用,杨春德教授等编著的《数学建模的认识与实践》也为本门课程的建设提供了素材。且《数学建模》已成为重庆市精品课程,“数学建模与数学实验”教学团队已获重庆市市级教学团队称号。
有鉴于此,我认为《数学建模与数学实验》已完全达到了重庆邮电大学重点课程的要求。
校外专家
(二)评价:
朱宁(全国大学生数学建模优秀指导教师,桂林电子科技大学教授)
全国大学生数学建模竞赛自90年代在我国开展以来,一直受到全国各高校的重视,把竞赛作为培养数学知识应用的一个平台。重庆邮电大学是较早参加这活动的高校,近几年,在竞赛中屡获佳绩,走在同类高校的前列,引起了广泛的重视。本人认为重庆邮电大学在数学建模赛成功的主要经验有如下几方面: 首先是有一支实力雄厚、敬业的师资队伍。《数学建模与数学实验》课程建设成员11名,其中有教授4人,副教授6人,4人具有博士学位,1人获全国大学生数学建模竞赛优秀指导教师称号。教学成果多,教学团队整体实力强,“数学建模与数学实验”教学团队已获重庆市市级教学团队称号。
其次《数学建模与数学实验》类课程形成了“三层次—两阶段”的教学和竞赛的课程改革方案,设计并探索了数学应用型人才培养理念,在教学模式和教学方法和评价方式等方面均有创新,形成了“教学-实践-竞赛” 的数学建模教学模式,形成了一套具有特色的加强数学模型思想的教学模式。
第三是注重校际间交流,吸取好的经验,完善教学过程。教师曾多次在国内外关于数学建模教学与应用会议上介绍经验,并先后在国内外核心期刊上发表论文数篇。每年参加赛区举办的数学实验课程和数学建模竞赛的教学经验交流会议。该课程建设已在西部地区起到了示范作用。
鉴于以上内容,个人认为《数学建模与数学实验》已达到了重庆邮电大学重点课程的要求。
校内同行评价
胡学刚(全国数学建模竞赛优秀指导教师,重庆邮电大学教务处副处长、教授)
《数学建模与数学实验》类课程先后为不同层次的学生开设了任选课、限选课和必修课。近年来,课程建设小组以《数学建模与数学实验》类课程为平台,以数学建模竞赛为契机,在工科数学类课程的教育教学改革中取得了突出成绩,主要表现在以下几个方面:
1.坚持数学建模类课程建设与工科数学教学改革相结合,数学建模类课程建设与数学建模竞赛相结合,理论教学与实验实践、课外活动相结合,将数学建模的思想融入到其它数学类课程的教学中,进一步深化工科数学类课程的教学改革。该课程建设特色鲜明,成效显著。
2.课题组老师热情指导学生开展数学建模活动,积极组织学生参加校内、国内及美国大学生数学建模竞赛。从最初的鼓励学生参赛,到现在同学们积极主动参赛;从最初的几个队参赛到现在的近百个队参赛,数学建模竞赛经历了一次次飞跃。经过多年的探索,课题组总结了一套成功的指导培训经验,使我校学生参加全国竞赛取得了优异成绩,近3年来,我校共有27个队获得国家级奖励,在重庆赛区位居前列,特别是2011年名列全国第二(公示中)。
3.师资队伍建设成效显著。近年来,课题组新增2位教师获得博士学位,1位教师博士即将毕业,教授由申报时的0人变为4人。队伍中现拥有全国模范教师、重庆市中青年骨干教师、重庆邮电大学优秀青年教师。他们多次在赛区组织的教练交流活动中介绍数学建模类课程程建设经验和竞赛经验,在重庆市乃至西部地区发挥了示范辐射作用。
4.课程建设成绩显著。在该门课程建设过程中,编著出版了《数学建模的认识与实践》一书,《数学建模》已成为重庆市精品课程,“数学建模与数学实验”已获重庆市市级教学团队称号,《数学建模理论与方法》于2011年成为重庆邮电大学立项建设教材。
有鉴于此,该课程是有较大影响的富有特色的课程,已具备了重庆邮电大学重点课程的条件。
学生评价
(一):
数学建模与数学实验这门课程是一门开放性和主动性的一门课程,它就是需要从现实生活、现实问题中抽象出数学模型,从而解决问题。这门课程融合了许多学科,对于学生来说,有机会广泛涉猎各种知识,这对于我们后续的发展是十分有好处的,因为目前在实际部门工作,也许不需要你对某一方面的有很深的知识,主要是遇到一个问题,能有解决的方法;再有就是对于继续深造的同学,也十分有益,因为通过广泛的知识储备,学生可以从中找到自己感兴趣的方向,继续深入的做下去,《数学建模与数学实验》这门课就为我们在这两方面打下了良好的基础。
同时,数学建模有利于培养学生的创造性思维能力,数学建模主要考查学生的数学思想方法,它是一种数学活动,而不单单像传统的数学练习题一样,做出来的答案是唯一的。相反,它可以有多种多样的答案,只要学生建立的模型是可行的,那它就是正确的。在学习这门课程的过程中,我也做过很多的实际题目,从那些过程中,我体会到的数学在实际生活中的应用,更重要的是培养了我们合作交流的方法、习惯,特别是促进学生的数学应用意识,提高了解决实际问题的能力。无论是数学研究还是数学学习,其目的之一就是将数学运用于社会,运用于现实,数学建模就重视培养学生的数学思维,加强数学应用意识,切实提高分析和解决实际问题的能力。
学习《数学建模与数学实验》是我大三的时候,朱伟老师将这门数学课讲得生动有趣,他没有介绍过于高深的理论,而是从实际应用出发。让我们对这门课程充满了兴趣,同时也对数学有了重新的认识,目前我正在进行硕士研究生阶段的学习,觉得那个时候学到的一些理论知识还有用,虽然那个时候没有过多的去深入研究那些知识,但现在当我遇到问题的时候,我知道有那样的一个理论存在,所以对于我来说就多了一些解决问题的方法。总之,在解决实际问题时,我们只有多了解一些方法,才能去掌握它,从而运用它,《数学建模与数学实验》就是一个连接理论与实际应用的桥梁。
(重庆邮电大学信息与计算科学专业,现西南财经大学统计学院硕士研究生 周黎)
校内学生(二)评价
大一的时候我就接触过数学建模,那是学校组织的数学建模竞赛,我们小组在比赛中获得了第三名,虽然是一个小小的第三名,当时还是给我很大的鼓舞,因为那时候大一能得奖好像只有两组,因此这学期一听说要开数模选修课,我就立马去报了名,抱着一点能学点东西的态度,认认真真的听完了前面大半的内容,后面由于很难坐倒好坐位,就只有自学了。
通过这门课的学习,我认识到了数模课多么的博大精深,虽然还是要靠一点小聪明,但主要还是要靠勤奋,因为数模涉及到太多的东西了,基本涉及到所有数学方面的知识,还有社会,科学等各方面的知识,要想能在这上面有所成就,只有靠平时的认真学习,打下牢实的基础。只有这样,才有可能在这上面有所发展。学习这门课,不管从学知识的角度,还是从学做学问的角度,对我而言,我都有很大的收获,衷心感谢各位数学组的老师在星期六不辞辛苦为我们上课。
(重庆邮电大学通信学院, 杨鹏)
校内学生(三)评价
从小到大,我对数学充满了爱好和兴趣,于是报名参加了数模学习辅导班。通过一个学期的数模学习,使自己学到了很多东西,不仅对数模的概念有了一定的了解,对数学建模的方法有了一定的掌握,同时也使自己加深了对数学知识的理解,能灵活运用数学解决一些实际吻题。数学建模是一种具有创造性的科学方法,它将现实问题简化,抽象为一个数学问题或者数学模型,然后采用恰当的数学方法求解,进而对现实问题进行定量分析和研究,最终达到解决实际问题的目的。随着计算机的运用和发展,数学建模成为高科技的一种“数学技术”,起着关键性的作用,作为计算机学员的一名学生,掌握新的技术和方法是必要的,是受益匪浅的。通过一个学期的学习,数模培养了我的洞察力,想象力,逻辑思维能力以及分析问题,解决问题的能力。在学习过程中,虽然碰到了很多的问题和困难,但是在老师的指点和教导下,使得很多问题都得到了解决,在这里要感谢辛勤教育我们的老师。虽然我没有去参加数模竞赛,但是我确实学到了很多东西,我相信这些我所学到的知识,对我的将来是有好处的。
(重庆邮电大学计算机学院:陈辉)
第三篇:数学建模教学大纲
数学建模教学大纲
(32学时)
一、课程内容简介
数学建模是研究如何将数学方法和计算机知识结合起来用于解决实际生活中存在问题的一门边缘交叉学科,数学建模是集经典数学、现代数学和实际问题为一体的一门新型课程,是应用数学解决实际问题的重要手段和途径。主要介绍数学建模的概述、初等模型、简单优化模型、微分方程模型、离散模型、线性规划模型、概率模型等模型的基本建模方法及求解方法。
二、教学目的及任务
数学建模是计算机类高职生继高等数学、线性代数之后进一步提高运用数学知识解决实际问题、基本技能,培育和训练综合能力所开设的一门新学科。通过具体实例引入使学生掌握数学建模基本思想、基本方法、基本类型。学会进行科学研究的一般过程,并能进入一个实际操作的状态。通过数学模型有关的概念、特征的学习和数学模型应用实例的介绍,培养学生双向翻译能力,数学推导计算和简化分析能力,熟练运用计算机能力;培养学生联想、洞察能力、综合分析能力;培养学生应用数学解决实际问题的能力。
三、本课程与其它课程的关系
在学习本课程前需要基本掌握下列课程内容:高等数学、线性代数、概率论与数理统计。由于本课程的学习,只要是使学生掌握数学知识,解决实际问题能力,这种能力提高有助其它专业课的学习。
四、本课程基本内容要求
以建立不同的数学模型作为教学项目载体,每个项目分解为若干个学习任务(学习情境),每个学习任务按照资讯、决策、计划、实施、检查、评估、拓展步骤进行教学组织和内容设计。教学内容按照教学做一体化的思路设计,实现实践教学与理论教学的相互渗透。教学内容
教学项目一:建立数学模型
学习学时:2 学习目标:(1)了解数学建模的历史和现状;开展数学建模的意义,熟悉数学模型的基本概念;数学模型的特点和分类;(2)掌握数学建模的方法及基本步骤的知识,并能用于指导全部课程的学习。(3)使学生正确地了解数学描写和数学建模的不同在于数学理论的思维特征。、教学内容:(1)数学建模的历史和现状(2)高职院校开设数学建模课的现实意义(3)数学模型的基本概念(4)数学模型的特点和分类(5)数学建模的方法及基本步骤。教学方法:(1)案例分析法(2)任务驱动法(3)演示法(4)分组讨论法
对学生要求:学生对数学建模有初步的了解,与较强地团结协作能力,具备计算机的基础和知识,具备上Internet 网查资料的能力,具备Office的基础知识和能力.教学项目二:初等数学建模
学习学时:2 学习目标:(1)掌握比例法,类比方法、图解法、定性分析方法建模的基本特点(2)能运用所学知识建立数学建模,并对模型进行综合分析‘
教学内容(1)初等函数建模法:基本初等函数数学模型;常用的经济函数模型(2)集合建模法:鸽笼原理;“奇偶效验”法;相识问题(3)比例与函数建模法:动物体型模型;双重玻璃的功效模型;席位分配模型。教学方法:(1)案例教学法(2)任务驱动教学法(3)探究式教学法
对学生要求:具有较扎实的初等数学功底;具有较强的理解能力;具有发现问题,探究问题的精神,具备计算机的基础知识,具备上Internet网查找资料的能力,具备office的基础知识和能力。
教学项目三:微分方程建模
学习学时:4 学习目标:(1)了解微分方程稳定性理论(2)熟悉微分模型的一般意义(3)掌握微分方程模型的建模建立思想求解方法(4)能够建立简单的微分方程模型解决实际问题。教学内容:(1)微分方程建模方法(2)熟悉微分方程建模案例:Malthus模型;Logistic模型;具有收获的单种群模型(3)经济增长模型;资金与劳动力的最佳分配;劳动生产率增长;(4)人口的预测和控制(5)微分方程稳定性理论简介 教学方法:(1)案例教学法(2)任务驱动教学法(3)探究式教学法
对学生要求:对微分方程知识有一定的了解;具有较强的理解能力、知识综合能力;具有发现问题,研究问题的精神;具有计算机的基础知识,具备上Internet网查找资料的能力,具备office的基础知识和能力。教学项目四:数学规划建模
学习学时:6 学习目标:(1)深刻理解数学规划模型的基本特点,理解模型的一般意义(2)较熟练的建立数学规划模型解决实际问题(3)能熟练的结合计算机软件求解数学规划模型。教学内容:(1)想行规划模型原理与案例:运输模型;食谱模型;河流污染与净化模型;合理下料模型(2)非线性规划模型原理与案例:投资决策模型;武器分配模型;防洪优化问题;森林救火费用最小模型(3)0-1规划模型原理与案例:饮料厂的生产与检修计划模型;
指派问题模型;投资决策问题模型 教学方法:(1)案例教学法(2)任务驱动教学法(3)探究式教学法
对学生要求:具有较扎实的线性规划知识功底;具有较强的理解能力、知识综合运用能力;具有发现问题,探究问题的精神;具有计算机的基础知识,具备上Internet网查找资料的能力,具备office的基础知识和能力。教学项目五:概率统计建模
教学学时:4 学习目标:(1)了解概率统计建模方法(2)熟悉案例模型
教学内容:报童卖报模型;随机存贮模型;商店进货策略模型。教学方法:(1)案例教学法(2)任务驱动教学法(3)探究式教学法
对学生要求:了解概率论数理统计、线性代数基础知识;具有较强的理解能力、知识综合运用能力;具有发现问题,探究问题的精神;具有计算机的基础知识,具备上Internet网查找资料的能力,具备office的基础知识和能力。教学项目六:层次分析建模
学习学时:4 学习目标: 1)了解层次分析法(2)深刻理解层次分析法建模的基本特征(3)熟练掌握层次分析法建模的典型案例级方法(4)能够运用层次分析法解决日常生活中简单的相关问题 教学内容:(1)层次分析法原理、步骤、特点(2)层次分析法案例:选拔干部模型;循环比赛的名次(3)效益的合理分配方法 教学方法(1)案例教学法(2)任务驱动教学法(3)探究式教学法
对学生要求:就有较扎实的线性代数知识功底;具有较强的理解能力、知识综合运用能力;具有发现问题,探究问题的精神;具备计算机的基础知识,具备上Internet网查找资料的能力,具备office的基础知识和能力 教学项目七:插值与拟合建模 学习目标:(1)了解插值、拟合的基本特点(2)熟练掌握插值与拟合的建模案例 教学内容:(1)插值方法与案例(2)拟合方法与案例 教学方法:(1)案例教学法(2)任务驱动教学法(3)探究式教学法
对学生要求:具有较强数据处理能力、理解能力、知识综合运用能力;具有发现问题,探究问题的精神;具备计算机的基础知识,具备上Internet网查找资料的能力,具备office的基础知识和能力 教学项目八:常用数学软件基础知识及其应用
学习目标:了解LINGO、MATLAB的基本知识能够及在建模中进行的应用 教学内容:(1)LINGO的基础知识(2)LINGO在建模中的应用案例(3)MATLAB的的基础知识(4)MATLAB在建模中的应用案例
六、本课程的教材和参考书
教 材:数学建模(第二版)徐全智主编,高等教育出版社 参考书:数学模型 姜启源 谢金星 叶俊 高等教育出版社
数学模型 杨启帆 浙江大学出版社
数学模型 任善强 高等教育出版社
数学模型与数学建模 刘来福 曾文艺 北京师范大学出版社
第四篇:数学建模实验小结
例1-1 >> r=2;V=4/3*pi*r^3 V =
33.5103 例2-1 计算s=...>> s=0;>> for n=1:100 s=s+1/n/n;end >> s s =
1.6350 例2-5 两个一元函数y=x3-x-1,y=x .2sin(5x)在区间-1 y=abs(x).^0.2.*sin(5*x);plot(x,y,':ro');hold off; 曲面图 >> xa=6:8;ya=1:4;>> [x,y]=meshgrid(xa,ya);>> z=x.^2+y.^2;>> mesh(x,y,z)>> [x,y,z] ans = 例2-6 二元函数图z=xexp(-x2-y2).xa=-2:0.2:2;ya=xa; [x,y]=meshgrid(xa,ya);z=x.*exp(-x.^2-y.^2);mesh(x,y,z);pause;surf(x,y,z);pause; contour(x,y,z,[0.1,0.1]);pause mesh(x,y,z); Page40 1.先在编辑器窗口写下列M函数,保存为ex2_1.m function [xbar,s]=ex2_1(x)n=length(x);xbar=sum(x)/n; s=sqrt((sum(x.^2)-n*xbar^2)/(n-1)); >> x=[81 70 65 51 76 66 90 87 61 77];>> [xbar,s]=ex2_1(x)xbar = 72.4000 s = 12.1124 2.s=log(1);n=0;while s<=100 n=n+1; s=s+log(1+n);end m=n 3.F(1)=1;F(2)=1;k=2;x=0;e=1e-8;a=(1+sqrt(5))/2;while abs(x-a)>e k=k+1;F(k)=F(k-1)+F(k-2);x=F(k)/F(k-1);end a,x,k m = a = 1.6180 x = 1.6180 k = 4.clear;tic;s=0;for i=1:1000000 s=s+sqrt(3)/2^i;end s,toc tic;s=0;i=1; while i<=1000000 s=s+sqrt(3)/2^i;i=i+1; end s,toc tic;s=0; i=1:1000000; s=sqrt(3)*sum(1./2.^i);s,toc s = 1.7321 Elapsed time is 2.038973 seconds.s = 1.7321 Elapsed time is 2.948968 seconds.s = 1.7321 Elapsed time is 0.453414 seconds 5.t=0:24; c=[15 14 14 14 14 15 16 18 20 22 23 25 28...31 32 31 29 27 25 24 22 20 18 17 16];plot(t,c) 6.(1)x=-2:0.1:2;y=x.^2.*sin(x.^2-x-2);plot(x,y)y=inline('x^2*sin(x^2-x-2)');fplot(y,[-2 2]) (2)t=linspace(0,2*pi,100); x=2*cos(t);y=3*sin(t);plot(x,y) (3)x=-3:0.1:3;y=x; [x,y]=meshgrid(x,y);z=x.^2+y.^2; surf(x,y,z) (4) x=-3:0.1:3;y=-3:0.1:13;[x,y]=meshgrid(x,y); z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6;surf(x,y,z) (5) t=0:0.01:2*pi; x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z) 7.x=linspace(0,pi,100); y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);plot(x,y1,x,y2,x,y3)%page41, ex7 x=-1.5:0.05:1.5; y=1.1*(x>1.1)+x.*(x<=1.1).*(x>=-1.1)-1.1*(x<-1.1);plot(x,y) Page59 1.>> a=[1 2 3];b=[2 4 3];>> a./b ans = 0.5000 0.5000 1.0000 >> a.b ans = >> a/b ans = 0.6552 >> ab ans = 0 0 0 0 0 0 0.6667 1.3333 1.0000 2.(1)>> a=[4 1-1;3 2-6;1-5 3];b=[9;-2;1];>> ab ans = 2.3830 1.4894 2.0213(2)>> a=[4-3 3;3 2-6;1-5 3],b=[-1;-2;1] a = -5 b = >> ab ans = -0.4706 -0.2941 0(3)>> a=[4 1;3 2;1-5],b=[1;1;1] a = -5 b = >> ab ans = 0.3311 -0.1219(4)>> a=[2 1-1 1;1 2 1-1;1 1 2 1],b=[1;2;3] a = -1 b = >> ab ans = 0 0 6.(1)>> a=[4 1-1;3 2-6;1-5 3];>> b=det(a),inv(a),[V,D]=eig(a)b = -94 ans = 0.2553 -0.0213 0.0426 0.1596 -0.1383 -0.2234 0.1809 -0.2234 -0.0532 V = 0.0185 -0.9009 -0.3066 -0.7693 -0.1240 -0.7248 -0.6386 -0.4158 0.6170 D = -3.0527 0 0 0 3.6760 0 0 0 8.3766(2)>> a=[1 1-1;0 2-1;-1 2 0];b=det(a),inv(a),[V,D]=eig(a)b = ans = 2.0000 -2.0000 1.0000 1.0000 -1.0000 1.0000 2.0000 -3.0000 2.0000 V = -0.5773 0.5774 + 0.0000i 0.57740.0000i 0.5773 + 0.0000i D = 1.0000 0 0 0 1.0000 + 0.0000i 0 0 0 1.00000.0000i -0.5773 0.5774 0.5774 -0.5774 0.57730.0000i >> det(V)ans = -5.0566e-028-5.1918e-017i %V的行列式接近0, 特征向量线性相关,不可对角化(3)>> a=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10];[V,D]=eig(a)V = 0.8304 0.0933 0.3963 0.3803 -0.5016 -0.3017 0.6149 0.5286 -0.2086 0.7603 -0.2716 0.5520 0.1237 -0.5676 -0.6254 0.5209 D = 0.0102 0 0 0 0 0.8431 0 0 0 0 3.8581 0 0 0 0 30.2887 >> inv(V)*a*V ans = 0.0102 0.0000 -0.0000 0.0000 0.0000 0.8431 -0.0000 -0.0000 -0.0000 0.0000 3.8581 -0.0000 -0.0000 -0.0000 0 30.2887 8 对称阵A为正定的充分必要条件是:A的特征值全为正。只有(3)对称, 且特征值全部大于零, 所以(3)是正定矩阵.例4.2用2次多项式拟合下列数据。>> clear;x=[0.1,0.2,0.15,0,-0.2,0.3];>> y=[0.95,0.84,0.86,1.06,1.50,0.72];>> p=polyfit(x,y,2)p = 1.7432 -1.6959 1.0850 得到二次拟合式:1.7432x^2-1.6959x+1.0850 >> xi=-0.2:0.01:0.3;>> yi=polyval(p,xi);plot(x,y,'o',xi,yi); 例4.3 求函数y=x*sin(x^2-x-1)在(-2,-0.1)内的零点。>> fun=inline('x*sin(x^2-x-1)','x')fun = Inline function: fun(x)= x*sin(x^2-x-1)>> fzero(fun,[-2,-0.1])??? Error using ==> fzero at 292 The function values at the interval endpoints must differ in sign.>> fplot(fun,[-2,-0.1]);grid on; >> [x,f,h]=fsolve(fun,-1.6),[x,f,h]=fsolve(fun,-0.6)Optimization terminated: first-order optimality is less than options.TolFun.x = -1.5956 f = 1.4909e-009 h = Optimization terminated: first-order optimality is less than options.TolFun.x = -0.6180 f = -3.3152e-012 h = 例4.4求下列方程组在原点附近的解 >> fun=inline('[4*x(1)-x(2)+exp(x(1))/10-1,-x(1)+4*x(2)+x(1)^2/8]','x');[x,f,h]=fsolve(fun,[0,0])Optimization terminated: first-order optimality is less than options.TolFun.x = 0.2326 0.0565 f = 1.0e-006 * 0.0908 0.1798 h = 例4.5 求二元函数f(x,y)=5-x^4-y^4+4*x*y在原点附近的极大值。(等价于求-f(x,y)的极小值) >> fun=inline('x(1)^4+x(2)^4-4*x(1)*x(2)-5');>> [x,g]=fminsearch(fun,[0,0])x = 1.0000 1.0000 g = -7.0000 例4.6 用Newton迭代法求下列方程的正根,要求精度为10的-6次 X^2-3x+e^x=2 >> fun=inline('x^2-3*x+exp(x)-2');>> fplot(fun,[0,2]);>> grid on; %M函数 newton.m function x =newton(fname,dfname,x0,e)if nargin<4,e=1e-4;end x=x0;x0=x+2*e;while abs(x0-x)>e x0=x;x=x0-feval(fname,x0)/feval(dfname,x0);end >> dfun=inline('2*x-3+exp(x)');format long;newton(fun,dfun,1.5,1e-6),format short ans = 1.*** 例4.7 用函数y=a*e^(b*x)拟合例4.2的数据。>> fun=inline('c(1)*exp(c(2)*x)','c','x');>> x=[0.1,0.2,0.15,0,-0.2,0.3];y=[0.95,0.84,0.86,1.06,1.50,0.72];>> c=lsqcurvefit(fun,[0,0],x,y)Optimization terminated: relative function value changing by less than OPTIONS.TolFun.c = 1.0997 -1.4923 PAGE 77 1.%Exercise 1(1)roots([1 1 1])%Exercise 1(2) roots([3 0-4 0 2-1])%Exercise 1(3)p=zeros(1,24); p([1 17 18 22])=[5-6 8-5];roots(p) %Exercise 1(4)p1=[2 3]; p2=conv(p1, p1);p3=conv(p1, p2); p3(end)=p3(end)-4;%原p3最后一个分量-4 roots(p3)2.%Exercise 2 fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x');fzero(fun,2)3.%Exercise 3 fun=inline('x^4-2^x');fplot(fun,[-2 2]);grid on; fzero(fun,-1),fzero(fun,1),fminbnd(fun,0.5,1.5) 4.%Exercise 4 fun=inline('x*sin(1/x)','x');fplot(fun, [-0.1 0.1]); x=zeros(1,10);for i=1:10, x(i)=fzero(fun,(i-0.5)*0.01);end;x=[x,-x] x = Columns 1 through 11 0.0050 0.0152 0.0245 0.0354 0.0455 0.0531 0.0637 0.0796 0.0796 0.1061 -0.0050 Columns 12 through 20 -0.0152 -0.0245 -0.0354 -0.0455 -0.0531 -0.0637 -0.0796 -0.0796 -0.1061 5.%Exercise 5 fun=inline('[9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);16*x(1)-x(1)^3-2*x(2)^2-16*x(3)^2]','x'); [a,b,c]=fsolve(fun,[0 0 0])6.%Exercise 6 fun=@(x)[x(1)-0.7*sin(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))];[a,b,c]=fsolve(fun,[0.5 0.5])7.%Exercise 7 clear;close;t=0:pi/100:2*pi; x1=2+sqrt(5)*cos(t);y1=3-2*x1+sqrt(5)*sin(t);x2=3+sqrt(2)*cos(t);y2=6*sin(t); plot(x1,y1,x2,y2);grid on;%作图发现4个解的大致位置,然后分别求解 y1=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.5,2])y2=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.8,-2])y3=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[3.5,-5])y4=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[4,-4]) 8.%Exercise 8(1)clear; fun=inline('x.^2.*sin(x.^2-x-2)');fplot(fun,[-2 2]);grid on;%作图观察 x(1)=-2; x(3)=fminbnd(fun,-1,-0.5);x(5)=fminbnd(fun,1,2); fun2=inline('-x.^2.*sin(x.^2-x-2)');x(2)=fminbnd(fun2,-2,-1);x(4)=fminbnd(fun2,-0.5,0.5);x(6)=2 feval(fun,x)x = -2.0000 -1.5326 -0.7315 -0.0000 1.5951 2.0000 ans = -3.0272 2.2364 -0.3582 -0.0000 -2.2080 0 %答案: 以上x(1)(3)(5)是局部极小,x(2)(4)(6)是局部极大,从最后一句知道x(1)全局最小,x(2)最大。 %Exercise 8(2)clear; fun=inline('3*x.^5-20*x.^3+10');fplot(fun,[-3 3]);grid on;%作图观察 x(1)=-3; x(3)=fminsearch(fun,2.5); fun2=inline('-(3*x.^5-20*x.^3+10)');x(2)=fminsearch(fun2,-2.5);x(4)=3;feval(fun,x)ans = -179 -54 199 %Exercise 8(3) fun=inline('abs(x^3-x^2-x-2)');fplot(fun,[0 3]);grid on;%作图观察 fminbnd(fun,1.5,2.5) fun2=inline('-abs(x^3-x^2-x-2)');fminbnd(fun2,0.5,1.5)ans = 2.0000 ans = 1.0000 9.%Exercise 9 close; x=-2:0.1:1;y=-7:0.1:1;[x,y]=meshgrid(x,y); z=y.^3/9+3*x.^2.*y+9*x.^2+y.^2+x.*y+9;mesh(x,y,z);grid on;%作图观察 fun=inline('x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9');x=fminsearch(fun,[0 0])%求极小值 fun2=inline('-(x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9)');x=fminsearch(fun2,[0-5])%求极大值 x = 0 0 x = -0.3333 -6.0000 10.clear;t=0:24; c=[15 14 14 14 14 15 16 18 20 22 23 25 28...31 32 31 29 27 25 24 22 20 18 17 16];p2=polyfit(t,c,2)p3=polyfit(t,c,3) fun=inline('a(1)*exp(a(2)*(t-14).^2)','a','t'); a=lsqcurvefit(fun,[0 0],t,c)%初值可以试探 f=feval(fun, a,t) norm(f-c)%拟合效果 plot(t,c,t,f)%作图检验 fun2=inline('b(1)*sin(pi/12*t+b(2))+20','b','t');%原题修改f(x)+20 b=lsqcurvefit(fun2,[0 0],t,c)figure f2=feval(fun2, b,t) norm(f2-c)%拟合效果 plot(t,c,t,f2)%作图检验 Page 94 chapter5 1.x=[0 4 10 12 15 22 28 34 40];y=[0 1 3 6 8 9 5 3 0];trapz(x,y)ans = 178.5000 2.>> x=[0 4 10 12 15 22 28 34 40];y=[0 1 3 6 8 9 5 3 0];diff(y)./diff(x) ans = 0.2500 0.3333 1.5000 0.6667 0.1429 -0.6667 -0.3333 -0.5000 3.xa=-1:0.1:1;ya=0:0.1:2;[x,y]=meshgrid(xa,ya);z=x.*exp(-x.^2-y.^3); [px,py] = gradient(z,xa,ya);Px 4.t=0:0.01:1.5;x=log(cos(t));y=cos(t)-t.*sin(t);dydx=gradient(y,x) [x_1,id]=min(abs(x-(-1)));%找最接近x=-1的点 dydx(id)5.(1)(2)fun=inline('exp(2*x).*cos(x).^3');quadl(fun,0,2*pi)(3)fun=@(x)x.*log(x.^4).*asin(1./x.^2);quadl(fun,1,3)(4)fun=@(x)sin(x)./x; quadl(fun,1e-10,1)%注意由于下限为0,被积函数没有意义,用很小的1e-10代替(5)(6)fun=inline('sqrt(1+r.^2.*sin(th))','r','th');dblquad(fun,0,1,0,2*pi)(7)首先建立84页函数dblquad2 clear; fun=@(x,y)1+x+y.^2;clo=@(x)-sqrt(2*x-x.^2);dup=@(x)sqrt(2*x-x.^2);dblquad2(fun,0,2,clo,dhi,100)%Exercise 6 t=linspace(0,2*pi,100);x=2*cos(t);y=3*sin(t); dx=gradient(x,t);dy=gradient(y,t);f=sqrt(dx.^2+dy.^2);trapz(t,f)10(1)(2) %先在程序编辑器,写下列函数,保存为ex5_10_2f function d=ex5_10_2f(fname,a,h0,e) h=h0;d=(feval(fname,a+h)-2*feval(fname,a)+feval(fname,a-h))/(h*h);d0=d+2*e; while abs(d-d0)>e d0=d;h0=h;h=h0/2; d=(feval(fname,a+h)-2*feval(fname,a)+feval(fname,a-h))/(h*h);end %再在指令窗口执行 fun=inline('x.^2*sin(x.^2-x-2)','x');d=ex5_10_2f(fun,1.4,0.1,1e-3)13.fun=inline('5400*v./(8.276*v.^2+2000)','v');quadl(fun,15,30) 《数学实验》教学大纲 课程编号:课程类型:公共基础课 课程名称:数学实验英文名称:Mathematical Experimentation学分:3适用对象:文理本科、专科 第一部分 大纲说明 一、课程的性质、目的和任务 数学实验课程系统地介绍数学实验中的一些常用方法及实例探索,通过课堂教学讨论,把一些数学概念直观而形象的显现出来。将形象思维与逻辑思维结合,并通过上机实验,将抽象的数学公式、定理通过实验得到验证和应用,调动学生学习数学的积极性,加强对学生的数学知识、软件知识、计算机知识和动手能力的培养。使学生了解数学实验的特性及实验的基本方法,并初步具备对实际问题如何实验的能力以及培养良好的思考习惯和归纳分析能力.二、课程的基本要求 掌握数学实验的主要步骤,学会使用相应数学软件解决数学实验问题。能对一些典型数学现象展开讨论。 三、本课程与相关课程的联系 本课程学生应具备高中数学知识基本知识和计算机操作基本技能.四、学时分配 本课程学分为3学分,建议开设54学时,第一学期讲授。 使用教材: 李尚志等编,数学实验(第二版)[M], 北京:高等教育出版社.2004 1、萧树铁、姜启源.数学实验[M].北京:高等教育出版社,2003.2、乐经良.数学实验[M].北京:高等教育出版社,2003.3、赵静.数学建模与数学实验[M].北京:高等教育出版社与施普林格出版社,2000.六、教学方法和手段建议 本课程以多媒体教学为主,板书讲解为辅,淡化数学的计算证明推导,强化使用计算机实验的能力。 七、课程考核方式 实验考核以学生的实验态度、掌握的实验理论、实际操作技能和实验报告等为主,各单项平时内容所占分数比例为:实验态度占10%、实验理论占15%;操作技能占50%;实验报告占25%。本课程进行平时作业和期末小论文的评判。 成绩评定方法:平时30%+期末70%.第二部分课程内容大纲 第一章 实验1平面几何初等函数(6学时) 一、本章的教学目的和要求 了解数学实验的基本方法和步骤,熟练掌握几何画板和Z+Z超级画板的使用。 二、教学内容 试验案例精选:初等几何著名定理的动态演示,初等函数的一阶和二阶导函数曲线;文本中的数学曲线编辑 重点:实验步骤 难点:参数选择.第二章 实验2一(多)元函数微积分(6学时) 一、本章的教学目的和要求 在EXCEL环境下,熟练掌握极限与导数的概念,定积分的核心概念,函数的做图 二、教学内容 试验案例精选:;试验:数学软件EXCEL; 重点:导数 积分.难点:$符号.第三章 实验3微分方程(6学时) 一、本章的教学目的和要求 掌握Mathematica软件的主要用法,能够动态显示微分曲线随初始条件变化的特点,会分析一般的方程数值解的相轨迹曲线的变化特点。 二、教学内容 试验案例精选:复杂的微分方程数值解和相轨线,方程的稳定性,极限环表现,混沌与孤立子表现。 重点:数值解 相轨线 难点:方程的稳定性分析 第四章 实验4线性代数、矩阵(6学时) 一、本章的教学目的和要求 综合使用EXCEL与MATHEMATICA软件进行矩阵讨论,病态矩阵的动态演示,广义逆矩阵的应用,迭代收敛的条件实验。线性规划初步(Lindo/Lingo软件初步) 二、教学内容 试验案例精选:矩阵的各种分解,线性方程组的解空间演示,矩阵特征根与特征向量动态变化,AHP应用,代数方程求根计算 重点:矩阵分解 难点:机器误差处理 第五章 实验5数论密码分析(6学时) 一、本章的教学目的和要求 综合使用EXCEL与MATHEMATICA软件进行数论问题与密码分析问题讨论。 二、教学内容 试验案例精选:数论基本定理实验证明,RSA公钥密码体系,椭圆密码体系(群结构分析),因数分解。 重点:RSA公钥 难点:群结构.第六章 实验6概率统计(6学时) 一、本章的教学目的和要求 综合使用EXCL和SPSS软件讨论随即现象,了解概率和统计中的基本问题及其实验观察。 二、教学内容 试验案例精选:伪随机数的产生,常用分布的动态计算,中心极限定理的演示,置信区间的置信度演示,回归分析的预测,假设检验中两类错误α与β的计量分析(样本量n)。 重点:数据模型的猜想 难点:α与β 第七章 实验7复变函数数理方程(6学时) 一、本章的教学目的和要求 学习MATLAB软件,用MATLAB重述以前各章的试验结果。讨论复变函数与偏微分方程概念与计算问题。 二、教学内容 试验案例精选:复变函数图像,留数定理,复函数映射问题,积分变换问题;数理方程数值解法及其动态显示 重点:复变函数 数理方程 难点:参数微调 第八章 分形混沌NP问题(6学时) 一、本章的教学目的和要求 综合使用各种数学软件编程讨论:分形现象与混沌现象,初步体验NP问题,发现计算机应用的NP特点 二、教学内容 试验案例精选:雪花曲线、分形集合计算机显示,洛伦兹吸引子,混沌分叉现象,十进制位的大整数因数分解NP难题(体验).重点:编程设计.难点:编程设计 第九章 实验9数学机械化证明(6学时) 一、本章的教学目的和要求 专题讨论我国科学家(吴文俊-高小山)在机械化证明中的工作(专题软件).二、教学内容 试验案例精选:吴序列方法,机器证明初步.重点:吴序列方法 难点:机器证明 执笔人:纪跃 审定人:魏兰阁 批准人:周和月 制定(修订)日期:2007.9第五篇:数学实验教学大纲