matlab实验题目.doc
MATLAB概率习题
数学实验（概率论）题目一．用MATLAB 计算随机变量的分布1．用MA TLAB 计算二项分布在一级品率为0.2的大批产品中，随机地抽取20个产品，求其中有2个一级品的概率。
1． 用MA TLAB 计算泊松分布用MATLAB 计算：保险公司售出某种寿险保单2500份.已知此项寿险每单需交保费120元,当被保人一年内死亡时,其家属可以从保险公司获得2万元的赔偿(即保额为2万元).若此类被保人一年内死亡的概率0.002,试求:(1)保险公司的此项寿险亏损的概率;(2)保险公司从此项寿险获利不少于10万元的概率; (3)获利不少于20万元的概率. 3．用MA TLAB 计算均匀分布乘客到车站候车时间ξ()0,6U ，计算()13P ξ<≤。
4．用MA TLAB 计算指数分布用MA TLAB 计算：某元件寿命ξ服从参数为λ（λ=11000-）的指数分布.3个这样的元件使用1000小时后，都没有损坏的概率是多少? 5。
用MATLAB 计算正态分布 某厂生产一种设备，其平均寿命为10年，标准差为2年.如该设备的寿命服从正态分布，求寿命不低于9年的设备占整批设备的比例？ 二．用MATLAB 计算随机变量的期望和方差 1．用MA TLAB 计算数学期望（1）用MATLAB 计算离散型随机变量的期望 1）。
一批产品中有一、二、三等品、等外品及废品5种，相应的概率分别为0.7、0.1、0.1、0.06及0.04,若其产值分别为6元、5.4元、5元、4元及0元.求产值的平均值 2）。
已知随机变量X 的分布列如下：{}kk X p 21== ,,2,1n k =计算.EX （2）用MATLAB 计算连续型随机变量的数学期望假定国际市场上对我国某种商品的年需求量是一个随机变量ξ（单位：吨），服从区间[],a b 上的均匀分布，其概率密度为： 1()0a x bx b aϕ⎧≤≤⎪=-⎨⎪⎩其它计算我国该种商品在国际市场上年销售量的期望.ξE .（3）用MATLAB 计算随机变量函数的数学期望假定国际市场每年对我国某种商品的需求量是随机变量X （单位：吨），服从[20,40]上的均匀分布，已知该商品每售出1吨，可获利3万美元，若销售不出去，则每吨要损失1万美元，如何组织货源，才可使收益最大？ 2． 用MA TLAB 计算方差（1）利用MATLAB 计算：设有甲、乙两种股票，今年的价格都是10元，一年后它们的试比较购买这两种股票时的投资风险.。
MATLAB实训实验
2015/2016学年下学期《信号与系统》实验报告班级:学号:学生姓名:指导教师:2016年3月8 日实验一 基本函数仿真实验项目: 基本函数仿真实验时间: 2016年 3 月 8 日 星期 二 第 34 节课 实验地点: 1501实验室 实验目的:1、 学习使用MATLAB 软件2、 学习MATLAB 中各种函数，并应用函数分析3、 对MATALB 的进一步的学习了解，熟练掌握MATALB 的各种操纵，学会使用MATALB 解决复杂的运算并学会用MATALB 解决平时学习4、 了解MATALB 的数值运算5、 了解MATALB 的基本函数和命令6、 学习掌握MATALB 有关命令 实验内容: 1、(1) 题目：应用MA TLAB 方法实现单位阶跃信号和矩形脉冲。
(2) 程序清单（源程序）解：对于阶跃函数，MATLAB 中有专门的stairs 绘图命令。
例如，实现)(t 和矩形脉冲的程序如下：t=-1:2; % 定义时间范围向量t x=(t>=0);subplot(1,2,1),stairs(t,x);axis([-1,2,-0.1,1.2]); grid on % 绘制单位阶跃信号波形 t=-1:0.001:1; % 定义时间范围向量t g=(t>=(-1/2))-(t>=(1/2));subplot(1,2,2),stairs(t,g);axis([-1,1,-0.1,1.2]); grid on % 绘制矩形脉冲波形(3) 运行结果（截图）00.20.40.60.8100.20.40.60.81图1 例1图（4）函数解析Subplot:使用方法：subplot （m,n,p ）或者subplot （m n p ）。
是将多个图画到一个平面上的工具。
其中，m 表示是图排成m 行，n 表示图排成n 列，也就是整个figure 中有n 个图是排成一行的，一共m 行，如果m=2就是表示2行图。
实验一 Matlab基本操作及运算
实验一Matlab基本操作及运算一、实验目的：1．熟悉MATLAB基本操作2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识：1.熟悉MATLAB环境熟悉MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器文件和搜索路径浏览器。
2.掌握MATLAB常用命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英文字母开头（3）长度不大于31个（4）区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。
MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算4. MATLAB的一维、二维数组的访问表6 子数组访问与赋值常用的相关指令格式5. MATLAB的基本运算表7 两种运算指令形式和实质内涵的异同表6.MATLAB的常用函数表8 标准数组生成函数表9 数组操作函数7.多项式运算poly——产生特征多项式系数向量roots——求多项式的根p=poly2str(c,‘x’)—（将特征多项式系数向量c转换为以习惯方式显示是多项式）conv, convs——多项式乘运算deconv——多项式除运算polyder(p)——求p的微分polyder(a, b)——求多项式a,b乘积的微分[p,q]=polyder(p1,p2)——求解多项式p1/p2微分的有理分式poly(p,A)——按数组运算规则求多项式p在自变量A的值polym(p,A)——按矩阵运算规则求多项式p在自变量A的值三、实验内容1、新建一个文件夹（自己的名字命名）2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。
方法如下：3、保存，关闭对话框（要求抓取自己实验的图，插入到自己的实验报告中）4、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）5、使用clc 、clear ，观察command window 、command history 和workspace等窗口的变化结果。
matlab简单编程21个题目及答案
1、设⎥⎦⎤⎢⎣⎡++=)1(sin35.0cos2xxxy，把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线。
第一题的matlab源程序：①考虑cos（x）为一个整体，然后乘以中括号里面的全部x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x).*(0.5+3*sin(x)./(1+x.^2)); %y的表达式plot(x,y)%画出图形图如下：②考虑对整体求解cos，先求x乘以括号中的部分x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x.*(0.5+3*sin(x)./(1+x.^2))); %y的表达式plot(x,y) %画出图形图如下：2、产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。
并求该矩阵全体数的平均值和均方差。
第二题的matlab源程序如下：R1=randn(8,6) %产生正态分布随机矩阵R1 =1.0933 -0.7697 1.5442 -0.1924 1.4193 0.21571.1093 0.3714 0.0859 0.8886 0.2916 -1.1658-0.8637 -0.2256 -1.4916 -0.7648 0.1978 -1.14800.0774 1.1174 -0.7423 -1.4023 1.5877 0.1049-1.2141 -1.0891 -1.0616 -1.4224 -0.8045 0.7223-1.1135 0.0326 2.3505 0.4882 0.6966 2.5855-0.0068 0.5525 -0.6156 -0.1774 0.8351 -0.66691.5326 1.1006 0.7481 -0.1961 -0.2437 0.1873aver=(sum(R1(1:end,1:end)))./8 %产生各行的平均值aver =0.0768 0.1363 0.1022 -0.3473 0.4975 0.1044a=std(R1(1:end,1:end)) %产生各行的均方差也就是标准差a =1.0819 0.8093 1.3456 0.8233 0.8079 1.2150aver1=(sum(R1(:)))./48 %全体数的平均值aver1 =0.0950b=std(R1(:)) %全体数的均方差即标准差b =1.01033、设x=rcost+3t,y=rsint+3,分别令r=2,3,4，画出参数t=0~10区间生成的x~y 曲线。
(完整word版)噪声干扰信号的Matlab仿真
雷达对抗实验报告实验题目:噪声干扰信号的Matlab仿真院系：电子与信息工程学院班级：姓名：学号:指导教师：实验时间： 2012 年 6 月噪声调幅、调频、调相信号的Matlab仿真一、实验目的通过实验，加深对噪声调幅、调频、调相信号的理解，加深对噪声调幅、调频、调相信号频谱分析的基本思想与实现方法的认识，并掌握Matlab对随机过程的仿真方法与其基本函数和语法的使用。
二、实验原理实验中要仿真的各种噪声的时域表达式及相应的频谱特性:1.射频噪声干扰窄带高斯过程:称为射频噪声干扰。
其中包络函数服从瑞利分布，相位函数服从［0，2］均匀分布，且与相互独立，载频为常数,且远大于的谱宽。
2.噪声调幅干扰广义平稳随机过程：称为噪声调幅干扰。
其中，调制噪声为零均值，方差为,在区间［—,分布的广义平稳随机过程，服从［0,2]均匀分布，且为与独立的随机变量,为常数。
噪声调幅信号的波形图,以及联合概率密度分布函数p（)以及各自的概率密度分布密度p(）存在下列关系：3.噪声调频干扰广义平稳随机过程：称为噪声调频干扰，其中调制噪声为零均值、广义平稳的随机过程，服从［0，2］均匀分布且与独立的随机变量,，噪声调频干扰中的调制噪声和噪声调频干扰信号的波形J（t）如下图示:4.噪声调相干扰广义平稳随机过程：称为噪声调频干扰,其中调制噪声为零均值、广义平稳的随机过程，服从［0,2］均匀分布且与独立的随机变量，，噪声调相干扰的功率谱如下图所示：三、实验内容利用Matlab仿真产生视频噪声：;射频噪声:；噪声调幅干扰：视频噪声，调制度m=0.1~1;噪声调频干扰：视频噪声；噪声调相干扰：视频噪声.等一系列干扰信号并分析特性。
四、实验思路与步骤1.产生一个高斯白噪声，2.利用Matlab自带的fir1函数产生一个低通滤波器，限制高斯白噪声的带宽，由此产生了视频噪声.3.利用产生的视频噪声，分别代入噪声调幅干扰的时域表达式，并且进行100次的积累后求平均值,由此画出噪声调幅干扰频域波形，对其进行快速傅里叶变换后，求出功率谱，由此画出噪声调幅干扰的功率谱波形。
MATLAB实验题答案
1、求以下变量的值，并在MATLAB 中验证。
( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 82、下列运算是否合法，为什么？如合法，结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 25 8 11 >> result4=b*d result4 =31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13-5 -8 7>> result6=a.*b result6 =2 8 -34 15 30 >> result7=a./b result7 =0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 =2 4 816 32 64 3、用MATLAB 求解下面的的方程组。
(1)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----017413231511222315921274321x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2)⎪⎪⎩⎪⎪⎨⎧=-++=--=-++=++56533332821w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6]>> B2=[1;8;3;5] >> x2=inv(A1)*B24、已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=1323151122231592127A (1)求矩阵A 的秩(rank) (2)求矩阵A 的行列式(determinant) (3)求矩阵A 的逆(inverse) (4)求矩阵A 的特征值及特征向量(eigenvalue and eigenvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3) >> a=det(A3) >> [V,D]=eig(A3) 5、109910101022222++++==---=∑n n y ,求y=？（运行format long g 命令后，查看y 的值） m1=0; for m=-10:10 m1=m1+2^m; end m1 m1 = 2047.9990234375 6、求分段函数的值。
MATLAB实验
MATLAB实验一：MATLAB语言基本概念实验实验目的：1. 熟悉MATLAB语言及使用环境；2.掌握MATLAB的常用命令；3．掌握MATLAB的工作空间的使用；4.掌握MATLAB的获得帮助的途径。
5.掌握科学计算的有关方法，熟悉MATLAB语言及其在科学计算中的运用；6.掌握MATLAB的命令运行方式和M文件运行方式；7.掌握矩阵在MATLAB中的运用。
实验方案分析及设计：本次实验主要目的是了解MATLAB的使用环境，以及常用的一些命令的使用；了解矩阵在MATLAB实验中的具体运用，以及相关的一些符号命令的使用。
实验器材：电脑一台，MATLAB软件实验步骤：打开MATLAB程序，将实验内容中的题目依次输入MATLAB中，运行得到并记录结果，最后再对所得结果进行验证。
实验内容及要求：1．熟悉MATLAB的菜单和快捷键的功能2．熟悉MATLAB的命令窗口的使用3．熟悉常用指令的使用format clc clear help lookfor who whos 4．熟悉命令历史窗口的使用5. 熟悉MATLAB工作空间的功能将工作空间中的变量保存为M文件，并提取该文件中的变量6.熟悉MATLAB获取帮助的途径将所有plot开头的函数列出来，并详细给出plotfis函数的使用方法1. 输入 A=[7 1 5;2 5 6;3 1 5]，B=[1 1 1; 2 2 2;3 3 3]，在命令窗口中执行下列表达式，掌握其含义：A(2, 3) A(:,2) A(3,:) A(:,1:2:3)A(:,3).*B(:,2) A(:,3)*B(2,:) A*BA.*BA^2 A.^2 B/A B./AA=[7 1 5;2 5 6;3 1 5]7 1 52 5 63 1 5>> B=[1 1 1; 2 2 2;3 3 3]1 1 12 2 23 3 3>> A(2, 3)6>> A(:,2)151>> A(3,:)3 1 5>> A(:,1:2:3)7 52 63 5>> A(:,3).*B(:,2)51215>> A(:,3)*B(2,:)10 10 1012 12 1210 10 10>> A*B24 24 2430 30 3020 20 20>> A.*B7 1 54 10 129 3 15>> A^266 17 6642 33 7038 13 46>> A.^249 1 254 25 369 1 25>> B/A0.1842 0.2105 -0.23680.3684 0.4211 -0.47370.5526 0.6316 -0.7105>> B./A0.1429 1.0000 0.20001.0000 0.4000 0.33331.0000 3.0000 0.60002．输入 C=1:2:20，则 C （i ）表示什么？其中 i=1,2,3, (10)1到19差为2，i 代表公差3. 试用 help 命令理解下面程序各指令的含义：cleart =0:0.001:2*pi;subplot(2,2,1);polar(t, 1+cos(t))subplot(2,2,2);plot(cos(t).^3,sin(t).^3)subplot(2,2,3);polar(t,abs(sin(t).*cos(t)))subplot(2,2,4);polar(t,(cos(2*t)).^0.5)4计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
(完整word版)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 onfor r=0:0.3:3.9x=[0.1];for i=2:150x(i)=r*sin(3.14*x(i-1));endpause(0.5)for i=101:150plot(r,x(i),'k.');endtext(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end加密迭代后clear;clf;axis([0,4,0,4]);hold onfor r=0:0.005:3.9x=[0.1];for i=2:150x(i)=r*sin(3.14*x(i-1));endpause(0.1)for i=101:150plot(r,x(i),'k.');endend运行后得到Feigenbaum图2.2实验题目二2.2.1实验问题某农夫有一个半径10米的圆形牛栏，长满了草。
他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？2.2.2问题分析如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。
问题要求区域ABCD等于圆ABC的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
前言一、上机实验目的上机实验的目的是提高学生对算法的理解程度，并掌握用实用工具进行数值计算的方法，通过实践环节理解数值分析的应用和研究方法。
二、实验基本内容本课程实验内容分为6个实验。
学生可以在课内机时先完成指导书中给出的程序或验证算法，理解所学的知识，在此基础上再编写其他应用程序。
指导书中的15个实验如下。
1．拉格朗日插值法。
2．最小二乘拟合。
3．数值积分。
4．范数计算和LU分解。
5．牛顿迭代法。
6．龙格库塔法。
三、实验任务与时间安排本课程是一门实践性很强的课程，除了在课内安排的实验外，鼓励同学在课外用相关技术进行编程练习。
具体实践课时安排如下：实验指导书实验1 拉格朗日插值法一、实验目的[1] 了解lagrange插值法的基本原理和方法；[2] 通过实例掌握用MA TLAB求插值的方法；[3] 编程实现lagrange插值二、实验内容1．在matlab中command窗口或新建M文件运行以下程序>> x0=[0.4 0.5 0.6 0.7 0.8];%输入节点的X值>> y0=[-0.916291 -0.693147 -0.510826 -0.356675 -0.223144];%输入节点的Y值>> T=interp1(x0,y0,0.54,'linear');%对插值节点使用线性插值并计算0.54处的函数值，改变红色处的参数可得到不同的插值结果，可选项为’cubic’,’ nearest’,’ spline’等。
>> x=0.2:0.05:1.0;%从0.2开始到1.0每隔0.05取一个待求点>> y=interp1(x0,y0,x,'linear');%用插值法求出待求点对应的函数值y>> plot(x0,y0,’r*’);%用红色星号标出插值点对>> hold on;plot(x,y);%蓝色线为所得的插值函数图形2．用C语言或matlab编写实现lagrange插值，要求：数据输入项（函数参数）为：插值节点及函数值，及待求点x的值输出为待求点x对应的函数值（程序流程图如下所示）三、实验步骤1．开启软件平台——MA TLAB，编程可以选用自己熟悉的软件；2．根据各种数值解法步骤编写M文件3．观察运行结果(数值或图形)；四、实验要求与任务每个同学独立完成编程，并在课后提交程序；同时验证matlab中相应函数的用法。
实验2 最小二乘拟合一、实验目的[1] 了解最小二乘拟合的基本原理和方法；[2] 掌握用MA TLAB作最小二乘多项式拟合和曲线拟合的方法；[3] 通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区别。
[4] 了解各种参数辨识的原理和方法；[5] 通过范例展现由机理分析确定模型结构，拟合方法辨识参数，误差分析等求解实际问题的过程；二、实验内容1．用MA TLAB 中的函数作一元函数的多项式拟合与曲线拟合； 2．用MA TLAB 中的函数作二元函数的最小二乘拟合，作出误差图； 3．针对预测和确定参数的实际问题，建立数学模型，并求解。
三、实验步骤1．开启软件平台——MA TLAB ，开启MA TLAB 编辑窗口； 2．根据各种数值解法步骤编写M 文件 3．观察运行结果(数值或图形)；四、实验要求与任务根据实验内容和步骤，完成以下具体实验旧车价格预测某年美国旧车价格的调查资料如下表，其中x i 表示轿车的使用年数，y i 表示相应的平均价格。
试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？多项式函数拟合：[a,s]=polyfit(xdata,ydata,n)其中n 表示多项式的最高阶数，xdata ，ydata 为将要拟合的数据，它是用数组的方式输入．输出参数a 为拟合多项式11++++=n n na x a x a y Λ的系数],,,[11+=n n a a a a Λ,s 为误差。
另：多项式在x 处的值y 可用下面程序计算．y=polyval(a,x)，对上面给出的数据做多项式拟合，可取不同的n 观测此时的误差，看取什么样的n 较好。
并在matlab 中画出拟合曲线(取拟合曲线上若干个点对，用plot 命令实现)实验3 数值积分一、实验目的[1] 了解数值积分的基本原理和方法；[2] 掌握用MATLAB求积分的方法；[3] 通过实例学习如何用几种方法求积分。
[4] 了解书上介绍的几种数值积分的不同原理和方法；二、实验内容1．用MA TLAB中的QUAD函数求出作业所示函数的数值积分；2．用MA TLAB中的QUADL函数求出作业中所示函数的数值积分。
三、实验步骤1．开启软件平台——MA TLAB，开启MA TLAB编辑窗口；2．根据各种数值解法步骤编写M文件3．观察运行结果(数值或图形)；四、实验要求与任务根据实验内容和步骤，完成以下具体实验根据提示计算积分dx e x ⎰-12实验提示：1．Q uad函数的使用方法，参数意义Q = QUAD(FUN,A,B)Fun为待求积的函数A,B为积分区间Quad函数是使用误差为10-6的递归自适应Simpson公式计算Fun函数在A，B区间上的积分值。
例子（摘自帮助）： 求dx x x ⎰--2035211）首先，定义myfun 函数为待求积函数function y = myfun(x)y = 1./(x.^3-2*x-5); 2）计算数值积分Q = quad(@myfun,0,2);2． Q uad8函数的使用方法，参数意义Q = QUAD8(FUN,A,B) Fun 为待求积的函数 A,B 为积分区间Quad8函数是使用误差为10-3的递归自适应Newton-Cotes 8公式计算Fun 函数在A ，B 区间上的积分值。
它是收敛速度比较高的一种方法，有些版本已经将此函数改为QUADL具体使用方法同quad3． Q UADL 函数的使用方法，参数意义Q = QUADL(FUN,A,B) Fun 为待求积的函数 A,B 为积分区间Quad8函数是使用自适应Lobatto 求积公式计算Fun 函数在A ，B 区间上的积分值。
它也是收敛速度比较高的一种方法，具体使用方法同quad注意，使用上面三个函数的时候A ，B ，Fun 函数的输出必须是以下三种类型的 float:，double ，single实验4 范数计算和LU 分解一、实验目的[1] 了解矩阵范数的定义和常见的几种矩阵范数； [2] 掌握用直接法解线性方程组的常见几种方法；[3] 通过实例学习如何使用直接法中的LU 分解求解线性方程组。
[4] 掌握判断矩阵，线性方程组本身性质的方法——条件数的计算； 二、实验内容1．用MA TLAB 中的NORM 函数求矩阵的范数； 2．用MA TLAB 中的Cond 计算矩阵条件数；3．用MA TLAB 中的Chol 函数计算Cholesky 分解（平方根），Lu 函数计算Doolittle 分解。
三、实验步骤1．开启软件平台——MA TLAB ，开启MA TLAB 编辑窗口； 2．根据各种数值解法步骤编写M 文件 3．观察运行结果(数值或图形)； 四、实验要求与任务根据实验内容和步骤，完成以下具体实验实验： 1． 用H=Hilb （n ）函数生成一个n 阶hilbert 矩阵H ，（n=4） 2． 计算生成的矩阵H 的范数 3． 计算H 的条件数 4．⎪⎪⎪⎭⎫⎝⎛=6858137576A 用平方根法分解A5．⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=764142321,461561552621C B ，B ，C 是否存在三角分解，若存在则用LU分解实现。
实验5 牛顿迭代法一、实验目的[1] 了解求解非线性方程的解的常见方法； [2] 掌握用二分法，迭代法解线性方程组；[3] 通过实例学习用牛顿法分解求解非线性方程。
[4] 掌握判断迭代法，牛顿法收敛的条件，和收敛阶的计算；二、实验内容1． 用牛顿迭代法求方程0634223=-+-x x x 在区间（1，5）之间的根。
三、实验步骤1．用C 语言或matlab 完成编程 2．观察运行结果实验运行结果； 四、实验要求与任务根据实验步骤，完成以下具体实验实验： 1． 按照以下流程图完成牛顿迭代法x = newton( fname,dfname,x0,e,N )2． 调用牛顿迭代法的函数计算实验内容中第一题。
实验6龙格库塔法一、实验目的[1] 了解matlab中提供的求解常微分方程解的常见方法——龙格库塔法的函数ODE45,ODE32的使用方法；二、实验内容求解⎪⎩⎪⎨⎧=-=1)0(2'x x x三、实验步骤1．在command 窗口中输入help ode32 或者help ode45 查看相关的函数调用格式和使用算法的情况2．用函数求解实验内容的原函数，并画相应的图。
四、实验要求与任务根据实验步骤，先观察帮助里面提供的例子，然后再求解实验内容。