matlab学习笔记

1).>> x=-2*pi:0.1:2*pi;>> y1=sin(x);>> y2=sin(x+pi/3)+2;>> y3=cos(x);>> plot(x,y1,'.-');>> hold on>> plot(x,y2,'.-')>> hold on>> plot(x,y3,'.-')>> plot(x,y1,'.-');>> hold on>> plot(x,y2,'.-')>> hold on>> plot(x,y3,'.-')>> h=legend('sin($x$)','sin($x+\frac{\pi}{3}$)','cos($x$)')h =159.0056Warning: Unable to interpret TeX string "\frac{\pi}{3}$)".>> set(h,'Interpreter','latex')>> xlabel('$x$','Interpreter','latex')>> ylabel('$y$','Interpreter','latex') %添加函数图标及纵横坐标轴名称，如下图显示>>2).常用画图命令：Surf （x,y,z ）、ezsurf(z)、ezmesh(z)%使用surf 画图，要求想x,y,z 是同维数的矩阵%使用ezsurf 和ezmesh 时，曲面函数必须能写成显函数的形式，如sin()/()Z xy xy =%可以使用@（x,y ）+函数表达式 定义匿名函数，如@(,)sin(*)/(*)z x y x y x y = %surf(x,y,z,x); %设置图形颜色变化方向为沿x 轴，此项可省，一般不要求。

图形属性： axis 属性Axis tight %将坐标轴减小到可能的最小尺寸 Axis equal %使x,y 轴的刻度设置相同 Axis square %使图形长宽相同 Axis off %不显示坐标轴Axis([x0 x1 y0 y1 z0 z1]); %手动改变轴的范围View （n,m ） %设定图形观察视角，n 为在水平方向测量的角度，m 为从水平面向上测量的仰俯角。

Subplot(n,m,x) %将图像窗口分为nxm 个子窗口，当前指的是第x 个，按行数。

当要对第y 个图做一定的设置时，先输入命令subplot(n,m,y),然后再输入要做的设置命令。

Legend （‘xxx ’） %给图形设置图标，也可指定图标的位置，默认为右上角。

Text （x,y,‘xxx ’） %给图像添加注释内容，注释位于图形中（x,y ）处。

Gtext(‘xxx ’); %给图像添加注释，位置随意，将十字光标点在哪儿就放哪儿。

delete(findobj('type','text')) %删除通过gtext()命令添加的注释内容 Close all %关闭所有图像窗口 Title （’xxx ‘） %给图形添加标题Material 属性，设置曲面的材料属性，值有dull,shiny,metallic ，该命令可结合用在surf,mesh,pcolor,fill,fill3等函数中。

Plotyy(x,y1,x,y2,’plot ’) %在同一个图中画出两个图形，公用一个横轴，纵轴不同，分列两边，且纵轴刻度设置可不同，自动的。

Semilogx 和semilogy 表示横轴和纵轴采用对数坐标画出x-y 的图形，用法为semilogx(x,y)和semilogy(x,y).Loglog(x,y) %横纵轴均采用对数形式画出x-y 图形Grid on %打开网格，系统默认为grid off ,即关闭网格。

A=Complex(x,y) %产生复数A=x+yi Real(A)=x,imag(A)=y.Fill （x,y,n ） %产生一个有x,y 坐标定义的填充颜色的图形，填充颜色根据n 对色图进行索引确定。

light('Color',选项一，'Style',选项二，'Position'，选项三)，选项一为表示光的颜色的矩阵，取RGB 三元组或相应的颜色字符。

选项二可取为'infinite'和'local'两个值，分别表示无穷远光和近光。

选项三为三维坐标点组成的向量形式[x,y,z]。

对远光，它表示光穿过该点射向原点；对于近光，它表示光源所在位置。

假如函数不包含任何参数，则采用缺省设置：白光、无穷远光、穿过(1，0，1)射向坐标原点。

一般后两个参数可省，默认。

--例： [x,y,z]=peaks; surf(x,y,z); shading interp;light('Color',[1 0 1],'Style','local','Position',[-4,-4,10])； %此命令表示在点[-4,-4,10]处有一处品红色光源例1.画出椭圆柱面图：3cos ,2sin ,,,x y z t t R αα=⎧⎪=⎨⎪=∈⎩编写的matlab 程序如下：@(,)3*cos();%y @(,)2*sin();@(,);(,,)x a t a a t a z a t t ezmesh x y z ===定义匿名函数%因为没能写成显函数的形式，所以调用命令函数ezmesh 时使用格式ezmesh(x,y,z)，当能写成显函数时可以直接调用格式ezmesh(z)。

