计算机仿真复习题.docx

1 .数字1.5e2,1.5e3中的哪个与1500相同？数字相同？答案1.5e3答案,2 .请指出如下5个变量名中，哪些是合法的？答案xyz_3 ABCDefgh答案abcd-2 xyz_33chan a 变量 ABCDefgh3 .在Matlab中，比1大的最小数是多少？答案1+eps4 .指出指令clear, elf, clc各有什么作用？答案clear:清除Matlab工作空间中的变量Clf:清除图形窗口中的图形Clc:清除当前指令行窗口中的内容5 .求矩阵行列式的指令是答案det6 . Matlab绘图中，用于多子图的指令是？答案subplot7 .指定绘图使用菱形符蓝色虚划线的字符号是？答案'd-b'8 .视觉控制指令View(az,el)中，az表示？答案方位角9 .用户想计算y=2sin(0.3 n )/1+5的值，那么应依次输入：2*sin(0.3*pi)/(l+sqrt(5)的值,10 .创建符号表达式 f(x)= sinx + x 答案»f=sym('sin(x)+x')12 .计算以上表达式在x =0, x = n /4, x =2 Ji计算以上表达式在>>x(3)处的值。ans=0.2答案答案» subs(f,x,0)答案»x(l 479) ans =0 ans=00.30.60.8» subs(f,x,pi/4)»x(5:end) ans =1.4925 ans=0.40.50.60.70.80.9» subs(f,x,2*pi)»x(find(x.0.4) ans =6.2832 ans=0.50.60.70.80.91.032 f13 .设x为符号变量，(x)= x4+ x2+l, 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');»件g ans =xA4+5*xA2+9+xA3+5*x » f*g ans =(xA4+xA2+l)*(xA3+4*xA2+5*x+8)兀14 .求 cos(4 x +3)dx 和6 cos(4 x +3) dx 求0>> syms x » f = cos(4*x+3);» int(f) ans =l/4*sin(4*x+3)» int(f,0,pi/6) ans =l/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)f f I »plot3(y,z,x)»gridon »fitle(' el;x')»xlabel('sin(t)')»ylabel('cos(t)')»zlabel()»text(0,0,0,'origin') x x<l 16.有一个函数 y=2xA2-l1.5x<103x+2x>=10编写一个M文件输入x,输出y.function y=y(x) ifx<l Y=x ;else if x>=l&x<10 Y=2*xA2-l ; else Y=3*x=2 end17.对于如下简单系统进行Simulink仿真。其中u(t)为系统输入，(t)为系统输出，t为系统时间。y (t )=2u (t), t >25lOu (t ),t W 25 y很显然，此系统服从简单系统的条件，为一简单系统。系统输出仅由系统当前时刻的输入决定。一、填空题1、 Matlab的命令提示符是（»）<,命令行以（回车）结束。2、 MATLAB的一个命令行如有多条命令，命令间以（逗号）分隔。3、 MATLAB的命令行以分号（；）结尾表示（不希望显示命令运算结果）。4、MATLAB的M文件由两种类型分别是（命令文件或脚本文件）和（函数文件）。5、命令文件中的变量是（全局变量），在命令文件执行后（保存）（填保存/不保存）于工作空间;函数文件中的变量一般是局部变量，调用执行后（不保存）（填保存/不保存）于工作空间。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、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)；D、矩阵(1)； E、矩阵中的行元素之间(4)； F、矩阵中的列之间(5)(1)方括号口(2)单引号''(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*y,a.l C.xy,al234 D.end,lbcx 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、计算三个多项式sl、s2和s3的乘积，则算式为CA.conv(s I,s2,s3)B.s 1*s2*s3C.conv(conv(s I,s2),s3)D.conv(s 1*s2*s3)13、运行以下命令：»x=l 23;456;»y=x+x*i »plot(y)则在图形窗口绘制C条曲线。A.3B.2C.6D.414、subplot(2,l,l)是指 BA.两行一列的上图B.两行一列的下图C.两列一行的左图D.两列一行的右图15、极坐标图是使用B来绘制的。A.原点和半径B.相角和距离C.纵横坐标值D.实部和虚部16、if结构的开始是“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=2/3；(C) y2=xA2；(D) ans=1/32 .在 MATLAB 命令窗口中键入命令 A=l 42;312;615;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的特征值为：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系统运行时，内存中有包括X和Y在内的多个变量(数据)，要删除所有变量(数据)，应该使用的命令是(A)(A) clear；(B) clc；(C) home；(D) clearXY5 .用赋值语句给定x数据，计算(7sin(3+2x)+ e21n3严对应的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=4123131324;y=hist(data,4), hist是用来计算data数组的分布情况，也就是分别统计1,2,3.4每个数出现的次数.其中，data中有3个1,2个33个3,2个4.所以输出是3232.结果是(B)(A) y=4123；(B) y=3232；(C) y=1324；(D) y=42117 .在 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=meshgrid(-2:2)返回数据中：(D)(A) x是行向量，y是列向量；(B) x是列向量，y是行向量；(C) x是行元素相同的矩阵；(D) x是列元素相同矩阵9.下面有关MATLAB函数的说法，哪一个是错误的(D)(A)函数文件的第一行必须由function开始，并有返回参数，函数名和输入参数；(B) MATLAB的函数可以有多个返回参数和多个输入参数；(C)如果函数文件内有多个函数，则只有第一个函数可以供外部调用；(D)在函数中可以用nargin检测用户调用函数时的输出参数个数10.将带小数的实数处理为整数称为取整，常用四种取整法则是：向正无穷大方向取整、向负无穷大方向取整、向零方向取整和四舍五入取整。MATLAB提供了如下四个取整函数，若a =-1.4,对a取整的结果是-1,则不应该选用下面哪个函数。(A)(A) floor；(B) round；(C) ceil；(D) fix一、单项选择题I.在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=1/(2+cos(x);ezplot(diff(F),结果是(D)(A)绘出函数F在0,2n 的图形；(B)绘出函数F在-2 n ,2 n 的图形；(C)绘函数F的导函数在0,2 n 的图形；(D)绘函数F的导函数在-2",2n 的图形3 .在 MATLAB 命令窗口中键入命令 B=8,l,6;3,5,7;4,9,2； B*B(:,2)。结果是(B)(A) ans=(B) ans=(C) ans=(D) ans=9167676767916767676791674 . MATLAB命令x =3:2:100将创建等差数列，该数列是(D)(A)以3为初值的98个数,；(B)以100为终值的98的个数；(C)以99为终值的97个数；(D)以3为初值的49个数。5 . MATLAB 语句 strcat(int2str(2008)；年是s；年)的功能是(C)(A)将数据2008转换为符号；(B)将数据2008与符号变量合并；(C)将几个符号变量合并为一个；(D)将符号变量转换为数值变量；9 . MATLAB 的语句 colormap(001)(D)(A)将三维网面图确定为红色；(B)将三维网面图确定为绿色；(C)将三维网面图确定为蓝色；(D)语句使用格式错误10 .设a,b,c表示三角形的三条边，表达式a+b<c|a+c<b|b+c<a,(B)(A)是三条边构成三角形的条件：(B)是三条边不构成三角形的条件；（C）构成三角形时逻辑值为真；（D）不构成三角形时逻辑值为假三、简答题：I、下面的变量名那些合法那些不合法.为什么？a. dogl合法b. Idog不合法，应以字母开头c. Do_you_know_the_way_to_sanJose 合法d. help不合法，应以字母开头e. Whafs up?不合法，和？不是合法字符。2、看下面的数组回答有关问题1.10.02.1-3.56.0-0.01.1-6.62.83.4array-2.10.10.3-0.41.3-1.45.10.01.10.0a. array 1的大小是多少？5*4b. arrayl(4,l)的值是多少？-1.4-1.10.0-c.arrayl（:,l:2）的大小和值为多少?大小为4*2,值为；-1.45.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 )这是一个数组或矩阵加法:10a =(b) 这是一个数组乘法：L21.T 0101(c )这是一个矩阵乘法：一121f =-25_(d )这是一个矩阵乘法：°=38一 O(e )操作非法，两数组形不同(f)数组与标量的加法：a +d (g )数组乘法:50105(h )矩阵乘法:105-"J o4、检测下面的MATLAB语句，是对是错？正确的，输出结果如何，错误的，错在哪里?4.1 if volts >125disp('WARNING: High voltage on line.');if volts <105dispC'WARNING: Low voltage on line/);elsedisp(tLine voltage is within tolerances.*);end不正确，第二个if应该用elseif语句当电压大于125,警告电压过高；当电压小于105,警告电压过低，否则（即电压大于或等于105且小于或等于125时，显示电压在允许范围内）。4.2 . color ='yellow1;switch( color);case Ted',disp('Stop now!*);case 'yellow',disp(Trepare to stop.*);case 'green,dispCTroceed through intersection.*);otherwise,disp(,Illegal color encountered.*);end正确，这是十字路口交通信号灯的交通规则要求，根据灯的颜色显示停车、等待停车或通行命令4、下面的语句用来判断一个人的体温是否处于危险状态(温度用的是华氏计量)。这些语句是否正确？如果不正确，指出错在那里？应该如何改正？if temp <97.5dispCTemperature below normal*);elseif temp >97.5disp(*Temperature normal');elseif temp >99.5disp(,Temperature slightly high*);elseif temp >103.0dispCTemperature dangerously high*);end正常体温的数值范围不正确，原因是elsfif的设置条件顺序不正确，而且if条件不完备。可以改成以下形式：if temp <97.5disp(*Temperature below normal');elseif temp >103.0disp(*Temperature dangerously high*);elseif temp >99.5disp(*Temperature slightly high*);elsedisp('Temperature normal*);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函数返回调用这个函数时所需要的实际输入参数的个数error (nargchk(l,2,nargin);%nargout函数：返回调用这个函数时所需要的实际输出参数的个数； nargchk函数：如果一个函数被调用时参数过多或过少，那么nargchk 函数将返回一个标准错误信息if nargin =2res = sqrt(x .A2+ y .人2);elseres = x;%允许输入函数不等于2end7、简要说明函数文件与命令文件的区别(1)函数文件的第一行必须包含function字符，一般以return结束，并且第一个注释行专用于解释函数功能；命令文件无此要求。(2)函数文件的第一行必须指定函数名、输入参数及输出参数，命令文件无此要求(即脚本文件没有输入参数和输出参数)。(3)一个函数文件可以含0个、1个或多个输入参数和返回值。(4)函数文件要在文件的开头定义函数名，如function yl,y2=fiinc(x,a,b,c),则该函数文件名必须存为func.m,而命令文件无此要求。(4)命令文件的变量在文件执行结束以后仍然保存在内存中而不会丢失(视为全局变量)，而函数文件的变量仅在函数运行期间有效(除非用global把变量说明成全局变量，否则函数文件中的变量均为局部变量)，当函数运行完毕后，这些变量也就消失了。8、矩阵运算和矩阵元素运算有何不同？举例说明。矩阵运算满足线性代数的矩阵运算法则；矩阵元素运算又称为数组运算，一般是两个矢量矩阵之间的对应元素进行相应运算。9、ones和eye两函数功能分别是什么？Ones:构建全部矩阵元素均为1的矩阵。Eye：构建单位矩阵，即除主对角线为1夕卜，其他元素均为0的矩阵；r兀2乃1而而9、用冒号表达式建立一个含有一百个元素的数组angles=Angles =(.01:.01:l)*pi10、已知 A(:,:)=0.70940.65510.95970.75470.16260.34040.27600.11900.58530.67970.49840.2238,求下列命令的执行结果(1) »size(A)%获取数组尺寸ans =43(2) »numel(A)%获取数组元素的个数ans =1211、MATLAB中的向量和矩阵有什么不同(1)向量是一维数组，包括行向量和列向量；只有一个元素的向量就是标量；(2)矩阵包括线性代数中的二维矩阵和多维矩阵。12、有几种建立矩阵的方法？各有什么优点？1、直接输入法，例：a=379;49101,输入方法方便简捷2、通过M文件建立矩阵，适用于建立尺寸较大的矩阵，并且易于修改3、由函数建立，如y=sin(x),可以用函数内部建立一些特殊的矩阵4、通过数据文件建立，该方法可以调用其他软件产生数据13、在sin(x)运算中，x是角度还是弧度？在sin(x)运算中，x是弧度，MATLAB规定所有的三角函数运算都是按弧度进行运算。14、四、程序分析题1、 写出以下命令的执行结果1.1 函数印per和lower分别把一个字符串中所有转化大定和小写。» result = upper(*This is test 1!') result =THIS IS TEST 1!» result = lower(*This is test 2V)result =this is test 2!1.2 写出以下命令的执行结果»a =l 23;456 a =1 23456» b = mat2str(a)b =123;4563、下面的语句产生的result的值是多少？(a)x=12+ 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 >95A95grade >86B86grade >76C76grade >66D66grade >0F4.1if grade >95.0disp(*The grade is A.');elseif grade >86.0disp(*The grade is B.');elseif grade >76.0disp(*The grade is CJ);elseif grade >66.0disp(*The grade is D/); elsedispC'The grade is F.');endif grade >95.0disp(*The grade is A.*);elseif grade >86.0disp(fThe grade is B.1);elseif grade >76.0disp(*The grade is CJ);elseif grade >66.0disp(*The grade is DJ);elsedisp('The grade is F/);endendendend5、说明下列程序段功能：根据输入自变量x的值,计算分段函数y的值,满足：丫4伙，x<。lg(+vl + x2)/2,x>0x=input(，请输入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,100);x3=linspace(0,4*pi,100);yl=sin(xl);y2=l+sin(x2);y3=2+sin(x3);x=xl;x2;x31;y=yi;y2;y3；plot(x,y)在一个图形窗口绘制出三条二维曲线，y 1=sinx,y2=1+sinx,y3=2+sinx,横轴x取值范围都是0100,步长分别为2冗、3 n和4n。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 Space*);xlabel(,X,);ylabel(,Y,);zlabel(,Z,);grid on;在X、Y、Z坐标系绘制带栅格的三维曲线，标题为Line in 3D Space。6.3 分析下列程序的功能。x,y=meshgrid(0:0.25:4*pi);z=sin(x+sin(y)-x/l 0;mesh(x,y,z);三维网格图axis(04*pi 04*pi -2.51);绘制三维曲面图z=sin(x+sin(y)-x/10o7、创建符号变量有几种方法？MATLAB提供了两种创建符号变量和表达式的函数：sym和symsosym用于创建一个符号变量或表达式，用法如x=sym(，x，)及f=sym(tx+y+z,), syms用于创建多个符号变量，用法如syms xyz。f=sym('x+y+z')相当于syms x y zf= x+y+z五、编程题I、骰子模拟。模拟掷骰子的情况在现实中非常有用。编写一个maltab程序模拟掷骰子，每次产生一个1到6之间的随机整数。(可以编写命令或函数文件)。1.1 命令文件dice.m1.2 函数文件fdice.m2、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示):price<200没有折扣200price<5003%折扣500price<10005%折扣1000price<25008%折扣2500price<500010%折扣5000 price14%折扣编程输入所售商品的价格，并求其实际销售价格。程序如下：price=input(,W输入商品价格')；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 num2ceH(10:24)%价格大于等于1000%但小于2500rate=8/100;case num2cell(25:49)%价格大于等于2500%但小于5000rate=10/100;otherwise%价格大于等于5000rate=14/100;endprice=price*(l-rate)%输出商品实际销售价格3、符号运算题(每题5分，共15分)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(<sinx4-exp(2*x)>);3.2 Solve the following system of equations by symbolic operation:x +2 xy -y =32x2 x-3x +5=0Pl=,xA2+2*x*y-y=3,;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*31 A(1/2)-29/64+19/64*i*31 A(1/2)4、编写MATLAB脚本文件，计算Fibonacci数列中的第一个大于10000的元素及n值。Fibonacci数列从1,1开始，由数列的最后两个元素之和生成新的元素，依次递推。clear, close allN=inpu(输入最大数值范围N=，)等待用户输入数据HU；i=l；%变量的初始化while f+f(i+l)<N %循环条件检验Ri+2尸f(i+l)+f(i); i=i+l;%求 fibonnaci 数的算式end5、编写MATLAB程序，以0.2为步长绘制正弦曲线sin (n t )和余弦曲线cos( n t )在0,2冗之间的图形，分别位于图形窗口的二个子图中，正弦曲线以红色虚线绘制，并且用星号标注每个节点，余弦曲线以红色实线绘制，用圆标注每个节点。6、在0WxW2兀区间内，绘制曲线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.2e05Xcos(4冗x)和y2=2e°5Xco§(冗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=find(abs(y 1-y2)<1 e-2);%查找yl与y2相等点(近似相等)的下标xl=x(k);%取yl与y2相等点的x坐标y3=0.2*exp(-0.5*xl).*cos(4*pi*xl);%求 yl 与 y2值相等点的 y 坐标plot(x,y 1,x,y2；k:xl ,y3；bp*);9、在QWxW2兀区间内，绘制曲线yl=2e05x和y2=cos(4nx),并给图形添加以下标注。(1)标题为：x from 0 to 2 Jr ；(2)X 轴标签为：Variable X, Y 轴标签为：Variable Y：(3)在坐标(0.8,1,5)处添加标注曲线yl=2eq5x,在坐标(2.5,1.1)处添加标注曲线y2Kos(4nx)；(4)添加图例Yl、 Y2.程序如下：x=0:pi/100:2*pi;yl=2*exp(-0.5*x);y2=cos(4*pi*x);plot(x,yl,x,y2)%加图形标题%加X轴说明%加Y轴说明%在指定位置添加图形说明%加图例titleCx from0to2pi*);xlabel('Variable X');ylabelfVariable Y*);text(0.8,l.5；曲线 y 1=2eA-0.5xt);text(2.5,L1；曲线 y2=cos(4!pix),);legend(fcyr； y2')10、求下列联立方程的解3x +4y -7z 12w =45x-7y+4z +2vv =4x +8-5卬=9-6x +5y-2z + l Ovv =4» a=34-7-12;5-742;180-55-65-210;» y=4494;» x=inv(a)*y'2.94471.5512-0.60171.270911、（共15分,每小题5分）利用函数int计算:(1) f 2%dx J(x2+l)(x2+x)»int(1/(xA2+1)*(xA2+x)ans =-1/4* log(xA2+1)-1/2*atan(x)+log(x)-1/2*log(x+1)(2) Jj(x + y)exydxdy» syms x y» int(int(x+y)*exp(-x*y)ans =-l/y*(-l/y*(-l/y*(-exp(-x*y)*x*y-exp(-x*y)+exp(-x*y)/y)-exp(-x*y)» syms x y» f=(x+y)*exp(-x*y);» F=int(int(f,x),y)F = x*Ei(l,x*y)+x*(l/x/y*exp(-x*y)-Ei(l,x*y)+l/x*exp(-x*y)» simple(F) simplify:exp(-x*y)*(x+y)/x/y(3)2=向。+3求» syms x y» y=sin(x)A4+cos(x)A4;» subs(y,x,10)ans =0.583312、(共10分)编制m文件，等待键盘输入，输入密码555,密码正确，显示输入密码正确，程序结束；否则提示，重新输入。a=input('input password please')；while(a=555)a=input('wrong password, input again：')； end13、（共14分）编写m函数文件求底面半径为r、高为h的圆柱体的体积。function ttr=input （'输入底面半径'）;h=input（'输入高度'）；v=pi*r*r*h;disp （v）;14、（共14分）从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和立方和。sum=0;sum3=0;n=0;a=input（，请输入数，）;while （a=0）sum=sum+a;sum3=sum34-aA3；n=n+l;a=input（'输入数：'）；endif n=0aver=sum/n;enddisp('平均数为:');disp(aver);dispC立方和为:9);disp(sum3);15、设x=rcost+3t, y=rsint+3,分别令r=2,3,4,画出参数t=010区间生成的 xy曲线。> > t=linspace(0,10);> > rl=2;> > xl=(rl*cos(t)+3*t);> > yl=r