matlab基础知识简介

1、MATLAB基 础知识一、MATLAB简介二、MATLAB基础知识1、命令窗口是用户与 MATLAB进行交互作业的主要场所，用户输入的 MATLAB交互命令均在命令窗口执行。例如：在 MATLAB命令窗口下键入 a=3 2 3;4 7 6;7 5 9 按回车键后，显示结果 a=3 2 34 7 67 5 92、求逆矩阵命令 格式： 变量=inv（ 参数 ） 例如：输入 b=inv（a） ，按回车后，显示 b= 1.3750 -0.1250 -0.37500.2500 0.2500 -0.2500 -1.2083 -0.0417 0.54173、MATLAB系统还具有保存历史纪录的功能，它将本次

2、启动 MATLAB系统 之后，用户输入的命令和创建的所有变量的值保存起来，用户通过方向 键可查找所需的命令。MATLAB提供了存储变量和删除变量的命令。 SAVE 文件名 变量名 1，变量名 2，功能：将命令中的变量保存在给出的文件中。 说明：（1）若过文件名省略，默认保存在 MATLAB.MA中T 。（2）若变量名省略，则保存所有的变量到指定的文件中。（ 3 ）若文件名和变量名都省略，则保存所有定义过的变量到 MATLAB.MA中T 。例如： SAVE AA.MAT a b c %将变量 a b c 保存在文件 AA.MAT中。SAVE BB.MAT % 将所有的变量保存到文件 BB.MAT

3、中。SAVE % 将所有变量保存到文件 MATLAB.MA中T clear 变量名 1 变量名 2， . 功能：删除指定的变量。 说明：若变量名表省略，表明删除当前工作空间中的所有变量。 例如： clear a b c % 删除变量 a b cclear % 删除当前工作空间的所有变量。 MATLAB还提供了一些命令，专门管理和控制命令窗口。例如： clc 格式： clc 功能：清除命令窗口。home格式：home功能：光标移动到左上角who 格式：who功能：查看当前的所有变量，只给出变量名。whos 格式：whos功能：查看当前的所有变量，给出变量的详细信息。信息同变量浏览器。clear

4、格式：clear 变量名 功能：删除后面列出的变量，如果变量名省略，则删除所有的变量。4、MATLAB中的常量MATLAB提供了整数、实数、复数和字符四种类型数据。对应的常量 类型也是这四种。实数在屏幕显时默认的小数位数为 4 位。可以用命令 改变实数的显示格式。命令 format格式： format 格式例如： format long : format short e输出实数为 16 位：5 位加指数format long e : 16位加指数format rat :有理数近似矩阵操作1、提取矩阵的元素 例如： A=1 2 3 3; 3 2 4 1; 3 4 5 6 执行 b=A(1,2)

5、，结果为： b= 2 执行 v=A(1,2,3,4) ，结果为： v= 3 34 12、矩阵的修改 例如：执行 A(1,2)=3 后，结果为： A= 1 3 3 33 2 4 13 4 5 6 作用：修改了第一行第二列的元素值。 执行 A(: ,3)= 后，结果为：A= 1 3 33 2 13 4 6 相当于删除了第三列。 接着执行 A(2, :)= 后，结果为： A= 1 3 33 4 6 相当于删除了第二行。说明： 只能删除矩阵的某一行或某一列， 不能删除某一个元素。3、有关矩阵的其他操作( 1) 提取矩阵的上三角矩阵和下三角矩阵A= 1 3 34 2 14 4 6执行 B=triu(A)

6、 后，结果为： B= 1 3 30 2 1执行 C=tril(A) 后， 结果为：C= 1 0 04 2 04 4 6(2) 提取矩阵的对角元素 执行 B=diag(A) 后，结果为：B= 1 2 6(3) 利用向量建立对角矩阵 执行 C=diag(B) 后，结果为：C= 1 0 00 2 00 0 6(4) 矩阵的秩、特征向量和特阵值A= 1 3 34 2 43 4 5% 求矩阵的秩% 求矩阵的特征值% 求矩阵的特征向量和对应的特征值组成的对角阵执行 rank(A) 后，结果为： ans=2执行 eig(A) 后， 结果为： ans=10.0000-2.0000-0.0000执行 v, d=

