matlab基础练习题(带答案)

1、For personal use only in study andresearch; not for commercial useMatlab基础练习题常量、变量、表达式1、2、3、4、MATLAB中，下面哪些变量名是合法的？(A) _num (B) num_ ( C) num- 在MATLAB中，要给出一个复数 z的模，B) abs(z)(C)的预定义特殊变量的是？none ( C) zero (D) exp 的内存工作区中，存放一个英文字符 需要占用2个字节。(错，都是(A) mod(z)下面属于MATLAB (A) eps( B)判断：在MATLAB 一个中文字符啊)(D) -num

2、应该使用( double(z)函数。(D)angle(z)a2个字节需要占用1个字节,)存放5、判断：MATLAB中，i和j都是虚数单位匚1，它们之间没有什么区别。( 对6、判断：MATLAB中，pi代表圆周率，它等于 3.14。(错，后面还有很多位小数7、在 MATLAB 中,若想计算的y1 = 2测(0严)值，那么应该在 MATLAB的指令窗中1+V59、输入的 MATLAB 在 MA TLAB 节。在 MA TLAB中,指令是仁2*sin(0.3*pi(1+sqrt(5) _。a = 1，b = i，贝U a占_8_个字节，b占_16_个字节，c占中,inf的含义是_无穷大_, nan的

3、含义是_非数(结果不定)数组1、2、X是一个一维数值数组, 的逆序排列输出，应该使用下面的(A) Xend:1( B) Xend:-1:1在MA TLAB中，A是一个字二维数组， 的命令是()。(A) class(A) ( B) sizeof(A) 在MATLAB 中，用指令 x=1:9 赋值为0,应该在指令窗中输入(A) x(2 7)=(0 0)(C) x(2,7)=0 0在 MATLAB中,现在要把数组X中的所有元素按原来次序)指令。(C) X (end:-1:1)( D) X(end:1)要获取A的行数和列数，应该使用的MATLAB(C) size(A)( D) isa(A)生成数组x。

4、现在要把x数组的第二和第七个元素都 )(B) x(2,7)=0,0(D) x(2 7)=(0 0)3、4、5、6、7、89、(A)b = -3-2 -1(B)b :=-2 -1 0 1(C)b = -5 -1 3(D)b :=-5 -2 1 4在 MA TLAB中，A =1:9，现在执行如下指令L1 = (A5)，则MATLAB的执行结果应该是L1 =111110 00 0 。(2)A.*B(3)A(:)(4)A(1,:)*B(:,3)(1) A=B11 11(2)A(:) = 1:2:24A 3)A(1:3:7)17211923在 MATLAB 中，依次执行以下指令

5、： clear;A=ones(3,4);A(:)=-6:5; 这时,若在指令窗中输入指令 b=A(:,2)，那么，MATLAB输出的结果应该是()在MATLAB中，要求在闭区间0,5上产生50个等距采样的一维数组b，请写出具体的 MATLAB 指令 _linspace(0,5,50) _。在 MATLAB 中，A = 0:1/2:2*pi，那么 sin(A) = _ 0 1 0 -1 0。在 MA TLAB 中，A=1,2,3;4,5,6;7,8,0，B=2,1,6;8,5,2;14,2,1。写出下面 MATLAB 语句执行的结果：(为节省篇幅，把矩阵写成mat2str的形式)0 0 0; 0

6、 1 0; 0 0 0 2 2 18;32 25 12;98 16 01 4 7 2 5 8 3 6 013在MATLAB中，写出下面 MATLAB语句执行的结果:(1) clear,A = on es(2,6)A =11111ans =1713(4)diag(diag(A)ans =1007(5) B = A(:,end:-1:1)B =21171392319151110、请编写一段matlab程序，完成以下功能：(1) 生成一个100行，200列的二维随机数组；A=ra nd(100,200);(2) 找出数组A中所有大于0.49且小于0.51的元素的单下标；ldx=find(A(:)0.

7、49 & A(:)0.51)(3) 数组A中满足(2)中的条件的元素有多少个？len gth(ldx)(4) 求出数组A中满足(2)中的条件的元素的和，并求出这些元素的平均值;sum(A(Idx), m=mea n( A(ldx)(5) 将(4)求出的平均值赋值给数组A中满足(1)中的条件的每个元素。A(ldx)=m;11、请编写一段matlab程序，完成以下功能：(1) 找出100到200之间的所有质数，将这些质数存放在一个行数组里;X=100:200;p=X(isprime(X)(2) 求出这些质数之和；sum(p)(3) 求出100到200之间的所有非质数之和(包括 100和200)。s

8、um(X(isprime(X)12、 y = 0.7 - 2cos2x sinx，编写一段 matlab 程序，要求如下(1+x)(1) 在0,2二区间，每隔0.01取一 x数值，计算出相应的 y的函数值;x=0:0.01:2*pi;y=(0.7+2*cos(x)./(1+x.A2).*s in(x)(2) 根据MATLAB计算出的数据，找出在0, 2二内该函数的极小值的坐标。my,idx=mi n(y) x(idx)数值计算1、在MATLAB中，A是一个10X 10数组，我们把该数组看成矩阵的话，则此矩阵的行列式值=_det(A)_,此矩阵的逆矩阵(如果存在的话) =_inv(A) _。(用

9、MATLAB 的函数表示)2、求解下列方程组：，3x + 4y-7z-12w = 45x - 7y + 4z + 2w = -3x 8z -5w = 9-6x 5y - 2z 10w = -8玄 + x2 + 3x3 _ x4 = 2 严 _X3 +X4 =1x1 x2 2x3 2x4 二 4x1 - x2 x3 - x4 = 0A=1 1 3 -1;0 1 -1 1;1 1 2 2;1 -1 1 -1; b=2;1;4;0;x=AbA=3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10; b=4;-3;9;-8;xyzw=AbC 33 c L(3x +x 丫 x

10、+0.5 3、 求有理分式 R =32的商多项式和余多项式(x +2x -2 丫5x +2x +1 )n=co nv(3 0 1 0,3 0 0 0.5);d=co nv(1 2 -2,5 2 0 1);q,r=dec onv(n,d)424、 一元多项式 p=2x -3x +4x，写出表示 p 的 MATLAB 语句 _p=2 0 -3 4 0_，求p=0的根的MATLAB语句是roots(p)_,求X = 43时p的数值的 MATLAB语句是 _ polyval(p,4.3)。程序设计1、M脚本文件和函数文件的有什么区别？最主要的差别是，函数有自己的独立工作区，和外界的信息交换通过输入输出

11、参数实现; 而脚本没有独立工作区，其直接使用宿主程序的工作区，并把结果留在宿主程序中。2、掌握下面程序设计中常用的函数及关键字：in putpausewarni ngnarginn argoutforwhileif elseifelse switch caseotherwiseendtrycatchbreakcon ti nueend633、分别用for和while循环语句计算 K = v 2i的程序。还请写出一种避免循环的计算程i =9序。(不要用 symsum函数)K=0;for i=0:63K=K+i;endK=0; i=0;while i= 10000，所以n = ceil(log(1

12、0000/0.06e-3)/log(2)10、编写MATLAB脚本输出“水仙花数”及水仙花数的个数。所谓“水仙花数”是一个3位数，其各位数字的立方和等于该数本身。例如15 13 53 33。c=0;for n=100:999m=fix(n/100) fix(mod(n,100)/10) mod(n,10);if n=sum(m.A3),fprin tf(%i = %iA3 + %jA3 + %iA3n, n,m)；c=c+1;endendfprintf(n 共 %i 个水仙花数 n,c)输出结果：153 = 1A3 + 5A3 + 3A3370 = 3A3 + 7A3 + 0A3371 = 3

13、A3 + 7A3 + 1A3407 = 4人3 + 0A3 + 7A3共4个水仙花数绘图1、MATLAB中，绘制三维曲面图的函数是()(A)surf( B) plot( C) subplot(D)plot32、MATLAB中，要绘制三维空间曲线，应该使用()函数。(A) polar( B) plot( C)subplot(D)plot33、在MATLAB中，A是一个1000行2列的二维数值数组，现在要把A的第一列数据作为横坐标，把 A的第二列数据作为纵坐标，画出一条曲线，试写出相应的MATLAB语句plot(A(:,1),A(:,2) 。4、MATLAB绘图指令中的_subplot_指令允许用

14、户在同一个图形窗里布置几个独立的子图。5、 编写一-段matlab程序，绘制出二元函数z二 一iy二维网线图，要求如下xy(1)x，y的取值范围为-9乞x乞9，-9 syms x limit(1+xF(1/x),0)ans =expsin 3x(2) 求函数y在X =0处的极限;tg5x syms x limit(si n(3*x)/ta n(5*x),0)ans =3/5(3)nx求函数y-在趋向正无穷处的极限;3x syms x n limit (n *x/3Ax,i nf)ans =0求函数ln2xy -x在趋向正无穷处的极限; syms x limit(log(x)A2/xA3, in

15、f)ans =07、求导数：1(1) 求函数y 2 _3x 3的50阶导数；x syms x y=1/xA2-3*x+3; diff(y,50)ans =1551118753287382280224243016469303211063259720016986112000000000000/xA52cta(2) 求函数y二as in be t 在t = b处的3阶导数；syms t a b cy=a*si n(b*exp(cAt)+tAa);simple(subs(diff(y,t,3),t,b)8、求不定积分：Ex2 dxa -xdx syms x in t(1/si n(xF3)ans =-

16、1/2/si n(x)A2*cos(x)+1/2*log(csc(x)-cot(x) syms x a in t(1/(aA2-xA2)ans =-1/2/a*log(a-x)+1/2/a*log(a+x) syms x in t(sqrt(xA2-3)-sqrt(xA2+3)/sqrt(xA4-9)ans =(xA4-9)A(1/2)/(xA2-3)A(1/2)/(xA2+3)A(1/2)*asi nh(1/3*3A(1/2)*x)-1/(xA2+3)A(1/2)*(xA4-9)A(1/2)/(xA2-3)A(1/2)*log(x+(xA2-3)A(1/2)9、求定积分及广义积分f Jx2

17、+adx-2sin2xcos2 xdx-312 x2221 十 x y dydx1 syms x a in t(sqrt(xA2+a),-2,2)ans =2*(4+a)A(1/2)+1/2*a*log(2+(4+a)A(1/2)-1/2*a*log(-2+(4+a)A(1/2) syms x in t(si n(x)A2*cos(x)A2,-pi,pi)ans =1/4*pi syms x y in t(i nt(xA2+yA2,y,1,xA2),1,2)ans =1006/10510、求下面的积分，给出50位精度的数值:sin2 x sin2 y dydx syms x y J=i nt(

18、i nt(si n(x)A2+si n(y)A2,y,1,xA2),1,2); vpa(J,50)ans =11、级数求和:n +迟 n()2.1540459589705316265997501755762001048498664176916_n、3n 1 z -1n 42k 1 壬 2 i x -1 z.7 2k 1 x 1 syms z n symsum(z-1)A n/(n A2*2A n),n ,1,i nf)ans =(1/2*z-1/2)*hypergeom(1, 1, 1,2, 2,1/2*z-1/2) syms z n symsum(3* n+1)*(z-1)A n,n ,1,

19、i nf)ans =(4*z-4)*(-1/(z-2)+3/4/(z-2)A2*(z-1) syms z n symsum( n*(_1)A( n+1)*zA n,n ,1,i nf)ans =z/(z+1)A2 syms x positive syms k simple(symsum(2/(2*k+1)*(x-1)/(x+1)A(2*k+1),k,0,i nf)ans =log(-(1+(xA2-2*x+1)/(xA2+2*x+1)A(1/2)/(-1+(xA2-2*x+1)/(xA2+2*x+1)A(1/2)仅供个人用于学习、研究；不得用于商业用途For personal use only in study and research; not for commercial use.Nur f u r den pers?nlichen f u r Studien, Forschung, zu kommerziellen Zwecken ve