1、数 学 实 验实验班级_学生姓名_学生学号_指导老师_华南农业大学理学院应用数学系2013-4-10数学实验 20131一、MABLAB 支持下的数学实验实验 (一) MABLAB 的基本命令、初等数值运算及图形班级专业： 姓名： 学号： 日期： 1、 实验目的1.了解Matlab基本界面，熟悉基本使用方法； 2. 应用Matlab 做基本的数学运算. 3.了解 MABLAB 的多项式运算4.掌握 MABLAB 的绘图功能二、实验准备阅读P916、P1721、P2124、P24303、 实验内容写出在命令窗口中的计算步骤和运行结果。1.计算 ；1.20(lnog)8e (log(pi)+log
2、(pi)/log(10)-exp(1.2)2/81ans =0.03482.在命令窗口中键入表达式 ，并求 时 的值。2ln3xyze4,2yxz syms x y z z=x2+exp(x+y)-y*log(x)-3z =x2+exp(x+y)-y*log(x)-3 x=2;y=4;eval(z)ans =401.65623已知多项式 ， ，求：323)(345xxf 131)(23xg（1 ） 的根； (2) 在闭区间-1,2上的最小值；)( syms x y z y=inline(3*x5-x4+2*x3+x2+3)y =Inline function:y(x) = 3*x5-x4+2*
3、x3+x2+3 fzero(y,0)ans =-0.8952数学实验 20132 z=inline(1/3*x3+x2-3*x-1)z =Inline function:z(x) = 1/3*x3+x2-3*x-1 fminbnd(z,-1,2)ans =1.00004输入变量 ,在工作空间中使用 who，whos，clear 命令，并用135.,2absave 命令将变量存入”D:exe01.mat”文件。 a=5.3a =5.3000 b=1 3; 2 5b =1 32 5 whoYour variables are:a b whosName Size Bytes Classa 1x1 8
4、 double arrayb 2x2 32 double arrayGrand total is 5 elements using 40 bytes clear 清除内存中在全部变量 save D:exe015、做出下列函数的图像：（1 ） ， （分别用 plot、fplot）)2sin()(2xxy 2xx= -2:0.01:2； y=x.2.*sin(x.2-x-2); plot(x,y) z=inline(x.2.*sin(x.2-x-2); fplot(z,-2,2)数学实验 20133(2) 在同一图形窗口中，画出四幅不同图形（用 subplot 命令）：22221 34sin(),
5、 cos(), 2sin, 2cosyxyxyxyx其中 0 x=linspace(0,2*pi,100); y1=x.2.*sin(x.2-2); y2=x.2.*cos(x.2-2); y3=2*sin(2*x); y4=2*cos(2*x); subplot(2,2,1),plot(x,y1),title(y1) subplot(2,2,2),plot(x,y2),title(y2) subplot(2,2,3),plot(x,y3),title(y3) subplot(2,2,4),plot(x,y4),title(y4)四、教师评语实验( 二) MABLAB 的 MATLAB 的程序
6、结构班级专业： 姓名： 学号： 日期： 一、实验目的1. MATLAB的程序结构，计算分段函数的函数值，2. 会使用MATLAB的程序结构解决一些问题，使用 M文件来编程。二、实验准备阅读数学实验第一章1.3节三、实验内容写出在命令窗口中的计算步骤和运行结果。1在程序编辑器中输入如下程序：x=4;y=2;z=5;sum=x+y+z,product=x*y*z结果是在哪个窗口?输出的结果是 ?Command窗口sum =11product =40数学实验 201342、 在 M文件编辑器中编写函数 ,并且利用该函数计算2345y=pol(x)1+xpoly(-3),poly(2),poly(5)
7、,poly(7)打开文件编辑器输入：function y=poly(x)y=1+x+x2+x3+x4+x5;保存为poly.m poly(7)ans = 19608 poly(5)ans = 3906 poly(2)ans = 63 poly(-3)ans = -1823. 定义下列分段函数定义分段函数 ，计算f(-7), f(1),f(5), 221()lnsiexxff(9)打开文件编辑器输入：function y=f(x)if x f(-7)ans = 43 f(1)ans = 8.2305 f(5)ans = 8.0396 f(9)ans = 9.99843. 函数 ，（使用2中的po
8、ly(x)函数）计()(2sin()+()polyxelyxpol算p(2)，p(4),p(6),p(10) 打开文件编辑器输入：function z=p(x)z=exp(poly(x)+2*sin(poly(x)+poly(x)2;保存为p.m p(2)ans = 2.2938e+027 p(4)数学实验 20135ans =Inf p(6)ans =Inf p(10)ans =Inf5一个皮球从100m高度自由落下，每次落地后反弹回原高度的一半开始再次下落，试编写程序，（1）输入某次数，给出皮球在此次的反弹高度和经过的总路程；（2）求出皮球经过的总路程.解答（1）在程序编辑窗口编写以下程序
9、，并以 zsq_xhw.m为名存入相应的子目录。function zsq=zsq_xhw(n)% 这是一个求皮球反弹的高度和经过总路程的程序% s是皮球经过的总路程，zsq_xhw 是弹起的高度，n是反弹次数s=0; zsq_xhw=100;for i=1: ns=s+ zsq_xhw;zsq_xhw= zsq_xhw/2;s=s+ zsq_xhw;enddisp(弹起高度：)disp(zsq_xhw)disp(经过路程：)disp(s)在MATLAB命令窗口中调用该M文件如下： zsq_xhw(3)弹起高度：12.5000经过路程：262.5000 zsq_xhw(10)弹起高度：0.097
10、7经过路程：299.7070（2 ）在程序编辑窗口编写以下程序，并以zsq_xhw1.m为名存入相应的子目录。% 这是求从100米高空自由落下的皮球到最后静止经过总路程的程序% 其中s是总路程，h是弹起高度，n是弹起次数s=0; h=100; n=0;while hepss=s+h;数学实验 20136h=h/2;s=s+h;n=n+1;enddisp(经过的总路程：)disp(s)disp(弹起的次数：)disp(n)在MATLAB命令窗口中调用该M文件如下： zsq_xhw1经过的总路程：300.0000弹起的次数： 59四、教师评语数学实验 20137实验( 三) MABLAB 的线性代
11、数运算班级专业： 姓名： 学号： 日期： 一、实验目的1、了解 MABLAB 的矩阵运算等差数列向量的创建, 矩阵的创建, 矩阵元素的表示2 矩阵的加法, 减法, 乘法 (点乘, 直接乘), 逆矩阵, 转置运算, 矩阵的行列式，矩阵的特征值和特征向量3.向量的相关性4.线性方程组解的判别5、学习使用 MATLAB 解线性方程组的一般方法。二、实验准备阅读数学实验第一章1.3节三、实验内容写出在命令窗口中的计算步骤和运行结果。1. 已知 , 输入 zeros(size(A), 看看结果分别是什么？6532A A=2 3;5 6A =2 35 6 zeros(size(A) 这个命令就是产生和 A
12、 同型的 0 矩阵ans =0 00 02.输入矩阵 , (1) 提取矩阵 A 的第 1 行所有列的元素; (2) 提取矩阵3512AA 的第 2 列所有行的元素; (3) 提取矩阵 A 中第 2 行到第 3 行第 2 列到第 3 列的元素. A=1 -2 2;3 0 5 ;1 5 3A =1 -2 23 0 51 5 3 A(1,:)ans =1 -2 2 A(:,2)ans =-205数学实验 20138 A(2:3,2:3)ans =0 55 33. 已知 , , 在MATLAB命令窗口中建立A、B矩阵并12A1304B对其进行以下操作：(1) 求 , , ; (2) 计算下列各式: ,
13、 , , ()Rank1-A2-*.AB.TA A=1 -2 2;3 0 5 ;1 5 3A =1 -2 23 0 51 5 3 B=1 3 4;2 0 -3;2 -1 1B =1 3 42 0 -32 -1 1 det(A)ans = 13 rank(A)ans = 3 inv(A)ans=-1.9231 1.2308 -0.7692-0.3077 0.0769 0.07691.1538 -0.5385 0.4615 2*A-B,A*B,A.*B,Aans =1 -7 04 0 130 11 5ans =1 1 1213 4 1717 0 -8ans = 1 -6 86 0 -152 -5
14、3数学实验 20139ans =1 3 1-2 0 52 5 34. , . 求 A+B, 2A-B, A*B 和 A.*B, A/B 和 AB, 7146A125043BA.B. A=1 2 6;2 4 4;1 -5 7A = 1 2 62 4 41 -5 7 B=1 3 4;3 0 5;2 -1 1B =1 3 43 0 52 -1 1 A+B,2*A-B,A*B,A.*B,A/B,AB,A.Bans =2 5 105 4 93 -6 8ans =1 1 81 8 30 -9 13ans =19 -3 2022 2 320 -4 -14ans =1 6 246 0 202 5 7ans =0.0714 1.5000 -1.78571.8571 -1.0000 1.5714-3.6429 5.5000 -5.9286ans =2.0714 -2.8571 0.7857-0.1607 0.6786 0.4821-0.1250 0.7500 0.3750
