西安交大数学实验报告
西安交大概率论上机实验报告-西安交通大学概率论实验报告
概率论与数理统计上机实验报告一、实验内容使用MATLAB 软件进行验证性实验，掌握用MATLAB 实现概率统计中的常见计算。
本次实验包括了对二维随机变量，各种分布函数及其图像以及频率直方图的考察。
1、列出常见分布的概率密度及分布函数的命令，并操作。
2、掷硬币150次，其中正面出现的概率为0.5，这150次中正面出现的次数记为X ，(1) 试计算45=X 的概率和45≤X 的概率；(2) 绘制分布函数图形和概率分布律图形。
3、用Matlab 软件生成服从二项分布的随机数，并验证泊松定理。
4、设22221),(y x e y x f +-=π是一个二维随机变量的联合概率密度函数，画出这一函数的联合概率密度图像。
5、来自某个总体的样本观察值如下，计算样本的样本均值、样本方差、画出频率直方图。
A=[16 25 19 20 25 33 24 23 20 24 25 17 15 21 22 26 15 23 2220 14 16 11 14 28 18 13 27 31 25 24 16 19 23 26 17 14 30 21 18 16 18 19 20 22 19 22 18 26 26 13 21 13 11 19 23 18 24 28 13 11 25 15 17 18 22 16 13 12 13 11 09 15 18 21 15 12 17 13 14 12 16 10 08 23 18 11 16 28 13 21 22 12 08 15 21 18 16 16 19 28 19 12 14 19 28 28 28 13 21 28 19 11 15 18 24 18 16 28 19 15 13 22 14 16 24 20 28 18 18 28 14 13 28 29 24 28 14 18 18 18 08 21 16 24 32 16 28 19 15 18 18 10 12 16 26 18 19 33 08 11 18 27 23 11 22 22 13 28 14 22 18 26 18 16 32 27 25 24 17 17 28 33 16 20 28 32 19 23 18 28 15 24 28 29 16 17 19 18] 6. 利用Matlab 软件模拟高尔顿板钉试验。
西安交通大学概率论实验报告-蒙特卡洛法
西安交通大学实验报告课程：概率论与数理统计实验日期：2013/12/22报告日期：2013/12/24专业班级：姓名：学号：实验内容：用蒙特卡洛方法估计积分值要求：（1）针对要估计的积分选择适当的概率分布设计蒙特卡洛方法；（2）利用计算机产生所选分布的随机数以估计积分值；（3）进行重复试验，通过计算样本均值以评价估计的无偏性；通过计算均方误差（针对第1类题）或样本方差（针对第2类题）以评价估计结果的精度。
目的：（1）能通过 MATLAB 或其他数学软件了解随机变量的概率密度、分布函数及其期望、方差、协方差等；（2）熟练使用 MATLAB 对样本进行基本统计，从而获取数据的基本信息；（3）能用 MATLAB 熟练进行样本的一元回归分析。
1用蒙特卡洛方法估计积分2sinx xdxπ⎰，2xe dx+∞⎰和22221x yx ye dxdy++≤⎰⎰的值，并将估计值与真值进行比较。
1)2sinx xdxπ⎰用区间为0-π/2的均匀分布产生；代码如下N=10000;x=unifrnd(0,pi/2,N,1); mean(x.*sin(x)*pi/2)计算出10次的数值计算出精确值：syms x ;int(x.*sin(x),0,pi/2)精确值为1；计算出均值：1.00158计算出均方误差：0.0000637580结论：这是一个计算积分的很好的近似，误差很小。
接下来考虑计算第二个积分：2)考虑2xe dx +∞⎰由对称性可以考虑正态分布N（0,1），代码如下：N=10000;x=normrnd(0,1,N,1)0.5*mean((sqrt(2.*pi)).*exp(-x.^2./2))求出均值为0.88598取0.8860计算出均方误差为：0.000018204说明误差允许范围内，可以用其作为积分的近似。
若考虑用参数为1的指数分布E（1）代码为：N=10000;x=exprnd(1,N,1)mean(exp(-x.^2./2+x))精确值为：0.8862计算出平均值为：1.25164计算出均方误差为：0.13356381和正态分布比相去甚远，效果不如正态分布3)22221x yx ye dxdy++≤⎰⎰利用代码计算出积分：N=10000;x=unifrnd(0,1,N,1) //已经转换为极坐标，r在[0,1]取值，取[0,1]均匀分布2*pi*mean(x.*exp(-x.^2))计算出十个值为：计算出平均值为：1.98397计算出均方误差为：0.000059其值与精确值非常接近，可以作为一个很好的近似第二类题：4) dx e x ⎰102用如下代码计算：N=10000;x=unifrnd(0,1,N,1) //[0,1]上的均匀分布mean(exp(x.^2))计算出平均值为：1.4619计算出标准偏差为：0.003304 ，说明波动性较小计算出均方误差为：0.000010其值与精确值非常接近，可以作为一个很好的近似5)用如下代码计算：N=10000;x=unifrnd(0,2,N,1) //转换为极坐标后取[0,2]的均匀分布4*pi*mean(x./sqrt(1+x.^2))计算出平均值为：7.76363计算出标准偏差为：0.015241，说明波动性较小计算出均方误差为：0.000217其值与精确值非常接近，可以作为一个很好的近似22x y x d y +≤⎰⎰。
西交大数字信号处理实验报告
数字信号处理实验报告学院：班级：姓名：学号：西安交通大学实验报告课程 数字信号处理 实验日期 年 月 日专业班号 交报告日期 年 月 日 姓名 学号 共 21 页 第 1 页 实验1 常见离散信号的MATLAB 产生和图形显示 一、实验内容1.编制程序产生上诉5种信号（长度可自行输入确定），并绘出其图形。
2.讨论复指数序列的性质。
二、实验结果及源代码1.单位抽样序列⎩⎨⎧=01)(n δ≠=n n 在MATLAB 中可以利用ZEROS()函数实现。
;1)1();,1(==x N zeros x如果)(n δ在时间轴上延迟了K 个单位，得到)(k n -δ即：⎩⎨⎧=-01)(k n δ0≠=n kn （1）单位抽样序列源程序： n1=-10; n2=10;k=0; %延时k 个单位 n=n1:n2;N=length(n);%N 为序列长度 nk=abs(k-n1)+1; x=zeros(1,N); x(nk)=1;stem(n,x,'fill');axis([n1,n2,0,1.1*max(x)]); title('单位脉冲序列'); xlabel('时间'); ylabel('幅度');实验结果：（2）延时后的单位脉冲序列源程序：n1=-10;n2=10;k=input('k='); %延时k个单位 n=n1:n2;N=length(n);%N为序列长度nk=abs(k-n1)+1;x=zeros(1,N);x(nk)=1;stem(n,x,'fill');axis([n1,n2,0,1.1*max(x)]); title('单位脉冲序列');xlabel('时间');ylabel('幅度');实验结果（延时k=5）：2．单位阶跃序列⎩⎨⎧01)(n u 00<≥n n在MATLAB 中可以利用ones()函数实现。
西安交大数学建模实验报告
数学建模实验报告1，存货问题（一）问题描述某企业对于某种材料的月需求量为随机变量，具有如下表概率分布：每次订货费为500元，每月每吨保管费为50元，每月每吨货物缺货费为1500元，每吨材料的购价为1000元。
该企业欲采用周期性盘点的),(S s 策略来控制库存量，求最佳的s ，S 值。
（注：),(S s 策略指的是若发现存货量少于s 时立即订货，将存货补充到S ，使得经济效益最佳。
）（二）问题分析随机产生每个月需求量的概率，取遍每一个S 和s 的值，将每种S ，s 的组合对应的每月平均花费保存在数组money 里，筛选数组，选出其中费用最小值，并求出对应的S 和s 。
模拟400个月的生产情况。
（三）程序代码clear;clc;need=0; remain=0; cost=0; mincostavg=inf; forsl=30:10:70 forsh=80:10:140 fornum=1:100000m=rand; if m<=0.1 need=50;elseif m<=0.3 need=60;elseif m<=0.45 need=70;elseif m<=0.7 need=80;elseif m<=0.75 need=90;elseif m<=0.85 need=100;elseif m<=0.95need=110;elseneed=120;endif remain<slcost=cost+(sh-remain)*1000+500;ifsh<needcost=cost+(need-sh)*1500;remain=0;elsecost=cost+(sh-need)*50;remain=sh-need;endelseif remain<needcost=cost+(need-remain)*1500;remain=0;elsecost=cost+(remain-need)*50;remain=remain-need;endendendcostavg=cost/100000;ifcostavg<mincostavgmincostavg=costavg;propersl=sl;propersh=sh;endfprintf('s=%d, S=%d\nMonthly average cost=%.1f\n',sl,sh,costavg);cost=0;endendfprintf('\nWhen s=%d, S=%d\nThe least monthly average cost=%.1f\n',propersl,propersh,mincostavg);（四）运行结果s=30, S=80Monthly average cost=85466.9s=30, S=90Monthly average cost=87007.6Monthly average cost=87114.2 s=30, S=110Monthly average cost=87951.0 s=30, S=120Monthly average cost=86778.9 s=30, S=130Monthly average cost=86411.8 s=30, S=140Monthly average cost=86374.8 s=40, S=80Monthly average cost=83707.2 s=40, S=90Monthly average cost=84026.6 s=40, S=100Monthly average cost=85089.1 s=40, S=110Monthly average cost=85386.0 s=40, S=120Monthly average cost=86294.0 s=40, S=130Monthly average cost=85148.0 s=40, S=140Monthly average cost=84992.9 s=50, S=80Monthly average cost=83693.0 s=50, S=90Monthly average cost=82548.0 s=50, S=100Monthly average cost=82730.9 s=50, S=110Monthly average cost=83873.1 s=50, S=120Monthly average cost=84029.5 s=50, S=130Monthly average cost=84908.4 s=50, S=140Monthly average cost=84134.1 s=60, S=80Monthly average cost=83615.9 s=60, S=90Monthly average cost=82503.9 s=60, S=100Monthly average cost=81677.0Monthly average cost=81905.5s=60, S=120Monthly average cost=82946.0s=60, S=130Monthly average cost=83449.2s=60, S=140Monthly average cost=83871.3s=70, S=80Monthly average cost=83522.6s=70, S=90Monthly average cost=82525.8s=70, S=100Monthly average cost=81627.9s=70, S=110Monthly average cost=81323.3s=70, S=120Monthly average cost=82005.5s=70, S=130Monthly average cost=82601.6s=70, S=140Monthly average cost=82858.3When s=70, S=110The least monthly average cost=81323.3（五）结果分析用计算机模拟的结果和用数学分析的结果有一定的差异，由于计算机模拟时一般情况都是要简化模型的，所以在一定程度上会有所差异，我们可以考虑能不能通过改进算法来消除该差异，但对于一般的生产要求亦可以满足。
数学实验报告
西安交通大学实验报告一、某棉纺厂的原棉需从仓库运送到各车间，各车间的原棉需求量，单位产品从各仓库运往各车间的运输费以及各仓库的库存如表所列，问如何安排运输任务使得总运费最小？问题分析：该题较为简单，只要根据表中数据确定不等式，找到上下限，在根据书上的已有例子，综合自己的判断，就可写出。
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;0,0,0,0,0,0,-1,0,0,0,0,0,1,0,0,0,0,0;0,0,0,0,0,-1,0,0,0,0,0,0,0,1,0,0,0,0];b=[-19;6;-3;0;0;0;0;0;0;0;0];[x,favl]=bintprog(c,a,b)favl=-favl;结果分析：有实验结果可知，连选前10门课才可达到学校的要求。
西安交通大学数字信号处理实验报告
数字信号处理实验报告班级：硕姓名：学号：实验1 常见离散信号的MATLAB 产生和图形显示实验目的：加深对常用离散信号的理解；实验内容：（1）单位抽样序列clc;x=zeros(1,11); x(1)=1; n=0:1:10;stem(n,x, 'fill'); title('单位抽样序列'); xlabel('n'); ylabel('x[n]')延迟5个单位：clc;x=zeros(1,11); x(6)=1; n=0:1:10;stem(n,x, 'fill'); title('单位抽样序列'); xlabel('n'); ylabel('x[n]')nx [n ]（2）单位阶跃序列clc;x=[zeros(1,5),ones(1,6)]; n=-5:1:5;stem(n,x,'fill'); title('单位阶跃序列'); xlabel('n'); ylabel('x[n]');nx [n ]（3）正弦序列clc; N=50; n=0:1:N-1; A=1; f=1; Fs=50; fai=pi;x=A*sin(2*pi*f*n/Fs+fai); stem(n,x,'fill'); title('正弦序列'); xlabel('n'); ylabel('x[n]'); axis([0 50 -1 1]);nx [n ]（4）复正弦序列clc; N=50; n=0:1:N-1; w=2*pi/50; x=exp(j*w*n); subplot(2,1,1); stem(n,real(x)); title('复正弦序列实部'); xlabel('n');ylabel('real(x[n])'); axis([0 50 -1 1]); subplot(2,1,2); stem(n,imag(x)); title('复正弦序列虚部'); xlabel('n');ylabel('imag(x[n])'); axis([0 50 -1 1]);nx [n ]（5）指数序列clc; N=10; n=0:1:N-1; a=0.5; x=a.^n;stem(n,x,'fill'); title('指数序列'); xlabel('n'); ylabel('x[n]'); axis([0 10 0 1]);nr e a l (x [n ])ni m a g (x [n ])（6）复指数序列性质讨论：0(j )()enx n σω+=将复指数表示成实部与虚部为00()e cos j sin n n x n n e n σσωω=+1.当σ=0时，它的实部和虚部都是正弦序列。
西安交通大学数学建模上机实验报告
问题一某大型制药厂销售部门为了找出某种注射药品销量与价钱之间的关系,通过市场调查搜集了过去30个销售周期的销量及销售价钱的数据,如表.按照这些数据至少成立两个数学模型, 作出图形,比较误差。
问题分析：该问题是通过已知的过去30个销售周期的销量及销售价钱的 数据，来寻觅一个最能反映该药销量与价钱之间的函数曲 线。
在数学上归结为最佳曲线拟合问题。
大体思想：曲线拟合问题的提法:已知一组二维数据，即平面上的n 个点),x i i y （ i=1,2,3.....n ，i x 互不相同，寻求一个函数)(f y x =，使)(x f 在某中准则下与所有数据点最为接近，即曲线拟合得最好。
最小二乘法是解决曲线拟合最常常利用的方式.大体思路：1122 ()()()()m m f x a r x a r x a r x =+++令其中rk(x) 是事前选定的一组函数，ak 是待定系数(k=1,2,…,m,m ＜n), 拟合准则是使n 个点(xi,yi) (i=1,2…,n)，与y=f(xi)的距离 的平方和最小，称最小二乘法准则。
一、系数的肯定22111 (,,)[()]n nm ii i i i J a a f x y δ====-∑∑记求m a a ,,1 使得使J 达到最小.0 (1,,)kJ k m a ∂==∂ 取得关于 m a a ,,1 的线性方程组：11111()[()]0 ()[()]0nmi k k i i i k n mm i k k i i i k r x a r x y r x a r x y ====⎧-=⎪⎪⎪⎨⎪⎪-=⎪⎩∑∑∑∑ 1 ,,().m a a f x 解出，即得散点图： 程序： x=[,,,,,,,,,,,,,,,,,,,,,,,,,,,,,]; y=[,,,,,,,,,,,,,,,,,,,,,,,,,,,,,]; plot(x,y,'r.')通过观察，结合实际情形。
西安交通大学数学实验报告(用MATLAB绘制二维、三维图形)
实验报告（二）完成人：L.W.Yohann注：本次实验主要学习了用MATLAB绘制二维、三维图形的基本命令、图形的标识与修饰以及用符号函数绘图，在学习完成后小组对52页的上机练习题进行了程序编辑和运行。
1.绘制数列变化趋势图.解：在编辑窗口输入：n=1:100;an=(1+1./n).^n;plot(n,an,'r*')grid并保存，命名为lab1；在命令窗口中输入lab1,得：2.绘制数列变化趋势图.解：在编辑窗口输入：n=1:0.1:50;an=n.^(1./n);plot(n,an,'r*')grid并保存，命名为lab2；在命令窗口中输入lab2,得：3.绘制函数在无定义点处的变化趋势.解：在编辑窗口输入：x=-10:0.05:10;y=sin(x)./x;plot(x,y,'r*')grid并保存，命名为lab3；在命令窗口中输入lab3,得：4.在同一坐标系中画出函数及其Taylor多项式的图像解：y=sinx在编辑窗口输入：syms xf=sin(x);T6=taylor(f,x);T8=taylor(f,x,'Order',8);T10=taylor(f,x,'Order',10);T12=taylor(f,x,'Order',12);fplot([T6 T8 T10 T12 f])xlim([-8 8])grid onlegend('approximation of sin(x) up to O(x^6)',...'approximation of sin(x) up to O(x^8)',...'approximation of sin(x) up to O(x^{10})',...'approximation of sin(x) up to O(x^{12})',...'sin(x)','Location','Best')title('Taylor Series Expansion')并保存，命名为lab4sin；在命令窗口中输入lab4sin,得：y=exp(x)在编辑窗口输入：syms xf=exp(x);T6=taylor(f,x);T8=taylor(f,x,'Order',8);T10=taylor(f,x,'Order',10);T12=taylor(f,x,'Order',12);fplot([T6 T8 T10 T12 f])xlim([-8 8])grid onlegend('approximation of exp(x) up to o(x^6)',...'approximation of exp(x) up to o(x^8)',...'approximation of exp(x) up to o(x^{10})',...'approximation of exp(x) up to o(x^{12})',...'exp(x)','Location','Best')title('Taylor Series Expansion')并保存，命名为lab4exp；在命令窗口中输入lab4exp,得：5.符号函数绘图.注：在matlab r2010b 和matlab r2019b中对绘制函数图像的输入方法有不同的要求，故此类题分两个版本来求解。
西安交通大学第二次汇编上机实验报告
西安交通大学第二次汇编上机实验报告
实验时间5.31姓名
学号
实验内容概述主要对第五章所学的循环了分支结构进行了实践。
第一题循环结构比大小，把最大值的最大值的相对位置存储在数组之后的两位第二题循环结构找两数组之间的相同数字第三题分支结构对数组里出现的元素进行计数实验步骤
一（1）对数据求源码（负数取反后加一）（2）取每个数的绝对值（3）从前向后两两比较，如果当前数字比max的记录值大，那么更新当前数字的max，并存储当前位置pos二（1）根据数组a的长度设置外循环次数（2）根据数组b 的长度设置外循环次数（3）外循环遍历a的每一个元素，将a中每一个元素与b中元素比较，如果相同，则终止当前内循环，将该值存储在c数组中三（1）定义一个分支结构，先判断当前值，如果相等，则跳转至对应子程序：使存储该数字出现次数的位置上加一（2）程序运行结束时，在dos窗口输出每个元素出现的次数，
实验中遇到的问题及其解决方法1、mov指令用错当时报错后，查找了mov 指令的课本内容，想起来存储单元之间不能直接赋值，已解决2、逻辑正确，数据错误当时确认逻辑无误之后，我的内心非常纠结，甚至出现了想要重装软件的想法，后来利用u和t指令，一步步找错，发现是比例因子出现错误，将inccx 改为addcx，2后解决实验结果
结果均正确
2021
年
5月
31日
建议无备注无。
大一上数学实验报告
西安交通大学实验报告第 1 页 共 12 页课程实 验 日 期 2012年12月25日系 别实 验 报 告 日 期 2012年12月 28日 专业班级_报 告 退 发( 订正 、 重做 ) 姓 名______学号___教 师 审 批 签 字实验I 一、实验问题单位球面被柱面所截以及单位球面里挖掉柱面1222=++z y x 222)21()21(=+-y x二、问题分析该问题应分解为单位球面被柱面所截和单位球面里挖掉柱面两部分考虑，两者的相同之处在于均需要绘制单位球面，区别在于前者需要将柱面显示出来，而后者需要将柱面隐去，且需要用到find 和nan 指令对图像实现精确绘制。
三、程序设计1.单位球面被柱面所截x1=-1:0.01:1; y1=-1:0.01:1;[X1,Y1]=meshgrid(x1,y1); Z1=sqrt(1-(X1.^2+Y1.^2)); i1=find(X1.^2+ Y1.^2> 1 ); Z1(i1)= nan; x2=0:0.01:1; z2=-1:0.01:1;[X2,Z2]=meshgrid(x2,z2); Y2=sqrt(-X2.^2+X2); mesh(X1,Y1,Z1) hold onmesh(X1,Y1,-Z1) mesh(X2,Y2,Z2) mesh(X2,-Y2,Z2)x=-1:0.01:1;y=-1:0.01:1;[X,Y]=meshgrid(x,y);Z=sqrt(1-(X.^2+Y.^2));i1=find(X.^2+ Y.^2> 1 );Z(i1)= nan;i2=find(X.^2-X+Y.^2<=0);Z(i2)=nan;mesh(X,Y,Z)hold onmesh(X,Y,-Z)四、问题求解结果与结论1.单位球面被柱面所截五、问题的进一步拓展与实验本实验中最显著的问题在于单位球面的绘制：步长大，运算快，但图像缺陷多；步长小，图像细腻，但运算太慢。
西安交通大学数学实验报告(MATLAB求解开普勒方程和方程求根)
实验报告（五）完成人：L.W.Yohann注：本次实验主要学习了用MATLAB求解开普勒方程和方程求根的问题，了解学习了用fzero命令、二分法、Newton迭代法、一般迭代法求解方程，以及学习了非线性方程组的求解问题，完成后，小组对第90页的上级练习题进行了程序编辑和运行。
1.绘图并观察函数零点的分布.解：在编辑窗口输入：f=inline('x-0.5*sin(x)-1');fplot(f,[0,1])grid存盘后运行得2. 利用fzero 命令求解方程.解：在编辑窗口输入：f=inline('x-0.5*sin(x)-1');c=fzero(f,[1,2])存盘后运行得c =1.49873. 用二分法求解方程.求解（1）方程x^2-2=0在（0，2）内的近似根；（2）圆x^2+y^2=2与曲线y=e^-x 的两个交点；（3）方程∫t 21+t 2x 0dt =12的近似根. （1）解：在编辑窗口输入：00.10.20.30.40.50.60.70.80.91-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3f=inline('x^2-2');x1=0;x2=2;while abs(x1-x2)>10^(-5)x3=(x1+x2)/2;if f(x3)==0break;elseif f(x1)*f(x3)>0x1=x3;else f(x2)*f(x3)>0;x2=x3;endendx0=x3存盘后运行得x0 =1.4142（2）解：在编辑窗口输入：f=inline('(x^2)*exp(2*x)+1-2*exp(2*x)');x1=0;x2=2;x5=-2;x6=0;while abs(x1-x2)>10^(-5)x3=(x1+x2)/2;if f(x3)==0break;elseif f(x1)*f(x3)>0x1=x3;else f(x2)*f(x3)>0;x2=x3;endendwhile abs(x5-x6)>10^(-5)x7=(x5+x6)/2;if f(x7)==0break;elseif f(x5)*f(x7)>0x5=x7;else f(x6)*f(x7)>0;x6=x7;endendx0=x3x4=x7存盘后运行得x0 =1.3922x4 =-0.3203（3）解：在编辑窗口输入：clear;clc;syms t xf1=(t^2)/(1+t^2);f2=int(f1,t,0,x);%¼ÆËã²»¶¨»ý·Öf=inline('x - atan(x)-0.5');x1=-5;x2=5;while abs(x1-x2)>10^(-5)x3=(x1+x2)/2;if f(x3)==0break;elseif f(x1)*f(x3)>0x1=x3;else f(x2)*f(x3)>0;x2=x3;endendx0=x3存盘后运行得x0 =1.47504.用Newton迭代法求解方程求解：x=0.5sinx+1的近似根；解：在编辑窗口输入：f=inline('x-0.5*sin(x)-1');df=inline('1-0.5*cos(x)');d2f=inline('0.5*sin(x)');a=1;b=2;dlt=1.0e-5;if f(a)*d2f(a)>0x0=a;elsex0=b;endm=min(abs(df(a)),abs(df(b)));k=0;while abs(f(x0))>m*dltk=k+1;x1=x0-f(x0)/df(x0);x0=x1;fprintf('k=%d x=%.5f\n',k,x0); end存盘后运行得k=1 x=1.54858k=2 x=1.49933k=3 x=1.498705.求解非线性方程组.试求非线性方程组{2x12−x1x2−5x1+1=0x1+3lgx1−x22=0的解，初值如下：（1）x0=[1.4,−1.5]（2）x0=[3.7,2.7]解：在编辑窗口输入：function f=group5(x)f=[2*x(1)^2-x(1)*x(2)-5*x(1)+1;x(1)+3*log10(x(1))-x(2)^2];(1)：输入：[f,fval]=fsolve('group2',[1.4,-1.5]) 运行得f =1.4589 -1.3968fval =1.0e-011 *0.0759-0.6178(2)：输入：[f,fval]=fsolve('group2',[3.7,2.7])运行得f =3.4874 2.2616fval =1.0e-006 *0.0059-0.20126.解决实际问题.为了在海岛I与某城市C之间铺设一条地下光缆，每千米光缆铺设成本在水下部分使C1万元，在地下部分使C2万元，为使得该光缆的总成本最低，光缆的转折点P（海岸线上）应该取在何处？如果实际测得海岛I与城市C之间的水平距离l=30km,海岛距海岸线垂直距离h1=15km，城市距海岸线垂直距离h=10km,C1=3000万元/km,C2=1500万元/km,求P点的坐标（误差<10−3km）.解：在编辑窗口输入：f=inline('(3000*x)/(x^2 + 225)^(1/2) + (750*(2*x - 60))/((x - 30)^2 + 100)^(1/2)'); x1=5;x2=10;while abs(x1-x2)>10^(-3)x3=(x1+x2)/2;if f(x3)==0break;elseif f(x1)*f(x3)>0x1=x3;else f(x2)*f(x3)>0;x2=x3;endendfprintf('x=%.5f',x3) 存盘后运行得x=7.69104>>。
西安交大概率论与数理统计实验报告
西安交大概率论与数理统计实验报告——蒙特卡洛算法计算积分姓名：学号：班级一、实验目的（1）能通过 MATLAB 或其他数学软件了解随机变量的概率密度、分布函数及其期望、方差、协方差等；（2）熟练使用 MATLAB 对样本进行基本统计，从而获取数据的基本信息；（3）能用 MATLAB 熟练进行样本的一元回归分析。
二、实验要求（1）针对要估计的积分选择适当的概率分布设计蒙特卡洛方法；（2）利用计算机产生所选分布的随机数以估计积分值；（3）进行重复试验，通过计算样本均值以评价估计的无偏性；通过计算均方误差（针对第1类题）或样本方差（针对第2类题）以评价估计结果的精度。
三、实验原理1. 蒙特卡洛法的思想简述当我们所求解问题是某种随机事件出现的概率，或者是某个随机变量的期望值时，通过某种“实验”的方法，以这种事件出现的频率估计这一随机事件的概率，或者得到这个随机变量的某些数字特征，并将其作为问题的解。
有一个例子我们可以比较直观地了解蒙特卡洛方法：假设我们要计算一个不规则图形的面积，那么图形的不规则程度和分析性计算（比如，积分）的复杂程度是成正比的。
蒙特卡洛方法是如下计算的：假想有一袋豆子，把豆子均匀地朝这个图形上撒，然后数这个图形之中有多少颗豆子，这个豆子的数目就是图形的面积。
当豆子越小，撒的越多的时候，结果就越精确。
在这里我们要假定豆子都在一个平面上，相互之间没有重叠。
2. 蒙特卡洛法与积分通常蒙特卡洛方法通过构造符合一定规则的随机数来解决数学上的各种问题。
对于那些由于计算过于复杂而难以得到解析解或者根本没有解析解的问题，蒙特卡洛方法是一种有效的求出数值解的方法。
一般蒙特卡洛方法在数学中最常见的应用就是蒙特卡洛积分。
非权重蒙特卡洛积分，也称确定性抽样，是对被积函数变量区间进行随机均匀抽样，然后对被抽样点的函数值求平均，从而可以得到函数积分的近似值。
此种方法的正确性是基于概率论的中心极限定理。
3. 本实验原理简述在本实验中，我们主要是计算积分值与误差比较。
数学实验实验报告
西安交通大学实验报告课程名称：数学实验实验名称：数值积分及软件实现学 院：航空航天学院 实验日期：2012年12月21日星期五班 级：力硕21 姓 名__汪振兴 _学号：21217020161．实验目的通过对实际问题的分析，运用数学方法解决问题，从而考察学生的数学素养、解决实际问题的动手能力，考评学生对MATLAB 软件的了解、掌握和运用能力2．实验任务与问题背景如下图所示的是某矿区的平面示意图：在选取自西向东的方向为x轴方向，自南向北的方向为y轴正向的情况下，测得该选定矿区边界上的一些数据如下表所示：矿区坐标数据8.612.021.835.442.945.661.669.8X 2.0(最西端)Y114.521.622.425.820.623.836.828.830.4Y294.576.891.5118.6124.2140.6172.8206.6175.5X120.5146.8168.5189.8220.4255.0286.6302.5350.1（最东端）Y152.898.655.569.423.451.946.852.866.6Y2145.6186.8352.4423.8386.7415.4358.2245.2198.9请你用多种方法估计该矿区的平面面积。
3.问题分析这是一道让我们求面积的题目，而且该矿区可以视为一个平面图形，并且题目给出了他的边界坐标，显然我们可以采用数值积分的方法解决这道问题。
但与一般求数值积分的题目不同的是，这个平面图形的边界并非由一条连续的光滑曲线（重点是有具体表达式），而是由一个个的点连成的折线。
但这并不影响我们对问题的求解，因为每相邻的两个点就确定了一条直线的直线方程。
比如说由（x(k),y(k)）,（x(k+1),y(k+1)）确定的直线方程为：Y={[y(k+1)-y(k)]/[x(k+1)-x(k)]}*[X-x(k)]+y(k)。
因此我们就可以用直线x=x(k),x=x(k+1)(k=1,2,3……)把整个矿区划分成一块块，每一块都可以用数值积分方法，比如矩形数值积分公式和梯形数值积分公式求出面积。
西安交大数学实验报告
数学实验报告作者：学号：班级：题目（一）：求函数 2x 11+的积分程序：clc;clear;n=100;x=0:1/n:1;sum=0;for i=1:n;sum=sum+(1/n).*((1./(1+x(i).^2))+1./(1+x(i+1).^2))*0.5; endsum结果:结果如图所示，为0.7854程序分析：本程序运用了梯形法求函数在区间[]10 的积分，将定义域分为100分，并取其微元的值。
我们知道，函数的实际积分为为π/4，即0.7853，与所运行的结果较类似。
题目（二）：用lagrange差值法求sin(49º)的近似值程序：在进行这个问题前，我们先编写一个lagrange函数文件function p=lagrange(x,y)L=length(x);A=ones(L);for j=2:LA(:,j)=A(:,j-1).*x';endX=inv(A)*y';for i=1:Lp(i)=X(L-i+1);end然后将数据导入clc;clear;x=[40 45 50 55 60];y=[0.64278 0.70710 0.76604 0.81915 0.86602];p=lagrange(x,y);t=49;u=polyval(p,t);u程序分析：由计算器可以算得sin(49º)的值为0.7547，与计算结果相同。
本题通过利用lagrange插值法，对原有数据利用多项式进行拟合，同时利用所求得的曲线对数值进行预测实验反思与结论：第一次做数学实验，第一次接触MATLAB软件，第一次进行简单编程，有许多错误和茫然指出，我在不断的尝试努力中进步着。
通过这次数学实验，我基本掌握了积分算法和插值法的运算同时还有各种赋值语句。
当然还有MATLAB软件中的一些特殊符号，对实验理论课上学的一些基本操作更加熟练。
最重要的是，很好的提高了我的数学逻辑思维，以及上机动手实验操作的能力，这对我以后的学习工作有很大帮助。
西安交通大学数学实验报告模板
二、
【实验题目】
煤矿的储量估计,下表给出了某露天煤矿在平面矩形区域(800m 600m)上,在纵横均匀的网 格交点处测得的煤层厚度(单位:m)(由于客观原因,有些点无法测量煤层厚度,这里用/标出),其 中的每个网格都为(10m 8m)的小矩形,试根据这些数据,来估算出该矩形区域煤矿的储藏量(体 积) A / / / 7.5 8.9 / / B / / 12 12.6 7.8 / / C 12.5 / 13.5 14.9 12.4 / 8.6 D 13.5 15.6 13.5 18.7 13.5 13.7 11.8 E 17.2 18.2 17.8 17.7 15.7 13.6 12.5 F / 13 16.9 17.5 17.6 16.5 11.3 G H 8.8 14.7 6.4 8.9 13.2 / 14.7 13 11.7 9.6 12.5 8.7 13.4 / I 8.0 9.2 / / 9.2 9.7 / J 13.0 11.7 / / 9.5 / / K / / / 6.5 8.6 / /
第三行: y=[12 13.5 13.5 17.8 16.9 13.2]; x=[80 160 240 320 400 480 ]; plot(x,y,'k.','markersize',15) axis([0 800 0 20]) grid; hold on t=0:80:800; p=polyfit(x,y,1); u=polyval(p,t) plot(t,u,'r--') u =12.4333 17.1190 13.0190 13.6048 14.1905 14.7762 15.3619 15.9476 16.5333
h
0
x 2dy ,
西安交通大学数学实验报告模板
成绩 西安交通大学实验报告课 程________概率论与数理统计__________________ 实验日期___2016.12.11________________________专业班号_物理51_____________________ 姓 名 _____________李淏淼_____________学 号_________2150900015_________________一、 实验问题1某大米生产厂将产品包装成1000克一袋出售，在众多因素的影响下包装封口后一袋的重量是随机变量，设其服从正态分布N(m ，)，其中σ已知，m 可以在包装时调整，出厂检验时精确地称量每袋重量，多余1000克的仍按1000克一袋出售，因而厂家吃亏；不足1000克的直接报废，这样厂方损失更大，问如何调整m 的值使得厂方损失最小？二、 问题分析（涉及的理论知识、数学建模与求解的方法等）设定x 为产品包装后的重量，依题意x 为一随机变量，且服从正态分布N ，概率密度函数为f （x ）当成品重量M 给定后，记：P 为x 大于等于M 的概率P ’为x 小于M 的概率故而有： P ＋P’＝1分析题意可知，厂方损失Y 由两部分组成：（1）x≥L 时，多余部分，重量为（x －L ）；（2）x<L 时，整袋报废，重量为x ；Y ＝()()()MM x M f x dx xf x dx ∞-∞-+⎰⎰＝m －MP生产N 袋大米报废总量为Nm －NMP成品袋数为NP则成品中，平均每袋损失的重量为J=mN MPN m M PN P-=- 求J 的最小值即可三、 程序设计1. 在MATLAB 中建立文件Jmin.m function J=Jmin(m)J=m/(1-normcdf( (1000-m),0,1));2. 在Matlab 的Medit 窗口建立文件figer.mfor m=1000:0.001:1020J=Jmin(m);plot(m,J)hold onend可得出函数图像根据图像，可知函数在该区间存在最小值3.在Matlab的Medit建立文件zuixiaozhi.mmin=1100;minm=0;for m=1000:0.001:1010J=Jmin(m);if J<=minmin=J;minm=m;endendminm,min运行程序得出结果为四、问题求解结果与结论m的值为1003.5时，厂方损失最小五、问题的进一步拓展与实验m的值为1003.5时，平均每袋的损失为多少？六、实验问题2设(X, Y)的联合分布律为求X与Y的协方差及相关系数。
西安交通大学概率论上机实验报告
西安交通大学一、试验目的概率论部分1.了解matlab软件的基本命令与操作；2.熟悉matlab用于描述性统计的基本菜单操作及命令；3.会用matlab求密度函数值、分布函数值、随机变量分布的上下侧分位数。
数理统计部分1.熟悉matlab进行参数估计、假设检验的基本命令与操作.2.掌握用matlab生成点估计量值的模拟方法3.会用matlab进行总体数学期望和方差的区间估计。
4.会用matlab进行单个、两个正态总体均值的假设检验。
5.会用matlab进行单个、两个正态总体方差的假设检验。
二、试验问题实验五、随机变量综合试验实验内容1. 产生(6),(10),F(6,10)和t(6)四种随机数，并画出相应的频率直方图；2. 在同一张图中画出了N(0,1)和t(6)随机数频率直方图，比较它们的异同；3. 写出计算上述四种分布的分布函数值和相应上侧分位点命令.实验七、对统计中参数估计进行计算机模拟验证实验内容：1.产生服从给定分布的随机数，模拟密度函数或概率分布；2.对分布包含的参数进行点估计，比较估计值与真值的误差；3. 对分布包含的参数进行区间估计，行区间估计，可信度。
三、实验源程序及结果实验5源程序：% 清空内存，清空输出屏幕clc;clear;% 首先是指数分布n = normpdf(-2::14,6);% 绘制频率直方图plot(-2::14,n,'color','r','linewidth',2);ylabel('概率密度');title('正态分布概率密度');% t分布h1 = figure;t = tpdf(-3::3,6);plot(-3::3,t,'color','g','linewidth',2);ylabel('对应频率');title('t分布频率密度');% F分布h2 = figure;f = fpdf(0::10,6,10);plot(0::10,f,'color','k','linewidth',2); ylabel('对应频率');title('F分布频率直方图');% 卡方分布h3 = figure;ka = chi2pdf(0::15,6);plot(0::15,ka,'color','y','linewidth',2); ylabel('对应频率');title('卡方分布频率直方图');% 再来绘图h4 = subplot(2,1,1);y1=normpdf(-10::10,0,1);plot(-10::10,y1,'color','b','linewidth',2); title('N(0,1)');h5 = subplot(2,1,2);t1 = tpdf(-10::10,6);plot(-10::10,t1,'color','r','linewidth',2);%上侧分位数norminv,0,1)tinv,6)chi2inv,6)finv,6,10)运行结果：正态分布T分布F分布N(0,1)和t(6)随机数频率直方图四种分布的分布函数值和相应上侧分位点实验7源程序：% 以正太分布为例% 清空内存，清空输出屏幕clc;clear;y=normrnd(10,1,10000,1);ymin=min(y);ymax=max(y);x=linspace(ymin,ymax,80);yy=hist(y,x);yy=yy/10000;bar(x,yy);grid;xlabel('(a)¸概率密度分布直方图 ');phat=mle(y,'distribution','norm','alpha',%对分布函数参数进行区间估计，并估计区间的可信度 [mu,sigma,m_ci,s_si]=normfit(y,运行结果：正态分布概率密度分布直方图得到估计参数m =σ =由上可知估计的m = ，而实际是 10。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验报告班级：姓名：学号：第一次实验代码：a = 1 + 3i;b=2 - i;a+ba-ba*ba/b运行结果：代码：x=-4.5/180*pi;y=7.6/180*pi;./,,,,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 ，其中、计算sin(abs(x)+y)/sqrt(abs(x+y))运行结果：（1）代码：x=0:0.01:2*piy1=sin(x);y2=cos(x);y3=exp(x);y4=log(x);plot(x,y1,x,y2,x,y3,x,y4);运行结果：（2）代码：x=0:0.01:3的图形。
下分别绘制）同一页面四个坐标系）同一坐标系下（、在（x y e y x y x y x ln ,,cos ,sin 214====y1=sin(x);y2=cos(x);y3=exp(x);y4=log(x);subplot(2,2,1)plot(x,y1)subplot(2,2,2)plot(x,y2)subplot(2,2,3)plot(x,y3)subplot(2,2,4)plot(x,y4)运行结果：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=randn(3,3)B=randn(3,2)AB=A*BC=B.^2sinB=sin(B) detA=det(A)Ani=A^(-1)X=A^(-1)*BA(2,:)=A(2,:)+1 D=ADet=det(D)运行结果：6、设y=(x2 + e x cosx+[x])/x,分别计算x=1, 3, 5, 7.4时y的值。
其中[x]表示x的取整函数。
代码：X=[1,3,5,7.4];Y=(X.^2+exp(X).*cos(X)+fix(X))./X;y1=Y(1)y2=Y(2)y3=Y(3)y4=Y(4)运行结果：7、已知某地区1-12月份的平均气温为5，8，13，20，25，28，31，33，27，24，19，10，请绘图表示。
代码：y=[5,8,13,20,25,28,31,33,27,24,19,10];x=[1,2,3,4,5,6,7,8,9,10,11,12];plot(x,y)xlabel('月份')ylabel('平均气温')grid on运行结果：5、画出半径为的圆的图形。
2代码：axis equalezplot('x^2+y^2=4')运行结果：第二次实验1.某大学本科生就业情况如下：458人考入研究生，60人签到外企或合资企业，184人签到国营大企业，87人签到私人企业，13人自主创业，画出饼图表示之；代码：x=[60 184 87 13 458];pie(x)pie(x,explode)name={'考入研究生458人','签到外企或合资企业60人','签到国营大企业184人','签到私人企业87人','自主创业13人'};pie(x,name)运行结果：2. 求极限x x ax a x )(lim +-∞→ 代码：syms a x ;y=((x-a)/(x+a))^x; limit(y,x,inf)运行结果：3.求极限x x x ln 10)(tan lim +→代码：syms x ;y=(tan(x))^(1/log(x)); limit(y,x,0,'right')运行结果：4.设)(',)11(lim )(2t f xt f x t x 求+=∞→ 代码：syms x t ;f=(1+1/x)^(2*t*x); f1=limit(f,x,inf); diff(f1,t)运行结果：5.展开多项式b a b a y 2)1()1(23++-++=； 代码：syms a b ;fc=(a+1)^3+(b-1)^2+2*b; fc1=expand(ps)运行结果：6.分解因式123245-++-=x x x x y 代码：syms x ;y=x^5-3*x^4+2*x^2+x-1; factor(y)运行结果：7.求方程0123=+-x x 的根代码：sym x ;f=sym('x^3-2*x+1');solve(f)运行结果：8.代码：syms x;y=sqrt(x+sqrt(x+sqrt(x)));diff(y,x) 运行结果：9. 求不定积分代码：syms x ; f=sym('1/(sqrt(2*x+3)+sqrt(2*x-1))');int(f)运行结果：10.代码：syms x ;f=sym('1/(x*sqrt(log(x)*(1-log(x))))');.的导数求函数x x x y ++=dx x x ⎰-++12321.cos 11⎰+dx xint(f,x,(exp)^(1/2),(exp)^(3/4))运行结果：10 求定积分代码：syms x ;f=sym('1/(x*sqrt(log(x)*(1-log(x))))'); int(f,x,(exp(1))^(1/2),(exp(1))^(3/4))运行结果：11.解方程组⎩⎨⎧=-=+1382yxyx代码：A=[2,1;1,-3];B=[8;1];X=A\B运行结果：12.求和∑=20121kk代码：syms k;A=symsum(1/(k^2),k,1,20)运行结果：13..156x2cos次泰勒多项式处的在求π=x代码：syms x;y=taylor(cos(2*x),x,1/6*pi,'Order',15).)ln1(ln43⎰-ee xxxdx运行结果：y =(2*3^(1/2)*(x - pi/6)^3)/3 - 3^(1/2)*(x - pi/6) - (2*3^(1/2)*(x - pi/6)^5)/15 + (4*3^(1/2)*(x - pi/6)^7)/315 - (2*3^(1/2)*(x - pi/6)^9)/2835 + (4*3^(1/2)*(x - pi/6)^11)/155925 - (4*3^(1/2)*(x - pi/6)^13)/6081075 - (x - pi/6)^2 + (x - pi/6)^4/3 - (2*(x - pi/6)^6)/45 + (x - pi/6)^8/315 - (2*(x - pi/6)^10)/14175 + (2*(x - pi/6)^12)/467775 - (4*(x - pi/6)^14)/42567525 + 1/214. 编写猜数游戏程序：首先由计算机随机产生一个[1,100] 之间的一个整数，然后由用户猜测所产生的这个数。
根据用户猜测的情况给出不同的提示，如果猜测的数大于产生的数，则显示“High” ，小于则显示“ Low ” ，等于则显示“You won！”，同时退出游戏。
用户最多有7 次机会。
代码：x=fix(rand(1)*100);for i=1:7y=input('ÇëÊäÈëÒ»¸öÕûÊý');if x==yfprintf('you won!')breakelseif y>xfprintf('High\n')elseif y<xfprintf('Low\n')endclear yend运行结果：第三次实验1.每门课程考试阅卷完毕,任课教师都要对各班的考试成绩进行统计,统计内容包括:全班人数,总得分,平均得分,不及格的人数及90分(包括90分)以上的人数.请编制程序解决这一问题,并自给一组数据验证程序的正确性.要求:使用者在提示下通过键盘输入学生成绩,计算机自动处理后,显示需要的结果.代码：x=input('请输入全班成绩');l=length(x);s=sum(x);a=mean(x);n=0;m=0;for i=1:lif x(i)>90n=n+1;elseif x(i)<60m=m+1;endendfprintf('学生总数=%i 总得分=%.1f 平均得分=%1.f 不及格人数=%i 90分以上人数=%i',l,s,a,m,n)运行结果：2.《数学实验》P67 练习1 第2题代码：a=input('请输入一个大于6的偶数');m=1;n=a-m;x=isprime(m);y=isprime(n);while x+y<2m=m+1;n=a-m;x=isprime(m);y=isprime(n);endfprintf('两个质数为a=%.i b=%.i',m,n)运行结果：3. 编写猜数游戏程序：首先由计算机随机产生一个 [1,100] 之间的一个整数，然后由用户猜测所产生的这个数。
根据用户猜测的情况给出不同的提示，如果猜测的数大于产生的数，则显示 “High” ，小于则显示 “ Low ” ，等于则显示 “You won ！”，同时退出游戏。
用户最多有 7 次机会。
代码：a=rand;b=fix(a*100);for n=1:1:7c=input('请输入一个1到100间的随机数c=');if c>bdisp ('high ')elseif c<bdisp('low ')else c==bdisp('you win!')endend运行结果：编写程序解决下列问题4.Fibonacci 数组元素满足Fibonacci 规则： );,2,1(,12 =+=++k a a a k k k .121==a a 其中（1）用循环语句形成前20个分量的数组；（2）求该数组中第一个大于10000的数.代码：(1)a(1)=1;a(2)=1;for n=1:18a(n+2)=a(n+1)+a(n);enda(2)a(1)=1;a(2)=1;for n=1:22a(n+2)=a(n+1)+a(n);endfor i=1:22while a(i)>10000fprint(‘第一个大于10000的数为：%.4f\n’,a(i)) returnendend运行结果：（1）（2）5.（1）编写一个function函数，调用该函数，可以求一组数的最大值及最小值.（2）通过键盘输入一组数，利用（1）编写的function函数，得出最大最小数.代码:(1)function maxandmin(a)max=a(1);n=length(a)for i=2:nif max<a(i)max=a(i);endendfprintf(' 这组数最大的为%.5f\n ',max);min=a(1);for i=2:nif min>a(i)min=a(i);endendfprintf('这组数最小的为\n',min)(2)x=input('请输入一组数');maxandmin(x)实验结果：。