数学实验报告题目

大学数学实验报告题目

实验报告一MATLAB基本操作实验题目

1.1 随机抽取1个班的《高等数学》课程成绩如下，并统计他们中的最高分、最低分以及他们的平均成绩。

60 84 83 69 39 60 79 88 75 82 80 80 95 73 60 74 60 67 71 83

63 60 60 76 65 72 76 90 98 77 76 86 60 61 97

1.2设A为2×3矩阵，A=

123

456

⎛⎫

⎪

⎝⎭

，试建立一个与矩阵A同样大小零矩阵，幺矩阵，单

位矩阵，0～1间均匀分布的随机矩阵，均值为0，方差为1的标准正态分布随机矩阵。1.3 建立随机矩阵。

(1) 在区间[10,100]内均匀分布的4阶随机矩阵；

(2) 均值为0.5、方差为0.7的4阶正态分布随机矩阵。

1.4 产生5阶随机方阵A，其元素为[5,70]区间的随机整数，然后判断A的元素是否能被3整除。

1.5 建立矩阵A=

510607

6060100

⎛⎫

⎪

⎝⎭

，然后找出大于4的元素的位置。

实验报告二微积分实验题目

教材62页：1（5），2（2），3（2），4（2），6（5），9（2）共6题

实验报告三函数作图题目

教材63页：8题

实验报告四线性代数实验题目

教材80页：1，3（2），4（3），5（3），6（3），7（2）

实验报告五MATLAB程序设计实验题目

1. 根据我国个人所得税计算方法,编制程序,要求:使用者在系统提示下通过键盘输入月工资薪金收入总数,计算机则在屏幕上显示个人所得税额,界面友好,方便使用.

个人所得税计算方法:

月个人所得税＝（月工资薪金收入－2000）*适用税率－速算扣除数

2求［2，999］中同时满足下列条件的数

（１）该数各位数字之和为奇数

（２）该数是素数

注：教材是汪晓银等. 数学软件与数学实验（第二版）. 科学出版社, 2010

高等数学实验报告一

高等数学数学实验报告实验人员：院（系）学号姓名实验地点：实验一实验题目：作出函数Y=ln(cosx^2+sinx) (-π/4, π/4)的函数图形和泰勒展开式图形，选取不同的x0和n，并进行比较。二、实验目的和意义熟悉Mathematica软件所具有的良好作图功能，并通过函数图形来认识函数，运用函数图形来观察分析函数的有关性态，熟悉泰勒多项式对函数的近似效果，并建立数形结合的思想三、程序设计 (一)确定x0=0,n=1,3,5……19; ●t = Table[Normal[Series[Log[Cos[x^2] + Sin[x]], {x, 0, i}]], {i, 1, 20, 2}];PrependTo[t, Log[Cos[x^2] + Sin[x]]];Plot[Evaluate[t], {x, -Pi/4, Pi/4}] (二)分别在x0=0.1,0.5,0.6,0.75,-0.1,-0.5,-0.75,0处进行泰勒展开，确定n=1,3,5……19； ●t = Table[Normal[Series[Log[Cos[x^2] + Sin[x]], {x, 0.1, i}]], {i, 1, 20, 2}];PrependTo[t, Log[Cos[x^2] + Sin[x]]];Plot[Evaluate[t], {x, -Pi/4, Pi/4}] ●t = Table[Normal[Series[Log[Cos[x^2] + Sin[x]], {x, 0.5, i}]], {i, 1, 20, 2}];PrependTo[t, Log[Cos[x^2] + Sin[x]]];Plot[Evaluate[t], {x,

《数学实验》实验报告

《数学实验》实验报告（ 2012 年 4 月 8 日）一、实验问题 1．（指派问题）考虑指定n个人完成n项任务（每人单独承担一项任务），使所需的总完成时间（成本）尽可能短. 已知某指派问题的有关数据（每人完成各任务所需的时间）如下表所示，试建模并求解该指派问题。 2．（二次指派问题）某公司指派n个员工到n个城市工作（每个城市单独一人），希望使所花费的总费用尽可能少。n个员工两两之间每个月通话的时间表示在下面的矩阵的上三角部分（因为通话的时间矩阵是对称的，没有必要写出下三角部分），n个城市两两之间通话费率表示在下面的矩阵的下三角部分（同样道理，因为通话的费率矩阵是对称的，没有必要写出上三角部分）. 试求解该二次指派问题。 3、金星第四章课后习题第1或3题任选一题。二、问题的分析（涉及的理论知识、数学建模与求解的方法等） 1)根据实际问题，建立数学优化模型 2)根据优化模型，利用LINGO 来求解模型。

