计算机仿真复习题.pdf

《计算机仿真复习题.pdf》由会员分享，可在线阅读，更多相关《计算机仿真复习题.pdf（41页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1.数字1.5e2,1.5e3中的哪个与1500相同？数字 相同？答 案1.5e3答案,2.请指出如下5个变量名中，哪些是合法的？答 案xyz_3 A B CD efgh答案abcd-2 xyz_3 3chan a 变量 ABCDefgh3.在M atlab中，比1大的最小数是多少？答案1+eps4.指出指令clear,elf,c lc各有什么作用？答 案clear:清 除M atlab工作空间中的变量Clf:清除图形窗口中的图形Clc:清除当前指令行窗口中的内容5.求矩阵行列式的指令是 答案det6.M atlab绘图中，用于多子图的指令是？答案subplot7.指定绘图使用菱形符蓝色虚划线

2、的字符号是？答案d-b 8.视觉控制指令View(az,el)中，a z表示？答案方位角9 .用户想计算y=2sin(0.3 n)/1+5的值，那么应依次输入：2*sin(0.3*pi)/(l+sqrt(5)的值，10.创建符号表达式 f(x)=sin x+x 答案 f=sym(sin(x)+x)12.计算以上表达式在x=0,x=/4,x=2”计算以上表达式在x(3)处 的 值。ans=0.2 答案 答案 subs(f,x,0)答案 x(l 4 7 9 )ans=0 ans=0 0.3 0.6 0.8 subs(f,x,pi/4)x(5:end)ans=1.49 25 ans=0.4 0.5

3、0.6 0.7 0.8 0.9 subs(f,x,2*pi)x(find(x.0.4)ans=6.28 32 ans=0.5 0.6 0.70.8 0.9 1.03 2 f13.设x为符号变量，(x)=x4+x2+1,g(x)=x+4x+5x+8,试进行如下运算。f(x)+g(x)f(x)X g(x)f=sym(xA4+xA2+1);g=sym(xA3+4*xA2+5*x+8);R g ans=xA4+5*xA2+9+xA3+5*x f*g ans=(xA4+xA2+1 )*(xA3+4*xA2+5*x+8)五14.求 cos(4 x+3)dx 和 6 cos(4 x+3)dx 求 0 sym

4、s x f=cos(4*x+3);int(f)ans=l/4*sin(4*x+3)int(f,0,pi/6)ans=1/8*3A(l/2)*cos(3)-3/8*sin(3)15.已 知 x 的 闭 区 间 0,10*pi 内,y=sin(x),z=cos(x),写 出 绘 制 如 图 示 的 三 维 螺 旋 线 的 指令 x=0:pi/50:1 O*pi y=sin(x);z=cos(x)ff1plot3(y,z,x)gridonfitle(c el;x)xlabel(sin(t)j y la b e l(cos(t)zlabel()text(0,0,0,origin)xxl16.有一个函数

5、 y=2xA2-l 1.5x=10编写一个M 文 件 输 入 x,输 出 y.function y=y(x)ifx=l&x2510u(t),t W 25 y很显然，此系统服从简单系统的条件，为一简单系统。系统输出仅由系统当前时刻的输入决定。一、填空题1、Matlab的命令提示符是（）o 命令行以（回车）结束。2、MATLAB的一个命令行如有多条命令，命令间以（逗号）分隔。3、MATLAB的命令行以分号（；）结尾表示（不希望显示命令运算结果）。4、MATLAB的M 文件由两种类型分别是（命令文件或脚本文件）和（函数文件）。5、命令文件中的变量是（全局变量），在命令文件执行后（保存）（填保存/不保

6、存）于工作空间；函数文件中的变量一般是局部变量，调用执行后（不保存）（填保存/不保存）于工作空间。6、函数文件保存时，文件名（必须是函数定义语句中的函数名）。7、（hold on）命令打开图形保持功能。8、（%）符号表示其后为程序注释。9、已知矩阵A,则矩阵A 的转置可表示为（AD,共枕转置可表示为（AD10、Matlab命 令 窗 口 中 的 符 号 表 示（）11、在 MATLAB中，图形的绘制必 须 在（图形）或（Figure）窗口中进行。二、选择题1、说明下列文件的后缀A、命令文件；B、Simulink仿真文件；C、数据文件；D、函数文件（1）.mdl（2）,m（3）.dat（4）2、

7、MATLAB是基于（函数）运行的软件，函数根据设计形式可以分为（内建函数或预定义函数）和（用户自定义函数）两类；根据功能又可以分为（通用函数）和专用函数。3、MATLAB的赋值符是（1）,相等符号是（2）(1)=(2)=(3)=4、A(:,j)表 示(1 )；A(i,:)表 示(3)；A(iJ)表 示(2)；A(i:i+m,:)表 示(4).(1)取 A 矩阵的第j 列全部元素；(2)取 A 矩阵第i 行、第j 列的元素。(3)A 矩阵第i 行的全部元素；(4)取 A 矩阵第ii+m行的全部元素；5、说明下列MATLAB数据对象的分隔符A、字符串(2)；B、向 量(3)C、函数的输出变量(1)

8、；D、矩 阵(1)；E、矩阵中的行元素之间(4)；F、矩阵中的列之间(5)(1)方括号口 单 引 号 (3)圆括号()(4)逗号，或空格(5)分号；6、命 令 subplot的用法和作用是(用来同时画出数个小图形于同一视窗中)，subplot(2,2,l)在一张图中同时显示(4)个子图。7、说明下列特殊变量的含义或取值A、inf()；B、pi()；C、NAN(不定项)；D、ans()；E、；F、()。8、在 MATLAB的命令窗口中执行(C)命令，将命令窗口的显示内容清空。A.clear B.clc C.echo off D.cd9、下列变量名中A 是合法的。A.char_l,i,j B.x*

9、y,a.l C.xy,al234 D.end,Ibex B。10、已知x 为一个向量，计 算 ln(x)的 MATLAB命令是A.A.ln(x)B.log(x)C.Ln(x)D.lglO(x)11、已 知 a=0:4,b=l:5,下面说法中错误的为DA.a+b B.a./b C.a*b D.a*b。12、计算三个多项式si,s 2 和 s 3 的乘积，则算式为CA.conv(s 1 ,s2,s3)B.s 1 *s2*s3 C.conv(conv(s I,s2),s3)D.conv(s 1 *s2*s3)13、运行以下命令：x=l 2 3;4 5 6;y=x+x*iplot(y)则在图形窗口绘制

10、C 条曲线。A.3B.2 C.6 D.414 subplot(2,l,l)是指 BA.两行-列的上图 B.两行一列的下图C.两列一行的左图D.两列一行的右图15、极坐标图是使用B 来绘制的。A.原点和半径 B.相角和距离 C.纵横坐标值 D.实部和虚部16、i f 结构的开始是“if”命令,结束是BA.End if B.end C.End一、单项选择题1.符号计算与一般数值计算有很大区别，它得到准确的符号表达式。在 MATLAB命 令 窗口中键入命令 syms x,yl=sqrt(x);y2=xA2；int(yl-y2,x,0,l),屏幕显示的结果是(D)(A)yl=xA(l/2)(B)ans

11、=2/3；(C)y2=xA2；(D)ans=1/32.在 MATLAB 命令窗口中键入命令 A=l 4 2;3 1 2;6 1 5;det(A(l:2,2:3).*A(l:2,2:3)o 结 果是(B)(A)ans=-143(B)ans=60(C)ans=-16(D)ans=-193,设 n 阶方阵A 的特征值为：入i(i=l,2,n),称 P(A)=max|入i|为矩阵Ai 的谱半径，则下列MATLAB求谱半径命令是(A)(A)max(abs(eig(A)；(B)abs(max(eig(A)；(C)max(norm(eig(A)；(D)norm(max(eig(A)4.MATLAB系统运行时

12、，内存中有包括X 和 Y 在内 的 多 个 变 量(数 据)，要删除所有 变量(数 据)，应该使用的命令是(A)(A)clear；(B)clc；(C)home；(D)clear X Y5.用赋值语句给定x 数据，计 算(7$皿3+2*)+62 1113)2对 应 的 MATLAB表达式是(A)(A)sqrt(7*sin(3+2*x)+exp(2)*log(3)(B)sqrt(7sin(3+2x)+exp(2)log(3)(C)sqr(7*sin(3+2*x)+eA2*log(3)(D)sqr(7sin(3+2x)+eA2 log(3)6.在 MATLAB命令窗口中输入命令data=4 1 2

13、3 1 3 1 3 2 4;y=hist(data,4)hist是用来计算data数组的分布情况，也就是分别统计1,2.3,4每个数出现的次数。其中，data中有3个1,2个2,3个3,2个4.所以输出是3 23 2.结果是(B)(A)y=4 1 2 3；(B)y=3 2 3 2；(C)y=1 3 24；(D)y=42 1 17.在 MATLAB 命令窗口中键入 A=magic(6);B=A(2:5,1:2:5)将得到矩阵 B,B 是(C)(A)2 行 5 列矩阵；(B)4 行两列矩阵；(C)4 行 3 列矩阵；(D)4 行 5 列矩阵8.MATLAB绘三维曲面需要构建网格数据，语句x,y=m

14、eshgrid(-2:2)返回数据中：(D)(A)x 是行向量，y 是列向量；(B)x 是列向量，y 是行向量；(C)x 是行元素相同的矩阵；(D)x 是列元素相同矩阵9.下面有关MATLAB函数的说法，哪一个是错误的(D)(A)函数文件的第一行必须由function开始，并有返回参数，函数名和输入参数；(B)MATLAB的函数可以有多个返回参数和多个输入参数；(C)如果函数文件内有多个函数，则只有第一个函数可以供外部调用；(D)在函数中可以用nargin检测用户调用函数时的输出参数个数1 0.将带小数的实数处理为整数称为取整，常用四种取整法则是：向正无穷大方向取整、向负无穷大方向取整、向零方

15、向取整和四舍五入取整。MATLAB提供了如下四个取整函数，若 a=-1.4,对 a 取整的结果是-1,则不应该选用下面哪个函数。(A)(A)floor；(B)round；(C)ceil；(D)fix一、单项选择题1.在 MATLAB命令窗口中，键入命令syms x；int(x*sin(x)结 果 是(A)(A)ans=sin(x)-x*cos(x)；(B)ans=cos(x)+x*sin(x)；(C)ans=sin(x)-cos(x)；(D)ans=-l/2*cos(x)*sin(x)+l/2*x2.在 MATLAB 命令窗口中，键入命令 syms x,F=l/(2+cos(x);ezplot

16、(difHF),结果是(D)(A)绘出函数F 在 0,2JI 的图形；(B)绘出函数F 在-2 冗，2k 的图形；(C)绘函数F 的导函数在 0,2 n 的图形；(D)绘函数F 的导函数在-2 n,2 n 的图形3.在 MATLAB 命令窗口中键入命令 B=8,l,6;3,5,7;4,9,2；B*B(:,2)。结果是(B)(A)ans=(B)ans=(C)ans=(D)ans=91 67 67 67 67 91 67 67 67 67 91 674.MATLAB命 令 x=3:2:100将创建等差数列，该数列是(D)(A)以 3 为初值的9 8 个数,；(B)以 100为终值的9 8 的个数；

17、(C)以 9 9 为终值的9 7 个数：(D)以 3 为初值的4 9 个数。5.MATLAB 语句 strcat(int2str(2008),年 是 s,年)的 功 能 是(C)(A)将数据2008转换为符号；(B)将数据2008与符号变量合并；(C)将几个符号变量合并为一个：(D)将符号变量转换为数值变量；9.MATLAB 的语句 colormap(0 0 1)(D)(A)将三维网面图确定为红色；(B)将三维网面图确定为绿色；(C)将三维网面图确定为蓝色；(D)语句使用格式错误1 0.设 a,b,c表示三角形的三条边，表达式a+bc|a+cb|b+c 2.1 0.1 0.3-0.4 1.3-

18、1.4 5.1 0.0 1.1 0.0a.array 1 的大小是多少？5*4b.arrayl(4,l)的值是多少？-1.4-1.1 0.0c.arrayl(:,l:2)的大小和值为多少?大小为4*2,值为-1.4 5.13、假设a,b,c和 d 的定义如下分别指出下列表达式的运算结果(a)a+b(b)a.*c(c)a*b(d)a*c(e)a+c(f)a+d(g)a.*d(h)a*d答案：(a)这是一个数组或矩阵加法：0 2-a+b=2 2_1 0a-(b)这是一个数组乘法：I?IJ 一 1 0-a.*b=0 1(c)这是一个矩阵乘法：一1 2a*b=-2 5_(d)这是一个矩阵乘法：3a*c

19、=一 8 O(e)操作非法，两数组形不同(f)数组与标量的加法：0 13d=52a+d=6756(g)数组乘法:a.*d=51005(h)矩阵乘法:4、检测下面的MATLAB语句，是对是错？正确的，输出结果如何，错误的，错在哪里？4.1 if volts 125disp(WARNING:High voltage on line.*);if volts 105disp(*WARNING:Low voltage on line/);elsedispfLine voltage is within tolerances.*);end不正确，第二个if应该用elseif语句当电压大于1 2 5,警告电压

20、过高；当电压小于1 0 5,警告电压过低，否则(即电压大于或等于105且小于或等于125时，显示电压在允许范围内)。4.2.color=yellow*;switch(color);case Ted,dispfStop now!1);case yellow；dispfPrepare to stop.*);case reen*,disp(fProceed through intersection.*);otherwise,dispClllegal color encountered.*);end正确，这是十字路口交通信号灯的交通规则要求，根据灯的颜色显示停车、等待停车或通行命令4、下面的语句用来判

21、断一个人的体温是否处于危险状态(温度用的是华氏计量)。这些语句是否正确？如果不正确，指出错在那里？应该如何改正？if temp 97.5dispCTemperature normal*);elseif temp 99.5dispCTemperature slightly high*);elseif temp 103.0disp(*Temperature dangerously high*);end正常体温的数值范围不正确，原因是elsflf的设置条件顺序不正确，而且i f条件不完备。可以改成以下形式：if temp 103.0disp(*Temperature dangerously high

22、*);elseif temp 99.5disp(*Temperature slightly high*);elsedisp(*Temperature normar);end6、请你确定下列两题中函数的调用是否正确。如果它是错误的，指出错误所在。6.1.不正确，test 1 函数有两个输入变量，调用时，必须给出两个实参out=test 1(6);function res=testl(x,y)res=sqrt(x.A2+y.A2);6.2.函数调用正确(该例说明：MATLAB函数允许改变输入输出参数的个数。)out=test2(12);function res=test2(x,y)%nargin函

23、数返回调用这个函数时所需要的实际输入参数的个数error(nargchk(l,2,nargin);%nargout函数：返回调用这个函数时所需要的实际输出参数的个数；nargchk函数：如果,个函数被调用时参数过多或过少，那 么 nargchk函数将返回一个标准错误信息if nargin=2res=sqrt(x.A2+y.A2);elseres=x;%允许输入函数不等于2end7、简要说明函数文件与命令文件的区别（1）函数文件的第一行必须包含function字符，一般以return结束，并且第一个注释行专用于解释函数功能；命令文件无此要求。（2）函数文件的第一行必须指定函数名、输入参数及输出参

24、数，命令文件无此要求（即脚本文件没有输入参数和输出参数）。（3）一个函数文件可以含0 个、1个或多个输入参数和返回值。（4）函数文件要在文件的开头定义函数名，如 function yl,y2=fimc（x,a,b,c）,则该函数文件名必须存为func.m,而命令文件无此要求。（4）命令文件的变量在文件执行结束以后仍然保存在内存中而不会丢失（视为全局变量），而函数文件的变量仅在函数运行期间有效（除非用global把变量说明成全局变量,否则函数文件中的变量均为局部变量），当函数运行完毕后，这些变量也就消失了。8、矩阵运算和矩阵元素运算有何不同？举例说明。矩阵运算满足线性代数的矩阵运算法则；矩阵元素

25、运算又称为数组运算，一般是两个矢量矩阵之间的对应元素进行相应运算。9、ones和 eye两函数功能分别是什么？Ones:构建全部矩阵元素均为1 的矩阵。Eye：构建单位矩阵，即除主对角线为1夕 卜，其他元素均为0 的矩阵；兀 2万 而丽,9、用冒号表达式建立一个含有一百个元素的数组angles=Angles=（.01:.01:l）*pi10、已知 A（:,:）=0.7094 0.6551 0.95970.7547 0.1626 0.34040.27 6 0 0.119 0 0.5 8 5 30.6 7 9 7 0.4 9 8 4 0.223 8,求下列命令的执行结果(1)size(A)%获取数

26、组尺寸ans=4 3(2)numel(A)%获取数组元素的个数ans=1211、MATLAB中的向量和矩阵有什么不同(1)向量是一维数组，包括行向量和列向量；只有一个元素的向量就是标量；(2)矩阵包括线性代数中的二维矩阵和多维矩阵。12、有几种建立矩阵的方法？各有什么优点？1、直接输入法，例：a=3 7 9;4 9 101,输入方法方便简捷2、通 过M文件建立矩阵，适用于建立尺寸较大的矩阵，并且易于修改3、由函数建立，如y=s i n(x),可以用函数内部建立一些特殊的矩阵4、通过数据文件建立，该方法可以调用其他软件产生数据13、在s i n(x)运算中，x是角度还是弧度？在s i n(x)运

27、算中，x是弧度，MATLAB规定所有的三角函数运算都是按弧度进行运算。14、四、程序分析题1、写出以下命令的执行结果1.1 函数up p e r和l o we r分别把一个字符串中所有转化大定和小写。result=upper(This is test 1!1)result=THIS IS TEST 1!result=lower(This is test 2!1)result=this is test 2!1.2 写出以下命令的执行结果 a=l 2 3;45 6a=1 2 34 5 6 b=mat2str(a)b=1 2 3;4 5 63、下面的语句产生的result的值是多少？(a)x=12+

28、i*5;y=5-i*13;result=x y;result=1,因为比较的是数的实部。(b)x=12+i*5;y=5-i*13;result=abs(x)abs(y);result=0,因为两个数的绝对值是确定的。(c)x=12+i*5;y=5-i*13;result=re al(x)-imag(y);result=254.说明下列程序段的功能：根据输入的数值分数，输出等级分数A 级、B 级和C 级，满足grade 95A95 2grade 86B86 2grade 76C76 2grade 66D66 2grade 0F4.1if grade 95.0disp(The grade is A

29、.);elseif grade 86.0disp(*The grade is BJ);elseif grade 76.0disp(*The grade is CJ);elseif grade 66.0disp(The grade is D.);elsedisp(*The grade is F.);end4.2if grade 95.0disp(*The grade is A.*);elseif grade 86.0disp(*The grade is B/);elseif grade 76.0disp(*The grade is C)；elseif grade 66.0disp(*The gr

30、ade is DJ);elsedisp(The grade is F.);endendendend5、说明下列程序段功能:根据输入自变量X的值,计算分段函数y的值,满足：y=(X +可 0 x=inputC请输入x的值:)；if x=0y=(x+sqrt(pi)/exp(2);elsey=log(x+sqrt(l+x*x)/2;end6、说明下列程序段功能6.1 分析下列程序绘制的曲线。x 1 =linspace(0,2*pi,l 00);x2=linspace(0,3*pi,l 00);x3=linspace(0,4*pi,100);yl=sin(xl);y2=l+sin(x2);y3=2+

31、sin(x3);x=xl;x2;x3;y=yl;y2;y3;plot(x,y)在一个图形窗口绘制出三条二维曲线，y 1 =sinx,y2=1 +sinx,y3=2+sinx,横轴x 取值范围都是0 7 0 0,步长分别为2 、3 兀和4 n。6.2 分析下列程序绘制的曲线。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 Space1);xlabel(X,);ylabel(Y,);zlabel(,Z,);grid on;在 X、Y、Z 坐标系绘制带栅格的三维曲线，

32、标题为Line in 3D Space。6.3 分析下列程序的功能。x,y=meshgrid(0:0.25:4*pi);z=sin(x+sin(y)-x/10;mesh(x,y,z);三维网格图axis(0 4*pi 0 4*pi-2.5 1);绘制三维曲面图z=sin(x+sin(y)-x/107、创建符号变量有几种方法？MATLAB提供了两种创建符号变量和表达式的函数：sym和 syms。sym用于创建一个符号变量或表达式，用法如x=sym(，x，)及 f=sym(1x+y+z,),syms用于创建多个符号变量，用法如syms x y z。f=sym(1x+y+z,)相当于syms x y

33、 zf=x+y+z五、编程题1、骰子模拟。模拟掷骰子的情况在现实中非常有用。编写一个maltab程序模拟掷骰子，每次产生一个1 到 6 之间的随机整数。(可以编写命令或函数文件)。1.1 命令文件dice.m1.2 函数文件fdice.m2、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示):price200没有折扣200price5003%折扣500price10005%折扣1000 price25008%折扣2500price500010%折扣5000price14%折扣编程输入所售商品的价格，并求其实际销售价格。程序如下：price=input(请输入商品价格)

34、；switch fix(price/100)case 0,1%价格小于 200rate=0;case 2,3,4%价格大于等于200但小于500rate=3/100;case num2cell(5:9)%价格大于等于500但%小于1000rate=5/100;case num2cell(10:24)%价格大于等于 1000%但小于2500rate=8/100;case num2cell(25:49)%价格大于等于 2500%但小于5000rate=10/100;otherwise%价格大于等于5000rate=14/100;endprice=price*(l-rate)%输 出商品实际销售价格

35、3、符号运 算 题(每 题 5 分，共 1 5 分)3.1 Calculate the third-order derivative(三阶导数)of following function by symbolic operation:sin x+e 2 xP=sinx+exp(2*x);DifT(p,x,3);ans=8*exp(2*x)或 p=sym(飞 inx+exp(2*x);3.2 Solve the following system of equations by symbolic operation:x 2+2 xy-y=32x2 x-3x+5=0Pl=,xA2+2*x*y y=3;

36、P2=,2*xA2-3*x+5=0,;x,y=solve(Pl,P2)x=3/4+l/4*i*31A(l/2)3/4-l/4*i*31A(l/2)y=-29/64-19/64*i*31A(l/2)-29/64+19/64*i*31A(l/2)4、编 写 MATLAB脚本文件，计 算 Fibonacci数列中的第一个 大 于 10000的元 素 及 n 值。Fibonacci数列从 1,1 开始，由数列的最后两个元素之和生成新的元素，依次递推。clear,close allN=inputC输入最大数值范围N=，)等待用户输入数据%变量的初始化while f(i)+f(i+l)N%循环条件检验好+

37、2尸f(i+l)+f(i);i=i+l;%求 fibonnaci 数的算式end5、编 写 MATLAB程序，以 0.2 为步长绘制正弦曲线sin(Ji t)和余弦曲线cos(兀 t)在 0,2冗 之间的图形，分别位于图形窗口的二个子图中，正弦曲线以红色虚线绘制，并且用星号标注每个节点，余弦曲线以红色实线绘制，用圆标注每个节点。6、在 0 x2兀区间内，绘制曲线y=2e-0.5xcos(4 n x)程序如下：x=0:pi/100:2*pi;y=2*exp(-0.5*x).*cos(4*pi*x);plot(x,y)7、采用图形保持，在同一坐标内绘制曲线yl=0.2e-5xcos(4 n x)和

38、 y2=2e5Xcos(冗 X)。程序如下：x=0:pi/100:2*pi;yl=0.2*exp(-0.5*x).*cos(4*pi*x);plot(x,yl)hold ony2=2*exp(-0.5*x).*cos(pi*x);plot(x,y2);hold off8、在同一坐标内，分别用不同线型和颜色绘制曲线yl=0.2e-0.5xcos(4兀 x)和 y2=2e-0.5xcos(n x),标记两曲线交叉点。程序如下：x=linspace(0,2*pi,1000);yl=0.2*exp(-0.5*x).*cos(4*pi*x);y2=2*exp(-0.5*x).*cos(pi*x);k=f

39、ind(abs(y 1 -y2)t=linspace(0,10);rl-2;xl=(rl*cos(t)+3*t);yl=rl*sin(t)+3;r2=3;x2=(r2*cos(t)+3*t);y2=r2*sin(t)+3;r3=4;x3=(r3*cos(t)+3*t);y3=r3*sin(t)+3;plot(xl,yl,r,x2,y2,b,x3,y3,m)16、设 f(x)=x-4x+3x2-2x+6(1)在 x=-2,8之间取100个点，画出曲线，看它有几个过零点。(提示:用 polyval函数)x=linspace(2,8,100);y=p o l yv a l(1 0-4 3-2 6,x

40、);(2)用r o o t s函数求此多项式的根。t=l 0-4 3-2 6p=r o o t s(t)17、绘制7=3-3+内在定义域x=-2,2,y=-2,2内的曲面。(利用me s h g r i d)x=-2:2;y=x;X,Y=me s h g r i d(x,y);Z=X 2*e x p(-(X 2+Y 2);me s h(X,Y,Z)4.求代数方程3x 5+4x+7x 3+2x+9x+12=0的所有根。(利用r o o t s函数)p=3 4 7 2 9 12;r o o t s(p)18、把1开五次方，并求其全部五个根。(提示：解(-1=0)p=l 0 0 0 0-1;r o

41、o t s (p)(-X2-4X-3)/2 -3 X-119、y=t _1X,画出在-3,3 上的曲线。(-/+4x-3)/2 1 A-aG，l)b=a(j,l);a(j,l)=a(j+l,l)；a(j+1,1)=b;endendenddispC排列后：)37、角度x=30 45 6 0 ,求X的正弦、余弦、正切和余切。x=30 45 60;xl=x/180*pi;sin(xl)ans=0.5000 0.7071 0.8660 cos(xl)ans=0.8660 0.7071 0.5000tan(xl)ans=0.5774 1.0000 1.7321 cot(xl)ans=1.7321 1.0

42、000 0.577438、用四舍五入的方法将数组2.4568 6.3982 3.9375 8.5042取整。b=2.4568 6.3982 3.9375 8.5042;round(b)ans=264939、基本编程题(每小题10分，共30分)1.Write a program to compute the following expressions10Z(4-耳)(%-弓)%=,，1 0；I 10 10 2J Z(X-吊 (%f)Y A=1 i=lHere,we suppose the variable x has existed in the workspace.for i=l:10for

43、j=l:10 xbar(i)=xbar(i)+x(i,j);endxbar(i)=xbar(i)/10;endfor i=l:10for j=l:10tl=0;t2=0;t3=0;for k=l:3tl=tl+(x(i,k)-xbar(i)*(x(j,k)-xbar(j);t2=t2+(x(i,k)-xbar(i)A2;t3=t3+(x(j,k)-xbar(j)A2;endr2(i,j)=tl/sqrt(t2*t3);endendr22.(1)Using plot()function to plot the curves ofy=sin(2x)+cos 3xe-2 J and z=sinxXi

44、n the range of XG-2,2,let their color are red and green respectively,and addthe grid to the figure.(2)Using fplot()function and ezplot()function to plot above-mentioned(上述的)curves respectively.(l)x=-2*pi:pi/100:2*pi;y=(sin(2*x)+cos(3*x).*exp(-2*x);z=sin(x)/x;plot(x,y,T,x,z,g)(2)fplot(sin(2*x)+cos(3*

45、x).*exp(-2*x),sin(x)/x,-2*pi 2*pi)3.Plot the 3D mesh figure and 3D surface figure of the function/(苍 )=9(17)%42-“2in the range of-44 and-4 4,respectively.x=-4:1/100:4;y4:1/100:4;z=9(1 -x)A2*exp(-x,A2/2-(y 9+l)A2)mesh(x,y,z);surf(x,y,z);得分四.综合编程题(每小题11分，共2 2分)1.Write a function program to compute the

46、 following functionJ|2 i|!,x=i(整数)凶 一14sin(2x)+”，其他and call this function program to compute the value of y=fi(-4)+fi(3)-fi(i4.37)inmain program or command window,there,i(|i shows the factorial(阶乘)of 小提示：对 x 取整的函数有 fix(x),floor(x),round(x)和 ceil(x).fuction y=fl(x)%fl.mif x-fix(x)=0y=prod(1 :abs(2*x)

47、;elsey=4*sin(2*x)+xA1.8；y=fl(-4)+fl(3)-fl(14.37)2.Write a program to compute the sample meanxbar=x=x,9n r=land the sample skewness(偏度)ske=-(x,-x)3.”,Tfor the following three case:(1)n and x=xt,x2,-,x exist in Workspace of Matlab.It is no any limit;(2)n and 工=入,刍,毛 are unknown variables.Please input

48、x by keyboard oneby one,and give a mark of end the data input.Examinees can use loop statementand length()function and input()function,but cannot use mean()function,std()function,and var()function directly(可以使用循环语句和 length。函数、input。函数等进行编程，但不能直接使用mean()、std()和 var()等函数).(l)s=0;for i=l:ns=s+x；endxbar=s/nm=0;for i=l:nm=m+(x(i)-xbar)A3;endske=m/n(2)x(1)=input(请输入 x);i=l;s=O;while(x(i)-a)s=s+x(i);i=i+l；x6=input(，请输入x，输a结束);endxbar=s/(i-l)m=0;for j=l:i-lm=m+(x(i)-xbar)A3;endske=m/(i-l)