当前位置：文档库 › 定积分的近似计算(数学实验报告matlab版)

定积分的近似计算(数学实验报告matlab版)

数学实验报告

针对不同

，在区间上取左端点，即取）中划分区间，在区间上取右端点，即取

）中划分区间，在区间上取中点，即取。

，

）．曲线上相应的点为

将曲线的每一段弧用过点，

来代替，这使得每个

，称此式为

MATLAB实验报告50059

实验一MATLAB操作基础实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。例1-1

例1-2 例1-3 例1-4

3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作：（1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。

5、访问mathworks公司的主页，查询有关MATLAB的产品信息。主要教学环节的组织：教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。思考题： 1、如何启动与退出MA TLAB集成环境？启动：（1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。（2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。（3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。退出：（1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。（2）在MA TLAB命令窗口中输入exit或quit命令。（3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？

计算方法上机实验报告

《计算方法》上机实验报告班级：XXXXXX 小组成员：XXXXXXX XXXXXXX XXXXXXX XXXXXXX 任课教师：XXX 二〇一八年五月二十五日

前言通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写。以下为本次上机实验报告，按照实验内容共分为六部分。实验一：一、实验名称及题目： Newton 迭代法例2.7(P38):应用Newton 迭代法求在附近的数值解，并使其满足 . 二、解题思路：设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交点的横坐标) (') (0001x f x f x x - =,称1x 为'x 的一次近似值,过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标) (') (1112x f x f x x - =称2x 为'x

的二次近似值,重复以上过程,得'x 的近似值序列{}n x ,把 ) (') (1n n n n x f x f x x - =+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。三、Matlab 程序代码： function newton_iteration(x0,tol) syms z %定义自变量 format long %定义精度 f=z*z*z-z-1; f1=diff(f);%求导 y=subs(f,z,x0); y1=subs(f1,z,x0);%向函数中代值 x1=x0-y/y1; k=1; while abs(x1-x0)>=tol x0=x1; y=subs(f,z,x0); y1=subs(f1,z,x0); x1=x0-y/y1;k=k+1; end x=double(x1) K 四、运行结果：实验二：

Matlab数学实验报告一

数学软件课程设计题目非线性方程求解班级数学081 姓名曹曼伦

实验目的：用二分法与Newton迭代法求解非线性方程的根；用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。编程实现二分法及Newton迭代法；学会使用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。实验内容：比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方） (1)在区间（0，1)内用二分法； (2)用迭代法x=(2-exp(x))/10，取初值x=0 。试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0，则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正，则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负，则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法； format long

syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果：（1）二分法： symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 （2）迭代法由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次，而用迭代法求得同样精度的解仅用5次，但由于迭代法一般只具有局部收敛性，因此通常不用二分法来求得非线性方程的精确解，而只用它求得根的一个近似解，再用收敛速度较快的迭代法求得其精确解。

数学实验“几种常见的求积分近似解的方法”实验报告(内含matlab程序)

西京学院数学软件实验任务书

实验二十一实验报告一、实验名称：Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。二、实验目的：进一步熟悉Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。三、实验要求：运用Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。四、实验原理： 1．Romberg 积分法：龙贝格积分法是用里查森外推算法来加快复合梯形求积公式的收敛速度，它的算法如下，其中()i m T 是通过一系列逼近原定积分的龙贝格分值. 计算(0)1[()()]2 b a T f a f b -= + 对1,2,3,k n = ，计算下列各步： 21()(1)1 111 1(21)()[()]222k k k k k j b a j b a T T f a ---=---=++∑

对1,2,,m k = 和,1,2,,1i k k k =-- ，计算111 441 m i i i m m m m T T T --+-=- 随着计算的步骤的增加，()i m T 越来越逼近积分()b a f x dx ?。 2．Gauss 型积分法：高斯积分公式的思想是用n 个不等距的节点123,,,n x x x x 对被积函数进行插值，然后对插值后的函数进行积分，其积分公式为： 1 1 1 ()()n k k k f x dx A f x -=≈∑? 如果积分区间不是[1,1]-，则需转换到此区间： 11()()222 b a b a b a b a f x dx f t dt ---+= +? ? 其中系数k A 、节点k x 与n 的关系如下表所示： 3．高斯-切比雪夫积分法：第一类切比雪夫积分形式为： 1 1 ()()n k k k f x dx A f x -=≈∑? 其中k A n π= ，21cos 2k k x n π-= 4．高斯-拉盖尔积分法：高斯-拉盖尔公式有两种形式： 1 ()()n x k k k e f x dx A f x +∞ -=≈∑?

