MATLAB初步学习

都是自己打的，有些知识直接复制的，但是代码都是运行调试过的，已经调了最全的基础知识了，基本都含括里面了，再高深的我也写不出来，有些地方没有代码因为太简单了，有些地方只有代码没有说明，就自己领悟吧！1矩阵及其基本运算1.1矩阵的表示1.1.1矩阵的生成1.矩阵的输入输入矩阵都是按行列，每一行之间的元素用逗号隔开，行与行之间用分号隔开。所有矩阵元素用[]括起来。2.复数矩阵输入在matlab中，输入一句命令后，后面加入分号表示结果不显示。1.1.2符号矩阵需要用到符号（symbol）矩阵定义函数sym，或者是符号定义函数syms，先定义符号变量，再输入符号矩阵。1．sym定义矩阵：还没研究出来2.用syms定义矩阵还没研究出来额，总体来说，编写m文件的时候，要定义函数的变量使用的是syms后面加上变量。3.Sym是把数值型转化为符号型转化之后是最接近原式数据的方式表达出来。1.1.3大矩阵的生成通过编写m文件输入大矩阵。size函数可以表示出矩阵的行列数。1.1.4多维数组有时候数据太大，需要拆分组装，cat函数A=cat[n,A1,A2,A3...]可以构造n维数组1.1.5特殊矩阵命令全零阵函数zeros格式B=zeros(n)%生成n×n全零阵B=zeros(m,n)%生成m×n全零阵B=zeros([mn])%生成m×n全零阵B=zeros(d1,d2,d3…)%生成d1×d2×d3×…全零阵或数组B=zeros([d1d2d3…])%生成d1×d2×d3×…全零阵或数组B=zeros(size(A))%生成与矩阵A相同大小的全零阵d1行，d2列，d3d4个命令单位阵函数eye格式Y=eye(n)%生成n×n单位阵Y=eye(m,n)%生成m×n单位阵Y=eye(size(A))%生成与矩阵A相同大小的单位阵命令全1阵函数ones格式Y=ones(n)%生成n×n全1阵Y=ones(m,n)%生成m×n全1阵Y=ones([mn])%生成m×n全1阵Y=ones(d1,d2,d3…)%生成d1×d2×d3×…全1阵或数组Y=ones([d1d2d3…])%生成d1×d2×d3×…全1阵或数组Y=ones(size(A))%生成与矩阵A相同大小的全1阵命令均匀分布随机矩阵函数rand格式Y=rand(n)%生成n×n随机矩阵，其元素在（0，1）内Y=rand(m,n)%生成m×n随机矩阵Y=rand([mn])%生成m×n随机矩阵Y=rand(m,n,p,…)%生成m×n×p×…随机矩阵或数组Y=rand([mnp…])%生成m×n×p×…随机矩阵或数组Y=rand(size(A))%生成与矩阵A相同大小的随机矩阵rand%无变量输入时只产生一个随机数命令正态分布随机矩阵函数randn格式Y=randn(n)%生成n×n正态分布随机矩阵Y=randn(m,n)%生成m×n正态分布随机矩阵Y=randn([mn])%生成m×n正态分布随机矩阵Y=randn(m,n,p,…)%生成m×n×p×…正态分布随机矩阵或数组Y=randn([mnp…])%生成m×n×p×…正态分布随机矩阵或数组Y=randn(size(A))%生成与矩阵A相同大小的正态分布随机矩阵randn%无变量输入时只产生一个正态分布随机数命令产生随机排列函数randperm格式p=randperm(n)%产生1～n之间整数的随机排列命令产生线性等分向量函数linspace格式y=linspace(a,b)%在(a,b)上产生100个线性等分点y=linspace(a,b,n)%在(a,b)上产生n个线性等分点命令产生对数等分向量函数logspace格式y=logspace(a,b)%在（）之间产生50个对数等分向量y=logspace(a,b,n)y=logspace(a,pi)命令计算矩阵中元素个数n=numel(a)%返回矩阵A的元素的个数命令产生以输入元素为对角线元素的矩阵函数blkdiag格式out=blkdiag(a,b,c,d,…)%产生以a,b,c,d,…为对角线元素的矩阵ba10,101.2矩阵运算1.2.1加减乘除向量叉乘函数cross格式C=cross(A,B)%若A、B为向量，则返回A与B的叉乘，即C=A×B，A、B必须是3个元素的向量；若A、B为矩阵，则返回一个3×n矩阵，其中的列是A与B对应列的叉积，A、B都是3×n矩阵。C=cross(A,B,dim)%在dim维数中给出向量A与B的叉积。A和B必须具有相同的维数，size(A,dim)和size(B,dim)必须是3。1.2.2集合运算交集（常用于检测数据相关性）函数intersect格式c=intersect(A,B,'rows')%A、B为相同列数的矩阵，返回元素相同的行。[c,ia,ib]=intersect(a,b)%c为a、b的公共元素，ia表示公共元素在a中的位置，ib表示公共元素在b中位置。检测元素是否在集合中函数ismember格式k=ismember(a,S)%当a中元素属于S时，k取1，否则，k取0。k=ismember(A,S,'rows')%A、S有相同的列，返回行相同k取1，不相同取0的列向量。1.2.3除法左除C=A\B是方程A*C=B的解1.2.4行列式按理说应该是0的，可能精度没调的原因吧！1.2.5矩阵的秩2绘图2.1二维图形2.1.1平面图形命令1plot功能线性二维图。plot(X1,Y1,X2,Y2,…)，其中Xi与Yi成对出现，plot(X1,Y1,X2,Y2,…)将分别按顺序取两数据Xi与Yi进行画图。说明参数LineSpec功能定义线的属性。Maltab允许用户对线条定义如下的特性：1．线型表2-1定义符---：-.线型实线（缺省值）划线点线点划线2．线条宽度指定线条的宽度，取值为整数（单位为像素点）3．颜色表2-2定义符R（red）G(green)b(blue)c(cyan)颜色红色绿色兰色青色定义符M(magenta)y(yellow)k(black)w(white)颜色品红黄色黑色白色4．标记类型表2-3定义符+o(字母)*.x标记类型加号小圆圈星号实点交叉号定义符d^v标记类型棱形向上三角形向下三角形向右三角形向左三角形定义符shP标记类型正方形正六角星正五角星051015202530354045-4-202468命令2fplot功能在指定的范围limits内画出一元函数y=f（x）的图形。注意：fplot采用自适应步长控制来画出函数function的示意图。fplot('tanh',[-22])命令3loglog功能双对数图形。命令4semilogx功能x轴对数图形。2.1.2极坐标命令1polar功能画极坐标图。t=0:.01:2*pi;polar(t,sin(3*t).*cos(2*t),'--r')命令2bar功能二维垂直条形图。用垂直条形显示向量或矩阵中的值。x=-2.9:0.2:2.9;bar(x,exp(x.*sin(x)))colormapgrayx=-2.9:0.2:2.9;命令5comet功能二维彗星图。t=0:.01:2*pi;x=exp(sin(2*t)).*(cos(t).^2/3);y=t.*(sin(t).^2);comet(x,y);命令6pie功能饼形图x=[130.52.52];explode=[01000];pie(x,explode)2.1.3二维图形注释命令1grid功能给二维或三维图形的坐标面增加分隔线。该命令会对当前坐标轴的Xgrid，Ygrid，Zgrid的属性有影响。用法gridon给当前的坐标轴增加分隔线。gridoff从当前的坐标轴中去掉分隔线。命令2gtext功能在当前二维图形中用鼠标放置文字。当光标进入图形窗口时，会变成一个大十字，表明系统正等待用户的动作。命令3title功能给当前轴加上标题。命令4xlabel、ylabel功能给x、y轴贴上标签2.2三维图形2.2.1三维等高线命令1contour功能曲面的等高线图用法contour(z)把矩阵z中的值作为一个二维函数的值，等高曲线是一个平面的曲线，平面的高度v是Matlab自动取的；contour(x,y,z)(x,y)是平面z=0上点的坐标矩阵，z为相应点的高度值矩阵。效果同上；contour(z,n)画出n条等高线；contour(x,y,z,n)画出n条等高线；contour(z,v)在指定的高度v上画出等高线；命令2pie3功能三维饼形图注意：命令pie3将x的每一个元素在所有元素的总和中所占的比例表达出来。若x中的分量和小于1（则所有元素小于1），则认为x中的值指明三维饼形图的每一部分的大小。x=[130.52.52]ex=[01000]pie3(x,ex)2.2.3曲面网格命令1mesh功能生成由X，Y和Z指定的网线面，由C指定的颜色的三维网格图。2.2.4其他命令axis功能坐标轴的刻度与外在显示用法axis([xminxmaxyminymax])设置当前坐标轴的x-轴与y-轴的范围。axisequal设置坐标轴的纵横比，使在每个方向的数据单位都相同。3符号运算3.1.基本运算命令1设置变量的精度函数digits格式digits(d)%设置当前的可变算术精度的位数为整数d位d=digits%返回当前的可变算术精度位数给ddigits%显示当前可变算术精度的位数命令2将符号转换为MATLAB的数值形式函数double格式R=double(S)命令3搜索符号表达式的最简形式函数simple格式r=simple(S)3.2基本运算3.2.1微积分命令1极限函数limit格式limit(F,x,a,'right')或limit(F,x,a,'left')%计算符号函数F的单侧极限：左极限x→a-或右极限x→a+。命令2导数（包括偏导数）函数diff格式diff(S,'v',n)%对表达式S中指定的符号变量v计算S的n阶导数。命令3符号函数的积分函数int格式R=int(S,v,a,b)%对表达式s中指定的符号变量v计算从a到b的定积分4数据分析4.1基本函数4.1.1三角函数函数sin、sinh功能正弦函数与双曲正弦函数格式Y=sin(X)函数asin、asinh功能反正弦函数与反双曲正弦函数格式Y=asin(X)函数cos、cosh功能余弦函数与双曲余弦函数格式Y=cos(X)函数acos、acosh功能反余弦函数与反双曲余弦函数格式Y=acos(X)函数tan、tanh功能正切函数与双曲正切函数格式Y=tan(X)函数atan、atanh功能反正切函数与反双曲正切函数格式Y=atan(X)函数cot、coth功能余切函数与双曲余切函数格式Y=cot(X)函数acot、acoth功能反余切函数与反双曲余切函数格式Y=acot(X)函数sec、sech功能正割函数与双曲正割函数格式Y=sec(X)函数asec、asech功能反正割函数与反双曲正割函数格式Y=asec(X)函数csc、csch功能余割函数与双曲余割函数格式Y=csc(X)函数acsc、acsch功能反余割函数与反双曲余割函数。格式Y=asec(X)函数atan2功能四象限的反正切函数格式P=atan2(Y,X)4.1.2其他函数函数fix功能朝零方向取整格式B=fix(A)函数roud功能朝最近的方向取整。格式Y=round(X)函数floor功能朝负无穷大方向取整格式B=floor(A)函数rem功能求作除法后的剩余数格式R=rem(X,Y)函数exp功能以e为底数的指数函数格式Y=exp(X)函数log功能自然对数，即以e为底数的对数。格式Y=log(X)函数ceil功能朝正无穷大方向取整格式B=ceil(A)函数log10功能常用对数，即以10为底数的对数。格式Y=log10(X)函数sort功能把输入参量中的元素按从小到大的方向重新排列格式B=sort(A)函数abs功能数值的绝对值与复数的幅值格式Y=abs(X)