5.实验五 MATLAB程序设计二 答案

实验五 MATLAB程序设计二

1. 编程将一维数组x中的所有元素逆序排列。如N=5，原来x为： x=[ 1 3 5 7 9 ]，而逆序处理后x为x=[ 9 7 5 3 1 ]。 解：方法一

x=input('输入一维数组x=')

x=fliplr(x)

方法二

x=input('输入一维数组x=')

x=rot90(rot90(x))

方法三

x=input('输入一维数组x=')

k=length(x);

for m=1:k

a(m)=x(k-m+1);

end

x=a

2. 生成一个5行5列方阵，使其每行每列及对角线的和均为60，并求该矩阵所有元素之和。

解：方法一

a=magic(5)-1;

b=a(:)';

s=0;

for m=b

s=s+m;

end

s

方法二：

a=magic(5)-1;

b=reshape(a,1,25)

s=0;

for m=b

s=s+m;

end

s

3. 生成一个10*10矩阵，矩阵元素为[0，100]的随机整数，求其中的最大和最小元素。

解：

方法一：

A=randint(10,10,[0,100])

M1=max(A);

M2=min(A);

M11=max(M1);

M22=min(M2);

disp(['最大元素是',num2str(M11)])

disp(['最小元素是',num2str(M22)])

方法二：

A=randint(10,10,[0,100])

M1=max(A(:));

M2=min(A(:));

disp(['最大元素是',num2str(M1)])

disp(['最小元素是',num2str(M2)])

4. 输入任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。

解：

a=input('输入任意整数a=')

while a~=1

if rem(a,2)==0

a=a/2

else

a=3*a+1

end

end

a

5. 一球从100米高度自由落下，每次落地后反跳回原高度的一半,再落下，求它在第10次落地时，共经过多少米？第10次反弹多高？ 解：

x=100;

for m=1:9

a=100/2^m;

x=x+2*a;

end

disp(['第10次落地时共经过',num2str(x),'米'])

disp(['第10次落地后反弹',num2str(100/2^10),'米'])

6. 找出1－10000范围内的完数。所谓完数即该数等于其真因数的和

如：28=1+2+4+7+14

解：

方法一：

clear;clc

k=0;

for m=1:10000

s=0;

for n=1:m/2

if rem(m,n)==0

s=s+n;

end

end

if s==m

k=k+1;

a(k)=m;

end

end

disp(['10000以内的完数有：

方法二：巧用“find”函数

clear;clc

for m=1:10000

s=0;

for n=1:m/2

if rem(m,n)==0

s=s+n; ',num2str(a)])

end

end

if s==m

a(m)=m;

end

end

find(a>0);

disp(['10000以内的完数有：',num2str(ans)])

7. 求【500，600】之间能被13整除的整数个数，以及第一个和最后一个这样的数。

解：

n=0;

for k=500:600

if rem(k,13)==0

n=n+1;

a(n)=k;

end

end

disp(['【500，600】之间能被13整除的整数有',num2str(n),'个'])

disp(['【500，600】之间能被13整除的第一个整数是',num2str(a(1))]) disp(['【500，600】之间能被13整除的最后一个整数是',num2str(a(n))])