厦门理工学院试卷matlab程序设计题目及答案
Matlab编程与应用习题和一些参考答案
Matlab编程与应用习题和一些参考答案Matlab 上机实验一、二3.求下列联立方程的解⎪⎪⎩⎪⎪⎨⎧=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\b4．设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。
>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];>> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];>> C1=A*B'>> C2=A'*B>> C3=A.*B>> inv(C1)>> inv(C2)>> inv(C3)5．设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。
>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。
并求该矩阵全体数的平均值和均方差。
（mean var ）a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)10.利用帮助查找limit 函数的用法，并自己编写，验证几个函数极限的例子。
大学matlab考试题及答案
大学matlab考试题及答案大学MATLAB考试题及答案一、选择题1. MATLAB的全称是什么？A. Matrix LaboratoryB. Microprocessor Application Tool SetC. Microsoft Advanced Tool SetD. Microprocessor Application Test System答案：A2. 在MATLAB中，以下哪个命令用于绘制三维图形？A. plotB. scatterC. surfD. bar答案：C3. MATLAB中用于求解线性方程组的函数是？A. solveB. linsolveC. linprogD. fsolve答案：A二、简答题1. 简述MATLAB的基本数据类型有哪些，并给出至少两个每种类型的示例。
答案：MATLAB的基本数据类型包括数值数组、字符数组和单元数组。
数值数组可以是向量、矩阵或多维数组。
例如，向量 `v = [1 2 3]`，矩阵 `M = [1 2; 3 4]`。
字符数组是由单引号或双引号括起来的字符序列，如 `C = 'Hello'`。
单元数组可以包含不同类型的数据，如`{1, 'text', [1; 2; 3]}`。
2. 描述如何在MATLAB中实现矩阵的转置和翻转。
答案：矩阵的转置可以通过 `'T'` 来实现，例如 `A'` 表示矩阵A 的转置。
矩阵的翻转可以通过 `flip` 函数实现，例如 `flip(A)` 可以翻转矩阵A的所有行和列，`flipud(A)` 仅翻转矩阵A的行，而`fliplr(A)` 仅翻转矩阵A的列。
三、编程题1. 编写一个MATLAB函数，该函数接受一个向量作为输入，并返回向量中所有元素的和以及平均值。
```matlabfunction [sumVal, avgVal] = calculateSumAndAverage(V)sumVal = sum(V);avgVal = mean(V);end```2. 设计一个MATLAB脚本来解决以下问题：给定一个3x3的矩阵，找出其中最大的元素，并显示其位置。
matlab期末考试试卷及参考答案
一、填空题(本题满分30分，每空3分)1.设有程序A=[2,0,1,0;1,0,1,5;0, 6 ,2, 3];B=[1,2,-1;1,3,1];C=A(1:3,[1,2]);D=size(C)- length(B(2,:));F= C.*B ’将下列命令的运行结果填在横线上0*、0D = [0 -1] ; F= [2 0;2 0;0 6] .2. A=[2,3,2,1;1,2,2,1; 1,2,3,-2];A1=sum(A>A(2,1));A(2,:)=[],A2=A.^2-4A1= ; A2= ;3.P=[1,2,5,4;2,2,1,3];max(P,[ ],1)= [2 2 5 4] ;mean(P,2)= [3 2] .二、(本题满分10分)试编写计算程序画出曲线与曲面图形.22,22,)1((2).);20(cos sin (1).223322≤≤-≤≤-++=≤≤⎪⎩⎪⎨⎧==--y x y x e x z t t y t x y xπ%1 %2t=0:0.1: 2*pi; [x,y]=meshgrid(-2:0.1:2);x= sin(t).^3; z=(x.^2+1).* exp(-x.^2- y.^2+)+x .^2.*y; y= cos (t).^3; mesh(x,y,z)plot(x,y)三、（本题满分12分）编写程序解决下列问题：（1）试产生100⨯4阶矩阵使其元素在0—100之间的随机整数,可将此矩阵理解为100名同学的四门课程(课程1, 课程2, 课程3, 课程4)的成绩.（2）计算出每门课程和每位同学的平均成绩;（3）统计出总的优秀人次(90分及以上为优秀);（4）统计出四门课程成绩至少有两门课程不及格的人数.解：(1) a=fix((100-0+1)*rand(100,4));(2) mean(a) mean(a,2)(3)sum(sum(a>=90,2)==4)(4)sum(sum(a<60,2)>=2)四、（本题满10分）).50,20(]2,0[)2()),(1:(;)1(,sin )1()(11==-=∑=+n n x fun y function kx k x f nk k 分别取内图形用子图画出该函数在程序的第一句为返回该函数并保存试编写函数式文件程序设π(1). function y=fun1(x,n) (2). x=0:0.1:2*pi;y=0; y1= fun1 (x,20);for k=1:n y2= fun1 (x,50);y=y+ ((-1)^(k+1)*sin(k*x))/k; subplot(2,1,1),plot(x,y1)end subplot(2,1,2),plot(x,y2)y五、(本题满分10分)试利用微分方程的数值解法写出下列微分方程的求解步骤和MATLAB 程序.0)0(,2)0(,3='==-'+''x x t x x t x .(t 的变化区间为[0,10]).第一步：先将高阶微分方程转化为一阶微分方程组：选择状态变量x x x x ='=21,，则原方程化为2)0(,0)0(.,32112211==⎩⎨⎧='++-='x x x x t x tx x第二步：建立函数文件function dxdt = fun2(t,x)dxdt = [-3*t*x(1)+ x(2)+t; x(1)];第三步：求解微分方程，命令如下：[t,x]=ode45(@fun2,[0,10],[0;2]);)(cos cos lim .3;,)(,,1ln sin .2;2.1:,sin 022x e x x x A dx x A dx dA x s e x x x A x e 、x x x x --⎥⎥⎦⎤⎢⎢⎣⎡+=+=→-⎰计算极限并计算生成符号矩阵的根求方程下列各题运用符号运算功能计算六 的解函数求解微分方程用1)0()0()0()0(,2.4)4(='''=''='==+x x x x e x x dsolve t六、参考程序：1． solve('exp(x)=x^2+2')2．syms x sA=[x*sin(x),log(x);exp(-x),1/(s+x)]diff(A)int(A)A*A3．syms xlimit((cos(x)-cos(x)^(1/2))/x/(exp(sin(x))-1))4．dsolve('D4x+x=2*exp(t)','x(0)=1','Dx(0)=1','D2x(0)=1','D3x(0)=1')⎪⎩⎪⎨⎧=-++=--++03320224:.2122212212221x x x x x x x x x fsolve 组求解下面的非线性方程函数以及符号运算功能用七解一：function y=fc(x)y(1)=4*x(1).^2+x(2).^2+2*x(1).*x(2)-x(2)-2y(2)=2*x(1).^2+x(2).^2+3*x(1).*x(2)-3y=[y(1) y(2)];x0=[1 1];fsolve(@fc,x0)解二: syms x1 x2[x1 x2]=solve(‘4*x1^2+ x2^2+2* x1* x2- x2-2=0’,’ 2*x1^2+ x2^2+3* x1* x2-3=0’,’ x1’,’ x2’)《MA TLAB程序设计》参考程序一、(30分)1．D =0 -1F =2 02 00 62．A1 =1 3 3 0A2 =0 5 0 -3-3 0 5 03．max(P,[ ],1)=2 2 5 4 mean(P,2)=324．x=155．b =0 1 4-1 0 1-4 -1 06．z =20107．285二、(10分)%1t=0:0.1: 2*pi;x= sin(t).^3;y= cos (t).^3;plot(x,y)%2[x,y]=meshgrid(-2:0.1:2);z=(x.^2+1).* exp(-x.^2- y.^2+)+x.^2.*y; mesh(x,y,z)三、(16分)syms xlimit((x-1)*log(x)/(exp(x-1)-1),x,1)%2solve('exp(x)-3*cos(x)+1=0')%3syms xM=[ exp(-x), x;sin(x),exp(-0.1*x)* cos(3*x)]diff(M,2)int(M,x,0,pi)%4 syms t xdsolve('Dx-6*x=exp(-t) ','x(0)=1')四、(12分)x=0.5:9.5;y=[ 0.16, 0.40, 0.66, 0.75, 0.88, 1.0, 1.1, 1.1, 1.2, 1.3];x1=1./x;y1=1./y;p=polyfit(x1,y1,1);a=p(2); b= p(1);ye=x./(a*x+b);plot(x,y,'ro', x, ye, 'b-')grid onxlabel('自变量x ')ylabel('因变量y ')title('因变量y与自变量x的拟合曲线图')五、(12分)%参考程序cj=round(normrnd(70,10,100,4));for i=1:100for j=1:4if cj(i,j)>100 cj(i,j)=100;elseif cj(i,j)<0 cj(i,j)=0;endendendkmean=round(mean(cj))disp('课程1 课程2 课程3 课程4')disp(sprintf(' %4d%7d%8d%8d ' ,kmean(1),kmean(2),kmean(3),kmean(4))) studentmean=round(mean(cj,2));studentnumber=[1:100]';[studentnumber, studentmean]c=sum(sum(cj>=90))d=sum(sum(cj<60,2)>=2)六、(10分)(1). function y=fun1(x,n)y=0;for k=1:ny=y+ ((-1)^k*sin(k*x))/k;endy(2). x=0:0.1:2*pi;y1= fun1 (x,20);y2= fun1 (x,50);subplot(211),plot(x,y1)subplot(212),plot(x,y2)七、(10分)第一步：先将高阶微分方程转化为一阶微分方程组：选择状态变量x x x x ='=21,，则原方程化为2)0(,0)0(.,32112211==⎩⎨⎧='++-='x x x x t x tx x第二步：建立函数文件function dxdt = fun2(t,x)dxdt = [-3*t*x(1)+ x(2)+t; x(1)];第三步：求解微分方程，命令如下：[t,x]=ode45(@fun2,[0,10],[0;2])补充：1.matlab 中清屏的命令是———clc2.matlab 中用于清除工作空间变量的命令是——clear3.在【0,5】上求函数f=（x-3）^2-1的最小值。
matlab程序设计例题及答案
matlab程序设计例题及答案1.编写程序：计算1/3+2/5+3/7+……+10/21法一： s=0;for i=1:10s=s+i/(2*i+1); end ss =法二：sum((1:10)./(3:2:21)) ans =2.编写程序：计算1~100中即能被3整除，又能被7整除的所有数之和。
s=0;for i=1:100if mod(i,3)==0&&mod(i,7)==0 s=s+i; end,end ss =2103.画出y=n！的图,阶乘的函数自己编写，禁用MATLAB 自带的阶乘函数。
x=1:10; for i=1:10try y(i)=y(i-1)*i; catch y(i)=1; end,end plot(x,y)106123456789104.一个数恰好等于它的因子之和，这个数就称为完数。
例如，6的因子为1，2,3，而6=1+2+3，因此6就是一个完数。
编程找出20XX以内的所有完数。
g=;for n=2:20XX s=0;for r=1:n-1if mod(n,r)==0 s=s+r; end endif s==ng=[g n]; end end gg =6 28 4965.编写一个函数，模拟numel函数的功能，函数中调用size函数。
function y=numelnumel(x) m=size(x); y=m(1)*m(2);numelnumel([1 2 3;4 5 6])ans =66. 编写一个函数，模拟length函数的功能，函数中调用size函数。
function y=lengthlength(x) m=size(x);y=max(m(1),m(2));lengthlength([1 2 3;4 5 6])ans =37.求矩阵rand的所有元素和及各行平均值，各列平均值。
s=rand(5);sum=sum(sum(s)) mean2=mean(s,2) mean1=mean(s)sum =mean2 =mean1 =8.编程判断1001,1003,1007，1009，1011为素数，若不是，输出其约数。
matlab试题及答案
matlab试题及答案# MATLAB试题及答案一、选择题1. MATLAB的基本数据单位是：A. 矩阵B. 向量C. 标量D. 数组答案：A2. 下列哪个命令可以用来绘制函数图形？A. `plot`B. `graph`C. `draw`D. `chart`答案：A3. MATLAB中，以下哪个是正确的矩阵转置操作？A. `transpose(A)`B. `A'`C. `A^T`D. `flip(A)`答案：B二、简答题1. 简述MATLAB中矩阵的基本操作。
答案：在MATLAB中，矩阵是最基本的数据结构，可以进行加、减、乘、除等基本运算。
矩阵的创建可以使用方括号`[]`，例如`A = [1 2;3 4]`。
矩阵的转置使用单引号`'`，例如`A'`。
矩阵的求逆使用`inv`函数，例如`inv(A)`。
2. MATLAB中如何实现循环结构？答案：MATLAB中实现循环结构主要有两种方式：`for`循环和`while`循环。
`for`循环用于已知迭代次数的情况，例如：```matlabfor i = 1:5disp(i);end````while`循环用于迭代次数未知的情况，例如：```matlabi = 1;while i <= 5disp(i);i = i + 1;end```三、计算题1. 给定矩阵A和B，请计算它们的乘积C，并求C的行列式。
A = [1 2; 3 4]B = [5 6; 7 8]答案：首先计算矩阵乘积C：```matlabC = A * B;```然后计算C的行列式：```matlabdetC = det(C);```结果为：```matlabC = [19 22; 43 50]detC = -16```2. 编写一个MATLAB函数，计算并返回一个向量的范数。
答案：```matlabfunction norm_value = vector_norm(v)norm_value = norm(v);end```四、编程题1. 编写一个MATLAB脚本，实现以下功能：- 随机生成一个3x3的矩阵。
Matlab考试卷试题库题集答案.doc
一、选择题共 10 小题二、填空题共 10 空三、看程序填空四、编程题 3 题填空题（计算部分）1、标点符号 ; _______可以使命令行不显示运算结果，%——用来表示该行为注释行。
2、x 为 0 ～4pi，步长为0.1pi的向量，使用命令_______创建。
x=0:0.1*pi:4*pi3、输入矩阵A= ，使用全下标方式用A(2,2) 取出元素“-5 ”，使用单下标方式用_______取出元素“-5 ”。
A(5)4、符号表达式sin(2*a+t)+m 中独立的符号变量为_______。
t5、M脚本文件和 M函数文件的主要区别是M脚本文件没有函数定义和件有函数定义 _______。
M函数文6.设 x 是一维数组， x 的倒数第 3 个元素表示为 _______；设 y 为二维数组，要删除 y 的第 34 行和 48 列，可使用命令; _______；x(_end-2_)y(34,:)=[]y(:,48)=[]7.将变量 x 以 Ascii 文本格式存储到文件 fname.txt，应使用命令 _________ _；save _x8.在 while 表达式 , 语句体 , End 循环语句中，表达式的值 __ __时表示循环条件为真，语句体将被执行，否则跳出该循环语句；非零9.要从键盘读入一个字符串并赋值给变量x，且给出提示“ Who is she?”，应使用命令 _________；x=input( ‘ Who is she? ’ , ’ s’ )_．设和B= 和 C= 均为 m*n 矩阵，且存在于 WorkSpace10A=中，要产生矩阵 D= ，可用命令 ________ _，计算可用命令 ________;D=(A-C)/B.^Cdet(inv(A ’ *B) 11.在 MATLAB 命令窗口中的“>>”标志为 MATLAB 的_______ 提示符，“│”标志为 _______提示符。
matlab开卷考试题及答案
matlab开卷考试题及答案MATLAB开卷考试题及答案一、选择题（每题2分，共20分）1. MATLAB中用于创建一个3x3的单位矩阵的命令是：A. `eye(3)`B. `ones(3)`C. `zeros(3)`D. `identity(3)`答案：A2. 下列哪个函数用于计算矩阵的逆？A. `inv()`B. `det()`C. `eig()`D. `trace()`答案：A3. MATLAB中，用于绘制二维图形的函数是：A. `plot()`B. `surf()`C. `mesh()`D. `bar()`答案：A4. 如何在MATLAB中计算向量的范数？A. `norm()`B. `abs()`C. `length()`D. `size()`答案：A5. MATLAB中，用于创建一个随机数矩阵的函数是：A. `rand()`B. `randi()`C. `randn()`D. `randperm()`答案：A6. 下列哪个命令用于在MATLAB中求解线性方程组？A. `solve()`B. `fsolve()`C. `fzero()`D. `linsolve()`答案：A7. MATLAB中，用于计算矩阵的行列式的函数是：A. `det()`B. `trace()`C. `eig()`D. `inv()`答案：A8. 如何在MATLAB中创建一个全1的3x3矩阵？A. `ones(3,3)`B. `eye(3,3)`C. `zeros(3,3)`D. `rand(3,3)`答案：A9. MATLAB中，用于计算矩阵的特征值的函数是：A. `eig()`B. `eigs()`C. `svd()`D. `qr()`答案：A10. 下列哪个函数用于计算矩阵的奇异值分解？A. `eig()`B. `svd()`C. `qr()`D. `lu()`答案：B二、填空题（每题3分，共30分）1. MATLAB中，用于创建一个3x3的零矩阵的命令是 `____`。
大学专业试题-Matlab期末考试及答案
大学专业试题-Matlab期末考试及答案一、填空题（一题3分，共15分）1、标点符号可以使命令行不显示运算结果。
答案：%2、x为0〜4，步长为0.1pi的向量，使用命令创建。
答案：x=0:0.1*pi:4。
3、输入矩阵A=[1 2 3;7 -5 8;9-1-7],使用全下标方式用A（2,2）取出元素“-5 ”，使用单下标方式用取出元素“-5 ”。
答案：A （5）。
4、在while表达式,语句体,End循环语句中，表达式的值__时表示循环条件为真，体语句将被执行，否则跳出该循环语句。
答案：非零。
5、要从键盘读入一个字符串并赋值给变量x,且给出提示“Hello"，应使用命令—答案：x=input（‘Hello,；s'）b二、选择题（一题3分，共15分）1、下列表示可以作为MATLAB的变量的是(D )A. abcd-2B. xyz_2#C. @hD. X_1_a2、如果x=1: 2 : 8,则U x(1)和x(4)分别是(B )A. 1, 8B. 1,7C. 2,8D. 2,73、指出下列错误的指令（B）。
A.symsa b;B. symsa, b;C.syms('a',’b');D. syms('a','b','positive');4、下列哪条指令是求矩阵的行列式的值(C)。
A .inv B. diag C. det D.eig5、清空Matlab工作空间内所有变量的指令是(C)。
A. clcB.clsC.clearD. clf三、程序分析题（一题10分，共20分）1、分析下面程序运行后s1,n1,s2,n2,m的值。
s1=0;s2=0;n1=0;n2=0；x=[1,-4,-8,3,10,-9,7,-3,10, 8,-5,-2,2,0 ]; m=length(x); fori=1:mif x(i)<0s1=s1+x(i); n1=n1+1; elses2=s2+x(i); n2=n2+1;end ends1,n1,s2,n2,m请回答s1,n1,s2,n2,m的结果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
8. 下面的程序执行后 array 的值为
for k=1:10 if k>6 break; else array(k) = k; end
end
(A) array = [1, 2, 3, 4, 5, 6]
(B) array = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
(C) array =6
3、 [ 2 5 ]
ห้องสมุดไป่ตู้
[5]
4、 ；
%
5、 if x<= 80 & x>60
6、 polyint(P,K)
7、 hold on
8、 save my_data
[k,r]= deconv (P,Q) grid on Load my_data
二、选择题
1、 A 2、 C 3、 B 4、 C 5、 C 6、 C 7、 D 8、 A 9、 C 10、 A
if x(k)==0 x(k)=k;
else x(k)=2*k+1;
end end disp(x);
3. 创建符号函数并求解，要求写出步骤和运行结果（7 分） (1)创建符号函数 f=ax2+bx+c (2)求 f=0 的解
4. 求解以下线性方程组，要求写出程序代码和运行结果（5 分）
2x1- 3x2+ x3+2x4=8 x1+3x2+ x4=6 x1- x2+ x3+8x4=1
三、写出程序的执行结果或写出给定要求的指令（总共 35 分）
1. 写出执行以下代码后C，D，E的值 (6分) A=[1,2,3;4:6;7:9]; C=[A;[10,11,12]], D=C(1:3,[2 3]) E=C(2,[1 2])
第3页共8页
2. 写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 (5分) x=[0,1,0,2,0,3,0,4]; for k=1:8
4. 清空 Matlab 工作空间内所有变量的指令是 (A) clc (B) cls (C) clear (D) clf
5. 用 round 函数四舍五入对数组[2.48 6.39 3.93 8.52]取整，结果为 (A) [2 6 3 8] (B) [2 6 4 8] (C) [2 6 4 9] (D) [3 7 4 9]
6. 已知 a=2:2:8, b=2:5，下面的运算表达式中，出错的为 (A) a' *b (B) a .*b (C) a*b (D) a-b
第2页共8页
7. 角度 x 30 45 60，计算其正弦函数的运算为
(A) SIN（deg2rad(x)） (B) SIN(x) (C) sin(x) (D) sin(deg2rad(x))
10. 求解方程 x4-4x3+12x-9 = 0 的所有解 (A)1.0000, 3.0000, 1.7321, -1.7321 (B)1.0000, 3.0000, 1.7321i, -1.7321i (C)1.0000i, 3.0000i, 1.7321, -1.7321 (D)-3.0000i, 3.0000i, 1.7321, -1.7321
(D) array =10.
9．i=2; a=2i; b=2*i; c=2*sqrt(-1); 程序执行后；a, b, c 的值分别是多少？
(A)a=4, b=4, c=2.0000i (B)a=4, b=2.0000i, c=2.0000i (C)a=2.0000i, b=4, c=2.0000i (D) a=2.0000i, b=2.0000i, c=2.0000i
一、填空题（每空 2 分，总共 30 分）
1. 已知 A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。
A & B 的运行结果 ans=
；
A ./ B 的运行结果 ans=
。
2. 产生 4 阶全 0 方阵的命令为
；产生 3 阶全 1 方阵的命令为
。
3. A=rand(2,5);b=size(A);c=length(A)；则 b 和 c 的值分别为
7x1+ x2-2x3+2x4=5
5．绘制函数曲线，要求写出程序代码(12 分) (1)在区间[0:2π]均匀的取 50 个点，构成向量π (2)在同一窗口绘制曲线 y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求 y1 曲线为红色点划线，标记点为圆 圈；y2 为蓝色虚线，标记点为星号
学号
栏 姓名 线
息
厦门理工学院试卷
2010－2011 学年 第 1 学期
课程名称
MATLAB 程序设计
专业 级 班级
试卷 A √ 卷别 B □
考试 闭卷 □ 方式 开卷 √
本试卷共 4 大题( 4 页)，满分 100 分，考试时间 120 分钟。
请在答题纸上作答，在试卷上作答无效。
上机考试说明： 1. 开考前可将准备程序拷到硬盘, 开考后不允许用移动盘，也不允许上网； 2. 领座考生试卷不同，开卷，可利用自己备用的书和其他资料，但不允许讨论，也不允 许借用其他考生的书和资料。 解答(指令行，答案等)全部用笔写在答题卷上。 更多 matlab 基础题目，请关注 算法源码云课堂。
四、使用 MATLAB 语言进行编程（15 分）
打印出所有的水仙花数。所谓“水仙花数”，是指一个三位数，其各位数字立方之和等于该数本身。
第4页共8页
答案：
一、填空题
1、
ans =
011
110
ans =
0 2.2500 2.0000
1.0000 0.6000
NaN
2、 zeros(4)
ones(3)
和
。
4. 标点符号 可以使命令行不显示运算结果， 用来表示该行为注释行。
5. 用 if 判断语句判断 80≥ x >60,在 MATLAB 中 if 语言后的判断应写
为
。
订
专业 08 级 级 班级
信
生
考
装
系 通信工程
电子与电气工程
第1页共8页
6. P, Q 分别是个多项式的系数矢量，求 P 对应的多项式的积分（对应的常数项为 K），使用的命令
x1 = 13/17
x2 = 320/187
x3 = 212/187
x4 = 19/187
5、 t = linspace(0, 2*pi, 50); plot(t,y1,'r-o' ,t,y2,'b--*');
四、使用 MATLAB 语言进行编程
for k = 100:999 a = fix(k/100); b = rem(fix(k/10),10);
是
；求 P/Q 的解，商和余数分别保存在 k 和 r，使用的命令是
；
7. 为了使两个 plot 的图形在同一个坐标显示，可以使用
用
命令为图形添加网格。
命令进行图形保持；可以使
8. MATLAB 的工作空间中有三个变量 v1, v2, v3，写出把它们保存到文件 my_data.mat 中的指
令
；写出把 my_data.mat 文件中的变量读取到 MATLAB 工作空间内的指
令
。
二、选择题（每空 2 分，总共 20 分）
1.下列哪个变量的定义是不合法的 (A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz
2. 下列哪条指令是求矩阵的行列式的值 (A) inv (B) diag (C) det (D) eig
3. 在循环结构中跳出循环，执行循环后面代码的命令为 (A) return (B) break (C) continue (D) keyboard
4、 >> syms a b c x >> syms x1 x2 x3 x4 >> clear >> syms x1 x2 x3 x4 >> f = 2*x1 + 3*x2 + x3 + 2*x4 -8; >> g = x1 + 3*x2 + x4 - 6; >> h = x1 - x2 + x3 + 8*x4 - 1; >> i = 7*x1 + x2 - 2*x3 + 2*x4 -5; >> [x1,x2,x3,x4] = solve(f,g,h,i)
7 17
第6页共8页
c = rem(k,10); if a.^3 + b.^3 + c.^3 == k
fprintf( “%u,\t\t”, k ) ; end end
第7页共8页
第8页共8页
三、写出程序的执行结果或写出给定要求的指令
1、 c=
0 + 2.0000i D=
23 56 89 E= 45
第5页共8页
2、 1 5 3 9 5 13
3、 >> syms a b c x >> s = a*x^2 + b*x + c; >> solve(s)
ans = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]