当前位置：文档之家› 数学实验1

数学实验1

重庆大学

学生实验报告

实验课程名称数学实验

开课实验室DS1422

学院计算机学院年级2013专业班物联网2班学生姓名卢锦晔学号20135195

开课时间2014 至2015 学年第2学期

总成绩

教师签名

数学与统计学院制

开课学院、实验室:DS1422实验时间：2015 年3月30 日

课程名称数学实验实验

项目

名

称

MATLAB软件入门实验项

目类型

验证演示综合设计其他

指导教师肖剑成

绩

实验目的

[1] 熟悉MATLAB 软件的用户环境； [2] 了解MATLAB 软件的一般目的命令； [3] 掌握MATLAB 数组操作与运算函数； [4] 掌握MATLAB 软件的基本绘图命令；

[5] 掌握MATLAB 语言的几种循环、条件和开关选择结构。

通过该实验的学习，使学生能灵活应用MATLAB 软件解决一些简单问题，能借助MATLAB 软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜想，发现进而证实其中的规律。

实验内容

1．MATLAB 软件的数组操作及运算练习； 2．直接使用MATLAB 软件进行作图练习； 3．用MATLAB 语言编写命令M-文件和函数M-文件。

基础实验

一、问题重述

1．设有分块矩阵??

???=????223

22333S O

R E

A ，其中E,R,O,S 分别为单

位阵、随机阵、零阵和对角阵，试通过数值计算验证

????

+=22S 0RS R E A 。

2．某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。

表1.1

货号 1 2 3 4 5 6 7 8 9 单件进价 7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30 单件售价 11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50 销量

568 1205 753 580 395 2104 1538 810 694

3．建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为153=13+53+33。

4．编写函数M-文件sq.m ：用迭代法求a

=x 的值。

求平方根的迭代公式为

)a (211n

n n x x x +=

迭代的终止条件为前后两次求出的x 的差的绝对值小于10-5。

5. 近景图将x 的取值范围局限于较小的区间内可以画出函数的近景图,用于显示函数的局部特性。

局部放大在绘图时,把x 的范围逐渐缩小,

可把函数的细节部分展现的很清楚.特别是观察极限问题时,这种方法比较便利.

远景图函数的远景图,是把x 的范围取得

比较大,使我们能够在大范围内观察函数图像.当研究x 趋向于∞时,这种方法给我们带来方便.

1）比较函数

3)(,)(,)(x x h x x x g x x f =+== 在x →0时函数

的性态。观察到什么现象？从观察到的现象，反映了什么结论。

2）在日常生活中我们有这样的经验：与幂函数相比，指数函数是急脾气，对数函数是慢性子。这就是说，

当x →∞时，再小的指数函数也比幂函数变化快，再大的对数函数也比幂函数变化慢。当x →∞时,比较

x y =与

y 1.1= 的大小.当x →∞时,比较

001

.0x y =与

y lg 1000= 的大

小.

）

在

同

一

个

坐

标

下

作

出

y 1=e x ,y 2=1+x,y 3=1+x+(1/2)x 2,y 4=

1+x+(1/2)x 2+(1/6)x 3这四条曲线的图形，要求在图

上加各种标注，观察到什么现象？发现有什么规律？同时用subplot 分别在不同的坐标系下作出这四条曲线，为每幅图形加上标题。

4)．作出下列曲面的3维图形，

)

sin(22y x z +π=；

5). 作出函数y=x 4-4x 3+3x+5 （x ∈[0,6]）的图形，用小红点标出其在[0,6]之间的最小值点，并在最小值点附近标出该最小值点的坐标值；

二、实验过程（一般应包括实验原理或问题分析，变量说明、程序、调试情况记录、图表等，实验结果及分析）

?源程序：

E=eye(3,3)

R=rand(3,2)

O=zeros(2,3)

a=[3 2]

S=diag(a)

A=[E R;O S]

B=[E R+R*S;O S^2]

B==A^2

disp('right')

else