MATLAB实验三-定积分的近似计算

实验三定积分的近似计算一、问题背景与实验目的利用牛顿—莱布尼兹公式虽然可以精确地计算定积分的值，但它仅适用于被积函数的原函数能用初等函数表达出来的情形．如果这点办不到或者不容易办到，这就有必要考虑近似计算的方法．在定积分的很多应用问题中，被积函数甚至没有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只能应用近似方法去计算相应的定积分．本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线法．对于定积分的近似数值计算，Matlab有专门函数可用．二、相关函数（命令）及简介 1．sum(a)：求数组a的和． 2．format long：长格式，即屏幕显示15位有效数字．（注：由于本实验要比较近似解法和精确求解间的误差，需要更高的精度）．3．double()：若输入的是字符则转化为相应的ASCII码；若输入的是整型数值则转化为相应的实型数值． 4．quad()：抛物线法求数值积分．格式： quad(fun,a,b) ，注意此处的fun是函数，并且为数值形式的，所以使用*、/、^等运算时要在其前加上小数点，即 .*、./、.^等．例：Q = quad('1./(x.^3-2*x-5)',0,2); 5．trapz()：梯形法求数值积分．格式：trapz(x,y) 其中x为带有步长的积分区间；y为数值形式的运算（相当于上面介绍的函数fun）例：计算 0sin()d x x π ? x=0:pi/100:pi;y=sin(x); trapz(x,y) 6．dblquad()：抛物线法求二重数值积分．格式：dblquad(fun,xmin,xmax,ymin,ymax)，fun可以用inline定义，也可以通过某个函数文件的句柄传递．例1：Q1 = dblquad(inline('y*sin(x)'), pi, 2*pi, 0, pi) 顺便计算下面的Q2，通过计算，比较Q1 与Q2结果（或加上手工验算），找出积分变量x、y的上下限的函数代入方法． Q2 = dblquad(inline('y*sin(x)'), 0, pi, pi, 2*pi)例2：Q3 = dblquad(@integrnd, pi, 2*pi, 0, pi) 这时必须存在一个函数文件integrnd.m：

matlab数学实验报告5

数学实验报告制作成员班级学号 2011年6月12日

培养容器温度变化率模型一、实验目的利用matlab软件估测培养容器温度变化率二、实验问题现在大棚技术越来越好，能够将温度控制在一定温度范围内。为利用这种优势，实验室现在需要培植某种适于在8.16℃到10.74℃下能够快速长大的甜菜品种。为达到实验所需温度，又尽可能地节约成本，研究所决定使用如下方式控制培养容器的温度：1，每天加热一次或两次，每次约两小时； 2，当温度降至8.16℃时，加热装置开始工作；当温度达到10.74℃时，加热装置停止工作。已知实验的时间是冬天，实验室为了其它实验的需要已经将实验室的温度大致稳定在0℃。下表记录的是该培养容器某一天的温度时间(h)温度（℃）时间（h）温度（℃）09.68 1.849.31 0.929.45 2.959.13 3.878.981 4.989.65 4.988.811 5.909.41 5.908.691 6.839.18 7.008.5217.938.92 7.938.3919.048.66 8.978.2219.968.43 9.89加热装置工作20.848.22 10.93加热装置工作22.02加热装置工作10.9510.8222.96加热装置工作12.0310.5023.8810.59 12.9510.2124.9910.35 13.889.9425.9110.18 三、建立数学模型 1，分析：由物理学中的傅利叶传热定律知温度变化率只取决于温度

