Matlab数学建模
➢M文件建立方法:1. 在Matlab中,点:File->New->M-file; 2. 在编辑窗口中输入程序内容;3. 点:File->Save,存盘,M 文件名必须与函数名一致。 ➢Matlab的应用程序也以M文件保存。
11
二、数组
➢1、创建简单的数组 ➢x=[a b c d e f ] 创建包含指定元素的行向量 ➢x=first:last ➢ 创建从first开始,加1计数,到last结束的行向量 ➢x=first:increment:last ➢ 创建从first开始,加increment计数,last结束的行向量 ➢x=linspace(first,last,n) ➢ 创建从first开始,到last结束,有n个元素的行向量 ➢x=logspace(first,last,n) ➢ 创建从开始,到结束,有n个元素的对数分隔行向量.
18
三、矩阵
➢(8)删除A的第j1~j2列,构成新矩阵:A(:, j1:j2)=[ ] ➢(9)将矩阵A和B拼接成新矩阵:[A B];[A;B] ➢3、矩阵的运算 ➢(1)标量-矩阵运算 ➢ 同标量-数组运算。 ➢(2)矩阵-矩阵运算 ➢ [1] 元素对元素的运算,同数组-数组运算。 ➢ [2]矩阵运算: ➢ 矩阵加法:A+B ➢ 矩阵乘法:A*B ➢ 方阵的行列式:det(A) ➢ 方阵的逆:inv(A) ➢ 方阵的特征值与特征向量:[V,D]=eig[A]
➢
{commands2}
➢
else if (expression3)
➢
{commands3}
➢
else if ……
➢
…………………………………
➢
else
➢
{commands}
➢
end
➢
end
➢
end
➢
……
➢
end
25
六、MATLAB作图——二维图形
➢1.曲线图
➢ Matlab作图是通过描点、连线来实现的,故在画一个曲线图 形之前,必须先取得该图形上的一系列的点的坐标(即横坐标 和纵坐标),然后将该点集的坐标传给Matlab函数画图.
➢例 在[0,2*pi]用红线画sin(x),用绿圈画cos(x).
➢解:x=linspace(0,2*pi,30);
➢ y=sin(x);
➢ z=cos(x);
➢ plot(x,y,'r',x,z,’g0') 1
0.8
0.6
0.4
0.2
0 -0.2
-0.4
-0.6
27
-0.8
六、MATLAB作图——二维图形
22
五、控制流
➢ ➢
例 对n=1,2,…,10,求xn= 2、While循环
sin
n
10
的值
➢与for循环以固定次数求一组命令相反,while循环以不定的次 数求一组语句的值.
➢
while expression
➢
{commands}
➢
end
➢只要在表达式(expression)里的所有元素为真,就执行while和 end语句之间的命令串{commands}.
15
二、数组
➢设:a=[a1,a2,…,an], b=[b1,b2,…,bn] ➢则:a+b= [a1+b1,a2+b2,…,an+bn] ➢ a.*b= [a1*b1,a2*b2,…,an*bn] ➢ a./b= [a1/b1,a2/b2,…,an/bn] ➢ a.\b=[b1/a1,b2/a2,…,bn/an] ➢ a.^b=[a1^b1,a2^b2,…,an^bn]
➢特殊矩阵的建立:
➢a=[ ] 产生一个空矩阵,当对一项操作无结果时,返回 空矩阵,空矩阵的大小为零.
17
三、矩阵
➢b=zeros(m,n) 产生一个m行、n列的零矩阵 ➢c=ones(m,n) 产生一个m行、n列的元素全为1的矩阵 ➢d=eye(m,n) 产生一个m行、n列的单位矩阵 ➢2、矩阵中元素的操作 ➢(1)矩阵A的第r行:A(r,:) ➢(2)矩阵A的第r列:A(:,r) ➢(3)依次提取矩阵A每1列,将A拉伸为一个列向量:A(:) ➢(4)取矩阵A的第i1~i2行、第j1~j2列构成新矩阵:A(i1:i2, j1:j2) ➢(5)以逆序提取A的第i1~i2行,构成新矩阵:A(i2:-1:i1,:) ➢(6)以逆序提取A的第j1~j2列,构成新矩阵:A(:, j2:-1:j1 ) ➢(7)删除A的第i1~i2行,构成新矩阵:A(i1:i2,:)=[ ]
例
设f
( x)
x2
1
x
1 ,
求f
(2),
f
(1)
2x x 1
➢ 先建立M文件fun1.m定义函数f(x),再在Matlab命令窗口 输入fun1(2),fun1(-1)即可。
24
五、控制流
➢2) 有三个或更多的选择的一般形式是:
➢
if (expression1)
➢
{commands1}
➢
else if (expression2)
➢例 设银行年利率为11.25%。将10000元钱存入银行,问多长 时间会连本带利翻一番?
23
五、控制流
➢3、If-Else-End结构
➢(1)有一个选择的一般形式是:
➢
if expression
➢
{commands}
➢
end
➢ 如果在表达式(expression)里的所有元素为真,就执行if 和end语句之间的命令串{commands}.
7
➢特殊变量表
特殊变量
ans pi eps flops inf NaN i,j nargin nargout realmin realmax
取值 用于结果的缺省变量名 圆周率 计算机的最小数,当和 1 相加就产生一个比 1 大的数 浮点运算数 无穷大,如 1/0 不定量,如 0/0
i=j= 1 所用函数的输入变量数目 所用函数的输出变量数目 最小可用正实数 最大可用正实数
14
二、数组
➢设:a=[a1,a2,…,an], c=标量 ➢则:a+c=[a1+c,a2+c,…,an+c] ➢ a.*c=[a1*c,a2*c,…,an*c] ➢ a./c= [a1/c,a2/c,…,an/c](右除) ➢ a.\c= [c/a1,c/a2,…,c/an] (左除) ➢ a.^c= [a1^c,a2^c,…,an^c] ➢ c.^a= [c^a1,c^a2,…,c^an] ➢(2)数组-数组运算 ➢ 当两个数组有相同维数时,加、减、乘、除、幂运算可按 元素对元素方式进行的,不同大小或维数的数组是不能进行运 算的.
➢2.符号函数(显函数、隐函数和参数方程)画图 ➢(1) ezplot ➢ ezplot(‘f(x)’,[a,b]) 表示在a<x<b绘制显函数f=f(x)的函数图 ➢ ezplot(‘f(x,y)’,[xmin,xmax,ymin,ymax]) ➢ 表示在xmin<x<xmax和 ymin<y<ymax绘制隐函数f(x,y)=0的函数图 ➢ezplot(‘x(t)’,’y(t)’,[tmin,tmax]) ➢ 表示在tmin<t<tmax绘制参数方程x=x(t),y=y(t)的函数图 ➢例 在[0,pi]上画y=cos(x)的图形 ➢解 输入命令 ➢ ezplot(‘sin(x)’,[0,pi])
12
二、数组
➢ 2、 数组元素的访问
➢(1)访问一个元素: x(i)表示访问数组x的第i个元素.
➢(2)访问一块元素: x(a :b :c)表示访问数组x的从第a个 元素开始,以步长为b到第c个元素(但不超过c),b可以为负 数,b缺损时为1.
➢(3)直接使用元素编址序号. x([a b c d]) 表示提取数组x的 第a、b、c、d个元素构成一个新的数组
➢这些结构经常包含大量的MATLAB命令,经常出现在 MATLAB程序中,而不是直接加在MATLAB提示符下.
➢1、for循环:允许一组命令以固定的和预定的次数重复
➢
for x=array
➢
{commands}
➢
end
➢ 在for和end语句之间的命令串{commands}按数组(array) 中的每一列执行一次. 在每一次迭代中,x被指定为数组的下一 列,即在第n次循环中,x=array(:,n)
➢(2)“%” 后面所有文字为注释.
➢(3) “...”表示续行.
9
一、变量与函数
➢3、数学函数
函数 名 称函 数 名 称
sin(x) cos(x) tan(x) abs(x) min(x) sqrt(x) log(x) sign(x)
正弦函数 余弦函数 正切函数 绝对值 最小值 开平方 自然对数 符号函数
命令为: PLOT(X,Y,S)
线型
•y 黄色 . 点 - 连线
•m 洋红 o 圈 : 短虚线
•c 蓝绿色 x x-符号 -. 长短线 r 红色 + 加号 -- 长虚线
X,Y是向量,分别表示点集的横坐标和纵坐标
26
六、MATLAB作图——二维图形
➢PLOT(X,Y)--画实线
➢PLOT(X,Y1,S1,X,Y2,S2,……,X,Yn,Sn) --将多条线画在一起
asin(x) acos(x) atan(x) max(x) sum(x) exp(x) log10 (x) fix(x)
反正弦函数 反余弦函数 反正切函数 最大值 元素的总和 以 e 为底的指数 以 10 为底的对数
取整
10
一、变量与函数