数学实验-实验指导书

1、数学实验2013年2月实验1 matlab基本特性与基本运算【实验目的】了解Matlab基本特性与基本运算【实验要求】1、熟悉MATLAB语言编程环境；2、熟悉MATLAB语言命令；3、熟悉Matlab基本运算命令【实验原理】MATLAB 是以复杂矩阵作为基本编程单元的一种程序设计语言。它提供了各种矩阵的运算与操作，并有较强的绘图功能。1.1 基本规则(1) 一般MATLAB命令格式为 输出参数1，输出参数2，=（命令名）（输入参数1，输入参数2，） 输出参数用方括号，输入参数用圆括号如果输出参数只有一个可不使用 括号。 (2) %后面的任意内容都将被忽略，而不作为命令执行，一般用于为代码加注

2、释。(3) 可用、键来重现已输入的数据或命令。用、键来移动光标进行修改。(4) 所有MATLAB命令都用小写字母。大写字母和小写字母分别表示不同的变量。(5) 常用预定义变量，如pi 、Inf、NaN、ans(6) 矩阵的输入要一行一行的进行，每行各元素用空格或“，”分开，每行用“；”分开。如 MATLAB书写格式为A=1 2 3 ；4 5 6 ；7 8 9 在MATLABZ中运行如下程序可得到A矩阵 a=1 2 3;4 5 6;7 8 9 a = 1 2 3 4 5 6 7 8 9(7) 需要显示命令的计算结果时，则语句后面不加“；”号，否则要加“；”号。 运行下面两种格式可以看出它们的区别

3、：a=1 2 3;4 5 6;7 8 9 a=1 2 3;4 5 6;7 8 9;a= 1 2 3 不显示结果 4 5 6 7 8 9(8) 当输入语句过长需要换行时，应加上“”后再回车，则可续行输入。1.2 文件管理常用命令(1) 帮助（HELP）命令 MATLAB有很多命令，因此很不容易记忆。使用HELP命令可以得到有关命令的屏幕帮助信息。如 在MATLAB环境下直接运行HELP命令就会在屏幕上给出MATLAB的基本命令以及相关的工具箱软件命令（请自行验证）。 而下面的格式可以得到指定命令的帮助信息，和详细说明。 help 命令名 (2) save 命令该命令存储定义的变量或演算结果，也可

4、以用来存储指定的变量。 命令格式为 save 文件名 变量名表 -append -ascii (3) what 命令 该命令可以在当目录下显示MATLAB文件和MAT数据文件(4) dir 命令 显示当前目录下的所有文件.(5) clear命令 如果输入 calear a b c ，则表示清除工作空间中指定变量a,b,c；如果仅仅输入calear命令，则清除整个工作空间。 1.3 矩阵运算 MATLAB具有强大的矩阵运算功能,这里我们只作简单介绍.(1) 在MATLAB中表示一个矢量要用方括号， 而列矢量的输入只需在行矢量输入格式基础上加转置符“”即可。 如 x=1 2 3;4 5 6 x =

5、 1 2 3 4 5 6 而 x=1 2 3;4 5 6' (加转置符) x = 1 4 2 5 3 6(2) 下面三条命令可以产生一个行矢量1 a=linspace(x,y,n) 2 a=logspace(x,y,n) 3 a=x:n:y第一条命令可以在线性空间产生一个值在x至y之间间隔点数为n的行矢量(一组数据)。其行矢量的起始值是x, 终值为y, 点数为n。第二条命令可以在对数空间产生一个值在10x至 10y之间等间隔的行矢量(一组数据)。第三条命令产生x至y步长为n的行矢量。但是, 三个命令之间存在差别，下面的例子可以说明这一点。 例1 x=logspace(0,5,6) x

6、= 1 10 100 1000 10000 100000 例2 x=linspace(0,10,11) x = 0 1 2 3 4 5 6 7 8 9 10 例3 x=0:1:10 x = 0 1 2 3 4 5 6 7 8 9 10 通过上面三个例子可以看出例一, 例二中n代表选取的点数。而在例三 中n则表示步长. 注意它们的区别。(3) 矩阵的加, 减, 乘, 除等, 和其它语言书写一样。 但要注意的是在运算 符前面加有“.”则表示是数组运算，为对应元素的相应操作。 例4 请注意下面的不同情况. A=1 2;3 4; b=2 2;2 2;c=a*b c=a.*b c = c = 6 6 2