差，与温度本身无关。因为培养容器最低温度和最高温度分别是：8.16℃和10.74℃；即最低温度差和最高温度差分别是：8.16℃和10.74℃。而且，16.8/74.10≈1.1467，约为1，故可以忽略温度对温度变化率的影响2，将温度变化率看成是时间的连续函数，为计算简单，不妨将温度变化率定义成单位时间温度变化的多少，即温度对时间连续变化的绝对值（温度是下降的），得到结果后再乘以一系数即可。四、问题求解和程序设计流程1）温度变化率的估计方法根据上表的数据，利用matlab 做出温度-时间散点图如下：下面计算温度变化率与时间的关系。由图选择将数据分三段，然后对每一段数据做如下处理：设某段数据为{(0x ,0y ),(1x ,1y ),(2x , 2y ),…,(n x ,n y )},相邻数据中点的平均温度变化率采取公式：温度变化率=（左端点的温度-右端点的温度）/区间长度算得即：v( 2 1i i x x ++)=(1+-i i y y )/(i i x x - +1). 每段首尾点的温度变化率采用下面的公式计算：v(0x )=(30y -41y +2y )/(2x -0x )v(n x )=(3n y -41+n y +2+n y )/(n x -2-n x )

MATLAB实验报告(1-4)

信号与系统MATLAB第一次实验报告一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。举例：计算一个函数并绘制出在对应区间上对应的值。

2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1：举例2：

数学实验报告

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

2019年4月15日

一、问题背景与提出图1是从百度地图中截取的山东省地图，试根据前面数值积分计算方法，计算山东省面积。图 1 二、实验目的 1、学会运用matlab 解决一些简单的数学应用问题。 2、学会运用matlab 建立数学模型。 3、学会运用一些常见的数值积分计算方法结算实际问题，并了解其实际意义，建立积分模型。三、实验原理与数学模型将积分区间 [a , b] n 等分，每个区间宽度均为h = (b - a) / n , h 称为积分步长。记 a = x 0 < x 1 < … < x k … < x n = b , 在小区间上用小矩形面积近似小曲边梯形的面积，若分别取左端点和右端点的函数值为小矩形的高，则分别得到两个曲边梯形的面积的近似公式： Ln = h ∑f (x k )n=1k=0 , h = b?a ?

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

MATLAB实验报告

数字信号处理及MATLAB 实验报告班级：学号：姓名：

4.7.2 例4,2 设x(n)是由两个正弦信号及白噪声的叠加，试用FFT文件对其作频谱分析。程序清单 %产生两个正弦加白噪声 N=256; f1=.1;f2=.2;fs=1; a1=5;a2=3; w=2*pi/fs; x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N); %应用FFT求频谱 subplot(2,2,1); plot(x(1:N/4)); title('原始信号'); f=-0.5:1/N:0.5-1/N; x=fft(x); y=ifft(x); subplot(2,2,2); plot(f,fftshift(abs(x))); title('频域信号'); subplot(2,2,3); plot(real(x(1:N/4))); title('时域信号');

例4.3 设x（n）为长度N=6的矩形序列，用MATLAB程序分析FFT取不同长度时x（n）频谱的变化。N=8,32,64,时x（n）的FFT MATLAB实现程序如下。 x=[1,1,1,1,1,1]; N=8; y1=fft(x,N); n=0:N-1; subplot(3,1,1);stem(n,abs(y1),'.k');axis([0,9,0,6]); N=32; y2=fft(x,N); n=0:N-1; subplot(3,1,2);stem(n,abs(y2),'.k');axis([0,40,0,6]); N=64; y3=fft(x,N); subplot(3,1,3);stem(n,abs(y3),'.k');axis([0,80,0,6]);

浅析Matlab数学实验报告

数学实验报告姓名: 班级: 学号: 第一次实验任务过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)

cos()sin(2=-=++y x y x y x ,其中、计算的图形。下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);

数学计算方法实验报告

数学计算方法实验报告习题二 2.估计用二分法求方程f(x)=x3+4x2-10=0在区间[1,2]内根的近似值，为使方程不超过10时所需的二分次数。f(x k) 程序过程： function two (tolerance) a=1;b=2;counter=0; while (abs(b-a)>tolerance) c=(a+b)/2; fa=a^3+4*a^2-10;

