2022年MATLAB实验题答案

1、MATLAB 试验题答 案 1,求以下变量的值 , 并在 MATLAB 中验 证； 1 a = 1 : 2 : 5 a = 31 22 22 40 49 13 -5 -8 7135 result6=a, *b result6 = 2 b = a , a , a ;a b = 28-3 111415 30 333 result7=a, /b result7 = 5551350,5000 0,5000 -3, 0000 3 c = a + b 2 , : c = 4,0000 1, 6667 1,2022 result8=a, c Attempt to reference field of non

2、-structure array, result9=a, b result9 = 4682,以下运算就是否合法 , 为什么？如合法 , 结果就是多少？ result2=a*b Error using *2,0000 2, 0000 -0, 3333 Inner matrix dimensions must agree , 0,2500 0,6000 0, 8333 result3=a+b result3 = result10=a, 2 result10 = 3621495811 16 25 36 result4=b*d result4 = result11=2, a result11 = 31

3、 22 22 24840 49 13 16 32 64 result5=b;c*d result5 = 3,用 MATLAB 求解下面的的方程 组； MATLAB 试验题答 案 7212x1 459, 2 92 10 y 10 2 n210 21 915 32x272211 5x3 1命令 n10 13213 x4 0, 求 y=？ 运行 format long g A=7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13 后 , 查瞧 y 的值 B=4 7 -1 0 m1=0; B=B for m=-10:10 x=invA*B m1=m1+2m; x y z 1en

4、d m1 2 x 2 y z w 82x y 3w 3m1 = 3x 3y 5z 6w 52047, 9990234375 A1=1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6 B2=1;8;3;5 6,求分段函数的值； x2=invA1*B2 y x 2 x 2 x 6x 0 x 5y 72125x 60 x 2x 1x 5A 915 322211 5用 if 语句实现 , 算出以下表中 x 对应的 4,已知 13213 值； 1 求矩阵 A 的秩 rank x=inputenter x=; 2 求矩阵 A 的行列式 determinant if x=0&x A3=7

5、 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13 else r=rankA3 y=x2-x-1; b=invA3 end a=detA3 y V,D=eigA3 第 2 页,共 10 页MATLAB 试验题答 案 7,分别用 if 与 switch 语句实现 , 将百分制 end case8 dispB case7 dispC case6 成果转换为成果等级 A, B, C, D, E； 其中 90100 分为 A,8089 分为 B,7079 分为 C,6069 分为 D,60 分以下为 E； 对超出百分制范畴的成果 , 给出错误提 示信息； if 结构程序 : x

6、=inputplease enter score=; if x=90&x=100 dispA dispD , elseif x=80 case0,1,2,3,4,5 dispB dispE elseif x=70 otherwise dispC disperror elseif x=60 end 8,摸索题 dispD elseif x=0 设计程序 , 完成成两位数的加,减,乘,除 四就运算 , 即产生两个两位随机整数 dispE 再输入一个运算符号 , 做相应的运算 , else 并显示相应的结果； disperror x=input 请输入运算符 end a=num2strfloorra

7、nd1*90+10; switch 结构程序 : ax=inputplease enter score=; b=num2strfloorrand1*90+10; switch fixx/10 case10,9 if x100 b if x=+ y=a+b; elseif x=- disperror else y=a-b; elseif x=* dispA 第 3 页,共 10 页MATLAB 试验题答 案 y=a*b; elseif x=/ y=a/b; else disperror end y x 2 61111, 1 2 22 3 2 n2 当 n 分别取 100,1000, 10000

8、时 , 求 x 的值分别就是多少？ a=input 请输入数值 n=0; for m=1:100 n=n+1/m*m; end n=6*n; x=sqrtn; 9,启动 MATLAB 后 , 点击 , 启MATLAB 的动 序编辑及调试器 程 Editor/Debugger, 编辑以下程序 , 点击储存程序 , 留意文 件名最好用英文字符；点击 Debug|Run 运行程序 , 在命令窗口查瞧运行结果 , 程序如有错误就改正； x 注 : 数论中一个好玩的题目 : 任意一个正整 m12k 10000 的 数 , 如为偶数 , 就用 2 除之 , 如为奇数 , 就与 3 相乘再加上 1；重复此过

9、程 , 最 11,编程求中意 k 终得到的结果为 1； 最小 mn=input 请输入 n 值 :; 值； sum=0; a=n; m=2; while n1 a=1; if remn,2=0 while a n=n/ 2; for i=1:m else sum=sum+2i; n=3*n+1; if sum10000 end a=0; a=a,n; end end end am=m+1; 10 , 根 据 end m第 4 页,共 10 页MATLAB 试验题答 案 12,摸索题 end 1, 0 已 知 y 与 t 的 函 数 关 14,编写一个函数求向量 x 中元素的平均 系 : 值,最

10、大值,最小值,均方根值； y t 1t 2 t / 2. 3 20 t / 3. .functtion / 20 . 求下面表格中与 t 对应的 y 值 m_x,max_x,min_x,rms_x=myfun2x t=input 请输入 t 值 : t 0, 2 0, 4 0,6 0, 8 sum=0; y for i=1:20 方均根值 Root Mean Square 的运算公式 b=factoriali; 为 : end sum=sum+ti/b; rms 1N2 xi Ni 1 sum=sum+1; 用下面数据测试您写的函数 : sum 1x=sin0:0 , 01:6*pi 13,编