7、 4 14 14 6 8 c=a/b c=a./b c = c = 1.0e+015 * 0.5000 1.0000 -2.2518 2.2518 1.5000 2.0000 -2.2518 2.25181.4 常见数学函数函 数 名数 学 计 算 功 能函 数 名数 学 计 算 功 能abs（x）实数的绝对值或复数的幅值floor（x）对x朝-方向取整acos（x）反余弦arcsingcd（m，n）求正整数m和n的最大公约数acosh（x）反双曲余弦arccoshimag（x）求复数x的虚部angle（x）在四象限内求复数 x 的相角lcm（m，n）求正整数m和n的最小公倍数asin（x）反

8、正弦arcsinlog（x）自然对数（以为底数）asinh（x）反双曲正弦arcsinhlog10（x）常用对数（以10为底数）atan（x）反正切arctanreal（x）求复数x的实部atan2（x,y）在四象限内求反正切rem（m，n）求正整数m和n的m/n之余数atanh（x）反双曲正切arctanhround（x）对x四舍五入到最接近的整数ceil（x）对x朝+方向取整sign（x）符号函数：求出x的符号conj（x）求复数x的共轭复数sin（x）正弦sincos（x）余弦cossinh（x）反双曲正弦sinhcosh（x）双曲余弦coshsqrt（x）求实数x的平方根：exp（x）

9、指数函数 tan（x）正切tanfix（x）对x朝原点方向取整tanh（x）双曲正切tanh例5：输入 x=-4.85 -2.3 -0.2 1.3 4.56 6.75，则： ceil(x)= -4 -2 0 2 5 7 fix(x) = -4 -2 0 1 4 6 floor(x) = -5 -3 -1 1 4 6 round(x) = -5 -2 0 1 5 71.4 关系及逻辑运算 (1) 关系运算允许常量(或矩阵中的元素与元素) 之间的比较. 如果比较结果为真,则答案为1, 否则为0. 常用的关系运算符如下:< 小于 , > 大于 , = 等于 <= 小于等于, >

10、;= 大于等于, = 不等于对于复数”=” 和”=” 既比较实部也比较虚部, 而其它运算仅比较实部。例6 比较 a=5:1:15 与 b=1 2 8 8 7 10 12 11 13 14 15 中的对应元素 是否相同,则运行下面程序a=5:1:15;b=1 2 8 8 7 10 12 11 13 14 15;y=a=b y = 0 0 0 1 0 1 0 0 1 1 1 (2) MATLAB有三个逻辑运算符. & 逻辑与 | 逻辑或 逻辑非在矩阵运算中,以上命令是将两个矩阵中的对应元素进行运算, 得到的结 果是具有同样元素的矩阵。例7 分别对下面两个矩阵进行逻辑运算. X=0 1;1

11、0; Y=0 0;1 0; X&Y ans = 0 0 1 0 X|Y ans = 0 1 1 0 X ans = 1 0 0 1 MATLAB在进行混合运算时, 其顺序为: 算术运算 关系比较 逻辑运算例8 计算 x&y+x>y 时, 程序首先执行x+y运算, 然后将相加结果与y比较, 最后将比较结果和x作与运算，运行结果如下。 x=0 1;1 0; y=0 0;1 0; x&y+x>y ans = 0 1 1 0【实验内容】熟悉MATLAB语言编程环境。在这一环境中，系统提供了许多编写，调试和执行MATLAB程序的便利工具。熟悉变量、MATLAB基本数学

12、函数使用，矩阵操作练习，主要包括矩阵的建立，算数运算，逻辑运算和关系运算。（一）熟悉Matlab的运行环境以及HELP命令，运行demo命令，观看matlab所具有的功能。（二）验证例题中所有操作。（三）完成以下各题，将步骤命令写入实验报告，并在机器上运行结果验证。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*B A.*B A2 A.2 B/A B./A2输入C=1:2:20，则C（i