fb=b^3+4*b^2-10; fc=c^3+4*c^2-10; if ((fa==0|fb==0)) disp(counter); elseif (fa*fc<0) b=c;counter=counter+1; elseif (fb*fc<0) a=c;counter=counter+1; elseif (fb==0) disp(counter); end end solution=(a+b)/2; disp(solution); disp(counter); 实验结果： 6.取x0=1.5,用牛顿迭代法求第三中的方程根.f(x)=x3+4x2-10=0的近似值（精确到||x k+1-x k|≦10-5,并将迭代次数与3题比较。程序过程： function six (g) a=1.5; fa=a^3+4*a^2-10;

ga=3*a^2+8*a; b=a-fa/ga; k=1; while(abs(b-a)>g) a=b; fa=a^3+4*a^2-10; ga=3*a^2+8*a; b=a-fa/ga; k=k+1; end format long; disp(a); disp(k); 实验结果：程序结果计算结果 8.用弦割法求方程f(x)=x3-3x2-x+9=0在区间[-2,-1]内的一个实根近似值x k,|f(x k)|≦10-5. 程序过程： function eight (t) a=-2; b=-1; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=1; while(abs(c-b)>t) a=b; b=c; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=k+1; end

《数学软件》实验报告-符号计算基础与符号微积分

实验报告课程名称：数学软件姓名：学院：专业：年级：学号：指导教师：职称：年月日

实验项目列表

附件三：实验报告（二）系：专业：年级：姓名学号：实验课程：实验室号：_ 实验设备号：实验时间：指导教师签字：成绩： 1. 实验项目名称：符号计算基础与符号微积分 2. 实验目的和要求 1.掌握定义符号对象的方法 2.掌握符号表达式的运算法则以及符号矩阵运算 3.掌握求符号函数极限及其导数的方法 4.掌握求符号函数定积分和不定积分的方法 3. 实验使用的主要仪器设备和软件方正商祺N260微机；MATLAB7. 0或以上版本 4. 实验的基本理论和方法（1）符号函数;sym(x)；syms a b …… （2）平方根：sqrt(x) （3）分解因式：factor（s）（4）符号表达式化简：simplify（s）（5）逆矩阵：inv(x) （6）下三角矩阵：tril（x）（7）矩阵行列式的值:det(x)

（8）符号函数求极限：limit （f ，x ，a ）；limit （f ，x ，a ，‘right ’）（9）符号函数求导：diff （f ，v ，n ）（10）符号函数求不定积分：int （f ，v ）（11）符号函数求定积分：int （f ，v ，a ，b ） 5. 实验内容与步骤（描述实验中应该做什么事情，如何做等，实验过程中记录发生的现象、中间结果、最终得到的结果，并进行分析说明）（包括：题目，写过程、答案）题目： 1. 已知x=6，y=5，利用符号表达式求 y x x z -++= 31。提示：定义符号常数)'5(')'6('sym y sym x ==，。 >> x=sym('6'); >> y=sym('5'); >> z=(x+1)/(sqrt(3+x)-sqrt(y)) z = 7/(3-5^(1/2)) 2. 分解因式：44y x - >> syms x y; >> A=x^4-y^4; >> factor(A) ans = (x-y)*(x+y)*(x^2+y^2) 3. 化简表达式 (1)2121sin cos cos sin ββββ- (2) 123842+++x x x (1) >> syms x y; >> f1=sin(x)*cos(y)-cos(x)*sin(y);

matlab实验报告

Matlab实验报告 ——定积分的近似计算学生姓名：学号：专业：数学与应用数学专业

数学实验报告实验序号:1001114030 日期：2012年10月20日班级应一姓名陈璐学号1001114030 实验名称：定积分的近似运算问题背景描述：利用牛顿—莱布尼茨公式虽然可以精确地计算定积分的值，但它仅适合于被积分函数的原函数能用初等函数表达出来的情形。如果这点办不到或不容易办到，这就有必要考虑近似计算的方法。在定积分的很多应用问题中，被积函数甚至没有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只能应用近似方法去计算相应的定积分。实验目的：本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线发。对于定积分的近似数值计算，Matlab有专门函数可用。实验原理与数学模型： 1.sum（a）：求数组a的和。 2.format long：长格式，即屏幕显示15位有效数字。 3.double（）：若输入的是字符则转化为相应的ASCII码；若输入的是整型数之则转化为相应的实型数值。 4.quad（）：抛物线法求数值积分。格式：quad（fun，a，b）。此处的fun是函数，并且