例2.绘制二元函数三维表面图：sin()xy z xy =程序如下：[,]([3:0.2:3]);%@(,)sin(*)/(*);()%(),x y meshgrid z x y x y x y ezmesh z ezsurf z =-=产生所定范围内的二维点阵或用所画图像颜色效果不同Meshgrid 的用法举例：两种情况3)积分命令：a).Int---符号积分法应用： int(s)符号表达式s 的不定积分.int(s,v)符号表达式s 关于变量v 的不定积分.int(s,a,b)符号表达式s 的定积分, a,b 分别为上﹑下限. int(s,v,a,b)符号表达式s 关于变量v 从 a 到b 的定积分. 当int 求不出符号解,会自动转求数值解.注:使用符号积分时需结合使用syms --符号定义命令，例如声明两个符号变量 x,y,matlab 语句格式为：Syms x,y;例3.求函数不定积分：sin()cos(2*)exp(.^2)x x x ++ Matlab 程序为：;sin()cos(2*x)exp(.^2);%Aint()%A x zsyms x A x x z A =++=函数对函数求关于的导数并赋给b).一重数值积分命令：quad(f,a,b)、quadl(f,a,b) 其中f 为积分函数，a 、b 为积分上下限。

例4.求数值积分12231x x++⎰Matlab 程序如下：@()(2*3)./(1.^2);(,0,1)f x x x I quad f =++=或用I=quadl(f,0,1),结果的精度不同c).二重数值积分命令：dblquad(f,a,b,c,d)其中a,b 为外重积分限,c,d 为内重积分限，f 为积分函数。

d).三重数值积分命令：triplequad(f,a,b,c,d,e,f)其中a,b 为外重积分限,c,d 为中间重积分限,e,f 为内重积分限，f 为积分函数。

4).几个特殊矩阵：eye(n):n 阶单位矩阵； Ones(n):n 阶1矩阵； Zeros(n):n 阶0矩阵；5).矩阵运算：乘法：格式 a*b; %注意和数组矩阵乘法区别 a.*b;注，乘以逆和标准除法等同右乘，后除可抵消对矩阵做的左乘效果。

求逆：格式 inv(a)和运算符向后除\，例inv(a)*c<=>a\c;注，c/a 表示''(\)'CB指数运算：k-kA A 和，其中k 为正数，分别表示A 连乘K 次，和？？？rot90(A,n):将矩阵A 逆时针旋转n 个90度 tril(A):保留A 的下三角并将上三角填充为0.unit8,unit16(A):将A 中数据转化为固定位数的无符号整数类型。

6)几个命令函数：find （）:查找命令；例：a(find(a==0))=inf ,表示把a 矩阵中的0元素全部替换为无穷。

Min():查找最小命令；Numel(a)，计算矩阵a的元素个数；[y n]=Max(x,[],2),查找最大元素及其位置,返回值y为最大值，n为对应位置7)Matlab与记事本的数据交互：A）.matlab读入记事本文件中的数据：如果记事本文件*.txt中内容全部为数值是，可直接调用命令A=load（‘*.txt’），由此matlab中会生成一个数值矩阵A。

如果文件中不全为数值，则不能用load命令，应该用[name,type,value]textread （‘*.txt’，‘%s %n %f’,2）,此格式代表文件中第一列为字符串，第二列为整数，第三列为浮点数，读取命令只读取文件中的前两行，生成三个列矩阵name type value,具体情况具体分析。

注：要求所读取的txt文件的位置必须在matlab中“当前目录”中包含，如当前目录为：D：\temp\work，则txt文件必须保存在该路径下，否则命令执行出错。

Fscanf():更强大的一个读取命令。

B）.将matlab数据写入记事本文件中：如果只想保存数值数据，则可调用save命令，格式为：save filename a b ...; 代表把变量a ,b ...等值存入名为filename的.mat格式的文件中，该文件的路径为“当前目录”。

且该文件只能用matlab打开，其他软件无法打开。

如果想把各种数据保存为.txt文件，则应该使用fprintf命令，格式为：Fid=fopen(‘file.txt’,wt); %打开一个文件，并返回一个指标fidFprintf(fid,’数据表1\n’); %可省Name=’sally’;types=1;x=3;y=4; %赋值一次Fprintf(fid,’%s type%n %f %f’,name,types,x,y); %写入一次Name=’jim’;types=2;x=1;y=3; %第二次赋值Fprintf(fid,’%s type%n %f %f’,name,types,x,y); %第二次写入Fclose(fid); %关闭文件则生成一个txt文件，内容为：数据表1Sally 3 4Jim 1 38)数据拟合：A).多项式拟合：指令拟合：p=polyfit(x,y,n); %返回系数矩阵p;Yi=Polyval(p,xi); %计算多项式的拟合函数值用命令poly2sym(p)输出拟合函数表达式图形窗口拟合：先画出散点图，plot（x,y,’r*’）;然后在图形窗口上点击Tools 下的basic fitting,打开拟合对话框进行操作。