13、）表示什么？其中i=1,2,3,10；3查找已创建变量的信息，删除无用的变量；4输入以下程序，察看结果，如果出错给出处理方法并验证。（1）a=magic(3);b=243; （2）a=magic(3);b=243;save myvar a b; save myvar a b ascii;clear; pause;clear;load myvar; load myvar ascii;whoa,b5求下列表达式的值。（1）（2），其中a=3.5,b=5,c=-9.8。6欲通过系统做一平面图，请查找相关的命令与函数，获取函数的帮助信息。实验2 matlab绘图【实验目的】学会用MATLAB语言绘图【

14、实验要求】1、熟悉matlab绘制二维图形2、熟悉matlab绘制三维图形3、了解特殊图形的绘制4、熟悉图形的标识与修饰方法【实验原理】1、 plot 该命令在线性坐标系下绘制y对应于x的轨迹。若其中之一为矩阵则该命令将对应于矩阵的行或者列绘制一簇曲线。该命令也可以在同一坐标轴下绘制多条曲线。例1 运行程序 t=0:1:10; x1=t; plot(x1)可在图形窗口得到下面图形: 例2 x=0:0.5:7; y=sin(x); y1=2*x; y2=3*x; plot(y,x,y1,x,y2,x)从上面的例子可以看出, 如果要在同一坐标轴下绘制多条曲线, 可用如下格式： plot(x1,y1

15、,x2,y2.xn,yn) % 对应于不同矢量时 其中x1,y1.xn,yn是对应的关系。2、 log log 在两个对数（标准对数）坐标间绘制图形。 3、 semilogx 或semilogy 在x轴或y轴为对数坐标, 而另一轴为线性轴的空间绘制图形。4、polar 在极坐标空间绘制图形。5、 grid 在图形中加栅格。6、 subplot (r,c,p) 该命令将屏幕分成r*c个窗口,p表示在第几个窗口. 同时绘制多个图形窗口时,绘图语句间应用逗号分开。 例3 下面的程序产生四个窗口,且分别绘制相应图形. t=0:0.01:2*pi; y=sin(t); x=cos(t); subplot

16、(221),plot(t,y),grid,subplot(222),plot(t,x),grid,subplot(223),polar(t,y),subplot(224),polar(t,x)从上面的程序可以看出,subplot语句中的r,c,p可以不用逗号分开。 上面程序执行后得到如下图形：7、 gtext 该命令用于鼠标定位的文字注释. 执行完绘图命令后再执行gtext(说明文字) 命令, 就可在屏幕上得到一个光标, 然后,用鼠标选择说明文字的位置。例4 t=0:1:10; x1=2*t; x2=0.5*t; plot(x1,t,x2,t); gtext('x=2t');

17、gtext('x=0.5t') 8、 如果需要同时在不同窗口绘制图形, 可用figure语句来完成。 例5 t=0:1:10; x1=2*t; x2=3*t; figure,plot(x1),figure,plot(x2) 9、有关绘图方面的还有： ² Title 给当前图像加标题。² xlabel 或 ylabel 给x或y轴加标注。² text 在当前图像上指定位置加注释。10、在同一座标系中绘制多个曲线时, 曲线的颜色可以自动确定，也可以自己选择。其使用格式如下:plot(x1,y1,选项1,x2,y2,选项2,xn,yn,选项n)11、三维

18、图形绘制函数plot3函数与plot函数的用法类似。还可以用mesh函数绘制三维网格图；用surf函数绘制三维曲面图。MATLAB提供了丰富的三维绘图函数，可以通过help命令查看。 例6 绘制三维曲线：t=0:pi/100:20*pi;x=sin(t) y=cos(t);z=t.*sin(t).*cos(t);plot3(x,y,z) title(line in 3-D space) xlabel(X) ;ylabel(Y);zlabel(z);grid on运行结果如下：【实验内容】（一）验证例题中的操作。（二）完成以下各题，将答案写入实验报告，并在机器上运行