三、计算过程、结论和结果分析 1. 模型： ij 44114 141 : 1,2,3,4 : 1234 1 i j a 0 i j x : i j model min 1 j=1,2,3,4.. 1 i=1,2,3,4ij ij ij i j ij i ij j m n a x a s t a ====?=???=????=??∑∑∑∑工人任务，，，第个人完成第项任务第个人不完成第项任务第个工人完成第项任务所用的时间 model : sets : m/1..4/; n/1..4/; link(m,n):a,x; endsets min =sum (link(i,j):x(i,j)*a(i,j)); for (m(i):sum (n(j):a(i,j))=1); for (n(j):sum (m(i):a(i,j))=1); data : x=15 18 21 24 19 23 22 18 26 18 16 19 19 21 23 17; enddata end 结果：Global optimal solution found. Objective value: 70.00000 Infeasibilities: 0.000000 Total solver iterations: 7 Variable Value Reduced Cost A( 1, 1) 0.000000 0.000000 A( 1, 2) 1.000000 0.000000 A( 1, 3) 0.000000 5.000000

数学实验报告1

数学实验报告实验一 1.题目：某车间有甲，乙，丙三台车床可用于加工三种零件，这三台车床可用于工作的最多时间分别为700，800和900，需要加工的三种零件的数量分别为300，400和500.不同车床加工不同的零件所用时间和费用如表所示，在完成任务的前提下，如何分配加工任务才能使加工费最少？

工时数分配表车床加工单位零件所需时数加工单位零件所需费用可用于工名称零件1 零件2 零件3 零件1 零件2 零件3 作的时数甲0.6 0.5 0.5 7 8 8 700 乙0.4 0.7 0.5 8 7 8 800 丙0.8 0.6 0.6 7 9 8 900 2.分析问题：此题考察用Matlab软件求线性规划问题。这是一个优约束的优化问题，其模型包括：甲生产零件1，零件2，零件3的个数分别为x1,x2,x3;以乙生产零件1，零件2，零件3的个数分别为x4,x5,x6;丙生产零件1，零件2，零件 3的个数分别为x7,x8,x9. 目标函数为： W=7x1+8x2+8x3+8x4+7x5+8x6+7x7+9x8+8x9 约束条件为： X1+x4+x7=300 X2+x5+x8=400 X3+x6+x9=500 0.6x1+0.5x2+0.5x3<=700 0.4x4+0.7x5+0.5x6<=800 0.8x7+0.6x8+0.6x9<=900 以及非负性约束x1,x2,x3,x4,x5,x6,x7,x8,x9都大于等于零。

3.建立模型： Min W=7x1+8x2+8x3+8x4+7x5+8x6+7x7+9x8+8x9; X1+x4+x7=300 X2+x5+x8=400 X3+x6+x9=500 0.6x1+0.5x2+0.5x3≤700 0.4x4+0.7x5+0.5x6≤800 0.8x7+0.6x8+0.6x9≤900 X1,x2,x3,x4,x5,x6,x7,x8,x9≥0 4.编写程序： c=[7,8,8,8,7,8,7,9,8] A=[0.6,0.5,0.5,0,0,0,0,0,0;0,0,0,0.4,0.7,0.5,0,0,0;0,0,0,0,0,0,0.8,0.6,0.6,] b=[700;800;900] aeq=[1,0,0,1,0,0,1,0,0;0,1,0,0,1,0,0,1,0;0,0,1,0,0,1,0,0,1] beq=[300;400;500] vlb=[0,0,0,0,0,0,0,0,0] vub=[] [x,minz]=linprog(c,A,b,aeq,beq,vlb,vub) 5.运行结果： c = 7 8 8 8 7 8 7 9 8

数学实验报告

