福建2003年高2第21单元第19节课_Matlab考试题库附加答案

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.wendangku.net/doc/8518790339.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试题解答

2．一单位反馈控制系统，若要求： （1）跟踪单位斜坡输入时系统的稳态误差为2。 （2）设该系统为三阶，其中一对复数闭环极点为-1±j1 求满足上述要求的开环传递函数。 利用MATALAB 完成以下操作： （1）求控制系统的单位阶跃响应、单位斜坡响应、单位加速度响应； （2）绘制控制系统的根轨迹，并判断系统的稳定性； （3）绘制控制系统的奈奎斯特图和Bode 图。 解：由于控制系统为单位反馈系统 则 H （s ）=1 由条件一可知，输入信号R(s)= 21s 时，稳态误差为： 2) (10)()(1)(==→+=s sG s s H s G s sR Ess 由条件二可知，系统为三阶且有。两个极点s1=-1+j1 ；s2=-1-j1 则可得，开环传递函数为 ) 222(1)(++=s s s s G （1）求控制系统的单位阶跃响应： G=tf([1],[1, 2, 2, 0]); Step(G) title('单位阶跃响应'); grid

(2)控制系统的单位斜坡响应：z0=0; p0=[0, -1+j,-1-j]; k0=1; [num0,den0]=zp2tf(z0,p0,k0); g0=tf(num0,den0) g=feedback(g0,1); num=g.num{1}; den=[g.den{1},0]; sys1=tf(num,den) step(sys1)

(3)控制系统的单位加速度响应：z0=0; p0=[0, -1+j,-1-j]; k0=1; [num0,den0]=zp2tf(z0,p0,k0); g0=tf(num0,den0) g=feedback(g0,1); num=g.num{1}; den=[g.den{1} 0 0]; sys2=tf(num,den) step(sys2)

matlab试卷合集附答案

一、填空题 1、求可逆矩阵A的逆矩阵的指令是 inv(A) 2、A=[0 2 3 0],B=[4 5 0 0],则X=A&B的显示结果为 0 1 0 0 3、求矩阵A的转置矩阵的指令是： A’ 4、x=[10 -20 30]; a=sign(x);显示结果为 1 -1 1 5、Notebook是Matlab的文字处理工具，启动Notebook有两种方法，即一 种是在 Word 中启动，另一种是在 Matlab 中启动。 6、Matlab7.3提供了非常强大的帮助系统，包括帮助命令、帮助 窗口、HTML格式帮助、pdf格式帮助和帮助演示等。 7、clc命令和clear命令的区别是 clc清除命令窗口所有显示内容； Clear清除所有显示的变量 . 8、命令窗口中的type命令的功能为显示文件的内容 . 9、matlab7.3提供了三种程序调试的工具，包括直接检测调试、专用调试命 令和使用 M文件编辑/调试器。 10、matlab除了最常用的主函数和子函数外，还提供了嵌套函数、私有函数、 重载函数和匿名函数。 11、在matlab7.3中，创建句柄图形对象的命令格式为____________ h_obj=funname(‘PropertyName’,PropertyValue,……)_ 12、计算5*e3*|cos175°|语句为 5*exp(3)*abs(cos(175*pi/180)). 13、创建随机矩阵的指令用rand,创建稀疏矩阵用sparse 14、为图形添加网格使用grid on，使用legend命令添加图例。 1、MATLAB常用操作界面包括命令窗口、工作空间窗口（浏览器）、 命令历史窗口、当前目录窗口、内存数组编辑器、M文件编辑/调试器、 帮助导航/浏览器、图形窗口等。 2、MATLABInf或inf表示无穷大、NaN或nan表示不是一个数、nargout表 示函数输出宗量数目。 3、MATLAB中逗号主要用作要显示计算结果的指令与其后指令的分隔； 用作输入量与输入量之间的分隔符；用作数组元素分隔符号。 4、工作空间浏览器主要用于内存变量的查阅、保存和编辑。 5、MATLAB实现将全下标转换为单下标的指令为Sub2ind、据单下标换算出全 下标的指令为Ind2sub。 6、二维数组的标识有“全下标”标识、“单下标”标识、“逻辑1”标识。 7、在一个元胞数组A中寻访第2行第3列元胞元素用A(2,3)；寻访数组第 2行第3列元胞中的内容用A{2,3}。 8、MATLAB中clf用于清除图形窗、clc用于清除指令窗中显示内容、clear 用于清除MATLAB工作空间中保存的变量。 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的值分别为和。

数学应用软件作业2Matlab作图

注意：上机作业文件夹以自己的姓名学号命名，文件夹中包括如下上机报告和Matlab程序。

5、用surf ，mesh 绘制曲面22 2z x y =+，]3,3[],3,3[-∈-∈y x 。 6、用polar 绘制阿基米德螺线r a θ=和三叶玫瑰线cos3r a θ=。（a=100） 7、在同一平面的两个窗口中分别画出心形线和马鞍面。 三. 上机方法与步骤 给出相应的问题分析及求解方法，并写出Matlab 程序。并有上机程序显示。 第1题：要在同一坐标系中分别画出四个函数的图形，可以利用plot 来将这四条曲线画在一起。 Matlab 程序： x=linspace(15,200,50); y1=(1+1./x).^x; y2=(1+1./x).^(x+1); y3=[1+1./(x+1)].^x; y4=2.7183; plot(x,y1,'r',x,y2,'g',x,y3,'b',x,y4,'co') 第2题：直接用ezplot 命令绘制函数的图形。 Matlab 程序： ezplot('（exp(x*y)-sin(x+y)）',[-3,3])

第3题：取a=1，直接用ezplot命令绘出这两条曲线。 Matlab程序： ezplot('[t-sin(t)]','[1-cos(t)]',[0,2*pi]) 第4题：利用subplot命令将这四个函数的图像画在两行两列的同一个图形上。Matlab程序： x=[0:0.1:5]; subplot(2,2,1); plot(x,sin(5*x)); grid on; title('plot-y=sin(5*x)'); subplot(2,2,2); plot(x,cos(3*x));

matlab练习题及答案

第1章MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ 答：MATLAB语言有一下特点： 1、起点高。 2、人机界面适合科技人员。 3、强大而简易的作图功能。 4、智能化程度高。 5、功能丰富，可扩展性强。 1.2 MATLAB系统由那些部分组成？ 答：MATLAB系统主要由开发环境、MATLAB数学函数库、MMATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即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 函数名”可以得到

matlab教程课后作业

【例1.3-5】图示复数i z i z 21,3421+=+=的和。 z1=4+3*i;z2=1+2i; z12=z1+z2 clf,hold on plot([0,z1,z12],'-b','LineWidth',3) plot([0,z12],'-r','LineWidth',3) plot([z1,z12],'ob','MarkerSize',8) hold off,grid on axis equal axis ([0,6,0,6]) text(3.5,2.3,'z1') text(5,4.5,'z2') text(2.5,3.5,'z12') xlabel('real') ylabel('image') shg z12 = a=-8; r_a=a^(1/3) p=[1,0,0,-a]; R=roots(p) MR=abs(R(1)); t=0:pi/20:2*pi; x=MR*sin(t); y=MR*cos(t); plot(x,y,'b:'),grid on hold on plot(R(2),'.','MarkerSize',30,'Color','r') plot(R([1,3]),'o','MarkerSize',15,'Color','b') axis([-3,3,-3,3]),axis square

hold off r_a = 1.0000 + 1.7321i R = -2.0000 1.0000 + 1.7321i 【例1.3-10】画出衰减振荡曲线t e y t 3sin 3-=，t 的取值范围是]4,0[π。 t=0:pi/50:4*pi; y=exp(-t/3).*sin(3*t); plot(t,y,'r','LineWidth',2) axis([0,4*pi,-1,1]) xlabel('t'),ylabel('y')

matlab考试题及答案

%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (iA(i)) y=A(i); end end disp('最小元素是：') disp(y) % 7、某系统由四个典型环节构成，试根据四阶龙格—库塔公式，求输出量y的阶跃响应，取仿真时间Tf=10s， %步长h=。

MATLAB期末考试试卷及其参考答案

MATLAB期末考试试卷及其参考答案 一、填空题（每空1分，20分） 1、MATLAB常用操作界面包括、工作空间窗口、 、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。 2、MATLAB中Inf或inf表示、NaN或nan表示、 nargout表示。 3、MATLAB中逗号主要用作； 用作输入量与输入量之间的分隔符；。 4、工作空间浏览器主要用于内存变量的、和。 5、MATLAB实现将全下标转换为单下标的指令为、据单下标换 算出全下标的指令为。 6、二维数组的标识有、、“逻辑1”标识。 7、在一个元胞数组A中寻访第2行第3列元胞元素用；寻访数组第2行第3列元胞中所存的内容用。 8、4、MATLAB中clf用于、clc用 于、 clear用于。 二、简答题（每题5分，共20分） 1、简述MA TLAB历史指令窗的主要作用。 2、简述空数组的功用。 3、简述MATLAB函数的基本结构。 4、简述绘制二维图形的一般步骤。 三、阅读程序并回答问题（每题4分，共28分） 1、写出下列指令运行结果。 A=zeros(2,4); A(:)=1:8;

s=[2 3 5]; A(s) Sa=[10 20 30]' A(s)=Sa 2、写出下列指令运行结果。 A=reshape(1:16,2,8) reshape(A,4,4) s=[1 3 6 8 9 11 14 16]; A(s)=0 3、写出下列指令运行结果。 A=[1,2;3,4]; B=[-1,-2;2;1]; S=3; A.*B A*B S.*A S*B 4、下面的函数主要完成什么功能？ function f=factor(n) if n<=1 f=1; else f=factor(n-1)*n; end 5、写出下列指令运行结果。 ch=‘ABc123d4e56Fg9’; subch=ch(1:5) revch=ch(end:-1:1)

MATLAB结课作业

4.10 上机操作步骤 1在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on n=1:1000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') 观察数列的散点图22,当n 趋于无穷大时,数列趋于 0 subplot(1,2,2) hold on grid on n=500:10000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') fplot('0.001',[500,10000]) fplot('-0.001',[500,10000]) axis([500,10000,-0.005,0.005]) 观察图23,当001.0=ε时,可以取N= 1000 ,当n>N 时有επε<< -2 co s n 1n . 图22 图23 2 在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on fplot('x.*x',[1,3])

