[工学]MATLAB复习.ppt

MATLAB复习 概述 帮助命令及区别 程序设计 变量命名 内存变量的管理 矩阵变量的建立 运算符 算术运算、关系运算、逻辑运算关系 函数 程序设计 函数文件 排序问题 数学家发现有这样一种数，它的约数（不含它本身）的和恰好等于它本身。如6=1+2+3，28=1+2+4+7+14。具有这种性质的第三个数是496，第四个数是8128……。人们把这样的数叫做完全数。编写程序输入一个数，判断是否为完全数。如果是完全数则输出“是完全数”，如果不是则输出“不是完全数”。 数值计算 矩阵的秩与特征值、特征向量 数据处理 曲线拟合 数值插值 3x1+ x2 - x3 = 3.6 x1+2x2+4x3 = 2.1 -x1+4x2+5x3 = -1.4 A=[3 1 -1;1 2 4;-1 4 5];b=[3.6;2.1;-1.4]; Ax=b 解为:x=A-1b x=inv(A)*b x = 1.4818 -0.4606 0.3848 A=[3,2, 11,1;7, 15,4,3; 5,10, 9,18]; 求每行每列的最大元素及矩阵的最大元素 y1=max(A,[],2); y2=max(A); y3=max(y1); A=[3,15,7;9,6,8]; [Y,I]=sort(A,2); Y = 3 7 15 6 8 9 I = 1 3 2 2 3 1 Y1=interp1(X,Y,X1,’method’); method插值方式:linear(线性插值),nearest(最近插值),spline(三次样条插值),cubic(三次多项式插值)默认插值是linear(线性插值). X=0:0.2:pi; Y=sin(X); y1=interp1(X,Y,pi/2); y2=interp1(X,Y,pi/2,nearest); x=0:0.1:3; y=sin(x); p=polyfit(x,y,8); y1=polyval(p,x); plot(x,y1,’k*’,x,y,’k-’); 符号运算 符号变量的创建 符号变量与数值变量的区别 微分、积分运算 泰勒级数 傅立叶变换 微分运算 diff(f,x,n) 求y=x^2*sin(x)的一阶二阶导数 程序:syms x; y=x^2*sin(x); diff(y,x); diff(y,x,2); 积分运算 int(f,x,a,b); 求解方程 solve(y,a) 代数方程求解 [x,y]=solve(‘x*y=2’,’x^2+y^2=4’); 微分方程和方程组求解 y=dsolve(‘Dy=sin(t)’); dsolve的默认变量是t 绘图 plot函数的用法 常用绘图选项 图形处理 窗口分割 条形图、直方图、饼图 例：用不同线形和颜色绘制曲线 及其包络线。 x=0:0.1:2*pi; y1=2*exp(-0.5*x); y2=2*exp(-0.5*x).*sin(2*pi*x); x1=0:0.5:6; y3=2*exp(-0.5*x1).*sin(2*pi*x1); y4=-2*exp(-0.5*x); plot(x,y1, g:,x,y4, g:, x,y2,b--,x1,y3,rp); 例： x=0:0.1:2*pi; y1=2*exp(-0.5*x); y2=2*exp(-0.5*x).*sin(2*pi*x); x1=[0:12]/2; y3=2*exp(-0.5*x1).*sin(2*pi*x1); y4=-2*exp(-0.5*x); plot(x,y1, g:,x,y4, g:, x,y2,b--,x1,y3,rp) title(曲线及其包络线); xlabel(变量x); ylabel(变量y); text(3.2,0.5,包络线); text(0.5,0.5,曲线y); text(1.4,0.15,离散数据点); legend(包络线,包络线,曲线y,离散数据点); 图形窗口的分割 subplot(m,n,p); [x,y]=meshgrid(-8:0.5:8); r=sqrt(x.^2+y.^2)+eps; z=sin(r)./r; mesh(x,y,z); bar(x,y) x为横坐标向量。y可以是向量或矩阵。y为向量时，每个元素对应一个竖条；y为m行n列时，将画出m组竖条，每组n条。 例