实验报告一、实验问题追击问题如图，在一边长为1的正方形跑道的四个顶点上各站有1人，他们同时开始以等速顺（1）四个人能否追到一起？（2）若能追到一起，则每个人跑过多少路程？（3）追到一起所需要的时间（设速率为1）？（4）如果四个人追逐的速度不一样，情况又如何呢？二、问题分析这是一个追及问题，模拟其追击过程就是将整个追击过程离散化，即以dt为时间间隔，每个人都一步步跳动，如此持续直到二者之间距离足够小为止。根据条件和追击规律，实时计算出四个人的位置：（以A点为坐标原点，边AD为x轴建立直角坐标系）刚开始四个人的坐标如右：A(0.00,0.00) B(0.00,1.00) C(1.00,1.00) D(1.00,0.00)，d=1.00。追逐开始后，第一个dt内，四个人各自在自己的初始方向上跑一步，则他们的坐标分别为A(0.00,0.01) B(0.01,1.00) C(1.00,0.99) D(0.99,0.00)，d=0.99，那么下一个dt内，A追B，A 沿AB单位向量的方向上跑K1*v*dt*e1的距离，则此时A的坐标为A+ K1*v*dt*e1，依此类推，B追C，B沿BC的方向上跑K2*v*dt*e2的距离，则此时B的坐标为B+K2*v*dt*e2, 同理，C追D，D追A，当第二个dt 时，记ABCD初始的坐标分别为（A1，A1），（B1，B1），（C1，C1）（D1，D1），则他们在这个dt（A2，A2）=（A1,A1）+K1v*dt*e11,e11为（B1,B1）-(A1,A1)的方向向量；（B2，B2）=(B1,B1)+K2*v*dt*e21,e21为（C1，C1）-（B1，B1）的方向向量，以此类推，可得到（C2,C2）和（D2，D2）的坐标。如此这般，随着n 的增加，四个人之间的距离逐渐缩小，直到距离小于程序所设定的精度为止。当dt很小时，这逐步的过程变可以模拟真实的情况，误差可以忽略。三、程序设计 clear;clc;clf; %清出内存变量，清理图型窗口 hold on %开启图形保持功能以便重复画点 axis([0 1.2 0 1.2]) %置坐标窗口 grid A=[0,0] ; %分别标记四人的初始位置B=[0,1]; C=[1,1]; D=[1,0]; C D

数学实验报告

西安交通大学实验报告一、某棉纺厂的原棉需从仓库运送到各车间，各车间的原棉需求量，单位产品从各仓库运往各车间的运输费以及各仓库的库存如表所列，问如何安排运输任务使得总运费最小？问题分析：该题较为简单，只要根据表中数据确定不等式，找到上下限，在根据书上的已有例子，综合自己的判断，就可写出。 f=[2,1,3,2,2,4,3,4,2]; A=[1,1,1,0,0,0,0,0,0;0,0,0,1,1,1,0,0,0;0,0,0,0,0,0,1,1,1]; b=[50;30;10]; aeq=[1,0,0,1,0,0,1,0,0;0,1,0,0,1,0,0,1,0;0,0,1,0,0,1,0,0,1] ; beq=[40,15,35]; vlb=[0,0,0,0,0,0,0,0,0]; vub=[];

[x,fval]=linprog(f,A,b,aeq,beq,vlb,vub) 结果分析：由运行结果可知，第一车间由1，2仓库分别运进10，20单位的原棉，第二车间由1仓库运进15单位的原棉，第三车间由1，3仓库分别运进25，10单位的原棉，即可使总运费最小。二、某校学生在大学三年级第一学期必须要选修的课程只有一门，可供限定选修的课程有8门，任意选修课程有10门，由于一些课程之间互有联系，所以可能在选修某门课程中必须同时选修其他课程，这18门课程的学分数和要求同时选修课程的相应信息如表：