观察函数图24, 当2x →时,2x y =的极限是 4 subplot(1,2,2) hold on grid on fplot('x.*x',[1.9,2.1]) fplot('4.001',[ 1.9,2.1]) fplot('3.999',[ 1.9,2.1]) axis([1.9997,2.0005,3.9989,4.0011]) % 调整显示图形的范围是该实验的重点 观察图25,当001.0=ε时, δ取 0.003 δ<-<2 0x 时,001.04<-y ? 图24 图25 3 在MatLab 的命令窗口输入: syms x limit((2.^x-log(2.^x)-1)./(1-cos(x)),x,0) 运行结果为 ans = log(2)^2 理论上用洛必达法则计算该极限: x x x cos 112ln 2lim 0x ---→= 1 4 在MatLab 的命令窗口输入如下命令序列: （1）syms x y=sqrt(x+2)*(3-x)^4/(x+1)^5 diff(y,x) ％求一阶导数 运行结果 ＝y'1/2/(x+2)^(1/2)*(3-x)^4/(x+1)^5-4*(x+2)^(1/2)*(3-x)^3/(x+1)^5-5*(x+2)^(1/2)*(3-x )^4/(x+1)^6 x=1; eval(y) ％求导数在x ＝1处的值 运行结果 1'=x y = 0.8660