11、写一个函数 , 运算下面函数的值 , 给出 标量 x 的值 , 调用该函数后 , 返回 y 的 2x=rand1,200, 得到的 x 为 200 个 0,1 值； 之间均匀分布的随机数； function y=myfun1x function m_x,max_x,min_x,rms_x=myfun2x %sin x, x 0求平均值 y x x, 0 x 3sum_x=sumx; %向量元素求与 x 6, x 3 m,n=sizex; %最好用 n=lengthx; m_x=sum_x/n; 选择一些数据测试您编写的函数； %求最大值 接受逐个比较方式 function y=myfun1x

12、if x1x2 if x0&x3 for k=3:n y=-x+6; 第 5 页,共 10 页MATLAB 试验题答 案 if max_xxk %可省略 min_x; max_x=xk; rms_x; %依据函数值行参次序输出结果 else 15 , 编 写 一 个 函 数 , 给 出 一 个 向 量 x x1, x2, xn , 生成如下范德 max_x=max_x; end 蒙矩阵； end function v=myvanderx %求最小值 111x1x2xnif x1xk min_x=xk; else n 1 x1 n 1 x2 n 1 xn 例如 :v=myvander2 3 4

13、5 得 v= 111123454916 25 827 64 125 function v=myvanderx min_x=min_x; %可省略 v1=vanderx; %生成范德蒙矩阵 end end %求均方根值 sum_x2=0; for k=1:n sum_x2=sum_x2+xk, 2; rms_x=sqrtsum_x2/n; end m_x; v2=v1; v=flipudv2; %实现矩阵上下翻转 16,摸索题 编写程序 , 用如下迭代公式求 a,a 的值 分别为 :3,17,113 ；迭代的终止条件为 xn 1xn 10 5, 迭 代 初 值 x0 , 迭代次数不超过 100

14、次； 分别对迭代结果与精确值进行比较 , 并 max_x; 第 6 页,共 10 页MATLAB 试验题答 案 统计迭代次数； x n1xn a t=0:0 ,1:4*pi; y1=t; 22xn y2=t0 , 5; 2 注 ：xa2 x2x 2 ax x a y2=sqrtt; xn x n 1 2x y3=4*pi*eps-0 2a2,1*t , *sint; 2 xn plott,y1,t,y2,t,y3 function x,n=sqrt_aa 18,编写程序 , 选择合适的步距 , 绘制下面函 x=1, 0; for k=1:100 m=x; x=x/2+a/2*x; if abs

15、x-m=10-5 break end end x; n=k; s=x-sqrta; if s=10-5 disp 正确 ; else disp 错误 ; end 17,在同一坐标系下绘制下面三个函数在 数在区间 -6,6 中的图象； sin x, x 0yx x, 0 x 3x 6, x 3 function y=fun1x y=; for x0=x if x00&x03 y=y,6-x0; end end x=-6:0, 1:6; y=fun1x; plotx,y 19,用 compass 函数画下面相量图 t 0,4 的图象； y 1t y 2t y 34e0 . 1 t sin t ua

16、 = 1; ub = 第 7 页,共 10 页MATLAB 试验题答 案 cos-2*pi/3+sin-2*pi/3*i ; 22,摸索题 uc=cos2*pi/3+sin2*pi/3*i; compassua,ub,uc,ua-ub,ub-uc,uc-ua ua=1; ob=cos-2*pi/ 3 ob=cos-2*pi/ 3+sin-2*pi/ 3*i; ub=cos-2*pi/ 3+sin-2*pi/ 3*i; 用对分法求解方程 2e x sin x 在 0,1 内的解 , 并验证 , 在程序中统计出对分 次数； 提 示 : 先 将 原 方 程 转 化 成 f x 2e x sin x

17、0的 形 式； uc=cos2*pi/ 3+sin2*pi/ 3*i; compassua,ub,uc,ua-ub,ub-uc,uc-ua 对 分 法 的 基 本 思 想 就 是 : 一 个 一 元 方 程 fx=0, 如 fx1*fx2 z=0:0, 1:4*pi; x=cosz; y=sinz; plot3x,y,z xm=x1+x2/2, 判定fx1 与 fx2 二 者 中 哪 一 个 与 fxm 异 号 , 如 fx1*fxm x0=0; x1=1; y=1; i=0; while y=0 i=i+1; x3=x0+x1/2; y=2*eps-x3-sinx3; y0=2*eps-x0

18、-sinx0; 21,用 mesh 或 surf 函数 , 绘制下面方程所 表示的三维空间曲面 ,x 与 y 的取值范 围设为 -3,3 ； x=-3:0,01:3; y=-3:0, 01:3; x,y=meshgridx,y; z=-x, *x/10+y , *y/10; y1=2*eps-x1-sinx1; meshx,y,z if y*y0 z=x, *x/10+y , *y/10; x0=x3; meshx,y,z elseif y*y1 plotx,y2 函数 : function y2=fzx y2=x,3-2*x ,2,*sinx+5*x ,*cosx+1, /x; 分析 : 依

19、据图像可知零点在 x=1 与 x=3 邻近 所以 End 23 , 求 代 数 方 程 fzerofz,1 ans = 1,5117 fzerofz,3 ans = 2,6095 25, 5 3x 4 4x 3 7 x 2 2 x 9x 12 0的 5 个根 , 并将其用星号 * 标记在复平面图上； 用 roots 与 plot 函数 a=3,4,7,2,9,12; y=rootsa ploty,* 求代数方程 5 x 1 0 的 5 个根 , 并将其用 星号 * 标记在复平面图上； 用 roots 与 plot 函数 程序 : a1=1,0,0,0,0,-1; y1=rootsa1 ploty1,* 已知 R=50 欧姆 ,U=4V, 二极管 D 正向电流与 电压的关系为 : 24,求下面函数在 0 ,5,4 区间内的过零点； I dI e s U d q 1 用 fzero 函 KT f x x3 2x 2 sin x 5x cosx