按学校规定，每个学生每学期选修的总学分不能少于21学分，因此，学生必须在上述18门课程中至少选修19学分学校同时还规定学生每学期选修任意选修课的学分不能少于3学分，也不能超过6学分，为了达到学校的要求，试为该学生确定一种选课方案。问题分析：本题是一道典型的0-1规划的问题，本体的难点在于，选了B一定要选A，但选了A却有选B，和不选B这两种方案，故不可采用以前普通的计算方式，考虑相减，即A-B>=0就可解决该问题。 c=[-5,-5,-4,-4,-3,-3,-3,-2,-3,-3,-3,-2,-2,-2,-1,-1,-1,-1]; a=[-5,-5,-4,-4,-3,-3,-3,-2,-3,-3,-3,-2,-2,-2,-1,-1,-1,-1; 0,0,0,0,0,0,0,0,3,3,3,2,2,2,1,1,1,1; 0,0,0,0,0,0,0,0,-3,-3,-3,-2,-2,-2,-1,-1,-1,-1; -1,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0; 0,-1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0; 0,0,0,0,0,0,0,-1,1,0,0,0,0,0,0,0,0,0; 0,0,0,0,0,-1,0,0,0,1,0,0,0,0,0,0,0,0; 0,0,0,-1,0,0,0,0,0,0,1,0,0,0,0,0,0,0; 0,0,0,0,-1,0,0,0,0,0,0,1,0,0,0,0,0,0;

工数实验报告上

工科数学分析数学实验学号：姓名：

高等数学数学实验报告实验人员：院（系）：学号：姓名：成绩_________ 实验一：观察数列的极限一、实验题目根据上面的实验步骤，通过作图，观察重要极限： e n =∞→n )n 1 + (1lim 二、实验目的和意义从点图上看出数列的收敛性，以及近似地观察出数列的收敛值；通过编程可以输出数列的任意多项值，以此来得到数列的收敛性。通过此实验对数列极限概念的理解形象化、具体化。三、计算公式四、程序设计五、程序运行结果六、结果的讨论和分析从点图可以看出，该数列是收敛的，并且收敛值在2.7左右，所以可以估计出e 的近似值为2.7

实验二：一元函数图形及其性态一、实验题目制作函数y=sin cx的图形动画，并观察参数c对函数图形的影响二、实验目的和意义通过作图形动画，观察参数c对函数性态（周期，最值，奇偶，凹凸）的影响，从而对函数的理解形象化、具体化。三、计算公式 sin(-x)=sin(x) sin(x+2π)=sin(x) sin(x+π)=-sin(x) 四、程序设计五、程序运行结果

六、结果的讨论和分析当参数|c|越大，函数的周期越小，并且符合T=2π/|c|; 参数c 的变化并不影响函数的最值，奇偶性（当p=0时，函数是既奇又偶函数），和凹凸性。参数c 的正负决定函数是在某一确定周期内的正负值。实验四定积分的近似计算一、实验题目分别用梯形法、抛物线法计算定积分dx x ? 2 2 sin π 二、实验目的和意义掌握积分公式背后的原理、方法三、计算公式梯形法: ])(2)()([)(11 ∑? -=-+++-≈n i b a n a b i a f b f a f n a b dx x f 抛物线法: )](2)(4)()([6)(1 1 2112∑∑? -==-+++-≈k i i k i i b a x x f f b f a f k a b dx x f 四、梯形法程序设计

数学实验报告4

实验报告4 实验名称数列与级数实验目的通过计算机图示的方法发现数列与级数的规律及其极限状态的性质。实验环境 Mathematica 4 实验内容 1. 分别取N=10，20，50，100，500，观察Fibonacci 数列的折线图。 2. 分别取N=2000，5000，10000，用直线去拟合N n F n n ,,2,1)),log(,( =的函数。 3. 分别取N=100，500，5000，演奏Fibonacci 数列的函数。 4. 分别取N=100，1000，5000，显示点列n i i i ,,2,1)),sin(,( =的函数。 5. 求级数∑ ∞ =11 n n α 的部分和。实验的基本理论和方法所谓一个无穷数列是指按一定顺序排列的一串数字 ,,,,21n a a a ，（1）而一个无穷级数则是用无穷项数字构成的和式 .211 ++++=∑∞ =n n n a a a a (2) 数列与级数有着密不可分的关系。给定一个无穷级数（2），它唯一确定了一个无穷数列

,,,21 S S 其中.,2,1,21 =+++=n a a a S n n 反过来，给定一个无穷数列（1），它也唯一地确定了一个无穷级数 ∑∞ =1 n n b ，这里.,2,1,,111 =-==-n a a b a b n n n 并且，无穷级数的和就是相应的无穷是咧的极限。因此，无穷数列与无穷级数是可以相互转化的。实验步骤 1. 用如下语句作图： FibShow[n_Integer]:= Module[ {t={},i}, For[i=1,i<=n,i++,AppendTo[t,{i,Fibonacci[i]}]]; ListPlot[t,PlotJoined-> True] ] FibShow[N] 2. 用如下语句计算： FibFit[n_Integer]:= Module[ {t={},i}, For[i=1,i<=n,i++,AppendTo[t,{i,Log[Fibonacci[i]]}]]; Fit[t,{1,x},x] ]

华工数学实验报告

华工数学实验报告篇一：华工数学实验报告微分方程《数学实验》报告学院：电子信息学院专业班级：信息工程电联班学号：姓名：实验名称：微分方程实验日期：XX/04/19 1．实验目的了解求微分方程解析解的方法了解求微分方程数值解的方法了解 dsolve,ode45 指令的使用方法 2．实验任务 1.用dsolve函数求解下列微分方程 ?y??(x)?y?(x)?2y(x)（2）? ?y(0)?1,y(0)?0? 2. 我辑私雷达发现，距离d处有一走私船正以匀速a 沿直线行驶，缉私舰立即以最大速度（匀速v）追赶。若用雷达进行跟踪，保持船的瞬时速度方向始终指向走私船，则辑私舰的运动轨迹是怎么的？是否能够追上走私船？如果能追上，需要多长时间？ M0 3．实验过程

3.1实验原理 dsolve(‘equation’,’condition’,’v’) （1） equation是方程式，condition是条件，v是自变量（缺省为t）（2）若不带条件，则解中带积分常数（3）如果没有显示解，则系统尝试给出隐式解（4）如果无隐式解，则返回空符号。以S0为原点建立坐标系。设缉私船出发的起点坐标为，根（x0,y0）据题意x02?y02?d2，经过时间t，走私船到达S(at,0)，缉私船到达M(x,y)，追赶时，缉私船总是向走私船所在的位置追赶，设在t+dt时刻，缉私船到达M'(x?dx,y?dy)，则M,M’,S三点一图2 dt时刻追击图由图可知，即 dy0?y? dxat?x(1) ?ydx?at?x dy(2) 此即缉私船的追辑模型。方程（2）两边对y求导，得 d2xdt?y2?a dydy(3) 又因为缉私船的速度恒为v，因此即 dy?dt?dy??dx?v2?????? ?dt??dt?22(4) (5)

数学实验报告3

实验目的：熟悉差分方程的求解，以及相关金融问题的数学建模方法。实验内容： 1、 2、小李夫妇曾经准备申请商业贷款10万元用于购置住房，每月还款880.66元，25年还清。房产商介绍的一家金融机构提出：贷款10万元，每半月还款440.33元, 22年还清, 不过由于中介费手续费等原因，贷款时要预付4000元。小李考虑，虽然预付费用不少，可是减少三年还款期意味着减少还款近3万2千元，而每月多跑一趟，那不算什么.这机构的条件似乎还是蛮优惠的。试通过计算两种贷款的利率水平，比较那种贷款更优惠。 3、试通过计算两种贷款的利率水平，比较那种贷款更优惠。试比较两种提前还款方式的优劣（附加）所谓提前还贷是指借款人在保证按月按额偿还个人住房贷款本息的基础上，提前偿还部分或全部购房借款的一种经济行为。每次提前还款后，相应冲减余贷款本金。银行根据尚未归还的贷款本金重新计算借款人的月均还款额，直至贷款本息全部还清。重新计算月还款金额有两种方式： A 、提前还款额冲抵最后月份的本金，每月的还款额度不变，还款时间缩短； B 、提前还款额冲抵本金后，将剩余的贷款重新计算月还款额减少，还款时间不变。例如，谢先生申请公积金贷款30万元，贷款期限为20年，在正常按月还了5年贷款后，谢先生决定提前还5万元本金，然后再继续按月还款。试比较两种提前还款方式的优劣？实验要求：撰写实验报告写出试验过程中所使用的Mathematica 程序或语句和计算结果第一题：第一年： In [1]≔A0=10000;k =1∗12;r =5.31100⁄12⁄ m =(A0∗(1+r )k ∗r )((1+r )k −1)⁄ A =m ∗k Out[2]=0.004425 Out[3]=857.496 Out[4]=10290.0 第二年和第三年：

MATLAB数学实验报告

Matlab 数学实验报告

一、实验目的通过以下四组实验,熟悉MATLAB的编程技巧,学会运用 MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x=λsin(πx),做出相应的 kk+1Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9

x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150 plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1));