matlab习题及答案

2. 用MATLAB 语句输入矩阵A 和B 3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来， 赋给B 矩阵，用magic(8)A =命令生成A 矩阵，用上述命令检验一下结果是不是正确。 4．用数值方法可以求出∑=++++++==63 63622284212i i S Λ，试不采用循环的 形式求出和式的数值解。由于数值方法是采用double 形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。 5．选择合适的步距绘制出下面的图形。 （1）)/1sin(t ，其中)1,1(-∈t ； （2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 6. 试绘制出二元函数2 2 2 2 )1(1)1(1),(y x y x y x f z +++ +-= =的三维图和三 视图 7. 试求出如下极限。 （1）x x x x 1)93(lim +∞ →； （2）1 1lim 0-+→→xy xy y x ； （3）2 2)()cos(1lim 2 2 220 0y x y x e y x y x +→→++- 8. 已知参数方程? ??-==t t t y t x sin cos cos ln ，试求出x y d d 和3 /2 2d d π=t x y 9. 假设?-=xy t t e y x f 0 d ),(2 ，试求2222 22y f y x f x f y x ??+???-?? 10. 试求出下面的极限。 （1）??????-++-+-+-∞→1)2(1 161141121lim 2222n n Λ； （2）)131211( lim 2 222π πππn n n n n n n ++++++++∞ →Λ 11. 试求出以下的曲线积分。 （1）?+l s y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=， )20(π≤≤t 。

MATLAB试卷

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.已知a=2:2:8, b=2:5，下面的运算表达式中，出错的 为（） (A) a'*b (B) a .*b (C) a*b (D) a-b 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. 下面的程序执行后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 (D) array =10. 7.下列关于脚本文件和函数文件的描述中不正确的是( ) A. 去掉函数文件第一行的定义行可转变成脚本文件; B. 函数文件可以在命令窗口直接运行; C. 脚本文件可以调用函数文件; D. 函数文件中的第一行必须以function开始;

8.对应MATLAB提供的绘制平面曲线、离散序列、三维曲线、三维网格曲线、等高线图的函数为( ) A. plot, stem, mesh, plot3, surf B. plot, stem, plot3, mesh, contour C. plot, hist, mesh, plot3, contour D. plot, hist, plot3, meshc, contour 9.设X=[0:0.5*pi:pi], Y=cos(X)且Z=sin(X)，那么Y.*Z, Y*Z’和cat(1,Y,Z)的结果分别为( ) A. [0 0 0] 0 [1 0 -1; 0 1 0] B. [0 0 0] 0 [1 0 -1 0 1 0] C. [0 0 0] [0 1 0; 0 0 0; 0 -1 0] [1 0 -1; 0 1 0] D. [0 1 0; 0 0 0; 0 -1 0] [0 0 0] [1 0 -1 0 1 0] 10.使用下列哪一条指令可以将图形窗体分割成二行三列，并且将第二行第二列的绘图区域设置为当前的绘图区域（） A.subplot(3,2,2) B.subplot(2,3,2) C.subplot(3 ,2,5) D.subplot(2,3,5)

有限元大作业matlab课程设计例子

有 限 元 大 作 业 程 序 设 计 学校：天津大学 院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\ 指导老师： 连续体平面问题的有限元程序分析 [题目]： 如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界 上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压 力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]： 由于连续平板的对称性， 只需要取其在第一象限的四分之一部分

