matlab实验答案

郑州航空工业管理学院

实验报告册

课程名称数值计算与仿真名

学号 ********* 名

姓名刘鑫豪名

郑州航空工业管理学院工业工程实验室IE LAB OF ZHENGZHOU INSTITUTE OF AERONAUTICAL INDUSTRY MANAGEMENT

实验一

日期 2013-10-29 地点 08B204 成绩

3、利用函数产生3X4阶单位矩阵和全部元素都是4.5的4X4阶常数矩阵

在命令窗口输入

eye(3,4)

4.5*ones(4）

4、利用函数产生5X5阶随机分布的矩阵和5X5阶正态分布的随机矩阵

在命令窗口输入：

A=rand(5)

B=randn(5)

A =

0.2760 0.4984 0.7513 0.9593 0.8407 0.6797 0.9597 0.2551 0.5472 0.2543 0.6551 0.3404 0.5060 0.1386 0.8143 0.1626 0.5853 0.6991 0.1493 0.2435 0.1190 0.2238 0.8909 0.2575 0.9293

B =

-0.7648 -0.1961 -0.8045 -1.1658 -0.6669 -1.4023 1.4193 0.6966 -1.1480 0.1873

-1.4224 0.2916 0.8351 0.1049 -0.0825 0.4882 0.1978 -0.2437 0.7223 -1.9330 -0.1774 1.5877 0.2157 2.5855 -0.4390

5、练习文件读写操作

第二章上机内容

(1)

1、 gcd 函数用于求两个整数的最大公约数。先用help 命令查看该函数的用法，然后利用该函数求15和35的最大公约数。

在命令窗口输入：

help gcd

syms g

g=gcd(15,35)

help gcd

syms g

g=gcd(15,35)

gcd

2 、已知矩阵 求下列表达式： 89513236711,2032185319A B --⎡⎤⎡⎤⎢⎥⎢⎥=-=⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦

(2)

clear

clc

n=0,s=0;

while n<=100

s=s+(2*n-1)^2 n=n+1

end

s

4、求下面分数序列前20项之和：2358 1235

clear

clc

s=1;

a=1

c=1

for i=1:20

b=a

a=a+c

c=b

s=s+a/b end

s

(2)

1、1行100列的Fibonacc数组a，

a(1)=a(2)=1,a(i)=a(i-1)+a(i-2),用for循环语句来寻求该数组中第一个大于10000的元素，并指出其位置i；要求编写M文件

clear

clc

a(1)=1;

a(2)=1;

for i=3:1000

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

if a(i)>1000

break

disp(i)

end

end

a(i)

i

2、创建一个4×3阶的服从0-1均匀分布的随机矩阵A，求出A中各行的平均

值B（列向量），将B补在A的右方构成4×4阶的矩阵C，并提取C的下三角矩阵。

A=rand(4,3);

B=mean(A');

C=[A,B];

tril(C)

3、根据麦克劳林公式可以得到：e≈1+1+1/2!+1/3!+…+1/n!，编写一段

程序, 求当n=10时e的近似值。

clear

clc

s=0;

for n=0:10

s=s+1/gamma(n);

end

s

4、编写一个函数，使其能够产生如下的分段函数：

function g=ff(x)

if x<=0

g=-2*abs(x)-1;

elseif x>0&x<1

g=x^2-1;

else

g=sqrt(x)-1;

End

对应问题运行结果

第2题

A = 1 0 -1

2 4 1

-2 0 5

B =

0 1 0

2 1 3

1 1 2

⎪

⎪

⎩

⎪⎪

⎨

⎧

≥

-

<

<

-

≤

-

=

1

,1

1

,1

,1

2

)

(2

x

x

x

x

x

x

x

f

实验二一

日期 2013-10-29 地点 08B204 成绩

第2题

第3题第4题

通过对MATLAB操作对软件有了更多的理解，例如在处理图像、图形方面，以及求解高等数学相关内容方面。对于曾经在学习高等数学时比较头疼的问题，在MATLAB软件环境下都可以轻而易举地解决。

同时，在对第五章的课后习题进行操作的时候，也了解到MATLAB软件在解方程，特别是微分方程方面，对我们也有很大的帮助。

实验三人口预测与数据拟合一

日期 2013-10-29 地点 08B204 成绩