Matab考试题库附加答案
一、辨析题（共2小题，每题 5 分，共计 10 分）
二、问答题（共3小题，每题 5 分，共计 15 分）
三、计算题（共7小题，每题 5 分，共计 35 分）
四、编程题（共4小题，每题 10 分，共计 40 分）
填空题（计算部分）
1、标点符号; _______可以使命令行不显示运算结果，
%——用来表示该行为注释行。
2、x为0～4pi，步长为0.1pi的向量，使用命令_______创建。
x=0:0.1*pi:4*pi
3、输入矩阵A=，使用全下标方式用A(2,2) 取出元素“-5”，使
用单下标方式用_______取出元素“-5”。
A(5)
4、符号表达式sin(2*a+t)+m中独立的符号变量为_______。
t
5、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 _x
8. 在while 表达式, 语句体, End 循环语句中，表达式的值__ __时表示循环条件为真，语句体将被执行，否则跳出该循环语句；
非零
9.要从键盘读入一个字符串并赋值给变量x，且给出提示“Who is she?”，应使用命令_________；
x=input(‘Who is she?’,’s’)_
10．设A=和B=和C=均为m*n矩阵，且存在于WorkSpace 中，要产生矩阵D=，可用命令________ _，
计算可用命令________;
D=(A-C)/B.^C
det(inv(A’*B)
11. 在MATLAB命令窗口中的“>>”标志为MATLAB的_______提示符，“│”标志为_______提示符。
命令行
输入
12.已知A=[1 2 3；4 5 0；7 8 9]；B=[1 0 3；1 5 0；0 1 2]；写出下列各指令运行的结果。
A+B；
A.*B；
A==B ；
ans= [2,2,6;5,10,0;7,9,11]
ans= [1,0,9;4,25,0;0,8,18]
ans= [1,0,1;0,1,1;0,0,0]
13.已知A是矩阵，求A的对角矩阵函数是_______，
求A的下三角矩阵函数是_______。
diag
tril
14.MATLAB的程序文件和Simulink模型文件的扩展名分别是_______、。
.m
.mdl
15.MATLAB最基本的绘图函数为_______。
plot()
16. A=[1,2,3;4,5,6]; A(:,[1,3])=[];A=__________________
[2；5]
17. fix(-1.5)=___ ________, round(-1.5)=__ _______________.
-1 -2
18.已知A=[1 3;4 6];C=[A,eye(size(A));A ,ones(size(A))] 则Ｃ＝
_____________.
[1 3 1 0；4 6 0 1；1 3 1 1；4 6 1 1]_ 19. A=[1,2,3;4,5,6]; A(4)=__________, A(3,2)=__________________
5_ 无答案
20.A=[1,2;3,4];B=[1,0;0,1]; A*B=__________, A.*B=____________.
[1 2；3 4]
[1 0；0 4]
21. A=[1 2 3;6 5 4]; reshape(A,3,2)=______________.
[1 5；6 3；2 4]
22. A=[1,2;3,1]; B=[1,0;0,1]; A~=B=___________.
[0 1；1 0]
23. A=[1,2,3;1,4,9;1,7,27]; triu(A,-1)=_____.提取矩阵上三角矩阵的函数
[1 2 3；1 4 9；0 7 27]
24. A=[4,15,-45,10,6;56,0,17,-45,0]; find(A>=10 & A<=20)=___ __________.
3 6 7
25.已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。
A &
B 的运行结果_______；
A ./
B 的运行结果_______。
ans=[0 1 1;1 1 0]
ans= [0 2.25 2;1 0.6 NaN]
26. 产生4阶全0方阵的命令为_______；
产生3阶全1方阵的命令为。
zeros(4)
ones(3)
27. A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为_______。
2 5和 5
28. 用if判断语句判断80≥ x >60,在MATLAB中if语言后的判断应写为_______。
if x<=80 & x>60
29.为了使两个pl ot的图形在同一个坐标显示，可以使用_______命令进行图形保持；可以使用_______命令为图形添加网格。
hold on
grid on
30. MATLAB的工作空间中有三个变量v1, v2, v3，写出把它们保存到文件my_data.mat中的指令_______；
save my_date
31. MATLAB提供了两种运行方式，是按_______和把_______文件copy到matlab里，回车。
F5
*.m
32．MATLAB的含义为_______。
矩阵实验室（Matrix Laboratory）33.MATLAB中，都包括哪些工具箱，举例不少于6个：_______
偏微分方程工具箱、最优化工具箱、插值运算工具箱、数理统计工具箱、信号处理工具箱、通信工具箱，滤波设计工具箱。
34、在MATLAB中，plot命令用于什么，。画曲线
35、size命令的功能，确定_______。阵列的长度
36、a=[1 2 2]，求a的转置_______。a=[1;2;2]
37、在线形图型格式的设置中，字符r表示什么颜色。红色
38. 在MATLAB中,矩阵的秩由函数_______来实现。rank
39．在MATLAB中,单位矩阵由函数_______来实现。eye
40. 在MATLAB中,二维饼图由函数_______来实现。pie
41. 在线形图型格式的设置中，字符b表示什么颜色。蓝色
42. MATLAB表达式3*2^3^2的结果是_______。192
43.条形图由函数_______来实现。bar
44. 倘若要是对x进行赋值，从5到25，间隔是0.1, x=5:0.1:25;
45. 倘若要是对x进行赋值，从2到20，中间间隔100个点,
x=linspace(2,20,100);
选择题（计算部分）
1. 下列变量中的最大数是()
A．eps B．realmin C．realmax D．-pi
2. 下列表示可以作为MATLAB的变量的是()
A．abcd- 2 B．xyz_2# C．@h D．X_1_a
3. 如果x=1: 2 : 8,则x(1)和x(4)分别是()
A．1，8 B．1, 7 C．2, 8 D．2, 7
4. 下列变量中比0大的最小数是()
A．i B．realmin C．realmax D．pi
5. 在MATLAB中下列数值的表示不正确的是（）
A．+99 B．1.3e-5 C．2-3*e^2 D．3-2*pi
c d b b b
6. 下列表示可以作为MATLAB的变量的是()
A．abcd-2 B．xyz_2a C．@h D．3hour
7. 如果x=1:-2:-8,则x(1)和x(5)分别是()
A．1，-8 B．-1,-7 C．-1,-8 D．1,-7
8. MATLAB表达式2*2^3^2的结果是()
A．128 B．4096 C. 262144 D．256
9. 下列哪个变量的定义是不合法的（）。
(A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz
10. 指出下列错误的指令（）。
(A) syms a b; (B) syms a, b; (C) syms(…a', …b'); (D) syms(…a','b','positive');
b d a A B
11.下列哪条指令是求矩阵的行列式的值（）。
(A) inv (B) diag (C) det (D) eig
inv是求逆矩阵
diag是矩阵对角元素的提取和创建对角阵
sqrt(x) ——求x的平方根
abs(x)——求x的绝对值
det(a)——求行列式的值
eig——计算矩阵A的特征值和特征向量的函数
clf; 用来清除图形的命令
12 .清空Matlab 工作空间内所有变量的指令是（）。
(A) clc (B) cls (C) clear (D) clf
13、下列变量名中————是合法的。
(A) char_1 ; (B) x*y ; (C) x\y ; (D) end
14.以下哪个说法是正确的（）。
(A)Matlab 进行数值计算的表达精度与其指令窗口中的数值显示精度相同。
(B)Matlab 指令窗口中显示的数值有效位数不能超过7 位。
(C)输入二维数值数组时，需要用到逗号和分号，它们可以在中文状态下输入。
(D)历史指令窗口所记录的内容与diary 指令所产生“日志”内容是不同的。
15. 产生四维单位矩阵的语句为().
A.ones(4)
B.eye(4)
C.zeros(4)
D.rand(4) >> eye(4)
ans =
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1 >> ones(4)
ans =
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1
>> zeros(4)
ans =
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0 >> rand(4)
ans =
0.8147 0.6324 0.9575 0.9572
0.9058 0.0975 0.9649 0.4854
0.1270 0.2785 0.1576 0.8003
0.9134 0.5469 0.9706 0.1419
C C a
D b
16.求解线性方程组Ax=b,当det(A)≠0时，方程的解是().
A. A\b
B.b/A
C.b\A
D.A/b
17. 运行如下程序后，输入９回车，命令窗口（command windows）显示的结果为()
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')));
elseif c>='0'& c<='9'
disp(abs(c)-abs('0'));
else
disp(c);
end
A. 9
B. 8
C. 7
D. 6
18. 运行下面的代码后，输入：１回车，２回车，１回车，则输出结果为（）a=input('a=?');
b=input('b=?');
c=input('c=?');
d=b*b-4*a*c;
x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];
disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);
Ａ. x1=num2str(1) x2=num2str(1)
B. x1=num2str(-1) x2=num2str(1)
C. x1= -1 x2=-1
D. x1=1 x2=-1
19. 清除工作空间（wordspace）的命令是（）
Ａ. clc; B. clear; C. clf; D.delete;
20. 运行如下程序后，输入９回车，命令窗口（command windows）显示的结果为( )
x=input('请输入x的值:');
if x==10
y=cos(x+1)+sqrt(x*x+1);
else
y=(3^(1/2))*sqrt(x+sqrt(x));
end
y
A. 9
B. 8
C. 7
D. 6
a a c
b d
21. 运行如下程序后, 命令窗口（command windows）显示的结果为( b)
s=0;
a=[12,13,14;15,16,17;18,19,20];
for k=a
s=s+k;
end
disp(s');
A. 144
B. 39 48 57
C.145
D. 45 48 51
22. 运行如下程序后, 命令窗口（command windows）显示的结果为(b )
k=0;
for n=100:200
if rem(n,21)~=0 %R=rem(X,Y),求余数函数,X,Y应该为正数
k=k+1;
continue
end
break;
end
A.105
B. 5
C.4
D.101
23. 图形窗口（Figure）显示网格的命令是（b）
A. axis on
B. grid on
C. box on
D. hold on
24. 已知函数文件如下，则factor(4)=(c )
function f=factor(n)
if n<=1
f=1;
else
f=factor(n-1)*n;
end
A. 4
B. 12
C. 24
D.48
25. 运行如下程序后, 命令窗口（command windows）显示的结果为( d)
A=[13,-56,78; 25,63,-735; 78,25,563; 1,0,-1];
y=max(max(A))
A. y=564
B.y=9
C.y=-735
D.y=563
26. 在图形指定位置加标注命令是（c）
A. title(x,y,?y=sin(x)?);
B. xlabel(x,y,?y=sin(x)?);
C. text(x,y,?y=sin(x)?);
D. legend(x,y,?y=sin(x)?); %添加图例的标注,
27.下列哪个函数为插值函数（b）
A. P=polyfit(X,Y,3)
B. Y1=interp1(X,Y,X1,'method')
C. [Y,I]=sort(A,dim)
D. R=corrcoef(X)
28．i=2; a=2i; b=2*i; c=2*sqrt(-1); 程序执行后；a, b, c的值分别是多少？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
29. 求解方程x4-4x3+12x-9 = 0 的所有解a
(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
30、在循环结构中跳出循环，但继续下次循环的命令为 。c (A) return; (B) break ; (C) continue ; (D) keyboard
31. 用round 函数四舍五入对数组[2.48 6.39 3.93 8.52]取整，结果为 c (A) [2 6 3 8] (B) [2 6 4 8] (C) [2 6 4 9] (D) [3 7 4 9]
32. 已知a=2:2:8, b=2:5，下面的运算表达式中，出错的为 c (A) a' *b (B) a .*b (C) a*b (D) a-b
辨析题：
1 MATLAB 有几种建立矩阵的方法？各有什么优点？。
2 MATLAB 提供编程语言，请辨析M 文件编程中的命令文件与函数文件两种形式。
3 辨析matlab 中的拟合与插值有什么区别
4. 判断下列一段程序的错误，并进行改正。
x=-2:0.2:2;
y=x.^2; plot(x,y)
str1 = num2str(min(x));% 数字转化为字符串 str2 = num2str(max(x));% 数字转化为字符串 out = [Value of f from 'str1' to 'str2']; xlabel(out);
out = ['Value of f from ' str1 ' to ' str2];
或out = ['Value of f from str1 to str2'];
5根据下面两个矩阵
??
??
?
?????=??????????=987654321,136782078451220124B A 执行下列的矩阵运算命令，并回答有关的问题
（1） A.*B 和 A*B 将分别给出什么结果，它们是否相同，为什么？ （2） 得出A.^B 、A/B 及A\B 的结果，并分别解释它们的物理意义。
问答题：
1 MATLAB系统包括哪些窗口，各个窗口功能是什么？
2 分析下面程序运行后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);
for i=1:m
if x(i)<0
s1=s1+x(i);
n1=n1+1;
else
s2=s2+x(i);
n2=n2+1;
end
end
s1,n1,s2,n2,m
请回答s1,n1,s2,n2,m的结果。
3 已建立一个简单的SIMULINK模型，如下图，用信号发生器产生一个仿真图，请对下图
的符号和功能进行解释。
3 A = [1 2; 3 4]; B = [1 2; 1 4];
写出下面几条指令的运行结果：
C = A*B
D = A+B
E = A.*B
答：
C =
3 10
7 22
D =
2 4
4 8
E =
1 4
3 16
4.mysum=0；i=1；
while(i<=100)
mysum=mysum+i；
i=i+1；
end
mysum
答：
mysum =
5050
5.在MATLAB窗口中，用哪两种方法启动Simulink？答：
在MATLAB命令窗口中输入“simulink”再回车，或者单击工具栏中的SIMULINK 图标
6. 请修改下面的程序，让他们没有for循环语句！
A=[1 2 3; 4 5 6; 7 8 9];
[r c]=size(A);
for i=1:1:r
for j=1:1:c
if (A(i,j)>8 | A(i,j)<2)
A(i,j)=0;
end
end
end
解：程序如下：a=[1 2 3;4 5 6;7 8 9];
b=(a<=8&a>=2);a=a.*b
嫁错恶灵进错门黑岩https://www.doczj.com/doc/494224609.html,/seshu/32687/index.html
7 编写一个函数，使其能够产生如下的分段函数：
并调用此函数，绘制x=[0,+2]范围内的f(x)*f(x+2)。
答：function y=f(x)
if x<=2
y=0.5*x;
else if x>6
y=0.5;
else y =1.5-0.25*x;
end
end
运行结果x=2
f(x)=1
x = 0:0.05:2;
y = diag(A2(x)'*A2(x+2));
plot(x,y);
xlabel('\bfx');
ylabel('\bfy');
0.2
0.4
0.6
0.8
1 1.2
1.4
1.6
1.8
2
00.050.10.150.20.250.30.350.40.45
0.5x
y
8（1） 将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。 ①
,x
（曲线图）；②f(x,y)=(-2<=x<=2,-4<=x<=4)
（曲
面图）。
答：subplot(2,2,2);
ezplot('(cos(x))^(1/2)',[-pi/2 pi/2]); ylabel('y'); subplot(2,2,3); x=-2:0.5:2; y=-4:1:4;
ezsurfc('x^2/2^2+y^2/4^2')
-1
01
0.5
1
x
(cos(x))1/2
y
x x 2/22+y 2/42
y
(2) A 是一个維度m×n的矩阵. 写一段程序, 算出A中有多少个零元素。答：>> A= input ('请输入一个矩阵')
[m,n]= size(A);
sig=0;
for i=1:m
for j=1:n
if A(i,j)==0
sig = sig+1;
end
end
end
请输入一个矩阵[0 1 2;1 0 2; 0 0 0]
A =
0 1 2
1 0 2
0 0 0
>>sig
sig =
5
9 向量A=[]. 写一段程序, 找出A中的最小元素。
答：A= input ('请输入一个向量')
[m,n]=size(A)
min =A(1,n);
for i=1:n
if A(1,i) min=A(1,i) end end 请输入一个向量[1 2 3 -5 2 0] A = 1 2 3 -5 2 0 min = -5 10 写出执行以下代码后C，D，E的值A=[1,2,3;4:6;7:9]; C=[A;[10,11,12]], D=C(1:3,[2 3]) E=C(2,[1 2]) C = 1 2 3 4 5 6 7 8 9 10 11 12 D = 2 3 5 6 8 9 E = 4 5 11 写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值。x=[0,1,0,2,0,3,0,4]; for k=1:8 if x(k)==0 x(k)=k; else x(k)=2*k+1; end end disp(x); 1 5 3 9 5 13 7 17 12 创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c (2)求f=0的解 syms a x b c; f=a*x^2+b*x+c; solve(f) ans = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2)) 13 求解以下线性方程组，要求写出程序代码和运行结果。 2x1-3x2+x3+2x4=8 x1+3x2+ x4=6 x1-x2+x3+8x4=1 7x1+x2-2x3+2x4=5 解：程序代码：a=[2 -3 1 2;1 3 0 1;1 -1 1 8;7 1 -2 2]; b=[8 6 1 5]'; ra=rank(a); rb=rank([a b]); det(a); xx=a\b 运行结果： xx = 3.1079 1.3278 7.9502 -1.0913 14. (1) 写出关系运算符中的等号、不等号、小于号、大于号、小于等于号和 大于等于号； (2) 写出逻辑操作中逻辑“与”、逻辑“或”及逻辑“非”的符号； (3) 并用语句表达“如果a等于b且c等于0就让d 等于3，否则如果 a大于b且c=1让d等于0，其余情况让d等于3”； 解：（1）== ~= < > <= >= （2）& | ~ （3）if a==b&c==0 d=3 else if a>b&c= =1 d=0 else d=3 end 15. 设A= 和B=是两个10行10列的矩阵（数组），试说 明命令A*B, A/B, A\B, A.*(B.^A), A./B, A.\B和sin(A).*exp(B)的涵义。 解：A,B的矩阵乘法，矩阵A乘以B的逆，矩阵B乘以A的逆，B的每个元素的对应A次幂乘以A对应的每个元素 A的每个元素除以B对应的每个元素，B的每个元素除以A对应的每个元素，A的每个元素的正弦与B对应的每个元素的指数求积， 16. 设A为14*17数组，B为一个行数大于3的数组，请给出 （1）删除A的第4、8、12三列的命令；（2）删除B的倒数第3行的命令； 解：(1)A(:,[4 8 12])=[] (2) B(end-2)=[] 编程简题 1 利用MATLAB数值运算，求解线性方程组(将程序保存为？？？.m文件) 2 求下列联立方程的解 3x+4y-7z-12w=4 5x-7y+4z+ 2w=-3 x +8z- 5w=9 -6x+5y-2z+10w=-8 求系数矩阵的秩；求出方程组的解。 解：（1） >> a=[3 4 -7 -12 5 -7 4 2 ; MATLAB数学实验练习题 ————————————————————————————————作者：————————————————————————————————日期： “MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) > ezplot('exp(x)-3*x^2') >> grid on ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2）1sin 02 x x - =至 少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3,x,0) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x) 1、求以下变量的值，并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A 实验二 基本矩阵操作 实验目的： ① 掌握matlab 变量和数据操作； ② 掌握matlab 矩阵的创立、拆分及特殊矩阵； ③ 掌握matlab 运算，掌握matlab 在矩阵分析中的应用，掌握稀疏矩阵的存储方式和创建方法； ④ 掌握矩阵的基本运算和结构体数组的创建。 ⑤ 了解字符串处理函数。 实验要求:给出程序和实验结果。 实验内容： 1、利用列向量()1,2,3,,10T 建立一个范得蒙矩阵A ，并利用位于矩阵A 的奇数行偶数列的元素建立一个新的矩阵B ，须保持这些元素的相对位置不变。 2、矩阵的基本运算与点运算的区别。 3、给出矩阵的两种存储方式的联系和区别，这两种存储方式在实际应用中主要应用于具有什么特点的矩阵？ 4*．将字符串'very good'转换为等值的整数。 5．按水平和竖直方向分别合并下述两个矩阵： 6． 分别删除第5题两个结果的第2行。 7． 分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。 8． 分别查看第5题两个结果的各方向长度。 9*． 分别判断pi 是否为字符串和浮点数。 10．分别将第5题两个结果均转换为2?9的矩阵。 11．计算第5题矩阵A 的转置。 12．分别计算第5题矩阵A 和B 的A +B 、A .* B 和A \B 。 13．判断第5题矩阵A 和B 中哪些元素值不小于4。 14*．分别用函数strcat()和矩阵合并符合并如下字符串：' The picture is '和' very good '。 15*．创建字符串数组，其中元素分别为‘Picture ’和'Pitch '。 16*．在第14题结果中查找字符串'e'。 17*．在第15题结果中匹配字符串'Pi'。 result5 = ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么如合法， 结果是多少 >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 58 11 >> result4=b*d result4 = 31 22 22 40 49 13 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 415 30 >> result7=a./b result7 = >> result8= Attempt to reference field of non-structure array. >> result9=a.\b result9 = >> result10=a92 result10 = 1 4 9 16 25 36 >> resultl 1=29a result11 = 2 4 8 16 32 64 >> result5=[b;c']*d 3、用MATLAB求解下面的的方程组。 1、求以下变量的值，并在MATLAB^验证。 1 2 x1 3 2 x2 11 5 x3 2 1 3 x4 >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 7 2 1 2 9 15 3 2 2 2 11 5 1 3 2 13 (1)求矩阵A的秩(rank) (2)求矩阵 A 的行列式(determinant) (3)求矩阵 A 的逆(inverse) (4)求矩阵 A 的特征值及特征向量 (eigenvalue and eigenvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3) >> a=det(A3) >> [V,D]=eig(A3) 10 n 10 查看y 的值) m1=0; for m=-10:10 m仁m1+2^m; end m1 m1 = 6、求分段函数的值。 用if 语句实现，算出下列表中x 对应的y 值。 x=input('enter x='); if x<0 y=x A2+x-6; elseif x>=0&&x<5 y=xA2-5*x+6; else y=xA2-x-1; end y 7、分别用if 和switch 语句实现，将百分 制成绩转换为成绩等级A、B、C、D、E。 其中90~1 00分为A，80~89 分为B，70~79 分为C，60~69 分为D，60 分以下为E。 对超出百分制范围的成绩，给出错误提示 信息。 if 结构程序： x=input('please enter score='); if x>=90&&x<=100 9 2 10 disp('A') 7 2 9 15 (1) 2 2 1 3 4 7 1 0 A 4、已知 2n 2 10 29 一、填空题：（每空2分）大约20空 （注：写出完成要求所需的命令语句或题目给定语句执行后输出的 结果） 1. Matlab 中用来表示不确定值的变量是 NaN ；默认变量i 表示 虚数单位 。 2. Matlab 中符号 ； （分号） 可以使命令行不显示运行的结果，表示该行为注释行的符号为 % （百分号） 。 3. Matlab 将工作空间所有变量保存在默认目录下的命令是 Save 。(命令复习) 4. 骰子会掷出1到6的随机数字，用函数rand 和取整函数(Important)模拟骰子某次投掷出的数字的命令为 ceil(6*rand) ，模拟两 个骰子某次投掷的数字组合之和情况的命令为 。 5. 生成一个含有10个元素的行向量x ，其元素在[-1,1]之间均匀分布的的命令为 x=2*rand(1,10)-1 ，使用向量x 的后5个元素作为对角线 元素生成对角阵A 的命令为 A=diag(x(end-4:end)) ，将A 的主对 角线下方的元素（不含主对角线）全部填充为2014的命令为 B=ones(5); 2014*B-triu(B)+A 6. 为函数22 110at t x +=分布定义内联函数命令是 ，定义匿名函数的命令是 7. 在循环结构中使用 break 语句可以终止当前循环。 8. Matlab 用于绘制统计直方图的命令是 hist ；用于绘制y 轴为对数坐标的单对数坐标图的命令为 semilog 。 二、解答题：大约七道题目 1. 在[0,π2]区间上绘制如下图所示的sin(x)与sin(x 2)的曲线，并如图为图像添加标题以及为坐标轴添加标注。 2.绘制上面曲面的图形（提示：该曲面由上下两部分组成：） 3.创建命令式的M 文件计算圆形的面积S ，提示用户输入圆形的半径r ，并输出圆形的面积。（注意输入输出都要对用户添加一定的提示语） %Function Circle 4.请利用matlab 的符号工具箱做下面的运算，并写出相应的命令。 (1) 求极限x e x L x x -+=→10)1(lim ； (2)求积分：dx x x ?∞+02sin ； Syms x; f=((1+x)^(1/x)-exp(1)/x; limit(f,x,0) 2222, [2,2], [1,1]z x y x y =+∈-∈- 复习题 1、写出3 2、i ｎv(A)表示Ａ的逆矩阵； 3、在命令窗口健入 clc,４、在命令窗口健入clea 5、在命令窗口健入６、x=-1:０.2:17、det(A)表示计算Ａ的行列式的值；8、三种插值方法：拉格朗日多项式插值,分段线性插值，三次样条插值。 9、若A=123456789?? ????????,则fliplr （Ａ）＝321654987?????????? Ａ-３=210123456--??????????A ．＾２=149162536496481?????????? tril(A)=100450789?????????? tri ｕ（Ａ，－1)=123456089??????????ｄiag(A ）＝100050009?????????? Ａ(:,2),＝2 58A(3,:)=369 10、nor ｍcd ｆ（1，1，2)=0.5%正态分布mu=1,s ｉgm ａ＝２,x ＝1处的概率 ｅ45(@ｆ,［ａ,b ］,x0),中参数的涵义是@fun 是求解方程的函数M 文 件,[a,b ］是输入向量即自变量的范围a 为初值,ｘ0为函数的初值，t 为输出指定的[ａ,b],x 为函数值 １5、写出下列命令的功能：te ｘt （1,2,‘ｙ=s ｉｎ(x)’ hold on 16fun ｃｔion 开头; 17 ,4) ３，４） 21、设x 是一向量，则）的功能是作出将Ｘ十等分的直方图 ２2、interp １(［1,2,3]，[3,４,5]，2.5) Ans=4.５ 23、建立一阶微分方程组? ??+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来） 二、写出运行结果: １、>＞ｅy ｅ（3，4)=1000 01000010 2、>>s ｉze([1,2,3])=1;3 3、设b=ro ｕnd （unifrnd(-５，5,１，4)),则=3 ５ 2 -5 >>[x,m]=ｍin(b)；x ＝－5；m=4 ，［x,n ］=sort(b ） －５ 2 3 ５ 4 3 １ 2 ｍea ｎ(ｂ）=１.2５，m ｅdian(b)＝2.5,range(ｂ）=10 ４、向量ｂ如上题,则 ＞>ａn ｙ(b),all(ｂ<2),ａｌl(b<6） Ans ＝１ 0 １ ５、＞>[5 6;7 8]>[7 8；5 6]=00 11 6、若1234B ??=???? ,则 7、>>ｄｉag(d ｉag （B ）)=10 04 8、＞>[4:-2:１].＊[-1,６]=-4 12 9、>>acos(0.5),a ｔan(1） ans ＝ 1.６5９8 ans= 第一套 华东交通大学2015-2016学年第一学期Matlab期末考试 一、填空题（4*6=24分） 1、在MATLAB命令窗口中的“>>”标志为MATLAB的提示符，“│”标志为提示符。 2、符号表达式sin(2*a+t)+m 中独立的符号变量为__ _。 3、在通常情况下，左除x=a\b是的解，右除x=b/a是的解，一般情况下，。 4、为了使两个plot的图形在同一个坐标显示，可以使用_____命令进行图形保持；可以使用_ __命令为图形添加网格。 5、倘若要是对x进行赋值，从5到25，间隔是0.1,为；倘若要是对x进行赋值，从2到20，中间间隔100个点,为 6、A=[1,2;3,4];B=[1,0;0,1]; A*B= _ , A.*B=_ _； 二、选择题（4*6=24分） 7、如果x=1: 2 : 8,则x(1)和x(4)分别是( ) A．1，8 B．1, 7 C．2, 8 D．2, 7 8、运行如下程序后，输入９回车，命令窗口（command windows）显示的结果为( ) 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'))); elseif c>='0'& c<='9' disp(abs(c)-abs('0')); else disp(c); end A. 9 B. 8 C. 7 D. 6 9、MATLAB表达式2*2^3^2的结果是() A．128 B．4096 C. 262144 D．256 10、在循环结构中跳出循环，执行循环后面代码的命令为( ) (A) return (B) break (C) continue (D) keyboard 11、在图形指定位置加标注命令是（） A. title(x,y,?y=sin(x)?); B. xlabel(x,y,?y=sin(x)?); C. text(x,y,?y=sin(x)?); D. legend(x,y,?y=sin(x)?); 12、下列哪个变量的定义是不合法的（）。 (A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz 三、程序题（52分） 13（10分）、在同一窗口绘制6个子图，已知t=0.1:0.1:2*pi，每个子图分别绘制（其中y1=sin(t),y2=cos(2t),y3=sin(t)+cos(t),y4=3t, y5=2ln(t),y6=et）； 14（10分）、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)： price<200 没有折扣 200≤price<500 3%折扣 500≤price<1000 5%折扣 1000≤price<2500 8%折扣 2500≤price<5000 10%折扣 一熟悉Matlab工作环境 1、熟悉Matlab的5个基本窗口 思考题： （1）变量如何声明，变量名须遵守什么规则、是否区分大小写。 答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则：?变量名必须以字母开头，只能由字母、数字或下划线组成。 ?变量名区分大小写。 ?变量名不能超过63个字符。 ?关键字不能作为变量名。 ?最好不要用特殊常量作为变量名。 （2）试说明分号、逗号、冒号的用法。 分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。 逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。 冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素。 （3）linspace（）称为“线性等分”函数，说明它的用法。 LINSPACE Linearly spaced vector. 线性等分函数 LINSPACE(X1, X2) generates a row vector of 100 linearly equally spaced points between X1 and X2. 以X1为首元素，X2为末元素平均生成100个元素的行向量。 LINSPACE(X1, X2, N) generates N points between X1 and X2. For N < 2, LINSPACE returns X2. 以X1为首元素，X2为末元素平均生成n个元素的行向量。如果n＜2，返回X2。 Class support for inputs X1,X2: float: double, single 数据类型：单精度、双精度浮点型。 （4）说明函数ones（）、zeros（）、eye（）的用法。 ones（）生成全1矩阵。 zeros（）生成全0矩阵。 eye（）生成单位矩阵。 2、Matlab的数值显示格式 1、求以下变量的值，并在MATLAB 中验证。 ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法，结果 是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 >> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB 求解下面的的方程组。 (1)????????????-=?????????????????????????----01741323151122231592127 4321x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2)???????=-++=--=-++=++56533332821 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知????????????----=1323151122231592127A (1)求矩阵A 的秩(rank) 实验二一维二维数组的创建和寻访 一、实验目的 1、掌握一维数组、二维数组创建和寻访的几种方法。 2、区别数组运算和矩阵运算的差别。 3、熟悉执行数组运算的常用数组操作函数。 4、掌握数组运算中的关系和逻辑操作及常用的关系、逻辑函数。 5、掌握“非数”、“空”数组在MA TLAB中的应用。 二、实验主要仪器与设备 装配有MA TLAB7.6软件的计算机 三、预习要求 做实验前必须认真复习第三章MATLAB的数值数组及向量化运算功能。 四、实验内容及实验步骤 1、一维数组的创建方法有哪几种？举例说明。 答：一维数组的创建方法有： ①递增/递减型一维数组的创建：冒号生成法：x=a：inc：b 线性（或对数）定点法：x=linspace(a,b,n)，x=logspace(a,b,n) ②逐个元素输入法：如x=[0.1，sin(pi/5)，-exp(-3)，-2*pi] ③运用MA TLAB函数生成法：例ones，rand等。 2、输入以下指令，并写出运行结果。本例演示：数组元素及子数组的各种标识和寻访格式；冒号的使用；end的作用。 A=zeros(2,6) %创建(2×6)的全零数组 A(:)=1:12 %赋值号左边：单下标寻访(2×6) 数组A的全部12个元素 %赋值号右边：拥有12个元素的一维数组 A(2,4) %双下标：A数组的第2行第4列元素 A(8) %单下标：数组A的第8个元素 A(: , [1,3]) %双下标：显示A的“第1列和第3列上全部行的元素” A([1, 2, 5, 6]') %单下标：把A数组第1,2,5,6个元素排成列向量 A(: , 4:end) %双下标：显示A的“从第4起到最后一列上全部行的元素” %在此end用于“列标识”，它表示“最后一列” A(2,1:2:5)=[-1, -3, -5] %把右边的3个数分别赋向A数组第2行的第1,3,5个元素位置 B=A([1, 2, 2, 2], [1, 3, 5]) %取A数组的1,3,5列的第1行元素作为B的第1行 %取A数组的1,3,5列的第2行分别作为B的第2,3,4行 L=A<3 %产生与A维数相同的“0，1”逻辑数组 A(L)=NaN %把逻辑1标识的位置上的元素赋为“非数” 运行结果： A = 0 0 0 0 0 0 0 0 0 0 0 0 数学建模 MATLAB上机实验练习题 1、给出一个系数矩阵A[234;541;132]，U=[123]，求出线性方程组的一个精确解。 2、给出两组数据x=[00.30.81.11.62.3]’y=[0.820.720.630.600.550.50]’，我们可以简单的 认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。 3、解线性方程 4、通过测量得到一组数据： 5、已知一组测量值 6、从某一个过程中通过测量得到： 分别采用多项式和指数函数进行曲线拟合。 7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近 插值法，双线性插值法，二重三次方插值法） 8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。分别使用不同的 颜色，线形和标识符。 9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额 占总销售额的比例，并分离第三种产品的切片。 X=19.322.151.6 34.270.382.4 61.482.990.8 50.554.959.1 29.436.347.0 10、对应时间矢量t，测得一组矢量y t00.30.8 1.1 1.6 2.3 y0.50.82 1.14 1.25 1.35 1.40 采用一个带有线性参数的指数函数进行拟合，y=a0+a1e-t+a2te-t，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。 11、请创建如图所示的结构数组（9分） 姓名编号指标 江明顺071023身高：176，体重：82 于越忠060134身高：168，体重：74 邓拓050839身高：182，体重：77 12、创建如图所示的元胞数组。（9分） 13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。 年份19901992199419961998200020022004200620082010 产量（万吨）75.99591.972105.711123.203131.669150.697179.323203.212226.505249.633256.344 14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。 time=[2.55.07.510.013.017.020.030.040.050.060.070.0] res=[0.290.560.771.051.361.522.002.272.813.053.253.56] 第1章MATLAB基础 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 答：在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 1.5 如何启动M文件编辑/调试器？ 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选 定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB 中有几种获得帮助的途径？ 答：（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。 （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 2.1 在MATLAB 中如何建立矩阵?? ? ? ??194375，并将其赋予变量a ？ 答：在Command Window 窗口输入操作: MATLAB期末考试题目 ------电气二班 一．选择题。 1．下列哪个变量的定义是不合法的（A） A. abcd-3 B. xyz_3 C. abcdef D. x3yz 2．在同一个图形窗口中画出三行两列的图形并指定当前可编辑图形为第3个图 的正确的 MATLAB命令是： A.subplot(3,2,1)B subplot(3,3,) C.subplot(3,2,3 ) D.subplot(1,3,2) 3.x=[1234],y=[1111],x.*y的结果是：(B) A．1234 B.出错信息 C.4 3 2 1 D.1111 4．角度x=[30 45 60]，计算其正弦函数的运算为（D） A. sin(deg2rad(x)) B. sin(x) C. sin(x) D. sin(deg2rad(x)) 5．i=2;a=2i;b=2*i;c=2*sqrt(-1);执行程序后，a,b,c的值分别是多少？（C） A.a=4,b=2,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 6. legend()命令是(A) A．在图中标注“图例” B.．在图中加上网格线C．在图中加上标题 D．打印图形 7.命令plot(x,y,′m:*′),其中字符串“m:o”中,第一个字符“m”(A) A. 表示曲线为m符号构成的线 B. 表示曲线颜色为紫色 C. 表示曲线x坐标的最大值为m变量的值 D. 表示曲线y坐标的最大值为m变量的值 8. 下列可作为matlab合法变量名的是(D) A、合计 B、123 C、@h D、xyz_2a 9. 使用语句t=0:7生成的是()个元素向量(A) A、8 B、7 C、6 D、5 10. 输入字符串或选项时，要用(D)括住 A、() B、【】 C、{} D、‘’ 11 . 清空 Matlab 工作空间内所有变量的指令是（A）。 A.clc B. cls C.clear D. clf 实验一：Matlab操作环境熟悉 一、实验目的 1．初步了解Matlab操作环境。 2．学习使用图形函数计算器命令funtool及其环境。 二、实验内容 熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算： 1．单函数运算操作。 求下列函数的符号导数 (1) y=sin(x); (2) y=(1+x)^3*(2-x); 求下列函数的符号积分 (1) y=cos(x); (2) y=1/(1+x^2); (3) y=1/sqrt(1-x^2); (4) y=(x-1)/(x+1)/(x+2); 求反函数 (1) y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2)); 代数式的化简 (1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4); (2) sin(x)^2+cos(x)^2; (3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x); 2．函数与参数的运算操作。 从y=x^2通过参数的选择去观察下列函数的图形变化 (1) y1=(x+1)^2 (2) y2=(x+2)^2 (3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6) y6=x^2/2 3．两个函数之间的操作 求和 (1) sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5 乘积 (1) exp(-x)*sin(x) Matlab上机实验答案 实验一MATLAB运算基础 1. 先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/200)/(1+exp(2)) z1 = 0.2375 >> x=[2 1+2i;-0.45 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = 0.7120 - 0.0253i 0.8968 + 0.3658i 0.2209 + 0.9343i 1.2041 - 0.0044i 2.9,,2.9, 3.0 >> a=-3.0:0.1:3.0; >> z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) (>> z33=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)可以验证z3==z33,是否都为1） z3 = Columns 1 through 5 0.7388 + 3.2020i 0.7696 + 3.2020i 0.7871 + 3.2020i 0.7920 + 3.2020i 0.7822 + 3.2020i Columns 6 through 10 0.7602 + 3.2020i 0.7254 + 3.2020i 0.6784 + 3.2020i 0.6206 + 3.2020i 0.5496 + 3.2020i Columns 11 through 20 0.4688 + 3.2020i 0.3780 + 3.2020i 0.2775 + 3.2020i 0.2080 + 3.2020i 0.0497 + 3.2020i Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w 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\b 4．设???? ??????------=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*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g) matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233，，1443678，i，，,,2357,,,,2335542，i,,,, A,1357B,,,2675342， i,,3239,,,,189543，，,,1894，， 再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4]; 1、完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。(2) 建立一个字符串向 量，删除其中的大写字母。 （1） m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 （2）>> ch='Maybe One Day' p=find(ch>='A'&ch<='Z') ch(p)=[] ch = Maybe One Day 2、自行产生一个5行5列的数组，分别得到最中间的三行三列矩阵、右下角2行2列矩阵， 奇数行矩阵、奇数列矩阵、奇数行奇数列矩阵。 >> t=rand(5)%生成矩阵 A=t(2:4,2:4)%中间三行散列矩阵 B=t(4:5,4:5)%右下角两行两列矩阵 C=t(1:2:end,:)%奇数行矩阵 D=t(:,1:2:end)%奇数列矩阵 E=t(1:2:end,1:2:end)%奇数行列矩阵 3、求方程组的根 syms x y z [X Y Z]=solve('x+4*y-3*z=2','2*x+5*y-z=11','x+6*y+z=12',x,y,z) 4、已知矩阵A=[1 2;3 4],运行指令B1=A.^(0.5), B2=A^(0.5), 可以观察到不同运算方法所得结果不同。（1）请分别写出根据B1, B2恢复原矩阵A的程序。（2）用指令检验所得的两个恢复矩阵是否相等(利用norm(…,’fro’)指令,误差矩阵F-范数，接近eps量级，认为实际相等)。 5、先运行clear,format long,rng('default'),A=rand(3,3)，然后根据A写出两个矩阵：一个对角 阵B，其相应元素由A的对角元素构成；另一个矩阵C，其对角元素全为0，而其余元素与对应的A阵元素相同（diag指令的使用）。 >> format long >> rand('twister',1) >> A=rand(3,3) A = 0.417022004702574 0.302332572631840 0.186260211377671MATLAB数学实验练习题
MATLAB实验题答案
matlab上机习题2
MATLAB实验题答案
《Matlab语言》期末考试试题样例
matlab数学实验复习题(有标准答案)
Matlab 期末考试题库(共12套卷)
河南城建学院MATLAB上机实验答案
MATLAB实验题答案
MATLAB实验题目及答案
MATLAB上机实验练习题答案
matlab基础练习题及答案讲解
MATLAB复习题.doc
Matlab实验指导书(含答案)汇总
Matlab上机实验答案 (1)
Matlab编程与应用习题和一些参考答案
matlab练习题和答案
MATLAB实验上机易错题汇总
相关主题
文本预览