19、结果验证。1. 在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)， t0，2 2. 绘制0，4区间上的x=10sint曲线，并要求： （1）线形为点划线、颜色为红色、数据点标记为加号； （2）坐标轴控制：显示范围、刻度线、比例、网络线 （3）标注控制：坐标轴名称、标题、相应文本。 3. 绘制不同插值算法（nearest、linear、cubic、spline）对peaks函数的插值效果图，要求使用subplot命令在一个图形窗口显示四幅图。4. 绘制三维圆柱螺线 ： ，要求给相应的坐标轴和标题附加标注，螺线为兰色虚线；并使用stem函数绘制三维火柴杆图。5

20、.利用图形显示命令作出下列函数的图形：（1），其中（2），其中(3) ，其中(4) ，其中(5) 6.制作平面振动动画（利用函数）7.我们选用球体和圆柱体，其中实验3 matlab程序设计【实验目的】了解Matlab编程原理和编程方法【实验要求】1、熟悉Matlab for-end、while-end循环结构2、数学Matkab if-end、if-else-end选择判断结构【实验原理】MATLAB提供三种决策或控制流结构： for循环、while循环、if-else-end结构. 这些结构经常包含大量的MATLAB命令，故经常出现在MATLAB程序中，而不是直接加在MATLAB提示符下.

21、1、for循环：允许一组命令以固定的和预定的次数重复 for x=array commands end 在for和end语句之间的命令串commands按数组（array）中的每一列执行一次. 在每一次迭代中，x被指定为数组的下一列，即在第n次循环中，x=array(：，n) 例 对n=1,2,10,求xn= 的值for n=1:10 x(n)=sin(n*pi/10);endx2、While循环与for循环以固定次数求一组命令相反，while循环以不定的次数求一组语句的值. while expression commands end 只要在表达式(expression)里的所有元素为真，就执

22、行while和end语句之间的命令串commands. 例 设银行年利率为11.25%。将10000元钱存入银行，问多长时间会连本带利翻一番？money=10000years=0while money<20000 years=years+1 money=money*(1+11.25/100)end3、If-Else-End结构（1）有一个选择的一般形式是： if expression commands end 如果在表达式(expression)里的所有元素为真，就执行if和end语句之间的命令串commands. 先建立M文件fun1.m定义函数f（x），再在Matlab命令窗口输入f

23、un1(2),fun1(-1)即可。function f=fun1(x)if x>1 f=x2+1endif x<=1 f=2*xend【实验内容】对以下问题,编写M文件:(1)用起泡法对10个数由小到大排序. 即将相邻两个数比较,将小的调到前头.(2)有一个矩阵,编程求出其最大值及其所处的位置.(3)编程求(4)一球从100米高度自由落下,每次落地后反跳回原高度的一半,再落下. 求它在第10次落地时,共经过多少米?第10次反弹有多高?(5)有一函数 ,写一程序,输入自变量的值,输出函数值.实验4 matlab基础数学实验高等数学极限【实验目的】1、熟悉函数极限的概念2、掌握各种求

