matlab期末考试试题答案

山东科技大学2014—2015学年第一学期

《计算机仿真》试卷（期中卷）

班级姓名学号

1.前寻已输入的命令的控制键为_方向键上键_，删除当前行全部内容的键为____Esc_____。

2.变量名是以字母开头，后接任意字母、数字或下划线的字符序列。P33

3.abs函数可以求绝对值、复数模、。

4.标点符号___；____可以使命令行不显示运算结果，_%___用来表示该行为注释行。

5.提取分子分母的指令为__numden_______

6.插值运算的指令为__interp_____，曲线拟合的指令为___polyfit___ P67

二．程序设计（共86分）注意：把以下试题的答案保存在电脑桌面上的文件夹里面，文件夹的名称为：班级学号姓名，如：自动化09-1班01号张三。

1. （12分）已知

352

a=918

467

??

??

??

??

??

，

（1）取a矩阵的1、3行，2~3列的数据赋予b1。

（2）建立与a矩阵同大小的全1矩阵b2。

（3）提取a矩阵的对角线元素给b4。

程序以m命令文件格式保存，名称为自己姓名的拼音字母1.m，如你的姓名为张三，则保存文件名称为zhangsan1.m，以下雷同。

a=[3 5 2;9 1 8;4 6 7]

b1=a([1 3],[2 3])

b2=ones(3,3)

b4=diag(a)

2. （12分）已知

152

a=346

798

??

??

??

??

??

，（1）提取a矩阵的上三角矩阵赋值为d1。

（2）把a矩阵的左右翻转赋值为d2。

（3）求a矩阵的秩d4。程序保存为姓名拼音字母2.m a=[1 5 2;3 4 6;7 9 8]

d1=triu(a)

d2=fliplr(a)

d4=rank(a)

3.（8分）数组的元素满足规则：a=a+a

k+2k k+1，)

,2,1

(

k；且a=1,a=2,

12

现要求该

数组中第一个大于1000的元素。用while语句或for循环完成写出程序，且求出k和a

k

的取值。程序保存为姓名拼音字母3.

a(1)=1;

a(2)=2;

i=2;

while a(i)<=10000

a(i+1)=a(i-1)+a(i);

i=i+1;

end

i

a(i)

4. （8分）在0~2π范围内产生100个数据点，用图形保持功能在同一坐标内以不

同颜色、线形和数据点绘制曲线

2

-0.05x

y1=2sin(3x),y2=3cos(2x),y3=3e sin(4x),给坐

标轴上加标注，给整个图形加上标题（自己的汉字姓名），在图形窗口添加文字说明，图例等，使图形做的尽量漂亮。

保存程序指令，不保存图形，程序保存为姓名拼音字母4.m

x=0:(2*pi)/100:2*pi;

y1=2*sin(3*x);

y2=3*cos(2*x);

y3=3*exp(-0.05*x.^2).*sin(4*x);

plot(x,y1,'r',x,y2,'b',x,y3,'g');

ylabel('y');

xlabel('x');

title('张三');

legend('y1','y2','y3',4)

5. （8分）在0~2π范围内间隔0.5产生数据点，用子图在左上部2、右上部1以条形图、(火柴）杆图绘制曲线0.5x

y=e sin(2x)，左下部3以饼图绘制[30 23 40 46 37]，右下部4绘制复向量图[2+i,5+5i,3-4i,-2-2i]。给每个图形加上标题（自己的汉字姓名）。

保存程序指令，不保存图形，程序保存为姓名拼音字母5.m

y=exp(0.5*x).*sin(2*x);

subplot(221);

bar(x,y);

title('张三');

subplot(222);

stem(x,y,'filled');

title('张三');

c=[30 23 40 46 37];

subplot(223);

explode=[0 0 0 1 0];

pie(c,explode);

title('张三');

d=[2+i,5+5*i,3-4*i,-2-2*i];

subplot(224);

compass(d);

title('张三');

6. （12分）53-210871-5f=-411490-16-6????????????

（1）求f 矩阵的各行最大值给f1。（2）求f 矩阵的各列元素累计和给f2。

（3）计算322x +3x +5x+6=0的全部根g2。

程序保存为姓名拼音字母6.m

f=[5 3 -2 10;8 7 1 -5;-4 11 4 9;0 -1 6 -6] f1=max(f')

f2=cumsum(f)

g2=[2 3 5 6]

roots(g2)

7. （12分）（1）23-51lim x x s x

+→= （2）3

3sin()s x dx π

=?

（3）111141223341920

s =++++???? 程序为姓名7.m

s1=sym('(x*x-5)/(x)')

limit(s1,'x',3,'right')

s2=sym('sin(x)')

int(s2,0,pi/3)

syms x k;

s4=symsum(1/(k*(k+1)),1,19)

8.（8分）将传递函数

4s+8

32

s+6s+5s

转换为零极点模型，求其单位阶跃响应。程序

保存为姓名拼音字母8.m

num=[4,8];

den=[1,6,5,0];

G1=tf(num,den)

G2=zpk(sys)

9.（5分）二阶仿真系统的图形如下图，在0.5s时刻输入单位阶跃信号，调节示波器观察清晰全部图形。加上标题（汉字姓名）保存图形为姓名拼音字母9.mdl文件。

设y=cos[0.5+((3sinx)/(1+x^2))]把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线；

解：>> x=linspace(0,2*pi,101);

y=cos(0.5+3.*sin(x)./(1+x.*x));

plot(x,y)

(5) 向量]a ,,a ,[a A n 11 . 写一段程序, 找出A 中的最小元素（10分） 答：A= input ('请输入一个向量')

[m,n]=sizeA

min =A(1,n);

for i=1:n

if A(1,i)

min=A(1,i)

end

end

打印出所有的水仙花数。所谓“水仙花数”，是指一个三位数，其各位数字立方之和等于该数本身。

for k = 100:999

a = fix(k/100);

b = rem(fix(k/10),10);

c = rem(k,10);

if a.^3 + b.^3 + c.^3 == k

fprintf( “%u,\t\t ”, k )

end

end

4 编写M 函数文件，利用for 循环或while 循环完成计算函数y=sum(n)= 的任务，并利用该函数计算 n=20时的和（将总程序保存为？？.m 文件）。

function sum=text4(n);

sum=0;

m=1;

while m<=n;

sum=sum+m^m

m=m+1;

end

当n=20时 用上述程序求得的和为： 1.0688e+026 G(S)=)()(S U S Y =5

1437105232+++++S S S S S MATLAB 程序代码如下：

b=[5,10,7];

z=roots(b)

a=[3,14,1,5];

p=roots(a)

k=3;

sys=zpk(z,p,k)