为数值形式，所以使用*、/、^等运算时要在其前加上小数点。 5.trapz():梯形法求数值积分。格式：trapz（x，y）。其中x为带有步长的积分区间;y为数值形式的运算。 6.fprintf（文件地址，格式，写入的变量）：把数据写入指定文件。 7.syms 变量1变量2……：定义变量为符号。 8.sym（'表达式'）：将表达式定义为符号。 9.int（f，v，a，b）：求f关于v积分，积分区间由a到b。 10.subs（f，'x'，a）：将a的值赋给符号表达式f中的x，并计算出值。若简单地使用subs （f），则将f的所有符号变量用可能的数值代入，并计算出值。实验所用软件及版本：Matlab 7.0.1

MATLAB数学实验报告

Matlab 数学实验报告

一、实验目的通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 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问题分析如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积，2a÷π×πx2=2aπ2，再求AB的面积，用扇形ABE的面积减去三角形ABE的面积即可。

离散数学实验报告()

《离散数学》实验报告专业网络工程班级姓名学号授课教师二 O 一六年十二月

目录实验一联结词的运算实验二根据矩阵的乘法求复合关系实验三利用warshall算法求关系的传递闭包实验四图的可达矩阵实现

实验一联结词的运算一．实验目的通过上机实验操作，将命题连接词运算融入到C语言的程序编写中，一方面加强对命题连接词运算的理解，另一方面通过编程实现命题连接词运算，帮助学生复习和锻炼C语言知识，将理论知识与实际操作结合，让学生更加容易理解和记忆命题连接词运算。二．实验原理 (1) 非运算, 符号: ,当P=T时，P为F, 当P=F时，P为T 。 (2) 合取, 符号: ∧ , 当且仅当P和Q的真值同为真，命题P∧Q的真值才为真；否则，P∧Q的真值为假。 (3) 析取, 符号: ∨ , 当且仅当P和Q的真值同为假，命题P∨Q的真值才为假；否则，P∨Q的真值为真。 (4) 异或, 符号: ▽ , 当且仅当P和Q的真值不同时，命题P▽Q的真值才为真；否则，P▽Q的真值为真。 (5) 蕴涵, 符号: →, 当且仅当P为T,Q为F时，命题P→Q的真值才为假；否则，P→Q 的真值为真。 (6) 等价, 符号: ?, 当且仅当P,Q的真值不同时，命题P?Q的真值才为假；否则，P→Q的真值为真。三．实验内容编写一个程序实现非运算、合取运算、析取运算、异或运算、蕴涵运算、等价运算。四．算法程序 #include void main() { printf("请输入P、Q的真值\n"); int a,b; scanf("%d%d",&a,&b); int c,d; if(a==1) c=0; else c=1; if(b==1) d=0;

定积分的近似计算

数学实验报告实验序号：4 日期：2012 年12 月13 日实验名称定积分的近似计算问题背景描述：利用牛顿—莱布尼兹公式虽然可以精确地计算定积分的值，但它仅适用于被积函数的原函数能用初等函数表达出来的情形．如果这点办不到或者不容易办到，这就有必要考虑近似计算的方法．在定积分的很多应用问题中，被积函数甚至没有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只能应用近似方法去计算相应的定积分．实验目的：本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线法。对于定积分的近似数值计算，Matlab有专门函数可用。

实验原理与数学模型： 1．矩形法根据定积分的定义，每一个积分和都可以看作是定积分的一个近似值，即在几何意义上，这是用一系列小矩形面积近似小曲边梯形的结果，所以把这个近似计算方法称为矩形法．不过，只有当积分区间被分割得很细时，矩形法才有一定的精确度．针对不同的取法，计算结果会有不同。（1）左点法：对等分区间，在区间上取左端点，即取。（2）右点法：同（1）中划分区间，在区间上取右端点，即取。（3）中点法：同（1）中划分区间，在区间上取中点，即取。2．梯形法等分区间，相应函数值为（）．