参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。 [ 用和单元信息文件DATA.OUT。 位移模式：用用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束 方程求解：针对半带宽刚度方程的Gauss消元法 输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：

（1） ID ： ID=2时为平面应变问题 （平面问题） ，LJK_ELE(I,1)，LJK_ELE(I,2)， X(I),Y(I)分别存放节点I 的x ，y 表示第I 个作用有节点载荷的节点x,y 方向的节点载荷数值 存放节点载荷向量，解方程后该矩 （2 READ_IN ： 读入数据 BAND_K ： 形成半带宽的整体刚度矩阵 FORM_KE ： 计算单元刚度矩阵 FORM_P ： 计算节点载荷 CAL_AREA ：计算单元面积 DO_BC ： 处理边界条件 CLA_DD ： 计算单元弹性矩阵 SOLVE ： 计算节点位移 CLA_BB ： 计算单元位移……应变关系矩阵 CAL_STS ：计算单元和节点应力 （3）文件管理： 源程序文件： chengxu.for 程序需读入的数据文件：

matlab习题及答案2

MATLAB 基本运算 1．在MATLAB 中如何建立矩阵?? ?? ??194375，并将其赋予变量a ？>>a=[573;491] 2．在进行算术运算时，数组运算和矩阵运算各有什么要求？ 进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 3．数组运算和矩阵运算的运算符有什么区别？ 在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 4．计算矩阵??????????897473535与??????????638976242之和。>>a=[535;374;798]; >>b=[242;679;836]; >>a+b ans = 7 779 1413151214 5．计算??????=572396a 与?? ????=864142b 的数组乘积。>>a=[693;275]; >>b=[241;468]; >>a.*b ans = 12 36384240 6．“左除”与“右除”有什么区别？ 在通常情况下，左除x=a\b 是a*x=b 的解，右除x=b/a 是x*a=b 的解，一般情况下，a\b ≠b/a 。

7．对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。>>A=[492;764;357]; >>B=[372628]’; >>X=A\B X = -0.5118 4.0427 1.3318 8．已知：???? ??????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 >>a=[123;456;789]; >>a.^2 ans = 1 4916 253649 6481 >>a^2 ans = 30 364266 81961021261509．[]7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。 相当于a=[11011]。 10．在sin(x )运算中，x 是角度还是弧度？ 在sin(x)运算中，x 是弧度，MATLAB 规定所有的三角函数运算都是按弧度进行运算。

MATLAB期末上机试题带答案

MATLAB 期末上机考试试题带答案版 姓名： 学号： 成绩： 1.请实现下图： 50 100 150 200 250 -1 -0.8-0.6-0.4-0.200.20.40.60.81x y y=sin(x) x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图： x=linspace(0,2*pi,100); y1=sin(x); subplot(2,2,1) plot(x,y1,'k--') grid on xlabel('x') ylabel('y') title('sin(x)') legend('y=sin(x)') y2=cos(x);

subplot(2,2,2) plot(x,y2,'r--') grid on xlabel('x') ylabel('y') title('cos(x)') legend('y=cos(x)') y3=tan(x); subplot(2,2,3) plot(x,y3,'k-') grid on xlabel('x') ylabel('y') title('tan(x)') legend('y=tan(x)') y4=cot(x); subplot(2,2,4) plot(x,y4) grid on xlabel('x') ylabel('y') title('cot(x)') legend('y=cot(x)') 3.解方程组： a=[3 2 1;1 -1 3;2 4 -4]； b=[7;6;-2] ； x=a\b 4.请实现下图： 2 4 6 8 10 12 14 -1-0.8-0.6-0.4-0.200.20.40.60.81x y '自己的姓名' sin(x)sin(2*x) x=linspace(0,4*pi,1000); y1=sin(x);

MATLAB期末上机试题带答案

MATLAB期末上机试题带答案

MATLAB 期末上机考试试题带答案版 姓名： 学号： 成绩： 1.请实现下图： x y y=sin(x) x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图：

x=linspace(0,2*pi,100); y1=sin(x); subplot(2,2,1) plot(x,y1,'k--') grid on xlabel('x') ylabel('y') title('sin(x)') legend('y=sin(x)') y2=cos(x); subplot(2,2,2) plot(x,y2,'r--') grid on xlabel('x') ylabel('y') title('cos(x)') legend('y=cos(x)') y3=tan(x); subplot(2,2,3) plot(x,y3,'k-') grid on xlabel('x')

ylabel('y') title('tan(x)') legend('y=tan(x)') y4=cot(x); subplot(2,2,4) plot(x,y4) grid on xlabel('x') ylabel('y') title('cot(x)') legend('y=cot(x)') 3.解方程组： a=[3 2 1;1 -1 3;2 4 -4]；b=[7;6;-2] ； x=a\b 4.请实现下图：

matlab结课论文

精选文库 山西大同大学matlab课程结课作业MATLAB程序应用 姓名： 课程序号： 2 班级： 学号： 2013年12月

1.实验内容：已知!123n n =?????L ，编写一个程序求满足100!10n ≤的 最大的n 值以及此时!n 的值。 function n n=2;m=1; while m<=10^100 m=m.*n;n=n+1; end m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 2.设)15113111191715131 1(22Λ+--++--+=π，试根据公式编出计算pi 的Mat lab 主程序文件，pi 的精度为0.00001。 程序： k=0;n=1;b=0;a=0; while abs((pi-a))>0.00001 a=2*sqrt(2)*k; k=( bcos( *pi/2)+sin(b*pi/2))/n+k; n=n+2; b=b+1; end a 输出a=3.141602572083633 ; a-pi= 9.918493839577991e-006 3.有两个矩阵A 和B 如下：????????????---=771175420132861-1A ，????????????------=0162310013125673B ， 将A 中所有等于-1的元素改为-2，将B 中所有小于0的元素改为1，然后将B 中等于0的元素的值改为A 的相应位置元素的值。请用Matlab 函数文件实现上述运算。 clear;

clc; A=[1 -1 6 8;2 3 -1 0;-2 4 5 7;1 -1 7 7]; B=[-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0]; C=A;A(A==-1)=-2;U=A; D=B;B(B<0)=1;V=B; A=C;B=D;[i,j]=find(B==0);A(i,j)=0;W=A; A=C;B=D; A,B,W,U,V %用函数文件实现矩阵中元素的变换。 %A、B为输入变量。 %U、V、W分别存放A、B中间变换结果。 ; 4.用matlab主程序文件产生动画：呈现一小圆（半径为1）在一大圆（半径为3）的圆周外部滚动的动画，要求连续滚动20周。 clea close;clc;r; axis([-6 6 -6 6],'equal','manual');hold on; ezplot('x^2+y^2-9'); h=ezplot('x^2+y^2-1'); x=get(h,'xdata'); y=get(h,'ydata'); for t=1:7200 set(h,'xdata',x+4*cosd(t),'ydata',y+4*sind(t)); drawnow; end

matlab练习题和答案

matlab练习题和答案 控制系统仿真实验 Matlab 部分实验结果 目录 实验一MATLAB基本操作..................................................................... .......................1 实验二Matlab 编程..................................................................... ...............................5 实验三Matlab 底层图形控制..................................................................... .....................6 实验四控制系统古典分析...................................................................... .......................12 实验五控制系统现代分析..................................................................... . (15) 实验六PID 控制器的设计..................................................................... ......................19 实验七系统状态空间设计...................................................................... .......................23 实验九直流双闭环调速系统仿真..................................................................... . (25) 实验一MATLAB基本操作

MATLAB大作业

M A T L A B大作业作业要求： （1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电 子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。 （2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果 分析、本课程学习体会、列出主要的参考文献。打印稿不要求MATLAB程序，但电 子稿要包含MATLAB程序。 （3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话 则长，无话则短。不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。 （4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。 （5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。 第一类：绘制图形。（B级） 问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契 螺旋线，如图所示。 问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三

角形，然后去掉中间的那一个小三角形。接下来对其余三个小三角形重复上述操作，如图所示。 问题三：其他分形曲线或图形。分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。这方面的资料很多（如），请分析构图原理并用MATLAB 实现。 问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。 点数 1 2 3 4 5 6 出现次数166 150 164 162 184 174 问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。 提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。 第二类：插值与拟合。（B级） 问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。注意，整个实验过程中从未减速。在一组时间点上测得汽车的速度如表所示。 （1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。 （2）绘制插值图形并标注样本点。

MATLAB试卷+答案

学习中心/函授站_ 姓名学号 西安电子科技大学网络与继续教育学院 2015学年上学期 《MATLAB与系统仿真》期末考试试题 （综合大作业） 考试说明： 1、大作业于2015年4月3日公布，2015年5月9日前在线提交； 2、考试必须独立完成，如发现抄袭、雷同、拷贝均按零分计。 3、程序设计题（三（8，10））要求写出完整的程序代码，并在matlab软件环境调试并运行通过，连同运行结果一并附上。 一、填空题（1’ ×25=25’） 1、Matlab的全称为MATrix LABoratory 。 2、在Matlab编辑器中运行程序的快捷键是：F5 。 3、Matlab的工作界面主要由以下五个部分组成，它们分别是：Command window(指令窗口、Current directory（当前路径浏览器）、Workspace(工作台窗口) 、Command history （指令历史窗口）和Launch pad（发布平台）。 4、在Matlab中inf表示：无穷大；clc表示：清除当前Command window区域内容；more表示：显示命令窗口输出一个整个萤幕的时间；who表示：返回当前工作空间中所有变量名；whos表示：返回当前工作空间中所有变量，以及它们的名字、尺寸、所占字节数、属性等信息。 5、在Matlab命令窗口中运行命令Simulink 可以打开Simulink模块库浏览器窗口。 6、求矩阵行列式的函数：det ；求矩阵特征值和特征向量的函数eig 。 7、Matlab预定义变量ans表示：当前未定义变量和变量名的答案；eps表示： 浮点相对精度（其值是从1.0到下一个最大浮点数之间的差值）；nargin表示：函 9、建立结构数组或转换结构数组的函数为：struct

Matlab程序设计(2016大作业)

Matlab程序设计 课程大作业 题目名称：_________________________________ 班级：_________________________________ 姓名：_________________________________ 学号：_________________________________ 课程教师：温海骏 学期：2015-2016学年第2学期 完成时间： MATLAB优化应用 §1 线性规划模型 一、线性规划问题： 问题1：生产计划问题 假设某厂计划生产甲、乙两种产品，现库存主要材料有A类3600公斤，B类2000公斤，C类3000公斤。每件甲产品需用材料A类9公斤，B类4公斤，C类3公斤。每件乙产品，需用材料A类4公斤，B类5公斤，C类10公斤。甲单位产品的利润70元，乙单位产品的利润120元。问如何安排生产，才能使该厂所获的利润最大。 问题2：投资问题 某公司有一批资金用于4个工程项目的投资，其投资各项目时所得的净收益(投入资金百分比)如下表：工程项目收益表 工程项目 A B C D 收益(%) 15 10

12 由于某种原因，决定用于项目A的投资不大于其他各项投资之和而用于项目B和C的投资要大于项目D的投资。试确定该公司收益最大的投资分配方案。 问题3：运输问题 有A、B、C三个食品加工厂，负责供给甲、乙、丙、丁四个市场。三个厂每天生产食品箱数上限如下表： 工厂 A B C 生产数 60 40 50 四个市场每天的需求量如下表： 市场 甲 乙 丙 丁 需求量 20 35 33 34 从各厂运到各市场的运输费(元/每箱)由下表给出： 收点 发点 市场 甲 乙 丙 丁 工 厂 A 2 1 3 2 B