disp('wrong’)

end

?运行结果：

E =

1 0 0

0 1 0

0 0 1

R =

0.2785 0.9649 0.5469 0.1576 0.9575 0.9706

O =

0 0 0 0 0 0

a =

3 2

S =

3 0

0 2

A =

1.0000 0 0 0.2785

0.9649

0 1.0000 0 0.5469

0.1576

0 0 1.0000 0.9575

0.9706

0 0 0 3.0000 0

0 0 0 0

2.0000

B =

1.0000 0 0 1.1140

2.8947

0 1.0000 0 2.1875

0.4728

0 0 1.0000 3.8300

2.9118

0 0 0 9.0000 0

0 0 0 0

4.0000

Right

分析：在该题中应该注意’*’与’.*’的区别，避免不必要的错误。题中涉及的矩阵函数有eye、

Rand、zeros、diag等，还有矩阵合并和结果输出的有关知识。

?源程序：

a=1:9

b=[7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30]

c=[11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50]

d=[568 1205 753 580 395 2104 1538 810 694]

x=c.*d

profit=(c-b).*d

f=sum(x)

g=sum(profit)

[income,number]=sort(x)

?运行结果：

?分析：此题最重要的部分是sort函数的应用，实现数组元素的排序和输出其在数组中对应的位置。另外，从中我还学会了利用数组对一系列的数据进行运算处理。

?计算程序

for a=1:1:9

for b=0:1:9

for c=0:1:9

A=a*100+b*10+c;

if A==a^3+b^3+c^3

disp(A)

end

?最终结果：

153 370 371 407

?分析：此题使我进一步熟悉了循环语句for、while 等的基本格式和应用方法，值得注意的是每一个循环语句或条件语句的出现，必定要对应的出现一个end，这一点很容易忽视。

?源程序：

function f=sq(a)

if a>=0

x=eps;

y=1/2*(x+a/x);

while (abs(x-y)>=10^(-5))

x=y

y=1/2*(x+a/x)

end

f=x

else

disp('their exist errors')

end

?运行结果

调用sq(16)，输出结果是4.0000

?分析：此题要求掌握利用循环语句进行迭代，实现开方方程的求解，另外x的取值有一定的讲究，在该题题设的情况下，很显然x>=0，但从迭代公式的形式可以得到，x不能等于0，这就要求x从一个较为接近0的数开始取，于是就得引用eps（计算机能辨别的

最小值）。

（1）

?源程序：

x=-1:0.0001:1;

y1=x;

y2=x.^3;

y3=y1+y2;

plot(x,y1,x,y2,x,y3)

?运行结果

分析

当x→0时，f(x)与g(x)很接近，而h(x)与前两个函数都不接近。

(2)

?源程序：

x=linspace(5000,8000,500);

y1=x.^10;

y2=1.1.^x;

Subplot(1,2,1),plot(x,y1),xlabel('x');ylabe l('y)'); grid;

title('y=x^1^0');

Subplot(1,2,2),plot(x,y2),xlabel('x');ylabe l('y)'); grid;

title('y=1.1^x');

?运行结果

分析：从上图可以看出来指数函数变化快

?源程序：

x=linspace(5000,8000,500);

y1=x.^0.001;

y2=1000.*log(x);

Subplot(1,2,1),plot(x,y1),xlabel('x');ylabe l('y)'); grid;

title('y=x^0.001');

Subplot(1,2,2),plot(x,y2),xlabel('x');ylabe l('y)'); grid;

title('y=1000.*log(x)');

?运行结果

分析：由以上函数图形可知对数函数变化比幂函数慢。

(3)

?源程序：

x=linspace(0,2.50);

y1=exp(x);

y2=1+x;

y3=1+x+0.5.*x.^2;

y4=1+x+0.5.*x.^2+1./6.*x.^3;

plot(x,y1,'b.'),gtext('y1=exp(x)');

hold on, plot(x,y2,'y-'),gtext('y2=1+x');

plot(x,y3,'g:'),gtext('y3=1+x+0.5.*x.^2');

plot(x,y4,'m--'),gtext('y4=1+x+0.5.*x^2+1./ 6.*x.^3');

hold off

?作图结果

(4)

?作图程序

x=-10:0.6:10;

y=x;

[X,Y]=meshgrid(x,y);

Z=sin(pi.*sqrt(X.^2+Y.^2));

surf(X,Y,Z)

xlabel('x')

ylabel('y')

zlabel('z')

title('z=sin(pi*sqrt(x.^2+y.^2)')

shading flat

?作图结果

?分析：在该题中，需要用到三维作图相关函数meshgrid、mesh、surf，注释函数xlabel、

ylabel、zlabel、title和开方函数sqrt，最后还涉及到了幂函数等。通过这道题，我充分运用并巩固了作图方面的基础知识。

(5)

?作图程序

x=0:0.1:6;

y=x.^4-4.*x.^3+3.*x+5;

plot(x,y)

[y1,x1]=min(y);

hold on

plot(x(x1),y1,'r.','MarkerSize',20)

hold off

a=['x=',num2str(x(x1))];

b=['y=',num2str(y1)];

min=char(a,b);

text(x(x1),y1+50,min)

?作图结果

6．自由发挥：自己提出问题，实验探索，广泛联想，发现规律，大胆猜想。比如函数cos(1/x)在x=0附近的振荡现象，有无规律可寻？

?问题

用两种方法在同一个坐标下作出y1=x2,y2=x3,y3=x4 y4=x5这四条曲线的图形，并要求用两种方法在图上加各种标注。

?程序（1）

x=linspace(-2,2,40);

y1=x.^2;

y2=x.^3;

y3=x.^4;

y4=x.^5;

plot(x,y1,'m')

hold on,plot(x,y2,'r');plot(x,y3,'g');plot(x,y4,'k');hold off

text(-0.6,0.4,'\leftarrowy1=x^2')

text(-0.9,0.2,'y2=x^3\rightarrow')

text(-0.9,-0.1,'y3=x^4\rightarrow')

text(-0.8,-0.6,'\leftarrowy4=x^5')

xlabel('x')

ylabel('y')

title('í?ê?')

?程序（2）

x=linspace(-1,1,20)

y1=x.^2;y2=x.^3;y3=x.^4 ;y4=x.^5;

数学实验答案-1

1.（1） [1 2 3 4;0 2 -1 1;1 -1 2 5;]+(1/2).*([2 1 4 10;0 -1 2 0;0 2 3 -2]) 2. A=[3 0 1;-1 2 1;3 4 2],B=[1 0 2;-1 1 1;2 1 1] X=(B+2*A)/2 3. A=[-4 -2 0 2 4;-3 -1 1 3 5] abs(A)>3 % 4. A=[-2 3 2 4;1 -2 3 2;3 2 3 4;0 4 -2 5] det(A),eig(A),rank(A),inv(A) 求计算机高手用matlab解决。 >> A=[-2,3,2,4;1,-2,3,2;3,2,3,4;0,4,-2,5] 求|A| >> abs(A) ans = （ 2 3 2 4 1 2 3 2 3 2 3 4 0 4 2 5 求r（A） >> rank(A) ans =

4 求A-1 《 >> A-1 ans = -3 2 1 3 0 -3 2 1 2 1 2 3 -1 3 -3 4 求特征值、特征向量 >> [V,D]=eig(A) %返回矩阵A的特征值矩阵D 与特征向量矩阵V , V = - + + - - + - + - + - + D = { + 0 0 0 0 - 0 0 0 0 + 0 0 0 0 - 将A的第2行与第3列联成一行赋给b >> b=[A(2,:),A(:,3)'] b = 《 1 - 2 3 2 2 3 3 -2

1． a=round(unifrnd(1,100)) i=7; while i>=0 i=i-1; b=input('请输入一个介于0到100的数字：'); if b==a ￥ disp('You won!'); break; else if b>a disp('High'); else if b