曲线上相应的点为（）将曲线的每一段弧用过点，的弦（线性函数）来代替，这使得每个上的曲边梯形成为真正的梯形，其面积为，．于是各个小梯形面积之和就是曲边梯形面积的近似值，，即，称此式为梯形公式。 3．抛物线法将积分区间作等分，分点依次为，，对应函数值为（），曲线上相应点为（）．现把区间上的曲线段用通过三点，，的抛物线

《数学实验》报告matlab-第五次作业

《数学实验》报告实验名称 matlab拟合与插值学院机械工程学院专业班级姓名学号

2011年 10月

一、【实验目的】掌握Matlab关于采用最小二乘法拟合曲线的方法。学会使用matlab求实际中得到数据的插值曲线。二、【实验任务】 P130第8、10、12题三、【实验程序】 P130第8题： x=[0.10,0.30,0.40,0.55,0.70,0.80,0.95]; y=[15,18,19,21,22.6,23.8,26]; p1=polyfit(x,y,1); p3=polyfit(x,y,3); p5=polyfit(x,y,5); disp('一阶拟合函数'),f1=poly2str(p1,'x') disp('三阶拟合函数'),f3=poly2str(p3,'x') disp('五阶拟合函数'),f5=poly2str(p5,'x') x1=0.1:0.0017:0.95; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1); plot(x,y,'rp',x1,y1,'--',x1,y3,'k-.',x1,y5); legend('拟合点','一次拟合','三次拟合','七次拟合') P130第10题 x=[10,15,20,25,30]; y=[25.2,29.8,31.2,31.7,29.4]; xi=10:.5:30; yi1=interp1(x,y,xi,'*nearest'); yi2=interp1(x,y,xi,'*linear'); yi3=interp1(x,y,xi,'*spline'); yi4=interp1(x,y,xi,'*cubic'); plot(x,y,'ro',xi,yi1,'--',xi,yi2,'-',xi,yi3,'k.-',xi,yi4,'m:') ,grid on

高等数学实验报告matlab

西安交通大学高等数学实验报告班级组员与学号 2013年

实验名称：学生成绩管理一、实验目的二、实验内容三、详细编程 clear for i=1:10 a{i}=89+i; b{i}=79+i; c{i}=69+i; d{i}=59+i; end c=[d,c]; Name=input('please input name:'); Score=input('please input score:'); n=length(Score); Rank=cell(1,n); S=struct('Name',Name,'Score',Score,'Rank',Rank); for i=1:n switch S(i).Score case 100 S(i).Rank='满分'; case a S(i).Rank='优秀'; case b S(i).Rank='良好'; case c S(i).Rank='及格'; otherwise S(i).Rank='不及格'; end end disp(['学生姓名 ','得分 ','等级']);

for i=1:n disp([S(i).Name,blanks(6),num2str(S(i).Score),blanks(6),S(i).Rank]); end s=0; for i=1:n s=S(i).Score+s; end averscore=s/n; t=S(1).Score; for i=1:(n-1) if(S(i).ScoreS(i+1).Score) m=S(i+1).Score; end end disp(['平均成绩']); disp([averscore]); disp(['最高分']); disp(t); disp(['最低分']); disp(m); 四、实验结果

定积分的近似计算(数学实验报告matlab版)

MATLAB实验报告50059

计算方法上机实验报告

Matlab数学实验报告一

数学实验“几种常见的求积分近似解的方法”实验报告(内含matlab程序)

MATLAB实验三-定积分的近似计算

matlab数学实验报告5

MATLAB实验报告(1-4)

数学实验报告

MATLAB实验报告

浅析Matlab数学实验报告

数学计算方法实验报告

《数学软件》实验报告-符号计算基础与符号微积分

matlab实验报告

MATLAB数学实验报告

离散数学实验报告()

定积分的近似计算

《数学实验》报告matlab-第五次作业

高等数学实验报告matlab

最新定积分的近似计算2