7、eig(A) 后， 结果为：v =-0.4243 0.5883 0.4243-0.5657 -0.7845 0.5657-0.7071 0.1961 -0.7071d =10.0000 000 -2.000000 0 -0.0000(5) 矩阵的算术运算和矩阵元素之间的运算 a= 1 23 2b= 2 1 1 2执行 c=a+b 后，结果为：c= 3 34 4执行 d=a-b 后，结果为：d= -1 12 0执行 e=a*b 后，结果为：执行 f=a.*b 后，结果为：f= 2 2 3 4执行 g=a.b 后，结果为：g= 2 1/21/3 1执行 h=a./b 后，结果为：h= 1/2 23

8、 1执行 j=a.b 后，结果为：j= 1 23 4执行 k=b.3 后，结果为：k= 8 1 1 8试验二 MATLAB程序设计1、关系运算符 小于 大于 = 大于等于 = = 等于 = 不等于 格式： A 关系运算符 B 例如： a=-1 2 4;5 4 8; b=0 1 5;5 1 2;则 c=abc=0 1 00 1 12、逻辑运算符 逻辑与： & , 格式： c=a&b 逻辑或： | ，格式： c=a|b 逻辑非： ，格式： c=a 例如： c=a&b c=0 1 11 1 1 c=a|b c=1 1 11 1 13、程序的控制结构选择结构：（ 1）if 语句格式 1： if end

9、格式 2： if elseend格式 3： if elseif elseif else格式 4：ifif elseelse if elseend（ 2）switch 语句 格式： switch 表达式case 表达式值 1语句组 1case 表达式值 2语句组 2case 表达式值 n语句组 notherwise语句组 n+1end 例 1 ： data=input( input the value of data); n=mod(data, 5);switch ncase 1 y=1;case 2y=2;case 3y=3;case 4 y=4 otherwise y=n;endy例 2 ：

10、var=input( input a number ); switch varcase 1disp( 1)case 2,3,4disp( 2 or 3 or 4 )case 5disp( 5 )otherwisedisp( something else )end(3) while 语句格式： while( 表达式 ) 语句组end例 3：使用 while 语句求 1+2+3+100 。i=1 ；sum=0; while(i=100) sum=sum+i: i=i+1;endstr= 计算结果为： num2str(sum);(4) for 循环语句格式： for 循环变量 =初始值：增量：终值

11、语句组； end例 4：使用 for 语句求 1+2+3+100 。sum=0;for i=1:100 sum=sum+i;endstr= 计算结果为： num2str(sum);disp(str)4、函数 M-文件 格式： function 输出参数表 = 函数名(形式参数) 语句组end说明：( 1)函数 M-文件独自保存为一个文件，文件名就是函数名。(2) 调用格式： 输出参数表 =函数名(实际参数)(3) 输出参数可以是多个变量，两两之间用逗号隔开，它们表示 要计算的项；4）形式参数是一组形式变量，本身没有任何意义，也不代表任何内容，只有在调用时赋予它实际值，才有意义。例 5：求向量的

12、平均值和标准差 function mean, stdev=stat（x） n=length（x ）;mean=sum（x）/n; stdev=sqrt（sum（x.2）/n-mean.2）; end 说明：函数 M-文件保存为 stat.m 。 length（x）: 求向量 x 的长度。 sum（x）: 求向量 x 中各个元素的总和，如果 x 是一矩阵，则按 列求和。sqrt（x）: 求 x 的平均值。 调用：重新打开一 m文件，输入： x=1 2 3;6 2 1;-1 -6 -4;mean, stdev=stat（ x ） ;然后保存为任一文件名即可。 （除函数名外）试验三函数极限与导数1、

13、函数的极限格式： limit（ 表达式，变量，常量， right 和 left ） 含义：当变量趋向于常量时， 表达式的左极限或右极限 。 说明：（ 1）当变量省略时，默认的变量是 x，没有 x 时，找字母表中离 x 最近的字母作为变量。（ 2）当常量省略时，默认为 0。 （3）在不说明是左极限还是右极限时，求表达式的极限。例如：syms x a t h limit(sin(x)/x); limit(x-2)/(x2-4),2);limit(1+x)/x)(2*x),x,inf); v=(1+a/x)x,exp(-x);limit(v,x,inf, left );2、导数与微分 求导公式 格式

