MATLAB复习题.docx

MATLAB复习题第一章MATLAB概述 (常用标点符号的功能见P9) 2、用“format”命令设置数据输出形式，将pi显示为3.14159265358979，将pi显示 3.1416e+000。(format对于数值显示格式的控制见 P11) 3、历史命令窗有何功能？ (P12) 4、在命令窗中使用命令显示当前目录(cd)，并将当前目录设置为“A:exe”。(cd A:exe) 1、标点符号可使命令行不显示运算结果，用来表示该行是注释行。 é12ù5、输入变量a=5.3，b=êú，使用save命令将变量存入A:exe0101.mat文件。 34ëû>> a=5.3; >> b = 1 2;3 4; >> save A:exe0101.mat a b 6、将A:exe加入搜索路径中 addpath A:exe 7、clear,clc,clf的功能 clear 从内存中清除变量和函数； clc 清除命令窗口 clf 清除当前图形窗口 8、下面说法正确吗？MATLAB的数值表达精度与其指令窗中的数据显示精度相同；MATLAB指令窗中显示数值有效位不超过7位。 (P10) 数值的显示精度并不是数值的存储精度；有效位数可达16位。 第二章 数值计算 1、下列变量名中_A_是合法变量。 变量的命名规则P31 A. char_1,i,j B. x*y , a.1 C. xy, a1234 D. end,1bcx 2、已知x为1个向量，计算其正弦函数的运算为_C_。 A. SIN(X) B. SIN(x) C. sin(x) D. sinx 3、已知x为1个向量，计算ln(x)的运算为_B_。 1 A. ln(x) B. log(x) C. Ln(x) D. log10(x) 4、当a=2.4，使用取整函数得出3，则该取整函数名为_ceil_。 关于数组的基本函数P50 A. fix B. round C. ceil D. floor 5、已知a=0:4,b=1:5,下面的运算表达式出错的为_D_。 A. a+b B. a./b C. a*b D. a*b é123ùú，4566、输入矩阵a=ê“3”，使用单下标方式 a(6) êú使用全下标方式 a(1,3) 取出元素êë789úû取出元素“8”，用 a(2,3,:) 取出后两行子矩阵块。 7、输入A为3*3的魔方阵，B为3*3的单位阵，由小矩阵构成3*6的大矩阵C和6*3的大矩阵D，将D矩阵的最后一行构成小矩阵E。 A = magic(3);B = eye(3); C = A B;D = A;B;E = D(end,:) 特殊矩阵的产生 P34 8、求矩阵A的转置矩阵，逆矩阵，矩阵的秩，矩阵的行列式值，矩阵的三次幂，矩阵的特征值和特征向量。 A；inv(A)；rank(A)；det(A); A3;v,d = eig(A) 常用矩阵运算函数P44 9、求解方程组 ì2x1-3x2+x3+2x4=8ïx+3x+x=6ï124 íx-x+x+8x=734ï12ïî7x1+x2-2x3+2x4=5>>A = 2 -3 1 2;1 3 0 4;1 -1 1 8;7 1 -2 2; >>b = 8;6;7;5; >>x = Ab 矩阵运算分为左除和右除，各表示什么含义，见P47。 10、2个多项式为a(x)=5x4+4x3+3x2+2x+1，b(x)=3x2+1，计算c(x)=a(x)*b(x)，并计算c(x) 的根。当x=2 时，计算c(x)的值，将b(x)/a(x)进行部分分式展开。 2 >>a = 5 4 3 2 1;b = 3 0 1; >>c = conv(a,b) >>roots(c) >>polyval(c,2) >>r,p,k = residue(b,a) 关于多项式的常用函数参见P62-66 11、x在0,20范围内，计算多项式y=5x4+4x3+3x2+2x+1的值，并根据x和y进行二阶、三阶和四阶拟合。 >>x = 0:0.5:20; >>p = 5 4 3 2 1; >>y = polyval(p,x) >>p1 = polyfit(x,y,2) >.p2 = polyfit(x,y,3) >.p3 = polyfit(x,y,4) 第三章 符号计算 1、确定下面各符号表达式中的自由符号变量： 确定自由符号变量的法则P92 1/(sin(t)+cos(w*t) w 2*i+a*j j z*exp(j*theta) z 2、已知表达式：f=1-sin2x,g=2x+1，计算当x=1时f的值，计算f与g的复合函数，g的逆函数 >>f = sym(1-sin(x)2);g = sym(2*x+1) >>subs(f,1) >>compose(f,g) >>finverse(g) 3、计算1*2+2*3+n*(n+1)的前3项的和 符号级数P103 symsum(k*(k+1),1,3) 2*a+theta theta 4、符号函数：f=x(-y)，分别对x进行定积分,y进行不定积分，对y定积分区间为(0,1). 3 符号积分 P102 >>f = sym(x(-y); >>int(f) >>int(f,y,0,1) dy5、求+ytanx=cosx的通解。 符号常微分方程 P109 dx>>dsolve(Dy+y*tan(x)=cos(x),x) 6、利用符号绘图函数绘制图形f(x)=sinxx，范围为1,10。 符号函数绘图指令P110 >>f = sym(sin(x)/x); >>ezplot(f,1,10) 第四章 MATLAB计算的可视化和GUI设计 1、常用的二维、三维绘图函数 2、图形句柄的概念，查询和设置的方法 3、GUI设计过程，以及常用的控件名称和含义。 4、z=xe-x2-y2，当x和y的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图。 >> x,y=meshgrid(-2:.2:2); >> z=x.*exp(-x.2-y.2); >> mesh(x,y,z) >> subplot(2,2,1), plot3(x,y,z) >> title('plot3 (x,y,z)') >> subplot(2,2,2), mesh(x,y,z) >> title('mesh (x,y,z)') >> subplot(2,2,3), surf(x,y,z) >> title('surf (x,y,z)') >> subplot(2,2,4), surf(x,y,z), shading interp >> title('surf (x,y,z), shading interp') 4 5、建立三个输入窗口的输入对话框，如图所示。 prompt='输入姓名','输入年龄','输入职业' title='信息登记' lines=1 1 1' def='张平','30','教师' answer=inputdlg(prompt,title,lines,def); 第五章 MATLAB程序设计 1、命令文件和函数文件的区别。 2、输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出 c=input('请输入一个字符','s'); if c>='A' & c<='Z' disp(setstr(abs(c)+abs('a')-abs('A'); elseif c>='a'& c<='z' disp(setstr(abs(c)- abs('a')+abs('A'); 5 elseif c>='0'& c<='9' disp(abs(c)-abs('0'); else disp(c); end 3、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price 来表示)： price<200 没有折扣 200price<500 3%折扣 500price<1000 5%折扣 1000price<2500 8%折扣 2500price<5000 10%折扣 5000price 14%折扣 输入所售商品的价格，求其实际销售价格 price=input('请输入商品价格'); switch fix(price/100) case 0,1 %价格小于200 rate=0; case 2,3,4 %价格大于等于200 但小于500 rate=3/100; case num2cell(5:9) %价格大于等于500但小于1000 rate=5/100; case num2cell(10:24) %价格大于等于1000但小于2500 rate=8/100; case num2cell(25:49) %价格大于等于2500但小于5000 rate=10/100; otherwise %价格大于等于5000 rate=14/100; end price=price*(1-rate) %输出商品实际销售价格 4、一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数 for m=100:999 m1=fix(m/100); %求m 的百位数字 59 m2=rem(fix(m/10),10); %求m 的十位数字 m3=rem(m,10); %求m 的个位数字 6 if m=m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 5、从键盘输入若干个数，当输入0 时结束输入，求这些数的平均值和它们之和 sum=0; cnt=0; val=input('Enter a number (end in 0):'); while (val=0) sum=sum+val; cnt=cnt+1; val=input('Enter a number (end in 0):'); end if (cnt > 0) sum mean=sum/cnt end 6、求100，200之间第一个能被21 整除的整数 for n=100:200 if rem(n,21)=0 continue end break end n 7、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6 是完数。求1,500 之间的全部完数 for m=1:500 s=0; for k=1:m/2 if rem(m,k)=0 s=s+k; end 7 end if m=s disp(m); end end 8、利用函数的递归调用，求n！ function f=exam5_14(n) if n<=1 f=1; else f=exam5_14(n-1)*n; %递归调用求(n-1)! end 9、nargin 用法示例 函数文件examp.m： function fout=charray(a,b,c) if nargin=1 fout=a; elseif nargin=2 fout=a+b; elseif nargin=3 fout=(a*b*c)/2; end 命令文件exam5_15.m： x=1:3; y=1;2;3; examp(x) examp(x,y') examp(x,y,3) 1010、编写M脚本，分别使用for和while循环语句计算sum=åii的程序。 i=1sum = 0; for n = 1:10 sum = nn + sum; 8 end sum n = 1; sum = 0; while n <= 10 sum = nn + sum; n = n + 1; end 11、编制一个解数论问题的函数文件：取任意整数，若是偶数，则除以2，否则乘3加1，重复此过程，直到整数变为1。 function c = collatz(n) % Collatz % Classic "3n+1" Plobrem from number theory c = n; while n > 1 if rem(n,2) = 0 n = n/2; else n= 3 * n + 1; end c = c n; end 12、编写M脚本文件，要求从键盘逐个输入数值(input)，然后判断输入的数是大于0还是小于0，并输出提示(disp)是正数，还是负数，同时记录输入的正数、负数的个数。当输入0时，中止M文件的运行；当输入第10个数字时，显示记录的正、负个数并终止程序。 n = 1;n1 = 0; n2 = 0; while n <= 10 a = input('请输入数据'); if a > 0 disp('正数') n1 = n1 + 1; else if a < 0 disp('负数') n2 = n2 + 1; else 9 break; end end n = n + 1; end disp('正数') n1 disp('负数') n2 13、编写M函数文件，实现分段绘制曲面，绘制每个曲线为1个子函数。 ì0.5457e-0.75y-3.75x-1.5xï22z(x,y)=í0.7575e-y-6xï0.5457e-0.75y2-3.75x2+1.5xîfunction ex0306 x = -2:0.1:2; y = -2:0.1:2; xx,yy = meshgrid(x,y); n1,m1 = size(xx); z = zeros(n1,m1); for n = 1:n1 for m = 1:m1 if xx(n,m) + yy(n,m) > 1 z (n,m) = z1(xx(n,m),yy(n,m); 22x+y>1-1<x+y£1 x+y£-1 elseif (xx(n,m) + yy(n,m) > -1) & (xx(n,m) + yy(n,m) < 1) z (n,m) = z2(xx(n,m),yy(n,m); elseif xx(n,m) + yy(n,m) < -1 z (n,m) = z3(xx(n,m),yy(n,m); end end end surf(xx,yy,z) function z = z1(x,y) 10 z = 0.5457 * exp(-0.75 * y2 - 3.75 * x2 - 1.5 * x); function z = z2(x,y) z = 0.7575 * exp(-y2- 6 * x2); function z = z3(x,y) z = 0.5457 * exp(-0.75 * y2 - 3.75 * x2 + 1.5 * x); 14、根据输入参数个数实现当没有输入参数时，显示信息；当有一个参数时，则以该参数为边长绘制正方形；当有两个参数时，以2个参数为长和宽绘制矩形。 function ex0307(varargin) if nargin = 0 disp('无输入参数'); elseif nargin = 1 r = varargin1; x = 0 0 r r; y = 0 r r 0; plot(x,y) elseif nargin = 2 r1 = varargin1; r2 = varargin2; x = 0 0 r1 r1; y = 0 r2 r2 0; plot(x,y) end 11