24、极限的方法3、会用Matlab命令求极限【实验要求】熟悉Matlab中求极限的命令limit【实验原理】例1：（1），（2）实验过程：（1）>> syms x a b>> limit(sin(a*x)/sin(b*x), x,0)运行结果：ans =a/b（2）>>syms x>> limit(1-cos(x)/(x*sin(x),x,0)运行结果：ans =1/2【实验内容】1、计算以下极限（1），（2），（3），（4）2、设计实验复利问题复利，即利滚利。不仅是一个经济问题，而且是一个古老又现代的经济社会问题。随着商品经济的发展，复利计算将日益

25、普遍，同时复利的期限将日益变短，即不仅用年息、月息，而且用旬息、日息、半日息表示利息率。现在我们已进入电子商务时代，允许储户随时存款或取款，如果一个储户连续不断存款和取款，结算本息的频率趋于无穷大，每次结算后将本息全部存入银行，这意味着银行不断地向储户支付利息，称为连续复利问题。 若银行一年活期年利率为0.06，那么储户存10万元的人民币，如果银行允许储户在一年内可任意次结算，在不计利息税的情况下，由于复利，显然这比一年结算一次要多，因为多次结算增加了复利。结算越频繁，获利越大。连续复利会造成总结算额无限增大吗？随着结算次数的无限增加，一年后该储户是否会成为百万富翁？实验5matlab基础数学

26、实验高等数学微分【实验目的】1.熟悉基本求导公式，掌握初等函数的求导方法.2.会求函数在给定点处的导数值.3.掌握隐函数求导的方法和步骤.4.掌握参量函数求一阶导数和二阶导数的方法和公式.【实验要求】1.熟悉，Matlab中的求导命令diff2.熟悉Matlab中解方程的命令solve和求导命令diff【实验原理】例1：求下列函数的导数（1），（2）例2.已知函数 ，求 ； 例3：设 ，求例4.已知 ，求实验步骤：1.（1）>> syms x>> y=exp(x)*(sin(x)+cos(x);>> diff(y)运行结果：ans =exp(x)*(sin(

27、x)+cos(x)+exp(x)*(cos(x)-sin(x) 即函数的导数为 （2）>> syms x>>y=log(x3+1)/(x2+1);>>diff(y)运行结果：ans =(3*x2/(x2+1)-2*(x3+1)/(x2+1)2*x)/(x3+1)*(x2+1)即函数的导数化简得 2、（1）>>syms x;>>f=1/x;>>f1=diff(f,x);>>ff=inline(f1);>>x=1;>>ff(1)运行结果：ans =-1>> x=-2;>&g

28、t; ff(-2)运行结果：ans = -0.25003、（1）解法一：>> syms x y;>> f=solve('x2+y2-R2=0',y);>> diff(f,x)运行结果：ans = -1/(-x2+R2)(1/2)*x 1/(-x2+R2)(1/2)*x即 或说明：对于能从方程中求出函数显示形式的题可以采用这种做法。解法二：>> syms x y R;>> f=x2+y2-R2;>> f1=diff(f,x);>> f2=diff(f,y);>> -f1/f2运行结果：

29、ans =-x/y 即 说明：对于不能从方程中解出函数显示形式的题要采用这种做法。4.(1)>> syms t;>> x=t2;>> y=4*t;>> f=diff(y,t);f1=diff(x,t);>> f2=f/f1运行结果：f2 =2/t 即 【实验内容】1、计算下列导数（1）（2）（3）已知函数，求，（4）设，求（5）已知，求2、设计实验最优价格问题某房地产公司拥有100套公寓当每套公寓的月租金为1000元时，公寓全部租出。当月租金每增加25元时，公寓就会少租出一套。1.请你为公司的月租金定价，使得公司的收益最大，并检验结论

30、2.若租出去的公寓每月每套平均花费20元维护费，又应该如何定价出租，才能使公司收益最大实验6 matlab基础数学实验高等数学积分【实验目的】1.掌握求函数的原函数的方法.2.熟悉基本积分公式和积分方法.3.掌握求函数定积分的方法.4.会求变上限函数的导数和带有变上限函数的极限.【实验要求】1.掌握Matlab中积分命令int.2.熟悉Matlab中求定积分的命令.【实验原理】例1、 求下列函数的一个原函数（1）（2）例2、 1.求下列定积分例3、 求变上限函数的导数例4、 求下列极限实验步骤：1.（1）>> syms x;>> f=1/x4;>> int(

31、f,x)运行结果：ans =-1/3/x3 即函数 的一个原函数为 .（2）>> syms x;>> f=exp(x)/(1+exp(x);>> int(f,x)运行结果：ans =log(1+exp(x) 即函数 的一个原函数为 .2、（1）>> syms x;>> f=sqrt(1-x2);>> int(f,x,0,1)运行结果：ans =1/4*pi 3、2.（1）>> syms t x;>> y=sin(t)/t;>> diff(int(y,t,0,x),x)运行结果：ans =

32、sin(x)/x 即 4、3.（1）>> syms x t;>> f=cos(t2);>> int(f,t,sin(x),0);>> f1=diff(int(f,t,sin(x),0),x)>> f2=f1/1>> limit(f2)运行结果：ans =-1【实验内容】1、求下列函数的一个原函数(1)、 (2) 2、求下列定积分（1）、（2）、3、求变上限函数的导数（1）、（2）、4、求极限5、现购买一栋别墅价值300万元，若首付50万元，以后分期付款，每年付款数目相同。10年付清，年利率为6%，按连续复利计算，问每年应付

33、款多少？（e -0.60.5448）实验7 matlab基础数学实验概率论与数理统计随机变量及其分布【实验目的】1会利用 Matlab 软件计算离散型随机变量的概率、连续型随机变量概率密度值, 以及产生离散型随机变量的概率分布(即分布律).2掌握常见分布的随机数产生的有关命令.【实验要求】1掌握常见分布的分布律和概率密度的产生命令，如binopdf，normpdf等2掌握常见分布的分布分布函数命令，如binocdf，normcdf等掌握常见分布的随机数产生命令，如binornd，normrnd等【实验原理】1.1 常见分布1. 事件 A 在每次试验中发生的概率是 0.3，计算（1）在 10 次

34、试验中 A 恰 好发生 6 次的概率；（2）生成事件A发生次数的概率分布；（3）在 10 次试验中 A 至少发生 6 次的概率；（4）设事件A发生次数为X，且X的分布函数为F(x)，求F(6.1)；又已知F（x）=0.345，求x。1（1）>> binopdf(6,10,0.3)运行结果为：ans = 0.0368（2）>>binopdf(0:10,10,0.3)运行结果为：ans = Columns 1 through 8 0.0282 0.1211 0.2335 0.2668 0.2001 0.1029 0.0368 0.0090 Columns 9 through

35、 11 0.0014 0.0001 0.0000 （3）>>binocdf(6,10,0.3)运行结果为：ans = 0.98944（4）>> binocdf(6.1,10,0.3)运行结果为：ans = 0.9894>> binoinv(0.345,10,0.3)运行结果为：ans = 2或者：>> icdf('bino',0.345,10,0.3)运行结果为：ans = 22 设随机变量 X服从参数是 3 的泊松分布， 求（1）概率 PX=6；（2）X的分布律前七项；（3）设X的分布函数为F(x)，求F(6.1)；又已知F（x

36、）=0.345，求x。（1）>> poisspdf(6,3)运行结果为：ans = 0.0504（2）>> poisspdf(0:6,3)运行结果为：ans = 0.0498 0.1494 0.2240 0.2240 0.1680 0.1008 0.0504（3）>> poisscdf(6.1,3)运行结果为：ans = 0.9665>> poissinv(0.345,3)运行结果为：ans = 2或者：>> icdf('poiss',0.345,3)运行结果为：ans = 23设随机变量 X服从区间2, 6上的均匀分

37、布，求（1）X=4 时的概率密度值； （2） PX5；（3）若PXx=0.345，求x。(1)>> unifpdf(4,2,6)运行结果为：ans = 0.2500 （2）>> unifcdf(5,2,6)运行结果为：ans =0.7500（3）>> unifinv(0.345,2,6)运行结果为：ans = 3.3800或者：>> icdf('unif',0.345,2,6)运行结果为：ans = 3.38004设随机变量 X 服从参数是 6 的指数分布， 求（1）X=0，1，2，3，4，5，6 时的概率密度值；（2）PX5；

38、（3）若PXx=0.345，求x。（1）>> exppdf(0:6,6)运行结果为：ans =0.1667 0.1411 0.1194 0.1011 0.0856 0.0724 0.0613(2) >> expcdf(5,6)运行结果为：ans = 0.5654（3）>> expinv(0.345,6)运行结果为：ans = 2.5387或者：>> icdf('exp',0.345,6)运行结果为：ans = 2.53875设随机变量 X 服从均值是 6，标准差是2 的正态分布，求：（1）X=3，4，5，6，7，8，9 时的概率密

39、度值；（2）X=3，4，5，6，7，8，9 时的分布函数值；（3）若PXx=0.345，求x；（4）求标准正态分布的上0.05分为点。（1）>> normpdf(3:9,6,2)运行结果为：ans =0.0648 0.1210 0.1760 0.1995 0.1760 0.1210 0.0648（2）>> normcdf(3:9,6,2)运行结果为：ans = 0.0668 0.1587 0.3085 0.5000 0.6915 0.8413 0.9332（3）>> norminv(0.345,6,2)运行结果为：ans = 5.2023或者：>>

40、; icdf('norm',0.345,6,2)运行结果为：ans = 5.2023（4）>> norminv(0.95,0,1)运行结果为：ans = 1.6449或者：>> icdf('norm',0.95,0,1)运行结果为：ans = 1.64491.2 随机数的生成1产生参数为 20, 概率为 0.25 的二项分布的随机数。(1) 产生 1 个随机数； (2) 产生 7个随机数； (3) 产生 21(要求 3 行7 列)个随机数。(1) >> binornd(20,0.25)运行结果为：ans = 4（2）>&

41、gt; binornd(20,0.25,1,7)运行结果为：ans = 4 7 3 7 7 5 7（3）>> binornd(20,0.25,3,7)运行结果为：ans = 4 9 2 0 6 5 7 5 6 3 4 4 5 2 2 5 2 0 5 6 5 或者：>> random('bino',20,0.25,3,7)运行结果为：ans = 2 9 5 1 3 4 4 3 4 2 5 3 7 1 5 4 9 6 7 6 42产生7个服从参数为6的泊松分布的随机数。>> poissrnd(6,1,7)运行结果为：ans = 3 10 5 9

42、3 4 6或者：>> random('poiss',6,1,7)运行结果为：ans = 6 8 8 6 7 13 73产生区间(1, 3)上的连续型均匀分布的随机数。 (1) 产生 6×6个随机数； (2) 产生 21(要求 3 行7列)个随机数。（1）>> unifrnd(1,3,6)运行结果为：ans = 2.0649 1.8077 1.4863 1.3516 2.2505 2.5248 2.4330 2.0971 1.3083 1.7207 2.0861 2.1521 1.3586 1.0975 2.9128 1.3776 1.8781

43、2.4953 1.6731 2.1055 2.8713 1.0024 1.5749 2.2911 1.3754 1.5496 2.6374 1.6328 2.0033 1.24641.6439 1.4830 2.4565 2.3992 2.5231 2.0088（2）>> unifrnd(1,3,3,7)运行结果为：ans = 1.6945 1.3963 2.3888 2.0646 2.8129 2.3429 1.1139 1.1843 2.3445 1.5136 1.5588 1.7854 2.6743 1.9006 1.2957 1.8630 1.0195 2.8925 1.0

44、497 2.9430 2.1649或者：>> random('unif',1,3,3,7)运行结果为：ans = 2.9571 1.9324 1.4606 2.1998 2.0276 1.9198 2.6108 2.6972 1.6513 2.1598 1.8969 1.8155 1.9018 2.4017 1.1013 2.2604 2.2063 1.0708 1.2161 2.1023 2.74454产生服从均值为4，均方差为1的正态分布的随机数。（1）产生7个随机数；(2) 产生 21(要求 3 行7列)个随机数。（1）>> normrnd(4,1

45、,1,7)运行结果为：ans =3.5674 2.3344 4.1253 4.2877 2.8535 5.1909 5.1892（2）>> normrnd(4,1,3,7)运行结果为：ans = 3.9624 3.8133 6.1832 5.0668 3.1677 4.7143 4.8580 4.3273 4.7258 3.8636 4.0593 4.2944 5.6236 5.2540 4.1746 3.4117 4.1139 3.9044 2.6638 3.3082 2.4063或者：>> random('norm',4,1,3,7)运行结果为：ans = 2.5590 4.6900 5.2902 2.7975 2.3959 5.4151 4.2193 4.5711 4.8156 4.6686 3.9802 4.2573 3.1949 3.0781 3.6001 4.7119 5.1908 3.8433 2.9435 4.5287 1.8293【实验内容】1、验证以上例题。2、请将教材上的部分题目用计算机实现。3、下图中每一个黑点表示钉在板上的一颗钉子，每排钉子等距排列，下一排的每个钉子恰在上一排两相邻钉子之间。假设有