14、： diff( 表达式，变量，阶数 ) 功能：它表示将“表达式”按照“变量”求“阶数”阶导数。 说明：若阶数省略，则表示求一阶导数。例如：syms xf=log(x);diff(f)diff(f,4) g=(x+exp(x)*sin(x)(1/2);diff(g) pretty(ans) diff(g,2) simplify(ans) pretty(ans)3、积分格式： int( 被积函数，积分变量，积分下限，积分上限 ) 功能：求“被积函数”对“积分变量”从“积分下限”到“积分上限” 的定积分。说明：当积分上下限省略时，求“被积函数”的不定积分。 例如：syms x a b int(x*l

15、og(x)int(1/(1+sqrt(1-x2) int(cos(log(x) eval(int(exp(-x2),x,-1,1) eval(int(exp(-x2),x,-inf,1) eval(int(exp(-x2),x,-inf,inf)试验四 微分方程、多元函数微积分一、解方程1、代数方程格式： solve(f,t)功能：对变量 t 解方程 f=0，t 缺省时默认为 x 或最接近字母 x的符号变 量。例如：求解一元二次方程 f=a*x2+b*x+c 的实根，执行 solve(f,x)结果为：ans=1/2/a*(-b+(b2-4*a*c)(1/2)1/2/a*(-b-(b2-4*a*

16、c)(1/2)又如：求解一元二次方程 g=4*x2+6*x-9 的实根，执行 solve(g,x);eval(ans)结果为： 0.9271 , -2.42712、微分方程(1) 常微分方程的同解 当微分方程中不含任何附加条件(初值条件)时，其解就是该方程 的通解。格式： dsolve( 微分方程 , 自变量 ) 说明：其中“微分方程”和“自变量”都按字符串形式给出，并用 单引号括住， 且当自变量为 t 时可省略。例如： syms x y z t a b dsolve( Dy-x*Dy=a*(y2+Dy) , x) 注： Dyy , D2yy ”, D3yy ”, 依次类推。dsolve( D

17、3y=exp(2*x)-cos(x) , x)(2) 常微分方程的特解格式：dsolve( 微分方程，初值条件 1，初值条件 2，自变量 )说明： y(x0)=y0y| x=x0=y0 ;Dy(x0)=y1 y | x=x0=y1 。 例如： syms x y z t a bdsolve( Dx=a*x*(b-x) , x(0)=x0 ) ; 自变量为 t dsolve( x*D2y+Dy=0, y(1)=1, Dy(1)=6 , x)(3) 常微分方程组求解 对于无初值条件的微分方程组，其通解的求法格式为： dsolve( 微分方程 1，微分方程 2, , 自变量 ) 对于带初值条件的微分方

18、程组，其特解的求解格式为：dsolve( 微分方程 1，微分方程 2, , 条件 1，条件 2，自变量 )例如：求f 3g sinx g f cosxf (2) 0 的通解和特解，其命令如下： f (0) 0 g(0) 0syms x f g hf1,g1=dsolve( D2f+3*g=sin(x), Df+Dg=cos(x) , x) f,g=dsolve( D2f+3*g=sin(x), Df+Dg=cos(x) , Df(2)=0, f(0)=0,g(0)=0 , x)、多元函数微积分1、多元函数极限 这里我们仅对极限存在的函数， 求沿坐标轴方向的极限。 ，即将求多元函数极限问题，化成

19、求多次单极限的问题。例如：syms x y z t a b limit(limit(x2+y2)/(sin(x)+cos(y),0),pi) limit(limit(limit(x+y*z+exp(sin(x*z)/(x+y/z),1),2),3) 2、多元函数偏导数 对多元函数中的某一变量求偏导数，等价于将其余变量看作常量， 仅对该变量求导数。syms x y z t u v diff(x3+2*x*y+y2-6,x) diff(diff(x3+2*x*y+y2-6,x),y)试验五 绘图和级数一、化简和代换simplify 利用各种恒等式化简代数式 expand 将乘积展开为和式 fact

20、or把多项式转换为乘积形式collet合并同类项horner 把多项式转换为嵌套表示形式 例如：合并同类项： collect(3*x3-0.5*x3+3*x2)结果为：ans=5/2*x3+3*x2 进行因式分解执行：factor(3*x3-0.5*x3+3*x2) 结果为：ans=1/2*x2*(5*x+6) 展开多项式执行：expand(x-a)*(x-b)结果为：ans=x2-x*b-a*x+a*b 把多项式转换成嵌套表示形式执行： horner(x3-6*x2+11*x-6)结果为：ans=-6+(11+(-6+x)*x)*x二、绘图命令matlab 提供了两个非常简便的绘图命令，分别

21、如下：1、格式： ezplot(f, xmin, xmax)功能：绘出以 x 为横坐标的曲线 f ，x 在xmin, xmax内取值，当 xmin 和 xmax缺省时，默认 xmin=-2*pi ，xmax=2*pi 。例如：描述函数 cos(2*x) 在区间 -2*pi, 2*pi 上的图形，执行命令 如下：ezplot(cos(2*x), -2*pi, 2*pi)2、格式： plot(x, y,格式串 , )功能：绘出以 x 为横坐标 y 为纵坐标的散点图或折线图。其中， x 和 y 为同维向量，“格式串”指定图形格式，包括颜色、线形等。例如：执行x=1 2 3 5 8;y=1 14 19

22、 22 37;plot(x, y, -*k )其中， -*k 表示黑色实线，点用“ * ”表示。 初等函数的图形： x=linspace(-pi,pi,60);z=linspace(0.1,2*pi,60); t=linspace(-1,1,60);u=linspace(-10,10,60);y1=x; y2=x.2; y3=x.3;y4=x.4; subplot(3,2,1);plot(x,y1, -,x,y2, : ,x,y3, -. ,x,y4, - - ); axis(-pi pi -2 5);text(-0.5,4, 幂函数 ) ；三、级数1、常数项级数 格式： symsum(一般项

23、，变量，起始，终止 ) 功能：用指定变量从“起始”到“终止”求级数和。当“终止”值取有 限值时，可求出级数指定项数的部分和。例如： symsum(1/n2,n,1,10);symsum(1/n,n,1,10); symsum(1+n)/(1+n2),n,1,10);当“终止”值为 inf 时，可求出级数的收敛和。 例如：symsum(1/n2,n,1,inf);symsum(1/n,n,1,inf); symsum(1+n)/(1+n2),n,1,inf);当“终止”值为变量，比如取 m时，可求得级数部分和数列的一般项， 该一般项以 m为变量，当 m 时，取得极限即为收敛和。例如：symsum

24、(1/n,n,1,m) symsum(1/n2,n,1,m);symsum(1+n)/(1+n2),n,1,m);2、幂级数幂级数的收敛半径：limannan 1r，求的收敛半径，用与求常数项级数收敛和相同的办法求幂级数的收敛和函数例如：syms x nlimit(n+1)/n,n,inf)symsum(-1)(n-1)*xn/n,n,1,inf) limit(2(n+1)*(n+1)/(2n*n),n,inf)symsum(x-1)n/(2n*n),n,1,inf)3、将函数展开成幂级数格式： taylor( 函数，阶数 n，x0 点，变量 x)功能： 数。说明：将给定“函数”在“ x0 点

25、”处按“变量 x”展开成 n-1 阶的泰勒级(1)阶数 n 省略，默认为 6；(2)x0 点省略，默认为 0；( 3)“变量 x”省略, 默认用 findsym(f) 查找默认变量。( 4)后面的三个参数顺序还可以打乱次序，但若 x0 点和阶数均为整数时，必须阶数在前， x0 点在后例如： taylor(exp(-x),x,8,7)表示将 e x在 x0=7 点处展开为七阶泰勒级数。taylor(exp(x)将e x在 x=0点展开成 5阶泰勒级数。2taylor(exp(-x2),x,n) %taylor(cos(x),8) % 注：参数终止给出一个整数，将e x2在 x=n点展开成 5阶泰

26、勒级数 .将 cos(x) 在 x=0 点展开成 7 阶泰勒级数 认为是阶数。四、傅立叶级数三角级数：a0(an cosnx bn sinnx)2 n 111其中， a0f(x)dx ，anf ( x)cos nxdx (n=0,1,2,31bnf ( x)sin nxdx (n=0,1,2,3, ) ，正弦级数： bn sin nx (n=0,1,2,3,)余弦级数： a20 an cosnx (n=0,1,2,3,)在 matlab 中，编写以下四个函数 mfourier(f), msin(f), mcos(f)和mlfourier(f, l) ，分别用于求函数 f 的傅立叶级数、正弦级数

27、、余弦级 数，以及周期为 2l 的傅立叶级数的系数。function a0, an, bn=mfourier(f) syms n xa0=int(f, -pi, pi)/pi ； an=int(f*cos(n*x), -pi, pi)/pi；bn=int(f*sin(n*x), -pi, pi)/pi；endfunction bn=msin(f) syms n x bn=int(f*sin(n*x), -pi, pi)/pi；endfunctiona0, an=mcos(f)syms n xa0=int(f, -l, l)/pi ； an=int(f*cos(n*x), -pi, pi)/pi

28、 ； endfunction a0, an, bn=mlfourier(f, l) syms n xa0=int(f, -pi, pi)/pi ； an=int(f*cos(n*pi*x/l), -l, l)/l；bn=int(f*sin(n*pi*x/l), -l, l)/l；end说明：以上四个函数分别保存在四个 M-文件中，以函数名作为文件名。试验六 方程组求解和数值微积分一、欠定方程组的求解若有方程组 AX=B，其中 X为 m个变量构成的列向量，一共有 n个方 程，则 A 为 n*m 矩阵。若 rank(A)m ，则方程有无穷多组解，称为欠定 方程组。若 b 中元素全为 0，称方程组为

29、其次线形方程组。1、齐次线性方程组的求解x1 2x2 3x3 4x4 0如方程组：x2 x3 x4 0x1 x3 2x4 0x1 3x2 4x3 5x4 0在 MATLAB中输入A=1 -2 3 -4;0 1 -1 1;-1 0 -1 2;1 -3 4 -5;a=null(A) % 返回矩阵 A 的零空间的标准正交基组成的矩阵 a。 a= 0.7888 -0.2110-0.4999 -0.2889-0.2110 -0.78880.2889 -0.4999因而，原方程的通解为： x10.7888 0.2110k2 0.2889 % 其中 k1 与 k2 为任意常数。0.78880.4999x2x

30、3k10.49990.21100.28892、欠定方程组的求解 从理论上说，欠定方程组有无穷多组解。这时，我们仍然可用广义 逆矩阵法和矩阵除法求解方程组。前者所求的解是所有解种范数最小的 一个，而后者求的解是所有解中含零个数最多的一个。他们求出的解是 所有解系中的一个特解。 要求其通解， 则先解其对应的齐次方程的通解， 然后再加上自身的一个特解即可。如：A=1 -2 3;0 1 -1;-1 0 -1;1 -3 4;b=4 -3 -4 1 ;x=pinv(A)*bx= -2.25491.21571.0392y=Abwaring: rank deficient, rank=2 tol=4.6151

31、e-015. y= 03.47063.2941B=null(A)B= 0.5774-0.5774-0.5774 故方程组的通解为x10.5774x2 k 0.5774 x3 0.57742.25491.21571.0392x1x2kx30.57740.57740.577403.47063.2941二、数值微积分1 、数值微分多项式求导法求微分若已知函数在某些接点处的值， (1) 只要用曲线拟合法得到一个函 数近似多项式(即用多项式近似代替函数) ，(2) 再对该多项式求微分， (3) 然后对微分后的多项式求值， 便可以方便的求出在拟合范围内的任一 点处的任意阶微分。几个命令：polyfit(x

32、, y, n)利用数据点向量 x 和 y 拟合一个 n 次多项式，返回多项式系数向量，系数由高次方朝低次方顺序排列polyder(P)返回多项式 P 求导后对系数向量。polyval(P, x0)计算多项式 P在点 x0处的值。 P是多项式系数。 根据上述几个函数，我们编写一个利用多项式方式求任何函数在任 一点处的任意阶导数的近似值。fpder(f, x0, n, m, d, k)f: 欲求导函数的符号表达式；x0: 固定点，即求函数在该点处的导数；n: 表示求导阶数，默认为一阶导数；m: 为多项式最高次数， 表示用 m次多项式近似模拟给定函数， 默认 为 6；d: 为半区间长度，即在 x0-

33、d, x0+d 区间内拟合多项式，默认值为 1；k: 为区间采样点数，即在闭区间 x0-d, x0+d上采 k+1 个点来拟合 多项式，默认值为 20；函数如下：function y=fpder(f, x0, n, m, d, k)if nargin6k=20;endif nargin5d=1.0;endif nargin4m=6;endif nargin= =2n=1;end x=x0-d:2*d/k:x0+d; y=eval(f); p=polyfit(x, y, m); pp=p;for n1=1:n pp=polyder(pp);end y=polyval(pp,x0);例如：fpde

34、r( sin(x) , pi/4, 3, 10, 0.5, 10)fpder( exp(-x) , 1)fpder( exp(-x.2) , 1, 5, 10, 0.5)2、数值积分 数值积分就是对那些不能求积分或不易求积分的函数提供的一种较 精确的积分方法。在 MATLAB中，有两个函数是专门求数值积分的函数， 分别是：quad( 函数名串 , a, b, tol, trace) 和 quad8( 函数名串 a, b, tol, trace)说明：函数名串是 matlab 内建函数的函数名或用户自定义函数的函数 名；a 为积分下限； b 为积分上限；tol 为数值方法允许误差；默认为 10-

35、3 ；trace 一般为零。quad8 ()和 quad()的用法完全相同。但在功能上有所区别： 前者用 newton-cotes 法求数值积分，后者用 simpson 法求数值积分。 例如： quad( sin ,0,pi) quad( sin ,0,pi,1e-6) quad( log ,0.1,1)再定义两个函数 f1(x) 和 f2(x) 如下：function y=f1(x)y=1./(1+x.4);function y=f2(x)y=log(1+tan(x);则有 quad( f1 ,0,1)quad( f2 ,0,pi/4)3、微分方程数值解一般情况下，除了简单的微分方程外，要找

36、出微分方程解的解析表 达式是极其困难的，有时甚至是不可能的。微分方程数值解法是能够计 算出解在若干个离散点上近似结果的一种通用方法。MATLAB提供了一些求解微分方程(组)的数值方法函数。ode23( )%低阶法解非刚性微分方程ode113( )%变序法解非刚性微分方程ode15s( )%变序法解刚性微分方程ode23s( )%高阶法解刚性微分方程以 ode23( ) 函数为例，其格式为ode( F ,Tspan, y0,OPTIONS) 说明：(1)F为函数名，即要将微分方程存入一个函数文件中，然后 调用。(2)Tspan为一向量 t0,tm 或 t0,t1,t2, ,tn ，前者是求解区

37、间，后者是指定要求解的点，必须按升序或降序排列。(3)y0 为初值条件，即 y|t t0 y0；(4) OPTIONS可以制定一系列可选项，包括允许绝对误差、允许 相对误差、是否作图等等。当该项省略时，使用软件默认值。( 5)当微分方程为方程组时， F 与 y0 均为向量， F 为微分方程 dy F(t, y)中的 F(t,y) 函数。dtdy t 2 y例 如 ： 求 微 分 方 程 dt t 在 区 间 1,1.1 的 数 值 解 ； 再 求 点y| t 1 1y| t 1 21,1.02,1.04,1.06,1.08,1.10 处的近似值。编写下述函数，并以 df1 为文件名存盘：fun

38、ction F=df1(t,y)F=t2+y/t ;然后分别用函数 ode23( ), ode45( ) 和 ode113( ) 在区间1,1.1 上求方 程的数值解。格式： T,Y= 函数(参数)说明：T位子变量值， Y为对应于 T的因变量值，若直接使用函数，可画 出数值解的图形。ode23(df1 ,1,1.1,0.5)T, Y=ode23( df1 ,1,1.1,0.5)T,Y=ode23( df1,1,1.02, 1.04,1.06,1.08,1.10,0.5)作业：求解下面四个微分方程在指定区间的函数值。(1) xy -ylny=0(2) 3x2+5x-5y =0,(3) y -xy =a(y 2+y)常染色体遗传模型随着人类的进化，人们为了揭示生命的奥妙，越来越重视遗传学定额研究，特别是遗传