end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图 2.2实验题目二 2.2.1实验问题某农夫有一个半径10米的圆形牛栏，长满了草。他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？问题分析 2.2.2

高数实验报告

高等数学数学实验报告实验一一、实验题目观察数列极限二、实验目的和意义通过作图观察数列极限：n趋向于无穷时，（1+1/n）^n 三、计算公式四、程序设计 data = Table[(1 + 1/i)^i, {i, 30}]; ListPlot[data, PlotRange -> {2, 3}, PlotStyle -> PointSize[0.018]] 五、程序运行结果六、结果的讨论和分析通过图像观察出数列趋向于重要极限e 实验二一、实验题目一元函数图形及其性态

二、实验目的和意义制作函数y=sincx的图形动画，并观察参数c对函数图形的影响三、计算公式请写出在程序中所需要的计算公式。比如定积分的数值计算中，如用梯形法计算的，请描述梯形法的公式。四、程序设计 Animate[Plot[Sin[c x], {x, 0, 10}, PlotRange -> {-1, 1}], {c, -1, 4, 1/3}] 五、程序运行结果 0.5 1.0 0.5 1.0

1.0 0.5 六、结果的讨论和分析通过图像观察出常数c 影响y=sincx 的周期和频率，函数周期为2Pi/c,频率为c/2Pi. 实验三一、实验题目泰勒公式与函数逼近二、实验目的和意义对y=cosx 分别在[-Pi,Pi],[-2Pi,2Pi]上进行n 阶泰勒展开三、计算公式请写出在程序中所需要的计算公式。比如定积分的数值计算中，如用梯形法计算的，请描述梯形法的公式。四、程序设计（1）t = Table[Normal[Series[Cos[x], {x, 0, i}]], {i, 0, 12, 2}]; PrependTo[t, Cos[x]]; Plot[Evaluate[t], {x, -Pi, Pi}] （2）For[i = 0, i <= 10, a = Normal[Series[Cos[x], {x, 0, i}]]; Plot[{a, Cos[x]}, {x, -Pi, Pi}, PlotStyle -> {RGBColor[0, 0, 1], RGBColor[1, 0, 0]}]; i = i + 2] (3)For[ =6, ≤16, =Normal[Series[Cos[ ],{ ,0, ,Cos[ ]},{ ,−2Pi,2Pi}, PlotStyle→{RGBC olor[0,0,1],RGBColor[1,0,0]}]; = +2] (4)tt[x0_]:=Normal[Series[Cos[ ],{ ,x0,6}]];gs0=tt[0];gs3=tt[3];gs6=tt[6];Plot[{Cos [ ],gs0,gs3,gs6},{ ,−3Pi,3Pi},PlotRange→{−2,2},PlotStyle→{RGBColor[0,0,1],RGB Color[1,0,1],RGBColor[1,0,0],RGBColor[0,1,0]}] (5) f[x_]:=Sin[x 2]; a=0;b=0.5Pi;m2=N[f''[0.0000635627]];dalta=10^(-4);n0=90; t[n_]:=(b-a)/n×((f[a]+f[b])/2+Sum[f[a+i×(b-a)/n],{i,1,n-1}]); Do[Print[n," ",N[t[n]]]; If[(b-a)^3/(12n^2)×m2

matlab数学实验报告

MATLAB 数学实验报告指导老师：班级：小组成员：时间：201_/_/_

Matlab第二次实验报告小组成员： 1题目：实验四;MATLAB选择结构与应用实验目的：掌握if选择结构与程序流程控制;重点掌握break;return;pause语句的应用.. 问题：问题1：验证“哥德巴赫猜想”;即：任何一个正偶数n>=6均可表示为两个质数的和..要求编制一个函数程序;输入一个正偶数;返回两个质数的和.. 问题分析：由用户输入一个大于6的偶数;由input语句实现..由if判断语句判断是否输入的数据符合条件..再引用质数判断函数来找出两个质数;再向屏幕输出两个质数即可.. 编程：function z1;z2=geden; n=input'please input n' if n<6

disp'data error'; return end if modn;2==0 for i=2:n/2 k=0; for j=2:sqrti if modi;j==0 k=k+1; end end for j=2:sqrtn-i if modn-i;j==0 k=k+1; end end if k==0 fprintf'two numbers are' fprintf'%.0f;%.0f';i;n-i break end end

end 结果分析：如上图;用户输入了大于6的偶数返回两个质数5和31;通过不断试验;即可验证哥德巴赫猜想.. 纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰;更快的解决问题.. 2题目：实验四;MATLAB选择结构与应用实验目的：用matlab联系生活实际;解决一些生活中常见的实际问题..

数学实验报告

实验过程记录（含基本步骤、主要程序清单及异常情况记录等）：第一题选择初速度v=0.6km/s,发射角a=45° X轴方向运动为x=cos a*v*t Y轴方向运动为y=sin a*v*t-1/2*g*t2 统一单位将0.6km/s化为600m/s 将数据代入利用函数做出运动轨迹，函数式为 8000 6000 4000 2000 5000100001500020000250003000035000 第二题确定速度为320m/s，求最佳角度使得轨迹与X轴交点为(10000,0) 先假定发射角为��/4 作图 ParametricPlot[{Cos[Pi/4]*320*t,Sin[Pi/4]*320*t-4.9*t ^2},{t,0,47},AspectRatio→Automatic] 2500 2000 1500 1000 500 200040006000800010000 进行调整角度调整为��/3.5 作图 ParametricPlot[{Cos[Pi/3.5]*320*t,Sin[Pi/3.5]*320*t-4 .9*t^2},{t,0,52},AspectRatio→Automatic] 3000 2500 2000 1500 1000 500 200040006000800010000

继续进行不断地调整，发现当发射角度为π/3.7时，落点十分接近(10000,0)点作图如下 2000 4000 6000 8000 10000 500 10001500200025003000 因此可以确定最合适的发射角就在π/3.7附近，此时可以利用FindRoot 函数找出准确值首先需要对已知式做等量变换： ∵X=cos a*v*t ∴t=x/(cos a*v) 将上式代入y=sina*v*t-1/2*g*t 2 中可得到 Y=tana*x-1/2*g*(x/(cosa*v))2 将y=0, x=10000, g=9.8, v=320代入利用FindRoot 函数求解a 的范围在π/3.7附近的a 的值：得出将这个值由弧度制化为360度制 a=53.4285° ∴最佳发射角为53.4285° 第三题由第二题的320m/s 起步进行研究 1.首先研究速度增大运用与第二题相似的研究方法，先大致计算符合要求的角度 (1)V=350m/s 时，最佳发射角为π/6.8: 2000 4000 6000 8000 10000 200 40060080010001200 (2)V=400m/s 时，最佳发射角为π/9.5: 2000 4000 6000 8000 10000 200 400600800

数学实验报告

《数学实验》报告题目：根据数值积分计算方法计算山东省面积学生姓名：学号：专业班级：机械工程17-1班 2019年4月15日

一、问题背景与提出图1是从百度地图中截取的山东省地图，试根据前面数值积分计算方法，计算山东省面积。图1 二、实验目的 1、学会运用matlab解决一些简单的数学应用问题。 2、学会运用matlab建立数学模型。 3、学会运用一些常见的数值积分计算方法结算实际问题，并了解其实际意义，建立积分模型。三、实验原理与数学模型将积分区间[a , b] n等分，每个区间宽度均为h = (b - a) / n , h称为积分步长。记a = x0< x1<…

, h = 如果将二者求平均值，则每个小区间上的小矩形变为小梯形，整个区间上的值变为：将山东省边界上的点反映在坐标化，运用梯形公式积分计算得山东省的面积。四、实验内容（要点） 1、将山东省的地图区域在matlab中画出。 2、在坐标系上运用积分方法将所求区域的面积求出。 3、通过比例尺将山东省的实际面积求出。五、实验过程记录（含基本步骤、主要程序清单及异常情况记录等） 1、在百度地图中标识出山东省的区域范围，标明对应的比例：图2 2、取出所截取图片中山东的边界的坐标，即将边界坐标化：

（1）运用imread函数和imshow函数导入山东省的区域图片。代码：运行结果：图3 （2）运用ginput函数，将边界的坐标点取出，即坐标化，并将x,y坐标分别存于x.txt和y.txt文本文件中。代码：

数学实验报告

高等数学数学实验报告实验人员：院〔系〕 ____**_______ 实验地点：计算机中心机房实验一空间曲线与曲面的绘制一、实验题目：〔实验习题1-2〕利用参数方程作图，做出由以下曲面所围成的立体图形： (1) x y x y x z =+--=2 222,1及*Oy 平面； (2) 01,=-+=y x xy z 及.0=z 二、实验目的和意义 1、利用数学软件Mathematica 绘制三维图形来观察空间曲线和空间曲面图形的特点，以加强几何的直观性。 2、学会用Mathematica 绘制空间立体图形。三、程序设计空间曲面的绘制作参数方程],[],,[,),(),() ,(max min max min v v v u u v u z z v u y y v u x x ∈∈⎪ ⎩⎪ ⎨⎧===所确定的曲面图形的 Mathematica 命令为： ParametricPlot3D[{*[u,v],y[u,v],z[u,v]},{u,umin,uma*}, {v,vmin,vma*},选项] (1) (2)

四、程序运行结果 (1) 〔2〕五、结果的讨论和分析 1、通过参数方程的方法做出的图形，可以比拟完整的显示出空间中的曲面和立体图形。 2、可以通过mathematica软件作出多重积分的积分区域，使积分能够较直观的被观察。 3、从(1)中的实验结果可以看出，所围成的立体图形是球面和圆柱面所围成的立体空间。 4、从(2)中的实验结果可以看出围成的立体图形的上面曲面的方程是z=，下底面的方程是z=0，右边的平面是0 xy x。 +y 1= - 实验一空间曲线与曲面的绘制一、实验题目：〔实验习题1-3〕观察二次曲面族kxy z+ x y =2 + 2的图形。特别注意确定k的这样一些值，当k经过这些值时，曲面从一种类型变成了另一种类型。

线性代数实验报告汇总

数学实验报告题目第一次实验题目一、实验目的 1MATLAB 的矩阵初等运算；．熟悉2 ．掌握求矩阵的秩、逆、化最简阶梯形的命令；3MABLAB 求解线性方程组．会用二、问题求解和程序设计流程344?221????????MATLABA1 B、，已知命令窗口中建立．，在 320B???50??3A????????112?153????矩阵并对其进行以下操作：(1) A 的行列式的值计算矩阵?)?Adet((2) 分别计算下列各式：、和、、、、 B?A.T112??B?BA?2A ABABAA ：解（1）编写程序如下： A=[4 -2 2;-3 0 5;1 5 3]; B=[1 3 4;-2 0 -3;2 -1 1]; a=det(A) 运行结果： a = -158 （2）编写程序如下： C=2*A-B D=A*B E=A.*B F=A/B G=A\B H=A*A K=A' 运行结果： C = 7 -7 0 -4 0 13 线性代数实验报告

0 11 5 D = 12 10 24 7 -14 -7 -3 0 -8 E = 4 -6 8 6 0 -15 2 -5 3 F = 0 0 2.0000 -2.7143 -8.0000 -8.1429 2.4286 3.0000 2.2857 G = 0.4873 0.4114 1.0000 0.3671 -0.4304 0 -0.1076 0.2468 0 H = 24 2 4 -7 31 9 -8 13 36 K = 4 -3 1 -2 0 5 2 5 3 2 MATLABrankinv 求下列矩阵的秩：中分别利用矩阵的初等变换及函数．在、函数线性代数实验报告