数学实验与数学软件实验报告一
mathematica 数学实验报告
数学实验报告实验一数学与统计学院信息与计算科学(1)班郝玉霞201171020107数学实验一一、实验名：微积分基础二、实验目的：学习使用Mathematica的一些基本功能来验证或观察得出微积分学的几个基本理论。
三、实验环境：学校机房，工具：计算机,软件：Mathematica。
四、实验的基本理论和方法：利用Mathematica作图来验证高中数学知识与大学数学内容。
五、实验的内容和步骤及结果内容一、验证定积分dttsx⎰=11与自然对数xb ln=是相等的。
步骤1、作积分dttsx⎰=11的图象；语句：S[x_]:=NIntegrate[1/t,{t,1,x}]Plot[S[x],{x,0.1,10}]实验结果如下：21图1dttsx⎰=11的图象步骤2、作自然对数xb ln=的图象语句：Plot[Log[x],{x,0.1,10}] 实验结果如下：2 1图2xb ln=的图象步骤3、在同一坐标系下作以上两函数的图象语句：Plot[{Log[x],S[x]},{x,0.1,10}] 实验结果如下：21图3dttsx⎰=11和xb ln=的图象内容二、观察级数与无穷乘积的一些基本规律。
（1）在同一坐标系里作出函数xy sin=和它的Taylor展开式的前几项构成的多项式函数3!3xxy-=，!5!353xxxy+-=，⋅⋅⋅的图象，观察这些多项式函数的图象向xy sin=的图像逼近的情况。
语句1：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,2]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}] 实验结果如下：64242图4x y sin =和它的二阶Taylor 展开式的图象语句2：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,3]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,1,1]}] 实验结果如下：642321图5x y sin =和它的三阶Taylor 展开式的图象语句3：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,4]},{x,-2Pi,2Pi},PlotStyle->{RGB[0,1,0]}] 实验结果如下：642321图6x y sin =和它的四阶Taylor 展开式的图象语句4：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}]Plot[{Sin[x],s[x,5]},{x,-2Pi,2Pi},PlotStyle->{RGB[1,0,0]}] 实验结果如下：642321图7x y sin =和它的五阶Taylor 展开式的图象语句5：s[x_,n_]:=Sum[(-1)^(k-1)x^(2k-1)/((2k-1)!),{k,1,n}] Plot[{Sin[x],s[x,2],s[x,3],s[x,4],s[x,5] },{x,-2Pi,2Pi}] 实验结果如下： 6422图8xy sin=和它的二、三、四、五阶Taylor展开式的图象（2）分别取n=10，20，100，画出函数xkkynk)12sin(1211--=∑=在区间[-3π，3π]上的图像，当n→∞时，这个函数趋向于什么函数？语句1：f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]Plot[f[x,10],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}]实验结果如下：6420.5图9 n=10时，xkkynk)12sin(1211--=∑=的图像语句2：f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]Plot[f[x,20],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}] 实验结果如下：6420.5图10 n=20时，xk k y nk )12sin(1211--=∑=的图像语句3：f[x_,n_]:=Sum[Sin[k*x]/k,{k,1,n,2}]Plot[f[x,100],{x,-2Pi,2Pi},PlotStyle->{RGB[0,0,1]}] 实验结果如下：6420.5图11 n=100时，xk k y nk )12sin(1211--=∑=的图像（3）分别取5,15,100,，在同一坐标系里作出函数x x f sin )(=与∏=-⋅=nk k x x x p 1222)1()(π在区间[-2π，2π]上的图像。
[数学软件及应用(Lingo)实验报告范文]lingo实验报告范文心得
![[数学软件及应用(Lingo)实验报告范文]lingo实验报告范文心得](https://img.taocdn.com/s3/m/b2e98c482f3f5727a5e9856a561252d380eb20b2.png)
[数学软件及应用(Lingo)实验报告范文]lingo实验报告范文心得2022~2022学年第二学期短学期《数学软件及应用(Lingo)》实验报告班级数学131班姓名张金库学号成绩实验名称奶制品的生产与销售方案的制定完成日期：2022年9月3日实验名称：奶制品的生产与销售方案的制定二、实验目的及任务了解并掌握LINGO的使用方法、功能与应用；学会利用LINGO去解决实际中的优化问题。
三、实验内容问题一奶制品加工厂用牛奶生产，两种奶制品，1桶牛奶可以在甲类设备上用12h加工成3kg，或者在乙类设备上用8h加工成4kg。
根据市场的需求，生产,全部能售出，且每千克获利24元，每千克获利16元。
现在现在加工场每天能的到50桶牛奶的供给，每天正式工人总的劳动时间为480h，并且甲类设备每天至多能加工100kg，乙类设备的加工能力没有限制。
为增加工厂的利益，开发奶制品的深加工技术：用2h和3元加工费，可将1kg加工成0.8kg高级奶制品，也可将1kg加工成0.75kg高级奶制品，每千克能获利44元，每千克能获利32元。
试为该工厂制订一个生产销售方案，使每天的净利润最大，并讨论以下问题：假设投资30元可以增加供给1桶牛奶，投资3元可以增加1h的劳动时间，应否做这些投资？假设每天投资150，可以赚回多少？每千克高级奶制品，的获利经常有10%的波动，对制订的生产销售方案有无影响？假设每千克获利下降10%，方案应该变化吗？假设公司已经签订了每天销售10kg的合同并且必须满足，该合同对公司的利润有什么影响？问题分析要求制定生产销售方案，决策变量可以先取作每天用多少桶牛奶生产，，再添上用多少千克加工，用多少千克加工，但是问题要分析，的获利对生产销售方案的影响，所以决策变量取作，，，每天的销售量更为方便。
目标函数是工厂每天的净利润——，，，的获利之和扣除深加工费用。
根本模型决策变量：设每天销售kg，kg，kg，kg，用kg加工，用kg加工。
数学的上机实验报告
实验题目：线性代数求解方程组一、实验目的1. 理解线性代数中方程组的求解方法。
2. 掌握利用计算机求解线性方程组的算法。
3. 熟悉数学软件（如MATLAB、Python等）在数学问题中的应用。
二、实验内容本次实验主要利用数学软件求解线性方程组。
线性方程组是线性代数中的一个基本问题，其求解方法有很多种，如高斯消元法、矩阵求逆法等。
本实验以高斯消元法为例，利用MATLAB软件求解线性方程组。
三、实验步骤1. 编写高斯消元法算法程序。
2. 输入方程组的系数矩阵和常数项。
3. 调用程序求解方程组。
4. 输出解向量。
四、实验代码及分析1. 高斯消元法算法程序```matlabfunction x = gaussElimination(A, b)[n, m] = size(A);assert(n == m, 'The matrix A must be square.');assert(n == length(b), 'The length of b must be equal to the number of rows in A.');% 初始化解向量x = zeros(n, 1);% 高斯消元for i = 1:n-1% 寻找最大元素[~, maxIdx] = max(abs(A(i:n, i)));maxIdx = maxIdx + i - 1;% 交换行A([i, maxIdx], :) = A([maxIdx, i], :);b([i, maxIdx]) = b([maxIdx, i]);% 消元for j = i+1:nfactor = A(j, i) / A(i, i);A(j, i:n) = A(j, i:n) - factor A(i, i:n); b(j) = b(j) - factor b(i);endend% 回代求解for i = n:-1:1x(i) = (b(i) - A(i, i+1:n) x(i+1:n)) / A(i, i); endend```2. 输入方程组的系数矩阵和常数项```matlabA = [2, 1, -1; 1, 2, 1; -1, 1, 2];b = [8; 5; 2];```3. 调用程序求解方程组```matlabx = gaussElimination(A, b);```4. 输出解向量```matlabdisp('解向量为：');disp(x);```五、实验结果与分析实验结果：```解向量为：2-13```实验分析：通过高斯消元法，我们成功求解了给定的线性方程组。
《数学实验》实验报告——用MATLAB绘制sin曲线
试验步骤（根据问题分析及试验目的所计划的试验步骤）：
（1）设计实验程序，程序如下：
（2）试运行实验程序，逐步修改程序。
试验过程（含详细试验步骤、程序清单及异常情况记录等）
（1）将编写的程序在软件中试运行；
（2）逐步修改指令，排除出现的各种异常情况；
（3）运行结果如下：
结果分析：
运行的结果能达到实验的目的和要求，通过在同一幅图上的（-pi，pi）区间，可绘制出用0.5的间隔的sinx的红色曲线，用0.1的间隔的sin(x+0.5)的绿色曲线，用0.01的间隔的sin(x+1)的蓝色曲线，并给图的x和y轴添加标注。
x1=-pi:0.5:pi
y1=sin(x1);
x2=-pi:0.1:pi
y2=sin(x2+0.5);
x3=-pi:0.01:pi
y3=sin(x3+1);
plot(x1,y1,'r',x2,y2,'g',x3,y3,'b')
xlabel('x1,x2,x3')
ylabel('sin(x1),sin(x2+0.5),sin(x3+1)')
有的工具绘制起来比较繁琐所达到的结果也不尽如人意所以就需要寻找一种比较容易操作及效果比较好的工具来完成sin曲线的绘实验报告
班级
****
学号
****
姓名
数学实验报告 (1)
（1）参数方程：z=2^2^/2^2^sin y x y x ++（-8<=x<=8,-8<=y<=8) （2）程序：[X,Y]=meshgrid(-8::8);r=sqrt(x.^2+y.^2)+eps;Z=sin(r)./r;Mesh(x,y,z)Axis square(3)程序的输出结果：3：球面,椭球面,双叶双曲面,单叶双曲面1球面: (4):参数方程：⎪⎩⎪⎨⎧===ϕθϕθϕcos *sin *sin *cos *sin *R z R y R x 0π<=θ<2* 0<=ϕ<π (5)程序：u=[0:pi/60:2*pi];v=[0:pi/60:pi];[U,V]=meshgrid(u,v);R=3;X=R*sin(v).*cos(u);Y=R*sin(v).*sin(u);Z=R*cos(v);Surf(x,y,z);axis equal;(3)程序输出结果：2椭球面： （1）参数方程：⎪⎩⎪⎨⎧===ϕθϕθϕcos *sin *sin *cos *sin *c z b y a x 0<=θ<2*π 0<=ϕ<=π （2）程序：ezsurf(‘3*sin(u)*cos(v) ,’3*sin(u)*sin(v)’,’1*cos(u)’,[0,pi,0,2*pi]);(3)程序的输出结果：3单叶双曲面： （1）参数方程：⎪⎩⎪⎨⎧===ϕθϕθϕtan sin *sec *cos *sec *z a y a x 0<=θ<2*π -π/2<ϕ<π/2 （2）程序：ezsurf(‘3*sec(u)*cos(v),’3*sec(u)*sin(v)’,’5*tan(u)’,[-pi/2,pi/2,0,2*pi]);axis auto（3）输出程序结果：4双叶双曲面： （1）参数方程：⎪⎩⎪⎨⎧===ϕθϕθϕsec *sin *tan *cos *tan *c z b y a x 0<=θ<2*π -π<ϕ<3*π/2，ϕ≠π/2（2）程序：ezsurf(‘3*tan(u)*cos(v)’,’3*tan(u)*sin(v)’,’5*sec(u)’,[-pi/2,3*pi/2,0,2*pi]);axis auto（4） （3）输出程序结果：抛物螺线： （1）参数方程：⎪⎩⎪⎨⎧===2^*sin **cos **t c z t t b y t t a x 0<T<+∞ （2）程序：ezplot3(‘2*t*cos(t)’,’2*t*sin(t)’,’t.^2/3’,[0,50]);（3）输出程序结果：（5）马鞍面： （1）参数方程：z=x^2/9-y^2/4 （-25<=x<=25,-25<=y<=25)（2）程序：[X,Y]=meshgrid(-25:1:25);Z=X.^2/9-Y.^2/4;Surf(X,Y,Z)Title(‘马鞍面’)grid off（3）输出程序结果：（6）黎曼函数：（1）程序：n=100;x=[];y=[];k=1;for q=2:nfor p=1:q-1if gcd(q,p)==1 %利用函数gcd(m,n)可求m和n的最大公约数x(k)=p/q;y(k)=1/q;k=k+1;endendendplot(x,y,’.b’); axis([0,1,0,1])（2）程序输出结果：。
数学实验实习报告
数学实验实习报告一、引言数学实验实习是数学专业学生在实践中提高数学建模能力、动手能力以及科学研究能力的重要环节。
本次实习报告旨在总结和分析实习过程中的实验内容、方法和结果，以及对实习的感悟和体会。
二、实验目的本次实习的目的是通过数学建模的方法，解决实际问题，培养学生的数学应用能力和创新思维。
具体实验目的如下：1. 掌握数学建模的基本原理和方法；2. 学习和运用数学软件和工具，如MATLAB、Mathematica等；3. 分析和解决实际问题，并给出科学合理的结论；4. 提升数据处理和实验报告撰写的能力。
三、实验内容本次实习的主题是“市场调研数据分析与预测”。
在实验过程中，我们使用了一系列数学模型和算法，对给定的市场调研数据进行了分析和预测，以期给公司提供决策支持。
具体的实验步骤如下：1. 数据收集：我们收集了与市场调研相关的数据，包括产品销售额、消费者满意度、竞争对手信息等。
2. 数据预处理：对收集到的数据进行清洗和整理，剔除异常值和缺失数据。
3. 数据分析：使用统计学和数据挖掘的方法，对数据进行分析和探索，包括描述统计、相关性分析、聚类分析等。
4. 模型构建：根据实际问题的要求，选择适当的数学模型建立预测模型，如线性回归、时间序列分析等。
5. 模型评估：对建立的模型进行评估，检验模型的准确性和稳定性，并提出改进意见。
6. 结果展示：根据模型分析结果，绘制相关图表，给出对市场趋势和销售预测的结论。
四、实验结果和讨论通过对市场调研数据的分析和预测，我们得到了以下结论：1. 市场趋势分析：根据历史数据和统计模型，预测市场的发展趋势，包括市场规模、增长率等。
2. 销售预测：通过建立销售预测模型，对未来一段时间内的销售额进行预测，为公司制定销售策略提供参考。
3. 消费者满意度分析：通过对消费者满意度调查数据的分析，找出关键因素和改进方向，提高产品竞争力。
4. 竞争对手分析：通过分析竞争对手的市场份额和策略，为公司制定竞争策略提供依据。
数学实验报告的总结(3篇)
第1篇一、实验背景随着科技的不断发展，数学实验在各个领域中的应用越来越广泛。
数学实验作为一种以计算机为工具，通过模拟、计算和验证等方法，对数学理论进行实践探索和研究的方法，已经成为数学研究的重要手段。
本次实验旨在通过数学实验，加深对数学理论的理解，提高数学应用能力，培养创新意识和团队协作精神。
二、实验目的1. 熟悉数学实验的基本方法，掌握数学实验的基本步骤。
2. 通过实验，加深对数学理论的理解，提高数学应用能力。
3. 培养创新意识和团队协作精神，提高自身综合素质。
三、实验内容本次实验主要包括以下内容：1. 实验一：线性方程组的求解通过编写程序，实现线性方程组的直接法、迭代法等求解方法，并对比分析各种方法的优缺点。
2. 实验二：矩阵运算实现矩阵的加法、减法、乘法、转置等基本运算，以及求逆矩阵、特征值和特征向量等高级运算。
3. 实验三：数值积分通过编写程序，实现定积分、变积分、高斯积分等数值积分方法，并分析各种方法的误差和适用范围。
4. 实验四：常微分方程的数值解法实现欧拉法、龙格-库塔法等常微分方程的数值解法，并对比分析各种方法的稳定性、精度和适用范围。
四、实验过程1. 确定实验内容，明确实验目的。
2. 设计实验方案，包括实验步骤、算法选择、数据准备等。
3. 编写实验程序，实现实验方案。
4. 运行实验程序，收集实验数据。
5. 分析实验数据，得出实验结论。
6. 撰写实验报告，总结实验过程和结果。
五、实验结果与分析1. 实验一：线性方程组的求解通过实验，验证了直接法和迭代法在求解线性方程组时的有效性。
直接法在求解大规模线性方程组时具有较好的性能，而迭代法在求解稀疏线性方程组时具有较好的性能。
2. 实验二：矩阵运算实验结果表明，矩阵运算的程序实现具有较高的精度和效率。
在实际应用中，可以根据具体需求选择合适的矩阵运算方法。
3. 实验三：数值积分通过实验，验证了各种数值积分方法的有效性。
高斯积分具有较高的精度，但在求解复杂函数时，需要调整积分区间和节点。
重庆大学-数学实验-1MATLAB软件入门
其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为153=1 计算程序
for a=1:1:9 for b=0:1:9 for c=0:1:9 A=a*100+b*10+c; if A==a^3+b^3+c^3 disp(A)
end end end end 最终结果： 153 370 371 407 分析：此题使我进一步熟悉了循环语句for、while等的基本格式和应用方法， 得注意的是每一个循环语句或条件语句的出现，必定要对应的出现一个end，这 点很容易忽视。 6．编写函数M-文件sq.m：用迭代法求的值。求平方根的迭代公式为 程序： function f=sq(a) a>=0 x=eps; y=1/2*(x+a/x); while (abs(x-y)>=10^(-5)) x=y y=1/2*(x+a/x) end f=x else disp('their exist errors') end 运行结果 调用sq(16)，输出结果是4.0000 分析：此题要求掌握利用循环语句进行迭代，实现开方方程的求解，另外x的取
分析在表达乘方的时候要用“y1=x.^2”的形式，在表达x分量的时候可以用 linspace和直接用 生成数组的形式。主要函数：plot：生成图片；hold on和old off用来在同一幅图 描图；text：用程序对图片加标注；gtext：用鼠标在图中加字符；title：用于对 加标题。
齐齐哈尔大学《数学实验》关于matlab软件设计结课报告
齐齐哈尔大学Matlab程序设计结课报告学院：理学院班级：信科152姓名：刘言学号：2015035000数学实验报告姓名：刘言班级:信科152 实验日期：2016.11.29实验名称:MATLAB实验汇总报告一、实验准备部分(一)实验基础：1.掌握基本的matlab编程基础，能够解决简单的数学问题2.GUI工具与图形用户界面3.数据可视化4.数学问题5.文件读取6.Matlab编译器7.本报告基于版本：matlab2016a进行编写，如遇不兼容问题，请放于此版本调试（二）实验基本理论与方法1、绘图函数polt实现一次函数图像、正多边形、二次函数的绘图功能:plot(x,y);绘图函数axis([-1.1,1.1,-1.1,1.1]);图像坐标2、添加GUI控件实现可视化操作：3、一元二次方程求解：一元二次方程的求根公式：将一元二次方程ax2＋bx＋c=0(a≠0)进行配方，当b2－4ac≥0时的根为5、矩阵基本计算：求转置、逆矩阵、加减运算、左除右除、秩等相关计算6、两个简单的数学游戏（三）实验设计（可以用流程图、功能图或其他方式表示，选择一种即可）1、程序流程图用户界面二、实验操作部分程序设计点击画图 用户界面菜单栏功能一元二次函数一元函数求解正比例函数多边形绘制MCC 编译为.EXE工具帮助二元函数求解矩阵求解数学游戏速算24智商测试（一）实验过程：1、创建GUI界面：button按钮、edit可编辑文本框、text静态文本框、Aitx显示图像、按钮组等一些列组件相互联系搭建成一个相对完美的用户操作系统，2、实现按钮功能化：例如：按钮控件的制作：在GUI界面新建一个‘按钮’，然后双击打开‘String’对文本进行编辑，即可改变按钮显示的文字，右击到按钮标签Tag的回调函数Callback，键入相应的代码或者.m文件，实现相应功能。
然后就是按钮的外观优化：字体颜色，背景色，字体大小等做细节处理。
数学实验报告 (1)
数学实验报告四实验项目名称MATLAB基础所属课程名称数学实验实验日期2012-10-10姓名（学号）周星（2010190135）成绩数学与计算科学学院数学实验室一、 实验目的1. 掌握使用plot 绘制二维图形；2. 掌握分段函数绘制；3. 掌握绘制图形的辅助操作二、 实验环境（使用软件）MATLAB V6.5三、 实验内容1. 设23sin (0.5)cos 1x y x x =++，在0~2x π=区间取等间隔101个点，绘制函数的曲线。
2. 在02x π≤≤区间内，绘制曲线0.52sin(2)x y e x π-=。
3. 生成10000×1的正态随机数矩阵，绘制直方图，要求30×1个长条。
4. 绘制曲线2cos(3)sin x t t t y t t ππ=⎧-≤≤⎨=⎩5. 已知21y x =，2cos(2)y x =，312y y y =⨯，[10,10]x ∈-完成下列操作：（1） 在同一坐标系下用不同颜色和线型绘制三条曲线，并在右上角给加入曲线说明；6. 绘制分段函数曲线04246()568218x x f x x x x ≤<≤<⎪⎪=⎨-≤<⎪⎪≥⎪⎩ 要求：（1）设置坐标轴范围为：横坐标范围为[0,10]，纵坐标范围为[0,2.5]； （2）给图形加上标题“分段函数曲线”； （3）给X,Y 轴分别添加说明“Variable X ”和”Variable Y ” （4）用鼠标在给分段曲线每段添加图形说明四、 实验解答1.解：x=linspace(0,2*pi,100);>> y=(0.5+3*sin(x)/(1+x.^2))*cos(x);>> plot(x,y,'o-')2．解：x=linspace(0,2*pi);>> y=2*exp(-0.5*x).*sin(2*pi*x); >> plot(x,y)>> hist(x,30);4.解：>> t=linspace(-pi,pi); >> x=t.*cos(3*t);>> y=t.*sin(t).^2;plot(x,y)>> y1=x.^2;>> y2=cos(2*x);>> y3=y1.*y2;>> plot(x,y1,'b-');>> hold on;>> plot(x,y2,'g:');>> plot(x,y3,'r-.');>> hold off;>> legend('x.^2','cos(2*x)','y1.*y2');x=linspace(-10,10);6.解：x=linspace(0,4);>> plot(x,sqrt(x));>> hold on;>> x=linspace(4,6);>> plot(x,2);>> x=linspace(6,8);>> plot(x,5-x/2);>> x=linspace(8,10);>> plot(x,1);>> hold off;>> xlabel('Variable X');>> ylabel('Variable Y ');>> title('分段函数曲线');>> axis([0,10,0,2.5]);>>gtext('y=sqrt(x)');gtext('y=2');gtext('y=5-x/2');gtext('y=1');。
数学实验报告之一——Rossler方程
数学实验报告之一——Rossler方程自动化学院2001级刘和松一 .实验目的掌握用MATLAB软件求解微分方程的基本方法，掌握求微分方程数值解的欧拉方法，了解龙格——库塔方法的思想；了解二次迭代过程中三种图形表示的含义，观察分枝与混沌现象，学会分析迭代对参数的敏感性。
Rossler方程组是非线性动力学中一个非常著名的方程，在理论和实际中都有非常重要的价值。
但是，由于它没有精确的解析解，我们只能通过数值方法对其进行求解。
二.试验内容——Rossler方程Rossler方程组当固定参数b=2，c=4时，试讨论随参数a由小到大变化（如a∈(0,0.65)）时方程解的变化情况，并且画出空间曲线图形，观察空间曲线是否形成混沌状。
三．问题分析对于Rossler方程组这样的非线性微分方程组，不可能求得其解析解。
利用MatLab可以用数值解法进行模拟。
在物理意义上，t一般表示时间，x，y，z通常表示几个随时间变化的物理量，例如几种生物种群的数量，或者是几个常见的物理量。
对于混沌现象，我们通常更关心变量的最终归宿。
因此，在研究时，我们有时舍去方程解的前几组数值。
四．MatLab求解建立函数M文件rossler.m，在其中用x(1)表示x，用x(2)表示y，用x(3)表示z.function r=rossler(t,x)global a;global b;global c;r=[-x(2)-x(3);x(1)+a*x(2);b+x(3)*(x(1)-c)];主程序如下：global a;global b;global c;b=2;c=4;t0=[0,200];for a=0:0.02:0.65[t,x]=ode45('rossler',t0,[0,0,0]);asubplot(1,2,1);plot(t,x(:,1),'r',t,x(:,2),'g',t,x(:,3),'b');title('x(红色),y(绿色),z(篮色)随t变化情况');xlabel('t');subplot(1,2,2);plot3(x(:,1),x(:,2),x(:,3))title('相图');xlabel('x');ylabel('y');zlabel('z');pauseend1． b=2,c=4a=0.00当a=0时，(x,y,z)收敛于（0,0.5,0.5）当a=0.05时，(x,y,z)仍然收敛，但收敛速度较小。
北科大Matlab_数学实验报告1~6次(全)
《数学实验》报告实验名称 Matlab 基础知识学院专业班级姓名学号2014年 6月一、【实验目的】1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。
2.掌握Matlab基本操作和常用命令。
3.了解Matlab常用函数，运算符和表达式。
4.掌握Matlab工作方式和M文件的相关知识。
5.学会Matlab中矩阵和数组的运算。
二、【实验任务】P16 第4题编写函数文件，计算1!nkk =∑，并求出当k=20时表达式的值。
P27第2题矩阵A=123456789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，B=468556322⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，计算A*B,A.*B,并比较两者的区别。
P27第3题已知矩阵A=5291⎡⎤⎢⎥⎣⎦，B=1292⎡⎤⎢⎥⎣⎦，做简单的关系运算A>B,A==B,A<B,并做逻辑运算(A==B)&(A<B),(A==B)&(A>B)。
P34 第1题用11114357π=-+-+……公式求π的近似值，直到某一项的绝对值小于-610为止。
三、【实验程序】P16 第4题function sum=jiecheng(n) sum=0;y=1;for k=1:nfor i=1:ky=y*i;endsum=sum+y;endsumP27第2题>>A=[1 2 3;4 5 6;7 8 9]>>B=[4 6 8;5 5 6;3 2 2]>>A*BP27第3题>> A=[5 2;9 1];B=[1 2;9 2];>>A>B>>A==B>>A<B>> (A==B)&(A<B)>> (A==B)&(A>B)P34 第1题t=1;pi=0;n=1;s=1;while abs(t)>=1e-6pi=pi+t;n=n+2;s=-s;t=s/n;endpi=4*pi;四、【实验结果】P16 第4题P27第2题两者的区别：A*B是按正规算法进行矩阵的计算, A.*B是对应元素相乘。
数学实验与数学软件
第一部分 基础实验数学实验与数学软件Mathematica1.1 实验目的了解高等数学实验的含义；初步掌握数学软件Mathematica的用法和基本功能。
1.2 实验内容1. 为什么要开设数学实验？长期以来，从小学到大学十几年，数学一直是我们学习的一门主课，老师所讲的、学生所练、所考的主要是定义叙述、定理证明、公式推算、计算方法、……，数学给我们的印象是，沿“定义→公理→定理→推论→证明”这么一条演绎道路进行的、一个十分严格的数学推理王国和一个充满美感的抽象世界。
然而，我们却不知道，也许也没有想过，这些如此严密、完整、美妙的结论是怎么来的？数学家是通过什么样的方式发现它们的？我们从这些可爱结论本身看不到数学家发现它们的艰辛，也体会不到数学家在发现它们之后的一种喜悦。
传统的数学教育过分强调了形式化的逻辑推导和形式化的结果，淡化了在科学突破上至关重要的实验、观察、归纳、猜想、验证和创新等能力的培养，以至使我们越来越不明白数学从何而来。
充满生机、充满活力、充满美感的数学成了内容多、负担重、枯燥乏味的公式、结论和习题的堆积。
计算机不仅使我们生活和工作的各个方面发生了可喜的变化，而且也使数学教育有了一种新的选择，这种新的选择所带来的模式正在发生一场变革，这场变革的特点就是“数学实验”。
与传统的课堂教学不同是，“数学实验”课的教学把教师的“教授→记忆→测试”的传统教学过程，变成“观察→直觉→探试→思考→归纳→猜想→证明”，将信息的单向交流变成多向交流。
“学习任何东西，最好的途径是自己去发现”，由我们自己动手，用我们熟悉的、喜欢“玩”的计算机去尝试发现一些问题，亲身感受在计算机的帮助下解决问题的酸甜苦辣，在自己的探索和实践中去体验那条被掩盖着的数学思维轨迹。
做，然后知不足！令即可完成。
打一个简单的比方，将我们设计的数学实验和大家熟悉的物理实验作一个比较，结果。
“思考题”是一些稍复杂的问题，是学生完成课堂实验后选做的实验，一般需要学生2.如何做数学实验？大家都做过物理实验、化学实验或生物实验，他们一般都有数目繁多的实验仪器、设备，或简单地说实验工具。
matlab实验一实验报告
matlab实验一实验报告实验一：Matlab实验报告引言：Matlab是一种强大的数学软件工具，广泛应用于科学计算、数据分析和工程设计等领域。
本实验旨在通过使用Matlab解决实际问题，探索其功能和应用。
一、实验目的本次实验的主要目的是熟悉Matlab的基本操作和常用函数，了解其在科学计算中的应用。
二、实验内容1. 数值计算在Matlab中，我们可以进行各种数值计算，包括基本的加减乘除运算，以及更复杂的矩阵运算和方程求解。
通过编写相应的代码，我们可以实现这些功能。
例如，我们可以使用Matlab计算两个矩阵的乘积，并输出结果。
代码如下：```matlabA = [1 2; 3 4];B = [5 6; 7 8];C = A * B;disp(C);```2. 数据可视化Matlab还提供了强大的数据可视化功能，可以将数据以图表的形式展示出来，更直观地观察数据的规律和趋势。
例如，我们可以使用Matlab绘制一个简单的折线图，来展示某个物体在不同时间下的位置变化。
代码如下：```matlabt = 0:0.1:10;x = sin(t);plot(t, x);xlabel('Time');ylabel('Position');title('Position vs. Time');```3. 图像处理Matlab还可以进行图像处理，包括图像的读取、处理和保存等操作。
我们可以通过Matlab对图像进行增强、滤波、分割等处理，以及进行图像的压缩和重建。
例如，我们可以使用Matlab读取一张图片，并对其进行灰度化处理。
代码如下：```matlabimg = imread('image.jpg');gray_img = rgb2gray(img);imshow(gray_img);```三、实验结果与分析在本次实验中，我们成功完成了数值计算、数据可视化和图像处理等任务。
Matlab数学实验报告
实验一 Matlab基本操作1．实验课程名称数学实验2．实验项目名称Matlab基本操作3．实验目的和要求了解Matlab的基本知识，熟悉其上机环境，掌握利用Matlab进行基本运算的方法。
4．实验内容和原理内容：三角形的面积的海伦公式为：area=)s-sa--)()(s(csb其中： s=(a+b+c)/2原理：将一般数学问题转化成对应的计算机模型并进行处理的能力。
了解Matlab的基本功能，会进行简单的操作。
5．主要仪器设备计算机与Windows 2000/XP系统；Matlab等软件。
6．操作方法与实验步骤步骤：（1）在M文件编辑窗口输入以下程序，并以文件名”area_helen.m”保存：a= input(‘a=‘) ； b= input(‘b=‘) ； c= input(‘c=‘) ；s= (a+b+c)/2;area=sqrt (s* (s-a) * (s-b) * (s-c))（2）在命令窗口输入文件名“area_helen”，按回车键，即可运行上面的程序，输入三边长，立即可得三角形面积（3）第二题在命令窗口输入b=6;a=3;c=a*b,d=c-2*b(4) 按回车键，即可运行上面的程序7．实验结果与分析<1> a=3; b=4; c=5;时，aera=6 当a为3，b为4，c为5时，s=6，aera=6<2> c= 18,d=6，a为3，b为6时，c=18，d=6实验二 Matlab的数值计算1．实验课程名称数学实验2．实验项目名称Matlab的数值计算3．实验目的和要求了解一些简单的矩阵、向量、数组和多项式的构造和运算方法实例，懂得编写简单的数值计算的Matlab程序。
熟悉一些Matlab的简单程序，会用Matlab的工具箱，懂得Matlab的安装和简单的使用。
4．实验内容和原理内容：从函数表：)1(),5.0(),2( ,0x 1x 021x 1x f(x) 32-⎪⎩⎪⎨⎧≤≤<>+=f f f x x求设)1(),2( ,1211)(2-⎩⎨⎧≤>+=f f x xx x x f 求设 原理：利用矩阵、向量、数组、和多项式的构造和运算方法，用常用的几种函数进行一般的数值问题求解。
数学实验报告范例1
本科学生综合性实验报告学号姓名成员：邓礼慧王慧莹闵文婷赵晓雪郝烨刘芳欣谢益华学院数学学院专业、班级B班实验课程名称利用Matlab数学软件绘图教师及职称黄炯老师（数学分析）开课学期2010 至2011 学年下学期填报时间2011 年06 月28 日云南师范大学教务处编印一、实验目的：熟悉MATLAB 软件关于常微分方程（组）的符号解命令；掌握用MATLAB 求常微分方程（组）解析解的方法。
二、实验原理：数学软件求常微分原理三、实验设备、材料：设备：计算机材料：matlab 软件 相关教材四、实验流程：Matlab 软件的下载→学习使用matlab 软件→进行实验→绘制基本函数图形→小组讨论→总结→制作报告五、实验过程：三、实验内容1.求下列方程（组）的通解，并作出解的图形：3(1)(1ln ln )(2)6130(3)25sin 2(4)69(1)233,(5)453,442.x x dy y y x y y y dx xy y y e x y y y x e dx d x y z dt dy x y z dt dz x y z dt '''=+- ++=''''''-+= -+=+⎧=-+⎪⎪⎪=-+ (6)⎨⎪⎪=-+⎪⎩,,.x y z dt dy x z dtdz x y dt ⎧=+⎪⎪⎪=+⎨⎪⎪=+⎪⎩解（1）：即（1）解为x/exp(C1*x)解（2）：即（2）结果为：C1*exp(-1/2*x)*sin(1/2*3^(1/2)*x)+C2*exp(-1/2*x)*cos(1/2*3^(1/2)*x) 解（3）：即（3）结果为：-4/25*exp(x)*cos(2*x)-3/25*exp(x)*sin(2*x)-3/2*x^2+C1*x+C2解（4）：即（4）的结果为：1/9*exp(3*x)*x+1/27*exp(3*x)-3/2*x^2+C1*x+C2解（5）：即（5）的结果为：x: [1x1 sym]y: [1x1 sym]z: [1x1 sym]六、实验总结：经过本次的绘图实验，我们小组学会了应用matlab绘图的基本方法，但还有很多复杂的复合函数图形不能画出，我们小组将继续进行学习和研究。
《数学实验》报告册(华南农业大学)
《数学实验》报告册(华南农业大学)
一、实验目的
本实验的目的是通过实践来深入学习高中数学中的函数、极限、导数等概念和运算，掌握一些应用技巧，加深对数学知识的理解和应用，提高数学思维能力和创新意识。
二、实验过程
本次实验共包括函数的变化趋势、导数与函数的性质、曲线的相关概念与方法、等差数列与等比数列等。
在实验中，我们用数学软件Geogebra进行数学模拟，还使用了一些工具在纸上进行计算和绘图。
1.函数的变化趋势
实验中，我们首先通过Geogebra绘制出一个函数图像，然后利用函数的导数等工具和知识，分析函数图像的运动变化趋势，并进行比较和总结。
2.导数与函数的性质
在这一部分中，我们通过对导数的定义和性质，结合具体的函数图像，来分析函数的性态变化趋势，并对函数的极值、最值、单调性、凸凹性等进行分析。
3.曲线的相关概念与方法
在这一部分中，我们通过曲线的方程和图像，来学习曲线的一些基本概念，如切线、法线、弧长、曲率半径等，同时还进行了一些曲线变换的操作，如平移、翻转、放缩等。
4.等差数列与等比数列
这一部分中，我们学习了等差数列和等比数列的基本概念和特点，掌握了求这些数列的和、通项公式等应用技巧。
三、实验收获
此外，我还学会了如何使用数学软件Geogebra，在使用过程中，我体会到了数学计算和思考的乐趣，并发现了数学工具的优越性，感受到了在数学实践中使用现代科技手段的重要性。
总之，本次实验使我更加熟练地掌握了高中数学中的一些基本概念和运算，提高了我的数学能力和综合素质，也拓展了我的视野和思考方式，使我更加自信地面对数学学习和应用。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
4010000000000000
c =
4003333333333333
是format short。
5.(1)请输入一个四阶方阵可逆A，并使用matlab语句求其行列式、逆矩阵；
(2)求A的元素的最大值及该最大值所处的行、列；
(3)将A中第2列3行元素赋值100；
(4)将A中第2列所有元素置为1；
(5)访问A中第5~12个元素并显示，体会矩阵的数据访问规则；
(6)输入一个4阶矩阵B，计算A+B，A*B，A.*B，
(7)使用语句C=[A,B]；D=reshape（C，8,4）得到D；
(8)使用语句E=[A;B];得到E.
通过以上语句的执行结果，掌握matlab中基本元素的赋值和访问
A=rand(4)
A2 =
-956/225 635/2302 1905/488 -1602/2779
1773/881 -1061/1413 -1367/784 1175/924
7946/1343 4055/3482 -3023/519 1403/1279
-1005/259 -2738/4331 2374/451 -1098/539
数学实验与数学软件
上机实验报告
实验一
课程名称：数学实验与数学软件
班级
日期
成绩评定
姓名
实验室
老师签名
学号
实验名称
MATLAB基础
所用软件
MATLAB
实
验
目
的
及
内
容
熟悉MATLAB的界面，熟悉和使用各个窗口的基本功能，了解基于矩阵的MATLAB运算的基本步骤，会利用帮助察看所需信息，熟悉常用的矩阵运算的函数和试用。
a =
2.0000e-01
b =
4
c =
2.4000e+00
>> format long g
>> a,b,c
a =
0.2
b =
4
c =
2.4
>> format bank
>> a,b,c
a =
0.20
b =
4.00
c =
2.40
>> format hex
>> a,b,c
a =
3fc999999999999a
A =
1 2 3 4
5 6 7 8
9 10 11 12
>> B=[-1.3,sqrt(3);(1+2)*4/5,sin(5);exp(2),6]
B =
-1.3000 1.7321
2.4000 -0.9589
7.3891 6.0000
3.求矩阵方程：设A、B满足关系式：AB＝2B+A,求B。
其中A=[3 0 1; 1 1 0; 0 1 4]。要求运用线性代数的矩阵方程的知识，使用matlab代码求解。
4.分别输入一个分数、整数、小数等，（如：a=1/9）,观察显示结果，并使用format函数控制数据的显示格式，如：分别输入format short、format long、format short e、format long g、format bank、format hex等，然后再在命令窗口中输入a，显示a的值的不同形式，并理解这些格式的含义。请问：在matlab中，默认的数据显示格式是什么？
5.(1)请输入一个四阶方阵可逆A，并使用matlab语句求其行列式、逆矩阵；
(2)求A的元素的最大值及该最大值所处的行、列；
(3)将A中第2列3行元素赋值100；
(4)将A中第2列所有元素置为1；
(5)访问A中第5~12个元素并显示，体会矩阵的数据访问规则；
(6)输入一个4阶矩阵B，计算A+B，A*B，A.*B，
A =
0.4218 0.6557 0.6787 0.6555
0.9157 0.0357 0.7577 0.1712
0.7922 0.8491 0.7431 0.7060
0.9595 0.9340 0.3922 0.0318
>> A1=det(A)
A1 =
-0.0961
>> format rat
>> A2=inv(A)
1.0000 1.0000 0.7577 0.1712
0.7922 1.0000 0.7431 0.7060
0.9595 1.0000 0.3922 0.0318
0.2769 0.6948 0.4387 0.1869
0.0462 0.3171 0.3816 0.4898
0.0971 0.9502 0.7655 0.4456
实
验
原
理
步
骤
、
1.在command window中输入如下指令，显示结果：
（1）x1=[1 2 4],x2=[1,2,1],x3=x1’
2）x1=3.4:6.7；x2=3.4:2:6.7；x3=2.6:-0.8:0
（3）x=linspace(0,1,5)
2.直接输入法创建简单矩阵，理解矩阵B的元素的构成
>> a=max(max(A))
a =
0.9595
>> [rows,cols]=find(a)
rows =
1
cols =
1
>> A(3,2)=100
A =
0.4218 0.6557 0.6787 0.6557
0.9157 0.0357 0.7577 0.1712
0.7922 100.0000 0.7431 0.7060
>> a=1/5,b=4,c=2.4
>> format long
>> a
a =
0.200000000000000
>> b,c
b =
4
c =
2.400000000000000
>> format short
>> a,b,c
a =
0.2000
b =
4
c =
2.4000
>> format short e
>> a,b,c
4.分别输入一个分数、整数、小数等，（如：a=1/9）,观察显示结果，并使用format函数控制数据的显示格式，如：分别输入format short、format long、format short e、format long g、format bank、format hex等，然后再在命令窗口中输入a，显示a的值的不同形式，并理解这些格式的含义。请问：在matlab中，默认的数据显示格式是什么？
A=[1 2 3 4; 5 6 7 8; 9 10 11 12]
B=[-1.3,sqrt(3);(1+2)*4/5,sin(5);exp(2),6]
3.求矩阵方程：设A、B满足关系式：AB＝2B+A,求B。
其中A=[3 0 1; 1 1 0; 0 1 4]。要求运用线性代数的矩阵方程的知识，使用matlab代码求解。
0.8235 0.0344 0.7952 0.6463
0.9595 0.9340 0.3922 0.0318
>> A(:,2)=1
A =
0.4218 1.0000 0.6787 0.6557
1.0000 1.0000 0.7577 0.1712
0.7922 1.0000 0.7431 0.7060
0.9595 1.0000 0.3922 0.0318
>> A(5:12)
x =
0 0.2500 0.5000 0.7500 1.0000
2.直接输入法创建简单矩阵，理解矩阵B的元素的构成
A=[1 2 3 4; 5 6 7 8; 9 10 11 12]
B=[-1.3,sqrt(3);(1+2)*4/5,sin(5);exp(2),6]
>> A=[1 2 3 4; 5 6 7 8; 9 10 11 12]
0.9595 0.3922 0.8235 0.7952
1.0000 0.6557 0.6948 0.1869
1.0000 0.1712 0.3171 0.4898
1.0000 0.7060 0.9502 0.4456
1.0000 0.0318 0.0344 0.6463
>> E=[A;B]
E =
0.4218 1.0000 0.6787 0.6557
（3）x=linspace(0,1,5)
>> x1=[1 2 4],x2=[1,2,1],x3=x1'
x1 =
1 2 4
2
4
>> x1=3.4:6.7;x2=3.4:2:6.7; x3=2.6:-0.8:0
x3 =
2.6000 1.8000 1.0000 0.2000
>> x=linspace(0,1,5)
(7)使用语句C=[A,B]；D=reshape（C，8,4）得到D；
(8)使用语句E=[A;B];得到E。
1.在command window中输入如下指令，显示结果：
（1）x1=[1 2 4],x2=[1,2,1],x3=x1’
2）x1=3.4:6.7；x2=3.4:2:6.7；x3=2.6:-0.8:0
>> A=[3 0 1; 1 1 0; 0 1 4];C=A-2;
>> format rat
>> inv(C)
ans =