MATLAB常用程序集合
MATLAB算术运算符、关系运算符、逻辑运算符、按位集合运算符
MATLAB算术运算符、关系运算符、逻辑运算符、按位集合运算符算术运算符MATLAB允许两种不同类型的算术运算 -•矩阵算术运算•数组算术运算矩阵算术运算与线性代数中定义的相同。
在一维和多维数组中，逐个元素执行数组运算。
矩阵运算符和数组运算符由句点符号(.)区分。
然而，由于对于矩阵和阵列的加减运算是相同的，因此对于这两种情况，运算符相同。
下表简要说明了算术运算符 -算术运算符示例代码的大小，除非是标量。
可以从任何大小的矩阵中减去标量。
矩阵乘法。
C = A * B是矩阵A和B的线性代数乘积。
更准确地说，执行公式：，对于非标量A和B，A的列数必须等于B的行数。
标量可以乘以任何大小的矩阵。
阵列乘法。
A .* B是数组A和B的逐个元素乘积。
A和B必须具有相同的关系运算符关系运算符也可以用于标量和非标量数据。
数组的关系运算符在两个数组之间执行逐个元素的比较，并返回相同大小的逻辑数组，如果为真，则元素设置为逻辑1(true)，如果为假，则元素设置为逻辑0(false)。
下表显示了MATLAB中可用的关系运算符：关系运算符示例代码操作符说明描述< 小于<= 小于或等于> 大于>= 大于或等于== 等于~= 不等于逻辑运算符MATLAB提供两种类型的逻辑运算符和函数：•逐元素 - 这些运算符对逻辑阵列的相应元素进行操作。
•短路 - 这些运算符在标量和逻辑表达式上运行。
元素逻辑运算符在逻辑数组上运行逐个元素。
符号＆，|和〜是逻辑数组运算符AND，OR和NOT。
短路逻辑运算符允许逻辑运算短路。
符号&&和||是逻辑短路运算符AND和OR。
逻辑运算符示例代码位运算按位运算符对位执行，并执行逐位运算。
＆，|和^的真值表如下 -假设A = 60和B = 13; 现在以二进制格式，它们将如下所示：A = 0011 1100B = 0000 1101-----------------A&B = 0000 1100A|B = 0011 1101A^B = 0011 0001~A = 1100 0011ShellMATLAB提供了诸如“按位与”和“按位或”以及“按位非”操作，移位操作等位操作的各种功能。
matlab集合运算
matlab集合运算Matlab是一种功能强大的数学计算和编程工具，它提供了丰富的集合运算函数，可以方便地对集合进行各种操作和计算。
本文将介绍一些常用的Matlab集合运算函数，并展示它们的用法和效果。
一、集合的定义和表示在Matlab中，集合可以用一维数组来表示，每个元素表示集合中的一个元素。
例如，我们可以用数组A=[1,2,3,4,5]表示一个包含5个元素的集合A。
二、集合的基本运算1. 并集运算并集运算是指将两个集合合并为一个集合，其中包含两个集合中的所有元素，且不重复。
在Matlab中，可以使用union函数实现并集运算。
例如，对于集合A和集合B，可以使用C=union(A,B)得到它们的并集C。
2. 交集运算交集运算是指找出两个集合中共同存在的元素，形成一个新的集合。
在Matlab中，可以使用intersect函数实现交集运算。
例如，对于集合A和集合B，可以使用C=intersect(A,B)得到它们的交集C。
3. 差集运算差集运算是指找出一个集合中存在，而另一个集合中不存在的元素，形成一个新的集合。
在Matlab中，可以使用setdiff函数实现差集运算。
例如，对于集合A和集合B，可以使用C=setdiff(A,B)得到A中存在而B中不存在的元素构成的集合C。
4. 补集运算补集运算是指找出一个集合中不属于另一个集合的所有元素，形成一个新的集合。
在Matlab中，可以使用setxor函数实现补集运算。
例如，对于集合A和集合B，可以使用C=setxor(A,B)得到A和B 的补集C。
三、集合的比较运算1. 子集判断子集判断是指判断一个集合是否是另一个集合的子集。
在Matlab 中，可以使用issubset函数进行子集判断。
例如，对于集合A和集合B，可以使用tf=issubset(A,B)判断A是否是B的子集，返回值tf为逻辑值。
2. 相等判断相等判断是指判断两个集合是否相等，即两个集合中的元素完全相同。
Matlab数值积分程序集合
Matlab数值积分程序集合[图书馆+网络收集]近来学习数值积分,手头积累了不少程序,也拿来和各位朋友分享一下。
主要是来自数值积分教材和网络，基本的原理也就不打算多说了，随便搜索一下就可以得到，那就开始上代码了，呵呵，非原创，但是全部验证过，有疑问可以给我e-mail：1 梯形数值积分的MATLAB主程序function T=rctrap(fun,a,b,m)%fun 函数，a 积分上限 b积分下限 m 递归次数n=1;h=b-a; T=zeros(1,m+1); x=a;T(1)=h*(feval(fun,a)+feval(fun,b))/2;for i=1:mh=h/2; n=2*n; s=0;for k=1:n/2x=a+h*(2*k-1); s=s+feval(fun,x);endT(i+1)=T(i)/2+h*s;endT=T(1:m);e.g运行后屏幕显示精确值F s，用rctrap计算的递归值T和T与精确值F s的绝对误差w T>> ) exp((-x^.2./2)./(sqrt(2*pi)))T=rctrap(fun,0,pi/2,14), syms tfi=int(exp((-t^2)/2)/(sqrt(2*pi)),t,0, pi/2);Fs= double(fi), wT= double(abs(fi-T))fun =@(x)exp((-x^.2./2)./(sqrt(2*pi)))T =Columns 1 through 71.4168 1.3578 1.3313 1.3195 1.3142 1.3119 1.3109Columns 8 through 141.3105 1.3103 1.3102 1.3102 1.3101 1.3101 1.3101Fs =0.4419wT =Columns 1 through 70.9749 0.9159 0.8894 0.8776 0.8723 0.8700 0.8690Columns 8 through 140.8686 0.8684 0.8683 0.8683 0.8683 0.8682 0.8682>>2 复合辛普森（Simpson）数值积分的MATLAB主程序function y=comsimpson(fun,a,b,n)% fun 函数 a 积分上限 b积分下限 n 分割小区间数z1=feval (fun,a)+ feval (fun,b);m=n/2;h=(b-a)/(2*m); x=a;z2=0; z3=0; x2=0; x3=0;for k=2:2:2*mx2=x+k*h; z2= z2+2*feval (fun,x2);for k=3:2:2*mx3=x+k*h; z3= z3+4*feval (fun,x3);endy=(z1+z2+z3)*h/3;由于Matlab自带了 quad就是这个算法所以比较少自己编3 龙贝格数值积分的MATLAB主程序function [RT,R,wugu,h]=romberg(fun,a,b, wucha,m)%fun被积函数 a，b积分上下限 wucha两次相邻迭代绝对差值 m 龙贝格积分表最大行数%RT 龙贝格积分表 R 数值积分结果 wucha 误差估计 h 最小步长n=1;h=b-a; wugu=1; x=a;k=0; RT=zeros(4,4);RT(1,1)=h*(feval(fun,a)+feval(fun,b))/2;while((wugu>wucha)&(k<m)|(k<4))k=k+1; h=h/2; s=0;for j=1:nx=a+h*(2*j-1); s=s+feval(fun,x);endRT(k+1,1)= RT(k,1)/2+h*s; n=2*n;for i=1:kRT(k+1,i+1)=((4^i)*RT(k+1,i)-RT(k,i))/(4^i-1);endwugu=abs(RT(k+1,k)-RT(k+1,k+1));endR=RT(k+1,k+1);>> F=inline('1./(1+x)'); [RT,R,wugu,h]=romberg(F,0,1.5,1.e-8,13) syms xfi=int(1/(1+x),x,0,1.5); Fs=double(fi),wR=double(abs(fi-R)), wR1= wR - wuguRT =1.0500 0 0 00 00.9536 0.9214 0 00 00.9260 0.9168 0.9165 00 00.9187 0.9163 0.9163 0.91630 00.9169 0.9163 0.9163 0.9163 0.91630.9164 0.9163 0.9163 0.9163 0.9163 0.9163R =0.9163wugu =2.9436e-011h =0.0469Fs =0.9163wR =9.8007e-011wR1 =6.8571e-011>>6 复合梯形法function [I,step] = CombineTraprl(f,a,b,eps)%f 被积函数%a,b 积分上下限%eps 精度%I 积分结果%step 积分的子区间数if(nargin ==3)eps=1.0e-4;endn=1;h=(b-a)/2;I1=0;I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h; while abs(I2-I1)>epsn=n+1;h=(b-a)/n;I1=I2;I2=0;for i=0:n-1x=a+h*i;x1=x+h;I2=I2+(h/2)*(subs(sym(f),findsym(sym(f)),x)+subs(sym(f ),findsym(sym(f)),x1));endendI=I2;step=n;7 辛普森法function [I,step] = IntSimpson(f,a,b,type,eps)%type = 1 辛普森公式%type = 2 辛普森3/8公式%type = 3 复合辛普森公式if(type==3 && nargin==4)eps=1.0e-4; %缺省精度为0.0001endI=0;switch typecase 1,I=((b-a)/6)*(subs(sym(f),findsym(sym(f)),a)+...4*subs(sym(f),findsym(sym(f)),(a+b)/2)+...subs(sym(f),findsym(sym(f)),b));step=1;case 2,I=((b-a)/8)*(subs(sym(f),findsym(sym(f)),a)+...3*subs(sym(f),findsym(sym(f)),(2*a+b)/3)+ ...3*subs(sym(f),findsym(sym(f)),(a+2*b)/3)+subs(sym( f),findsym(sym(f)),b));step=1;case 3,n=2;h=(b-a)/2;I1=0;I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym (sym(f)),b))/h;while abs(I2-I1)>epsn=n+1;h=(b-a)/n;I1=I2;I2=0;for i=0:n-1x=a+h*i;x1=x+h;I2=I2+(h/6)*(subs(sym(f),findsym(sym(f )),x)+...4*subs(sym(f),findsym(sym(f)), (x+x1)/2)+...subs(sym(f),findsym(sym(f)),x1 ));endendI=I2;step=n;end8 牛顿－科茨法function I = NewtonCotes(f,a,b,type)%type = 1 科茨公式%type = 2 牛顿－科茨六点公式%type = 3 牛顿－科茨七点公式I=0;switch typecase 1,I=((b-a)/90)*(7*subs(sym(f),findsym(sym(f)),a)+...32*subs(sym(f),findsym(sym(f)),(3*a+b)/4)+...12*subs(sym(f),findsym(sym(f)),(a+b)/2)+...32*subs(sym(f),findsym(sym(f)),(a+3*b)/4)+7*su bs(sym(f),findsym(sym(f)),b));case 2,I=((b-a)/288)*(19*subs(sym(f),findsym(sym(f)),a)+...75*subs(sym(f),findsym(sym(f)),(4*a+b)/5)+...50*subs(sym(f),findsym(sym(f)),(3*a+2*b)/5)+...50*subs(sym(f),findsym(sym(f)),(2*a+3*b)/5)+...75*subs(sym(f),findsym(sym(f)),(a+4*b)/5)+19*s ubs(sym(f),findsym(sym(f)),b));case 3,I=((b-a)/840)*(41*subs(sym(f),findsym(sym(f)),a)+...216*subs(sym(f),findsym(sym(f)),(5*a+b)/6)+...27*subs(sym(f),findsym(sym(f)),(2*a+b)/3)+...272*subs(sym(f),findsym(sym(f)),(a+b)/2)+...27*subs(sym(f),findsym(sym(f)),(a+2*b)/3)+...216*subs(sym(f),findsym(sym(f)),(a+5*b)/6)+41* subs(sym(f),findsym(sym(f)),b));end9 高斯公式function I = IntGauss(f,a,b,n,AK,XK)if(n<5 && nargin == 4)AK = 0;XK = 0;elseXK1=((b-a)/2)*XK+((a+b)/2);I=((b-a)/2)*sum(AK.*subs(sym(f),findsym(f),XK1));endta = (b-a)/2;tb = (a+b)/2;switch ncase 0,I=2*ta*subs(sym(f),findsym(sym(f)),tb);case 1,I=ta*(subs(sym(f),findsym(sym(f)),ta*0.5773503+tb)+...subs(sym(f),findsym(sym(f)),-ta*0.5773503+tb));case 2,I=ta*(0.55555556*subs(sym(f),findsym(sym(f)),ta*0.7745 967+tb)+...0.55555556*subs(sym(f),findsym(sym(f)),-ta*0.7 745967+tb)+...0.88888889*subs(sym(f),findsym(sym(f)),tb));case 3,I=ta*(0.3478548*subs(sym(f),findsym(sym(f)),ta*0.86113 63+tb)+...0.3478548*subs(sym(f),findsym(sym(f)),-ta*0.86 11363+tb)+...0.6521452*subs(sym(f),findsym(sym(f)),ta*0.339 8810+tb)...+0.6521452*subs(sym(f),findsym(sym(f)),-ta*0.3 398810+tb));case 4,I=ta*(0.2369269*subs(sym(f),findsym(sym(f)),ta*0.90617 93+tb)+...0.2369269*subs(sym(f),findsym(sym(f)),-ta*0.90 61793+tb)+...0.4786287*subs(sym(f),findsym(sym(f)),ta*0.538 4693+tb)...+0.4786287*subs(sym(f),findsym(sym(f)),-ta*0.5 384693+tb)+...0.5688889*subs(sym(f),findsym(sym(f)),tb)); end10 区间逐次分半梯形法function q=DblTraprl(f,a,A,b,B,m,n)if(m==1 && n==1) %梯形公式q=((B-b)*(A-a)/4)*(subs(sym(f),findsym(sym(f)),{a,b})+...subs(sym(f),findsym(sym(f)),{a,B})+...subs(sym(f),findsym(sym(f)),{A,b})+...subs(sym(f),findsym(sym(f)),{A,B}));else %复合梯形公式 C=4*ones(n+1,m+1);C(1,:)=2;C(:,1)=2;C(n+1,:)=2;C(:,m+1)=2;C(1,1)=1;C(1,m+1)=1;C(n+1,1)=1;C(n+1,m+1)=1; %C矩阵endF=zeros(n+1,m+1);q=0;for i=0:nfor j=0:mx=a+i*(A-a)/n;y=b+j*(B-b)/m;F(i+1,j+1)=subs(sym(f),findsym(sym(f)),{x,y});q=q+F(i+1,j+1)*C(i+1,j+1);endendq=((B-b)*(A-a)/4/m/n)*q;11 区间逐次分半布尔法function [I,step] = DDBuer(f,a,b,eps)if(nargin==3)eps=1.0e-4;end;n=1;h=b-a;I1=0;I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h/ 2;tol=1;while tol>epsn=n+1;h=h/4;I1=I2;I2=0;for i=0:(4^(n-1)-1)x=a+h*4*i;x1=x+h;x2=x1+h;x3=x2+h;x4=x3+h;I2=I2+(2*h/45)*(7*subs(sym(f),findsym(sym(f)),x)+...32*subs(sym(f),findsym(sym(f)),x1)+...12*subs(sym(f),findsym(sym(f)),x2)+...32*subs(sym(f),findsym(sym(f)),x3)+...7*subs(sym(f),findsym(sym(f)),x4));endtol=abs(I2-I1);endI=I2;step=n;10 自适应求积法function I=SmartSimpson(f,a,b,eps)if(nargin==3)eps=1.0e-4;end;e=5*eps;I=SubSmartSimpson(f,a,b,e);function q=SubSmartSimpson(f,a,b,eps)QA=IntSimpson(f,a,b,1,eps);QLeft=IntSimpson(f,a,(a+b)/2,1,eps);QRight=IntSimpson(f,(a+b)/2,b,1,eps);if(abs(QLeft+QRight-QA)<=eps)q=QA;elseq=SubSmartSimpson(f,a,(a+b)/2,eps)+SubSmartSimpson(f,(a+b)/2,b ,eps);ende.g>> SmartSimpson('x*sin(x)',0,1)ans =0.3011>>。
matlab交集和并集
matlab交集和并集Matlab是一种高级计算机语言和环境，常用于科学计算、数据分析和工程应用。
其中，交集和并集是两个常用的集合操作，可以在Matlab中通过多种方式进行计算和操作。
本文将介绍Matlab中交集和并集的定义、符号表示以及使用示例。
一、交集在数学中，交集指的是两个或多个集合中共有的元素构成的新集合。
在Matlab中，可以使用“&”符号来表示交集操作。
示例1：计算两个向量的交集```matlabA = [1 2 3 4];B = [3 4 5 6];C = A & B;disp(C);```在上述示例中，A和B分别代表两个向量，将A和B进行交集操作后，将结果赋值给向量C。
执行程序后，控制台会输出交集结果C。
在上述示例中，输出结果为[0 0 1 1]，表示在A和B中，下标为3和4的元素是共有的。
示例2：计算两个矩阵的交集```matlabA = [1 2 3; 4 5 6; 7 8 9];B = [4 5 6; 1 2 3; 7 8 9];C = A & B;disp(C);```在上述示例中，A和B分别代表两个矩阵，将A和B进行交集操作后，将结果赋值给矩阵C。
执行程序后，控制台会输出交集结果C。
在上述示例中，输出结果是一个逻辑矩阵，表示在A和B中，对应位置上的元素是否相同。
二、并集在数学中，并集指的是两个或多个集合中的所有元素组合而成的新集合。
在Matlab中，可以使用“|”符号来表示并集操作。
示例1：计算两个向量的并集```matlabA = [1 2 3 4];B = [3 4 5 6];C = A | B;disp(C);```在上述示例中，A和B分别代表两个向量，将A和B进行并集操作后，将结果赋值给向量C。
执行程序后，控制台会输出并集结果C。
在上述示例中，输出结果为[1 1 1 1 1 1]，表示将A和B中所有的元素放在一个向量中。
示例2：计算两个矩阵的并集```matlabA = [1 2 3; 4 5 6; 7 8 9];B = [4 5 6; 1 2 3; 7 8 9];C = A | B;disp(C);```在上述示例中，A和B分别代表两个矩阵，将A和B进行并集操作后，将结果赋值给矩阵C。
Matlab常用命令集合
Matlab常用命令集合1.常用命令-->管理命令和函数2.常用命令-->管理变量和工作区(输入输出、内存管理等)3.常用命令-->管理命令控制窗口(command窗口)4.常用命令-->文件和工作环境5.常用命令-->启动和退出Matlab6.语言结构和调试命令-->程序设计7.语言结构和调试命令-->流程控制8.语言结构和调试命令-->交互输入9.语言结构和调试命令-->面向对象编程10.语言结构和调试命令-->程序调试11.语言结构和调试命令-->lasterr,lastwarn12.创建图形用户界面-->对话框13.矩阵函数－数值线性代数-->矩阵分析14.矩阵函数－数值线性代数-->线性方程15.矩阵函数－数值线性代数-->特征值和奇异值16.矩阵函数－数值线性代数-->矩阵函数17.矩阵函数－数值线性代数-->低级函数18.创建图形用户界面-->用户界面对象19.创建图形用户界面-->其他函数20.创建图形用户界面-->其他函数21.图象可视化函数-->三维绘图函数22.图象可视化函数-->绘制标注和网络23.图象可视化函数-->体数据可视化24.图象可视化函数-->表面、网格和轮廓绘制25.图象可视化函数-->域生成26.图象可视化函数-->专门图形绘制27.图象可视化函数-->视觉控制28.图象可视化函数-->颜色操作29.图象可视化函数-->打印函数30.图象可视化函数-->图形图象处理31.双重函数和非线性数值方法32.多项式和插值函数-->多项式33.多项式和插值函数-->数据插值34.多项式命令补遗35.数组函数-->单元数组函数36.数组函数-->多维数组函数37.位操作函数、结构和对象函数38.底层文件输入输出函数39.底层文件输入输出函数-->特殊函数40.字符串函数-->常用函数41.字符串函数-->字符串操作42.字符串函数-->字符串和数值的转换43.字符串函数-->基转换44.声音处理函数-->常用声音函数45.声音处理函数-->特殊声音函数46.声音处理函数-->WAV声音文件47.稀疏矩阵函数-->基本稀疏矩阵48.稀疏矩阵函数-->满阵和稀疏矩阵的转换49.稀疏矩阵函数-->稀疏矩阵的非零元素操作50.稀疏矩阵函数-->稀疏矩阵的可视化51.稀疏矩阵函数-->排序算法52.稀疏矩阵函数-->范数、条件数和秩53.稀疏矩阵函数-->线性方程的稀疏系统54.稀疏矩阵函数-->线性方程的稀疏系统55.稀疏矩阵函数-->杂项函数56.数据分析和傅立叶变换-->基本运算57.数据分析和傅立叶变换-->有限差分58.数据分析和傅立叶变换-->相关59.数据分析和傅立叶变换-->滤波和卷积60.数据分析和傅立叶变换-->傅立叶变换61.数据分析和傅立叶变换-->向量函数62.数学函数和坐标变换-->基本数学函数63.数学函数和坐标变换-->特殊函数64.数学函数和坐标变换-->坐标变换65.矩阵和矩阵操作基础-->矩阵和数组基础66.矩阵和矩阵操作基础-->特殊变量和常量67.矩阵和矩阵操作基础-->时间和日期68.矩阵和矩阵操作基础-->矩阵操作69.矩阵和矩阵操作基础-->特殊函数矩阵70.运算符和逻辑函数-->算数运算符71.运算符和逻辑函数-->关系运算符72.运算符和逻辑函数-->特殊函数73.运算符和逻辑函数-->逻辑函数74.simulink的命令集75.rtw相关的命令1．常用命令-->管理命令和函数addpath 添加目录到MATLAB搜索路径doc 在Web浏览器上现实HTML文档help 显示Matlab命令和M文件的在线帮助helpwinhelpdesk helplookfor 在基于Matlab搜索路径的所有M文件中搜索关键字partialpath 部分路径名path 所有关于路径名的处理pathtool 一个不错的窗口路径处理界面rmpath 删除搜索路径中指定目录type 显示指定文件的内容ver 版本信息version 版本号web 打开web页what 列出当前目录吓所有的M文件 Mat文件和 Mex文件whatsnew 显示readme文件which 显示文件位置（返回）2．常用命令-->管理变量和工作区(输入输出、内存管理等) clear 从内存中删除disp 显示文本或数组内容length 数组长度(最长维数)load 重新载入变量(从磁盘上)mlock 锁定文件，防止文件被错误删除munlock 解锁文件openvar 在数组编辑器中打开变量pack 整理内存空间save 保存变量到文件 8*)size 数组维数who 列出内存变量whos 列出内存变量,同时显示变量维数workspace 显示工作空间窗口（返回）3．常用命令-->管理命令控制窗口(command窗口)clc 清空命令窗口echo 禁止或允许显示执行过程format 设置输出显示格式home 光标移动到命令窗口左上角more 设置命令窗口页输出格式（返回）4．常用命令-->文件和工作环境cd 改变工作目录copyfile 复制文件delete 删除文件和图形对象diary 把命令窗口的人机交互保存到文件dir 显示目录edit 编辑文本文件fileparts 返回文件的各个部分fullfile 使用指定部分建立文件全名inmem 返回内存(伪代码区)的matlab函数名ls 在unix系统中列出目录(win中亦可)matlabroot 根目录mkdir 新建目录open 打开文件pwd 显示当前目录tempdir 返回系统临时目录的名字tempname 随机给出一个临时字符串(可用作文件名)! 直接调用操作系统command命令（返回）5．常用命令-->启动和退出Matlabmatlabrc Matlab的启动M文件exit 退出Matlabquit 退出Matlabstartup 运行matlab启动文件（返回）6．语言结构和调试命令-->程序设计builtin 从可重载方法中调用内置函数eval 执行包含可执行表达式的字符串evalc 计算并返回表达式的值evalin 执行某个工作空间中的包含表达式的字符串feval 执行函数(从函数名或函数句柄)function 函数头global 定义全局变量nargchk 检查输入参数数目persistent 定义常量script 作为脚本的M文件（返回）7．语言结构和调试命令-->流程控制break 停止执行循环case case语句,switch语句一部分switch switch_exprcase case_exprstatmentscaseotherwiseendcatch try-catch语句一部分，捕捉程序else if条件语句一部分elseif if条件语句一部分end 终止for、while、switch、try和if语句error 显示错误信息for 循环语句一(确定次数)if if条件语句other switch语句一部分return 返回到调用函数switch 开关语句try try程序块warning 类似于disp，但可被禁止while 循环语句(次数不确定)（返回）8．语言结构和调试命令-->交互输入input 交互输入keyboard 在M文件中遇到keyboard时将在命令窗口产生交互直到输入return命令menu 为用户输入产生一个选择菜单pause 暂停（返回）9．语言结构和调试命令-->面向对象编程class 创建一个对象或者返回一个对象类double 转换为双精度inferiorto 亚类关系inline 创建一个内联函数int8,int16,int32 转换到符号整数isa 检查是否为所给类的对象loadobj load函数用户定义扩展saveobj save函数用户定义扩展single 转换为单精度superiorto 超类关系uint8,uint16,uint32 转换到无符号整数（返回）10语言结构和调试命令-->程序调试dbclear 断点清除dbcont 重新开始执行dbdown 改变当前工作空间dbmex 调试Mex文件dbquit 退出调试模式dbstack 显示函数调用堆栈dbstatus 列出所有断点dbstep 从断点处执行dbstop 设置断点dbtype 列出带行号的M文件内容dbup 改变当前工作空间（返回）11．语言结构和调试命令-->lasterr,lastwarnlasterr 返回Matlab中产生的最后一个异常信息。
MATLAB操作命令大全
MATLAB操作命令大全1.基本操作：- clear: 清除工作区中的所有变量。
- clc: 清除命令窗口的内容。
- close all: 关闭所有图形窗口。
- help function-name: 显示与函数相关的帮助文档。
- who: 显示当前工作区中的所有变量。
- save file-name: 保存当前工作区中的所有变量到指定的文件。
- load file-name: 从文件中加载变量到当前工作区。
2.变量操作：-=:赋值操作符，将右边的值赋给左边的变量。
-+:加法操作符。
--:减法操作符。
-*:乘法操作符。
-/:除法操作符。
-^:幂运算操作符。
- sqrt(x): 计算 x 的平方根。
- abs(x): 计算 x 的绝对值。
- max(x): 返回 x 中的最大值。
- min(x): 返回 x 中的最小值。
- sum(x): 计算 x 中所有元素的和。
3.数组操作：- zeros(m, n): 创建一个 m 行 n 列的全零数组。
- ones(m, n): 创建一个 m 行 n 列的全一数组。
- eye(n): 创建一个 n 行 n 列的单位矩阵。
- size(x): 返回 x 的维度。
- length(x): 返回 x 的长度。
- reshape(x, m, n): 将 x 重新排列为一个 m 行 n 列矩阵。
- transpose(x): 将 x 的行和列互换。
4.控制流程：- if-else: 条件语句，根据条件执行不同的代码块。
- for loop: 循环语句，执行指定次数的代码块。
- while loop: 循环语句，根据条件反复执行代码块。
- break: 在循环中使用，用来跳出当前循环。
- continue: 在循环中使用，用来跳过当前循环的剩余部分。
5.统计分析：- mean(x): 计算 x 的平均值。
- median(x): 计算 x 的中位数。
- std(x): 计算 x 的标准差。
MATLAB命令与程序清单
matlab命令与程序清单(matlab commands and functions list)A aabs 绝对值、模、字符的ASCII码值acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割align 启动图形对象几何位置排列工具all 所有元素非零为真angle 相角ans 表达式计算结果的缺省变量名any 所有元素非全零为真area 面域图argnames 函数M文件宗量名asec 反正割asech 反双曲正割asin 反正弦asinh 反双曲正弦assignin 向变量赋值atan 反正切atan2 四象限反正切atanh 反双曲正切autumn 红黄调秋色图阵axes 创建轴对象的低层指令axis 控制轴刻度和风格的高层指令B bbar 二维直方图bar3 三维直方图bar3h 三维水平直方图barh 二维水平直方图base2dec X进制转换为十进制bin2dec 二进制转换为十进制blanks 创建空格串bone 蓝色调黑白色图阵box 框状坐标轴break while 或for 环中断指令brighten 亮度控制C ccapture （3版以前）捕获当前图形cart2pol 直角坐标变为极或柱坐标cart2sph 直角坐标变为球坐标cat 串接成高维数组caxis 色标尺刻度cd 指定当前目录cdedit 启动用户菜单、控件回调函数设计工具cdf2rdf 复数特征值对角阵转为实数块对角阵ceil 向正无穷取整cell 创建元胞数组cell2struct 元胞数组转换为构架数组celldisp 显示元胞数组内容cellplot 元胞数组内部结构图示char 把数值、符号、内联类转换为字符对象chi2cdf 分布累计概率函数chi2inv 分布逆累计概率函数chi2pdf 分布概率密度函数chi2rnd 分布随机数发生器chol Cholesky分解clabel 等位线标识cla 清除当前轴class 获知对象类别或创建对象clc 清除指令窗clear 清除内存变量和函数clf 清除图对象clock 时钟colorcube 三浓淡多彩交叉色图矩阵colordef 设置色彩缺省值colormap 色图colspace 列空间的基close 关闭指定窗口colperm 列排序置换向量comet 彗星状轨迹图comet3 三维彗星轨迹图compass 射线图compose 求复合函数cond （逆）条件数condeig 计算特征值、特征向量同时给出条件数condest 范–1条件数估计conj 复数共轭contour 等位线contourf 填色等位线contour3 三维等位线contourslice 四维切片等位线图conv 多项式乘、卷积cool 青紫调冷色图copper 古铜调色图cos 余弦cosh 双曲余弦cot 余切coth 双曲余切cplxpair 复数共轭成对排列csc 余割csch 双曲余割cumsum 元素累计和cumtrapz 累计梯形积分cylinder 创建圆柱D ddblquad 二重数值积分deal 分配宗量deblank 删去串尾部的空格符dec2base 十进制转换为X进制dec2bin 十进制转换为二进制dec2hex 十进制转换为十六进制deconv 多项式除、解卷delaunay Delaunay 三角剖分del2 离散Laplacian差分demo Matlab演示det 行列式diag 矩阵对角元素提取、创建对角阵diary Matlab指令窗文本内容记录diff 数值差分、符号微分digits 符号计算中设置符号数值的精度dir 目录列表disp 显示数组display 显示对象内容的重载函数dlinmod 离散系统的线性化模型dmperm 矩阵Dulmage-Mendelsohn 分解dos 执行DOS 指令并返回结果double 把其他类型对象转换为双精度数值drawnow 更新事件队列强迫Matlab 刷新屏幕dsolve 符号计算解微分方程E eecho M文件被执行指令的显示edit 启动M文件编辑器eig 求特征值和特征向量eigs 求指定的几个特征值end 控制流FOR等结构体的结尾元素下标eps 浮点相对精度error 显示出错信息并中断执行errortrap错误发生后程序是否继续执行的控制erf 误差函数erfc 误差补函数erfcx 刻度误差补函数erfinv 逆误差函数errorbar 带误差限的曲线图etreeplot画消去树eval 串演算指令evalin 跨空间串演算指令exist 检查变量或函数是否已定义exit 退出Matlab环境exp 指数函数expand 符号计算中的展开操作expint 指数积分函数expm 常用矩阵指数函数expm1 Pade法求矩阵指数expm2 Taylor法求矩阵指数expm3 特征值分解法求矩阵指数eye 单位阵ezcontour 画等位线的简捷指令ezcontourf 画填色等位线的简捷指令ezgraph3 画表面图的通用简捷指令ezmesh 画网线图的简捷指令ezmeshc 画带等位线的网线图的简捷指令ezplot 画二维曲线的简捷指令ezplot3 画三维曲线的简捷指令ezpolar 画极坐标图的简捷指令ezsurf 画表面图的简捷指令ezsurfc 画带等位线的表面图的简捷指令F ffactor 符号计算的因式分解feather 羽毛图feedback 反馈连接feval 执行由串指定的函数fft 离散Fourier变换fft2 二维离散Fourier变换fftn 高维离散Fourier变换fftshift 直流分量对中的谱fieldnames 构架域名figure 创建图形窗fill3 三维多边形填色图find 寻找非零元素下标findobj 寻找具有指定属性的对象图柄findstr 寻找短串的起始字符下标findsym 机器确定内存中的符号变量finverse 符号计算中求反函数fix 向零取整flag 红白蓝黑交错色图阵fliplr 矩阵的左右翻转flipud 矩阵的上下翻转flipdim 矩阵沿指定维翻转floor 向负无穷取整flops 浮点运算次数flow Matlab提供的演示数据fmin 求单变量非线性函数极小值点（旧版）fminbnd 求单变量非线性函数极小值点fmins 单纯形法求多变量函数极小值点（旧版）fminunc 拟牛顿法求多变量函数极小值点fminsearch 单纯形法求多变量函数极小值点fnder 对样条函数求导fnint 利用样条函数求积分fnval 计算样条函数区间内任意一点的值fnplt 绘制样条函数图形fopen 打开外部文件for 构成for环用format 设置输出格式fourier Fourier 变换fplot 返函绘图指令fprintf 设置显示格式fread 从文件读二进制数据fsolve 求多元函数的零点full 把稀疏矩阵转换为非稀疏阵funm 计算一般矩阵函数funtool 函数计算器图形用户界面fzero 求单变量非线性函数的零点G ggamma 函数gammainc 不完全函数gammaln 函数的对数gca 获得当前轴句柄gcbo 获得正执行“回调”的对象句柄gcf 获得当前图对象句柄gco 获得当前对象句柄geomean 几何平均值get 获知对象属性getfield 获知构架数组的域getframe 获取影片的帧画面ginput 从图形窗获取数据global 定义全局变量gplot 依图论法则画图gradient 近似梯度gray 黑白灰度grid 画分格线griddata 规则化数据和曲面拟合gtext 由鼠标放置注释文字guide 启动图形用户界面交互设计工具H hharmmean 调和平均值help 在线帮助helpwin 交互式在线帮助helpdesk 打开超文本形式用户指南hex2dec 十六进制转换为十进制hex2num 十六进制转换为浮点数hidden 透视和消隐开关hilb Hilbert矩阵hist 频数计算或频数直方图histc 端点定位频数直方图histfit 带正态拟合的频数直方图hold 当前图上重画的切换开关horner 分解成嵌套形式hot 黑红黄白色图hsv 饱和色图I iif-else-elseif 条件分支结构ifft 离散Fourier反变换ifft2 二维离散Fourier反变换ifftn 高维离散Fourier反变换ifftshift 直流分量对中的谱的反操作ifourier Fourier反变换i, j 缺省的“虚单元”变量ilaplace Laplace反变换imag 复数虚部image 显示图象imagesc 显示亮度图象imfinfo 获取图形文件信息imread 从文件读取图象imwrite 把图象写成文件ind2sub 单下标转变为多下标inf 无穷大info MathWorks公司网点地址inline 构造内联函数对象inmem 列出内存中的函数名input 提示用户输入inputname 输入宗量名int 符号积分int2str 把整数数组转换为串数组interp1 一维插值interp2 二维插值interp3 三维插值interpn N维插值interpft 利用FFT插值intro Matlab自带的入门引导inv 求矩阵逆invhilb Hilbert矩阵的准确逆ipermute 广义反转置isa 检测是否给定类的对象ischar 若是字符串则为真isequal 若两数组相同则为真isempty 若是空阵则为真isfinite 若全部元素都有限则为真isfield 若是构架域则为真isglobal 若是全局变量则为真ishandle 若是图形句柄则为真ishold 若当前图形处于保留状态则为真isieee 若计算机执行IEEE规则则为真isinf 若是无穷数据则为真isletter 若是英文字母则为真islogical 若是逻辑数组则为真ismember 检查是否属于指定集isnan 若是非数则为真isnumeric 若是数值数组则为真isobject 若是对象则为真isprime 若是质数则为真isreal 若是实数则为真isspace 若是空格则为真issparse 若是稀疏矩阵则为真isstruct 若是构架则为真isstudent 若是Matlab学生版则为真iztrans 符号计算Z反变换J j , K kjacobian 符号计算中求Jacobian 矩阵jet 蓝头红尾饱和色jordan 符号计算中获得Jordan标准型keyboard 键盘获得控制权kron Kronecker乘法规则产生的数组L llaplace Laplace变换lasterr 显示最新出错信息lastwarn 显示最新警告信息leastsq 解非线性最小二乘问题（旧版）legend 图形图例lighting 照明模式line 创建线对象lines 采用plot 画线色linmod 获连续系统的线性化模型linmod2 获连续系统的线性化精良模型linspace 线性等分向量ln 矩阵自然对数load 从MAT文件读取变量log 自然对数log10 常用对数log2 底为2的对数loglog 双对数刻度图形logm 矩阵对数logspace对数分度向量lookfor 按关键字搜索M文件lower 转换为小写字母lsqnonlin 解非线性最小二乘问题lu LU分解M mmad 平均绝对值偏差magic 魔方阵maple 运作Maple格式指令mat2str 把数值数组转换成输入形态串数组material 材料反射模式max 找向量中最大元素mbuild 产生EXE文件编译环境的预设置指令mcc 创建MEX或EXE文件的编译指令mean 求向量元素的平均值median 求中位数menuedit 启动设计用户菜单的交互式编辑工具mesh 网线图meshz 垂帘网线图meshgrid 产生“格点”矩阵methods 获知对指定类定义的所有方法函数mex 产生MEX文件编译环境的预设置指令mfunlis 能被mfun计算的MAPLE经典函数列表mhelp 引出Maple的在线帮助min 找向量中最小元素mkdir 创建目录mkpp 逐段多项式数据的明晰化mod 模运算more 指令窗中内容的分页显示movie 放映影片动画moviein 影片帧画面的内存预置mtaylor 符号计算多变量Taylor级数展开N nndims 求数组维数NaN 非数（预定义）变量nargchk 输入宗量数验证nargin 函数输入宗量数nargout 函数输出宗量数ndgrid 产生高维格点矩阵newplot 准备新的缺省图、轴nextpow2 取最接近的较大2次幂nnz 矩阵的非零元素总数nonzeros 矩阵的非零元素norm 矩阵或向量范数normcdf 正态分布累计概率密度函数normest 估计矩阵2范数norminv 正态分布逆累计概率密度函数normpdf 正态分布概率密度函数normrnd 正态随机数发生器notebook 启动Matlab和Word的集成环境null 零空间num2str 把非整数数组转换为串numden 获取最小公分母和相应的分子表达式nzmax 指定存放非零元素所需内存O oode1 非Stiff 微分方程变步长解算器ode15s Stiff 微分方程变步长解算器ode23t 适度Stiff 微分方程解算器ode23tb Stiff 微分方程解算器ode45 非Stiff 微分方程变步长解算器odefile ODE 文件模板odeget 获知ODE 选项设置参数odephas2 ODE 输出函数的二维相平面图odephas3 ODE 输出函数的三维相空间图odeplot ODE 输出函数的时间轨迹图odeprint 在Matlab指令窗显示结果odeset 创建或改写ODE选项构架参数值ones 全1数组optimset 创建或改写优化泛函指令的选项参数值orient 设定图形的排放方式orth 值空间正交化P ppack 收集Matlab内存碎块扩大内存pagedlg 调出图形排版对话框patch 创建块对象path 设置Matlab搜索路径的指令pathtool 搜索路径管理器pause 暂停pcode 创建预解译P码文件pcolor 伪彩图peaks Matlab提供的典型三维曲面permute 广义转置pi （预定义变量）圆周率pie 二维饼图pie3 三维饼图pink 粉红色图矩阵pinv 伪逆plot 平面线图plot3 三维线图plotmatrix 矩阵的散点图plotyy 双纵坐标图poissinv 泊松分布逆累计概率分布函数poissrnd 泊松分布随机数发生器pol2cart 极或柱坐标变为直角坐标polar 极坐标图poly 矩阵的特征多项式、根集对应的多项式poly2str 以习惯方式显示多项式poly2sym 双精度多项式系数转变为向量符号多项式polyder 多项式导数polyfit 数据的多项式拟合polyval 计算多项式的值polyvalm 计算矩阵多项式pow2 2的幂ppval 计算分段多项式pretty 以习惯方式显示符号表达式print 打印图形或SIMULINK模型printsys 以习惯方式显示有理分式prism 光谱色图矩阵procread 向MAPLE输送计算程序profile 函数文件性能评估器propedit 图形对象属性编辑器pwd 显示当前工作目录Q qquad 低阶法计算数值积分quad8 高阶法计算数值积分quit 推出Matlab 环境quiver 二维方向箭头图quiver3 三维方向箭头图R rrand 产生均匀分布随机数randn 产生正态分布随机数randperm 随机置换向量range 样本极差rank 矩阵的秩rats 有理输出rcond 矩阵倒条件数估计real 复数的实部reallog 在实数域内计算自然对数realpow 在实数域内计算乘方realsqrt 在实数域内计算平方根realmax 最大正浮点数realmin 最小正浮点数rectangle 画“长方框”rem 求余数repmat 铺放模块数组reshape 改变数组维数、大小residue 部分分式展开return 返回ribbon 把二维曲线画成三维彩带图rmfield 删去构架的域roots 求多项式的根rose 数扇形图rot90 矩阵旋转90度rotate 指定的原点和方向旋转rotate3d 启动三维图形视角的交互设置功能round 向最近整数圆整rref 简化矩阵为梯形形式rsf2csf 实数块对角阵转为复数特征值对角阵rsums Riemann和S ssave 把内存变量保存为文件scatter 散点图scatter3 三维散点图sec 正割sech 双曲正割semilogx X轴对数刻度坐标图semilogy Y轴对数刻度坐标图series 串联连接set 设置图形对象属性setfield 设置构架数组的域setstr 将ASCII码转换为字符的旧版指令sign 根据符号取值函数signum 符号计算中的符号取值函数sim 运行SIMULINK模型simget 获取SIMULINK模型设置的仿真参数simple 寻找最短形式的符号解simplify 符号计算中进行简化操作simset 对SIMULINK模型的仿真参数进行设置simulink 启动SIMULINK模块库浏览器sin 正弦sinh 双曲正弦size 矩阵的大小slice 立体切片图solve 求代数方程的符号解spalloc 为非零元素配置内存sparse 创建稀疏矩阵spconvert 把外部数据转换为稀疏矩阵spdiags 稀疏对角阵spfun 求非零元素的函数值sph2cart 球坐标变为直角坐标sphere 产生球面spinmap 色图彩色的周期变化spline 样条插值spones 用1置换非零元素sprandsym 稀疏随机对称阵sprank 结构秩spring 紫黄调春色图sprintf 把格式数据写成串spy 画稀疏结构图sqrt 平方根sqrtm 方根矩阵squeeze 删去大小为1的“孤维”sscanf 按指定格式读串stairs 阶梯图std 标准差stem 二维杆图step 阶跃响应指令str2double 串转换为双精度值str2mat 创建多行串数组str2num 串转换为数strcat 接成长串strcmp 串比较strjust 串对齐strmatch 搜索指定串strncmp 串中前若干字符比较strrep 串替换strtok 寻找第一间隔符前的内容struct 创建构架数组struct2cell 把构架转换为元胞数组strvcat 创建多行串数组sub2ind 多下标转换为单下标subexpr 通过子表达式重写符号对象subplot 创建子图subs 符号计算中的符号变量置换subspace 两子空间夹角sum 元素和summer 绿黄调夏色图superiorto 设定优先级surf 三维着色表面图surface 创建面对象surfc 带等位线的表面图surfl 带光照的三维表面图surfnorm 空间表面的法线svd 奇异值分解svds 求指定的若干奇异值switch-case-otherwise 多分支结构sym2poly 符号多项式转变为双精度多项式系数向量symmmd 对称最小度排序symrcm 反向Cuthill-McKee排序syms 创建多个符号对象T ttan 正切tanh 双曲正切taylortool 进行Taylor逼近分析的交互界面text 文字注释tf 创建传递函数对象tic 启动计时器title 图名toc 关闭计时器trapz 梯形法数值积分treelayout 展开树、林treeplot 画树图tril 下三角阵trim 求系统平衡点trimesh 不规则格点网线图trisurf 不规则格点表面图triu 上三角阵try-catch 控制流中的Try-catch结构type 显示M文件U uuicontextmenu创建现场菜单uicontrol 创建用户控件uimenu 创建用户菜单unmkpp 逐段多项式数据的反明晰化unwrap 自然态相角upper 转换为大写字母V vvar 方差varargin 变长度输入宗量varargout 变长度输出宗量vectorize 使串表达式或内联函数适于数组运算ver 版本信息的获取view 三维图形的视角控制voronoi V oronoi多边形vpa 任意精度（符号类）数值W wwarning 显示警告信息what 列出当前目录上的文件whatsnew 显示Matlab中Readme文件的内容which 确定函数、文件的位置while 控制流中的While环结构white 全白色图矩阵whitebg 指定轴的背景色who 列出内存中的变量名whos 列出内存中变量的详细信息winter 蓝绿调冬色图workspace 启动内存浏览器X x , Y y , Z zxlabel X轴名xor 或非逻辑yesinput 智能输入指令ylabel Y轴名zeros 全零数组zlabel Z轴名zoom 图形的变焦放大和缩小ztrans 符号计算Z变换。
【精品】MATLAB程序应用大全精编版
1.3 MATLAB 第六版外观
3. Current Directory Browser 当前目录浏览器窗口默认地出现在MATLAB界面左上方后台，用来设置当前目录 ，可以随时显示当前目录下的MATLAB文件的信息，包括文件名、文件类型、最 后修改时间和文件说明等。并可以复制、编辑和运行M文件，以及加载MAT类型 的数据文件。 4. Workspace Browser 工作空间浏览器窗口在MATLAB界面左上方后台，用来显示MATLAB工作空间中 的变量名、数据结构、类型、大小和在内存中的字节数。在该窗口中，还可以对 工作空间中的变量进行观察、编辑、提取和保存。
SIMULINK库浏览器：SIMULINK是一个仿真工具，它不是MATLAB的一个部分 ，所有如果您没有安装SIMULINK，就不会有此图标。
在线帮助：此图标打开MATLAB的在线帮助窗口，您可以一路单击所要的类别或 命令，直到得到所需要的帮助信息。
中北大学 信息与通信工程学院
1.3 MATLAB 第六版外观
中北大学 信息与通信工程学院
2 初探MATLAB
2.1 使用变量与基本运算 2.2 向量与矩阵的处理 2.3 常用数学函数 2.4 程序流程控制
2.5 M文件 2.6 搜索路径 2.7 工作空间与变量
本章介绍 MATLAB 的基本环境和操作，如果您是MATLAB的初学者，建议您先 熟悉本章各节的内容，能够帮助您很快进入状态。当然，如果您是MATLAB老手 ，可以直接跳到后面各章节。
[A; 1 1] %在A中加入第4行
ans=
1
3
4
6
7
9
1
1
%取出矩阵的一部分形成一个新矩阵
B=A(1:2,1:2)
MATLAB简单程序大全
MATLAB简单程序大全求特征值特征向量A=[2 3 4;1 5 9;8 5 2]det(A)A'rank(A)inv(A)rref(A)eig(A)%求特征值和特征向量卫星运行问题h=200,H=51000,R=6378;a=(h+H+2*R)/2;c=(H-h)/2;b=(a^2-c^2)^(1/2);e=c/a;f=sqrt(1-exp(2).*cos(t)^2);l=int(f,t,0,pi/2)L=4*a.*l动态玫瑰线n=3;N=10000;theta=2*pi*(0:N)/N;r=cos(n*theta);x=r.*cos(theta);y=r.*sin(theta);comet(x,y)二重积分syms x yf=x^2*sin(y);int(int(f,x,0,1),y,0,pi)ezmesh(f,[0,1,0,pi])函数画图syms x;f=exp(-0.2*x)*sin(0.5*x);ezplot(f,[0,8*pi])玫瑰线theta=0:0.01:2*pi;r=cos(3*theta);polar(theta,r,'r')求x^2+y^2=1和x^2+z^2=1所围成的体积syms x y z Rr=1;Z=sqrt(1-x^2);y0=Z;V=8*int(int(Z,y,0,y0),x,0,1)求导数及图像f='1/(5+4*cos(x))';subplot(1,2,1);ezplot(f)f1=diff(f)subplot(1,2,2);ezplot(f1)绕x轴旋转t=(0:20)*pi/10;r=exp(-.2*t).*sin(.5*t);theta=t;x=t'*ones(size(t));y=r'*cos(theta);z=r'*sin(theta);mesh(x,y,z)colormap([0 0 0])某年是否闰年year=input('input year:=');n1=year/4;n2=year/100;n3=year/400;if n1==fix(n1)&n2~=fix(n2)disp('是闰年')elseif n1==fix(n1)&n3==fix(n3) disp('是闰年')elsedisp('不是闰年')End玫瑰线的绘制theta=0:0.001:2*pi;rho=2*cos(3*theta);figure(1);polar(theta,rho)x=rho.*cos(theta);y=rho.*sin(theta);figure(2);comet(x,y)相遇问题function k=moto(A,B)if nargin==0,A=0;B=100;end va=10;vb=8;vc=60;f=1;k=0;while(B-A)>0.2if f==1tk=(B-A)/(vb+vc);elsetk=(B-A)/(vc+va);endA=A+va*tk;B=B-vb*tk;f=-f;k=k+1;End数学实验syms x; x=1:100;f(x)=exp(-0.2*x).*sin(0.5*x); syms x; x=1:100;f(x)=exp(-0.2*x).*sin(0.5*x); plot(x,f(x))syms x; x=0:0.1:2*pi;f(x)=exp(-0.2*x)*sin(0.5*x);plot(x,f(x))syms x; x=0:0.1:2*pi;f(x)=exp(-0.2*x).*sin(0.5*x);plot(x,f(x))syms x; x=0:0.1:2*pi;f(x)=exp(-0.2*x).*sin(0.5*x); zplot(x,f(x))clcf=sym('sin(x)/x');limit(f)f1=sym('(1+a/x)^x')limit(f1)f1=sym('(1+a/x)^x'),x=inf;limit(f1)clcf1=sym('(1+1/x)^x');limit(f1,'inf')f1=sym('(1+1/x)^x');limit(f1,'inf')clcf1=sym('(1+1/x)^x');limit(f1,'x',inf)Clc相遇问题A=0;B=100;va=10;vb=8;vc=60;f=1;k=0;plot(A,0,'ro',B,0,'go'),hold on while(B-A)>0.2if f==1tk=(B-A)/(vb+vc);elsetk=(B-A)/(vc+va);endA=A+va*tk;B=B-vb*tk;plot(A,0,'R.',B,0,'g.'),pause(1)f=-f;k=k+1;endk,tk,A,B数学实验f=[0 9 13.5;0.1 0 0;0 0.2 0]x=[100;100;100]y=xfor n=1:3x=f*xy=[y,x]endyfigure(1),bar(y(1,:))figure(2),bar(y(2,:))figure(3),bar(y(3,:))[p,d]=eig(f)a=p(:,1)b=a/sum(a)x1=b*300y1=x1for n=1:3x1=f*x1y1=[y1,x1]endy1figure(4),bar(y1(1,:))figure(5),bar(y1(2,:))figure(6),bar(y1(3,:))杨辉三角function Y=yang(n)if nargin==0,n=3;endY=eye(n);Y(:,1)=ones(n,1);for k=3:nY(k,2:k-1)=Y(k-1,1:k-2)+Y(k-1,2:k-1);End实例2 寻找水仙花数一、实验内容数153有一个特殊的性质,即153=1^3+5^3+3^3,那么有没有其他这样的数呢？我们怎样利用MATLAB程序设计快速寻求这些水仙花数呢？二、实验分析把任一数各位数字的立方和求出,再求所得的各位数字的立方和,依次进行下去,总可以得到一个水仙花数.如把76按这样的运算进行下去:333333333333+=++=++=+++= 76559,559979,9791801,1801514, 333333333333++=++=++=++= 514190,190730,730370,370370结果得到一个水仙花数370。
Matlab100个实例程序
程序代码：(代码标记[code]...[/code] ) 1-32是：图形应用篇33-66是：界面设计篇67-84是：图形处理篇85-100是：数值分析篇实例1：三角函数曲线（1）function shili01h0=figure('toolbar','none',...'position',[198****0300],...'name','实例01');h1=axes('parent',h0,...'visible','off');x=-pi:0.05:pi;y=sin(x);plot(x,y);xlabel('自变量X');ylabel('函数值Y');title('SIN( )函数曲线');grid on实例2：三角函数曲线（2）function shili02h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例02');x=-pi:0.05:pi;y=sin(x)+cos(x);plot(x,y,'-*r','linewidth',1);grid onxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例3：图形的叠加function shili03h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例03');x=-pi:0.05:pi;y1=sin(x);y2=cos(x);plot(x,y1,...'-*r',...x,y2,...'--og');grid onxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例4：双y轴图形的绘制function shili04h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例04');x=0:900;a=1000;b=0.005;y1=2*x;y2=cos(b*x);[haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot'); axes(haxes(1))ylabel('semilog plot');axes(haxes(2))ylabel('linear plot');实例5：单个轴窗口显示多个图形function shili05h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例05');t=0:pi/10:2*pi;[x,y]=meshgrid(t);subplot(2,2,1)plot(sin(t),cos(t))axis equalsubplot(2,2,2)z=sin(x)-cos(y);plot(t,z)axis([0 2*pi -2 2])subplot(2,2,3)h=sin(x)+cos(y);plot(t,h)axis([0 2*pi -2 2])subplot(2,2,4)g=(sin(x).^2)-(cos(y).^2);plot(t,g)axis([0 2*pi -1 1])实例6：图形标注function shili06h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例06');t=0:pi/10:2*pi;h=plot(t,sin(t));xlabel('t=0到2\pi','fontsize',16);ylabel('sin(t)','fontsize',16);title('\it{从0to2\pi 的正弦曲线}','fontsize',16) x=get(h,'xdata');y=get(h,'ydata');imin=find(min(y)==y);imax=find(max(y)==y);text(x(imin),y(imin),...['\leftarrow最小值=',num2str(y(imin))],... 'fontsize',16)text(x(imax),y(imax),...['\leftarrow最大值=',num2str(y(imax))],...'fontsize',16)实例7：条形图形function shili07h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例07');tiao1=[562 548 224 545 41 445 745 512];tiao2=[47 48 57 58 54 52 65 48];t=0:7;bar(t,tiao1)xlabel('X轴');ylabel('TIAO1值');h1=gca;h2=axes('position',get(h1,'position'));plot(t,tiao2,'linewidth',3)set(h2,'yaxislocation','right','color','none','xticklabel',[])实例8：区域图形function shili08h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例08');x=91:95;profits1=[88 75 84 93 77];profits2=[51 64 54 56 68];profits3=[42 54 34 25 24];profits4=[26 38 18 15 4];area(x,profits1,'facecolor',[0.5 0.9 0.6],...'edgecolor','b',...'linewidth',3)hold onarea(x,profits2,'facecolor',[0.9 0.85 0.7],...'edgecolor','y',...'linewidth',3)hold onarea(x,profits3,'facecolor',[0.3 0.6 0.7],... 'edgecolor','r',...'linewidth',3)hold onarea(x,profits4,'facecolor',[0.6 0.5 0.9],... 'edgecolor','m',...'linewidth',3)hold offset(gca,'xtick',[91:95])set(gca,'layer','top')gtext('\leftarrow第一季度销量')gtext('\leftarrow第二季度销量')gtext('\leftarrow第三季度销量')gtext('\leftarrow第四季度销量')xlabel('年','fontsize',16);ylabel('销售量','fontsize',16);实例9：饼图的绘制function shili09h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例09');t=[54 21 35;68 54 35;45 25 12;48 68 45;68 54 69];x=sum(t);h=pie(x);textobjs=findobj(h,'type','text');str1=get(textobjs,{'string'});val1=get(textobjs,{'extent'});oldext=cat(1,val1{:});names={'商品一：';'商品二：';'商品三：'}; str2=strcat(names,str1);set(textobjs,{'string'},str2)val2=get(textobjs,{'extent'});newext=cat(1,val2{:});offset=sign(oldext(:,1)).*(newext(:,3)-oldext(:,3))/2; pos=get(textobjs,{'position'});textpos=cat(1,pos{:});textpos(:,1)=textpos(:,1)+offset;set(textobjs,{'position'},num2cell(textpos,[3,2]))实例10：阶梯图function shili10h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例10');a=0.01;b=0.5;t=0:10;f=exp(-a*t).*sin(b*t);stairs(t,f)hold onplot(t,f,':*')hold offglabel='函数e^{-(\alpha*t)}sin\beta*t的阶梯图'; gtext(glabel,'fontsize',16)xlabel('t=0:10','fontsize',16)axis([0 10 -1.2 1.2])实例11：枝干图function shili11h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例11');x=0:pi/20:2*pi;y1=sin(x);y2=cos(x);h1=stem(x,y1+y2);hold onh2=plot(x,y1,'^r',x,y2,'*g');h3=[h1(1);h2];legend(h3,'y1+y2','y1=sin(x)','y2=cos(x)') xlabel('自变量X');ylabel('函数值Y');title('正弦函数与余弦函数的线性组合');实例12：罗盘图function shili12h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例12');winddirection=[54 24 65 84256 12 235 62125 324 34 254];windpower=[2 5 5 36 8 12 76 14 10 8];rdirection=winddirection*pi/180;[x,y]=pol2cart(rdirection,windpower); compass(x,y);desc={'风向和风力','北京气象台','10月1日0：00到','10月1日12：00'};gtext(desc)实例13：轮廓图function shili13h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例13');[th,r]=meshgrid((0:10:360)*pi/180,0:0.05:1); [x,y]=pol2cart(th,r);z=x+i*y;f=(z.^4-1).^(0.25);contour(x,y,abs(f),20)xlabel('实部','fontsize',16);ylabel('虚部','fontsize',16);h=polar([0 2*pi],[0 1]);delete(h)hold oncontour(x,y,abs(f),20)实例14：交互式图形function shili14h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例14');axis([0 10 0 10]);hold onx=[];y=[];n=0;disp('单击鼠标左键点取需要的点'); disp('单击鼠标右键点取最后一个点'); but=1;while but==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('单击鼠标左键点取下一个点'); x(n,1)=xi;y(n,1)=yi;endt=1:n;ts=1:0.1:n;xs=spline(t,x,ts);ys=spline(t,y,ts);plot(xs,ys,'r-');hold off实例15：变换的傅立叶函数曲线function shili15h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例15');axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren')h=uicontrol('style','slider','position',... [100 10 500 20],'min',1,'max',20) for j=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0 0 1 1]);movie(m,30)实例16：劳伦兹非线形方程的无序活动function shili15h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例15');axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren')h=uicontrol('style','slider','position',... [100 10 500 20],'min',1,'max',20) for j=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0 0 1 1]);movie(m,30)实例17：填充图function shili17h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例17');t=(1:2:15)*pi/8;x=sin(t);y=cos(t);fill(x,y,'r')axis square offtext(0,0,'STOP',...'color',[1 1 1],...'fontsize',50,...'horizontalalignment','center')实例18：条形图和阶梯形图function shili18h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例18');subplot(2,2,1)x=-3:0.2:3;y=exp(-x.*x);bar(x,y)title('2-D Bar Chart')subplot(2,2,2)x=-3:0.2:3;y=exp(-x.*x);bar3(x,y,'r')title('3-D Bar Chart')subplot(2,2,3)x=-3:0.2:3;y=exp(-x.*x);stairs(x,y)title('Stair Chart')subplot(2,2,4)x=-3:0.2:3;y=exp(-x.*x);barh(x,y)title('Horizontal Bar Chart')实例19：三维曲线图function shili19h0=figure('toolbar','none',...'position',[200 150 450 400],... 'name','实例19');subplot(2,1,1)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,y1,z1,x,y2,z2,x,y3,z3) grid onxlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure1:3-D Plot')subplot(2,1,2)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,z1,y1,x,z2,y2,x,z3,y3) grid onxlabel('X轴');zlabel('Z轴');title('Figure2:3-D Plot')实例20：图形的隐藏属性function shili20h0=figure('toolbar','none',...'position',[200 150 450 300],... 'name','实例20');subplot(1,2,1)[x,y,z]=sphere(10);mesh(x,y,z)axis offtitle('Figure1:Opaque')hidden onsubplot(1,2,2)[x,y,z]=sphere(10);mesh(x,y,z)axis offtitle('Figure2:Transparent') hidden off实例21：PEAKS函数曲线function shili21h0=figure('toolbar','none',...'position',[200 100 450 450],... 'name','实例21');[x,y,z]=peaks(30);subplot(2,1,1)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfc(x,y,z)xlabel('X轴');ylabel('Y轴');title('Figure1:surfc函数形成的曲面')subplot(2,1,2)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfl(x,y,z)xlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure2:surfl函数形成的曲面')实例22：片状图function shili22h0=figure('toolbar','none',...'position',[200 150 550 350],...'name','实例22');subplot(1,2,1)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trimesh(t,x,y,z)hidden offtitle('Figure1:Triangular Surface Plot');subplot(1,2,2)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trisurf(t,x,y,z)title('Figure1:Triangular Surface Plot');实例23：视角的调整function shili23h0=figure('toolbar','none',...'position',[200 150 450 350],... 'name','实例23');x=-5:0.5:5;[x,y]=meshgrid(x);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;subplot(2,2,1)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure1')view(-37.5,30)subplot(2,2,2)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure2')view(-37.5+90,30)subplot(2,2,3)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure3')view(-37.5,60)subplot(2,2,4)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure4')view(180,0)实例24：向量场的绘制function shili24h0=figure('toolbar','none',...'position',[200 150 450 350],... 'name','实例24');subplot(2,2,1)z=peaks;ribbon(z)title('Figure1')subplot(2,2,2)[x,y,z]=peaks(15);[dx,dy]=gradient(z,0.5,0.5); contour(x,y,z,10)hold onquiver(x,y,dx,dy)hold offtitle('Figure2')subplot(2,2,3)[x,y,z]=peaks(15);[nx,ny,nz]=surfnorm(x,y,z);surf(x,y,z)hold onquiver3(x,y,z,nx,ny,nz)hold offtitle('Figure3')subplot(2,2,4)x=rand(3,5);y=rand(3,5);z=rand(3,5);c=rand(3,5);fill3(x,y,z,c)grid ontitle('Figure4')实例25：灯光定位function shili25h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例25');vert=[1 1 1;1 2 1;2 2 1;2 1 1;1 1 2;12 2;2 2 2;2 1 2];fac=[1 2 3 4;2 6 7 3;4 3 7 8;15 8 4;1 2 6 5;5 6 7 8];grid offsphere(36)h=findobj('type','surface');set(h,'facelighting','phong',...'facecolor',...'interp',...'edgecolor',[0.4 0.4 0.4],...'backfacelighting',...'lit')hold onpatch('faces',fac,'vertices',vert,... 'facecolor','y');light('position',[1 3 2]);light('position',[-3 -1 3]); material shinyaxis vis3d offhold off实例26：柱状图function shili26h0=figure('toolbar','none',...'position',[200 50 450 450],...'name','实例26');subplot(2,1,1)x=[5 2 18 7 39 8 65 5 54 3 2];bar(x)xlabel('X轴');ylabel('Y轴');title('第一子图');subplot(2,1,2)y=[5 2 18 7 39 8 65 5 54 3 2];barh(y)xlabel('X轴');ylabel('Y轴');title('第二子图');实例27：设置照明方式function shili27h0=figure('toolbar','none',...'position',[200 150 450 350],... 'name','实例27');subplot(2,2,1)sphereshading flatcamlight leftcamlight rightlighting flatcolorbaraxis offtitle('Figure1')subplot(2,2,2)sphereshading flatcamlight leftcamlight rightlighting gouraudcolorbaraxis offtitle('Figure2')subplot(2,2,3)sphereshading interpcamlight rightcamlight leftlighting phongcolorbaraxis offtitle('Figure3')subplot(2,2,4)sphereshading flatcamlight leftcamlight rightlighting nonecolorbaraxis offtitle('Figure4')实例28：羽状图function shili28h0=figure('toolbar','none',...'position',[200 150 450 350],... 'name','实例28');subplot(2,1,1)alpha=90:-10:0;r=ones(size(alpha));m=alpha*pi/180;n=r*10;[u,v]=pol2cart(m,n);feather(u,v)title('羽状图')axis([0 20 0 10])subplot(2,1,2)t=0:0.5:10;x=0.05+i;y=exp(-x*t);feather(y)title('复数矩阵的羽状图')实例29：立体透视（1）function shili29h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例29');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2); grid onfor i=-2:0.5:2;h1=surf(linspace(-2,2,20),...linspace(-2,2,20),...zeros(20)+i);rotate(h1,[1 -1 1],30)dx=get(h1,'xdata');dy=get(h1,'ydata');dz=get(h1,'zdata');delete(h1)slice(x,y,z,v,[-2 2],2,-2)hold onslice(x,y,z,v,dx,dy,dz)hold offaxis tightview(-5,10)drawnowend实例30：立体透视（2）function shili30h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例30');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2); [dx,dy,dz]=cylinder;slice(x,y,z,v,[-2 2],2,-2)for i=-2:0.2:2h=surface(dx+i,dy,dz);rotate(h,[1 0 0],90)xp=get(h,'xdata');yp=get(h,'ydata');zp=get(h,'zdata');delete(h)hold onhs=slice(x,y,z,v,xp,yp,zp);axis tightxlim([-3 3])view(-10,35)drawnowdelete(hs)hold offend实例31：表面图形function shili31h0=figure('toolbar','none',...'position',[200 150 550 250],...'name','实例31');subplot(1,2,1)x=rand(100,1)*16-8;y=rand(100,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),33); ylin=linspace(min(y),max(y),33); [X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic'); mesh(X,Y,Z)axis tighthold onplot3(x,y,z,'.','Markersize',20)subplot(1,2,2)k=5;n=2^k-1;theta=pi*(-n:2:n)/n;phi=(pi/2)*(-n:2:n)'/n;X=cos(phi)*cos(theta);Y=cos(phi)*sin(theta);Z=sin(phi)*ones(size(theta)); colormap([0 0 0;1 1 1])C=hadamard(2^k);surf(X,Y,Z,C)axis square实例32：沿曲线移动的小球h0=figure('toolbar','none',...'position',[198****8468],... 'name','实例32');h1=axes('parent',h0,...'position',[0.15 0.45 0.7 0.5],... 'visible','on');t=0:pi/24:4*pi;y=sin(t);plot(t,y,'b')n=length(t);h=line('color',[0 0.5 0.5],...'linestyle','.',...'markersize',25,...'erasemode','xor');k1=uicontrol('parent',h0,...'style','pushbutton',...'position',[80 100 50 30],...'string','开始',...'callback',[...'i=1;',...'k=1;,',...'m=0;,',...'while 1,',...'if k==0,',...'break,',...'end,',...'if k~=0,',...'set(h,''xdata'',t(i),''ydata'',y(i)),',...'drawnow;,',...'i=i+1;,',...'if i>n,',...'m=m+1;,',...'i=1;,',...'end,',...'end,',...'end']);k2=uicontrol('parent',h0,...'style','pushbutton',...'position',[180 100 50 30],...'string','停止',...'callback',[...'k=0;,',...'set(e1,''string'',m),',...'p=get(h,''xdata'');,',...'q=get(h,''ydata'');,',...'set(e2,''string'',p);,',...'set(e3,''string'',q)']);k3=uicontrol('parent',h0,...'style','pushbutton',...'position',[280 100 50 30],... 'string','关闭',...'callback','close');e1=uicontrol('parent',h0,...'style','edit',...'position',[60 30 60 20]);t1=uicontrol('parent',h0,...'style','text',...'string','循环次数',...'position',[60 50 60 20]);e2=uicontrol('parent',h0,...'style','edit',...'position',[180 30 50 20]);t2=uicontrol('parent',h0,...'style','text',...'string','终点的X坐标值',...'position',[155 50 100 20]);e3=uicontrol('parent',h0,...'style','edit',...'position',[300 30 50 20]);t3=uicontrol('parent',h0,...'style','text',...'string','终点的Y坐标值',...'position',[275 50 100 20]);实例33：曲线转换按钮h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例33');x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid on'if i==1,',...'i=0;,',...'y=cos(x);,',...'delete(h),',...'set(hm,''string'',''正弦函数''),',...'h=plot(x,y);,',...'grid on,',...'else if i==0,',...'i=1;,',...'y=sin(x);,',...'set(hm,''string'',''余弦函数''),',...'delete(h),',...'h=plot(x,y);,',...'grid on,',...'end,',...'end'];hm=uicontrol(gcf,'style','pushbutton',... 'string','余弦函数',...'callback',huidiao);i=1;set(hm,'position',[250 20 60 20]);set(gca,'position',[0.2 0.2 0.6 0.6]) title('按钮的使用')hold on实例34：栅格控制按钮h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例34');x=0:0.5:2*pi;y=sin(x);plot(x,y)huidiao1=[...'set(h_toggle2,''value'',0),',...'grid on,',...];'set(h_toggle1,''value'',0),',...'grid off,',...];h_toggle1=uicontrol(gcf,'style','togglebutton',... 'string','grid on',...'value',0,...'position',[20 45 50 20],...'callback',huidiao1);h_toggle2=uicontrol(gcf,'style','togglebutton',... 'string','grid off',...'value',0,...'position',[20 20 50 20],...'callback',huidiao2);set(gca,'position',[0.2 0.2 0.6 0.6])title('开关按钮的使用')实例35：编辑框的使用h0=figure('toolbar','none',...'position',[200 150 350 250],...'name','实例35');f='Please input the letter';huidiao1=[...'g=upper(f);,',...'set(h2_edit,''string'',g),',...];huidiao2=[...'g=lower(f);,',...'set(h2_edit,''string'',g),',...];h1_edit=uicontrol(gcf,'style','edit',...'position',[100 200 100 50],...'HorizontalAlignment','left',...'string','Please input the letter',...'callback','f=get(h1_edit,''string'');',...'background','w',...'max',5,...'min',1);h2_edit=uicontrol(gcf,'style','edit',...'HorizontalAlignment','left',...'position',[100 100 100 50],...'background','w',...'max',5,...'min',1);h1_button=uicontrol(gcf,'style','pushbutton',... 'string','小写变大写',...'position',[100 45 100 20],...'callback',huidiao1);h2_button=uicontrol(gcf,'style','pushbutton',... 'string','大写变小写',...'position',[100 20 100 20],...'callback',huidiao2);实例36：弹出式菜单h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例36');x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid onhm=uicontrol(gcf,'style','popupmenu',...'string',...'sin(x)|cos(x)|sin(x)+cos(x)|exp(-sin(x))',... 'position',[250 20 50 20]);set(hm,'value',1)huidiao=[...'v=get(hm,''value'');,',...'switch v,',...'case 1,',...'delete(h),',...'y=sin(x);,',...'h=plot(x,y);,',...'grid on,',...'case 2,',...'delete(h),',...'y=cos(x);,',...'h=plot(x,y);,',...'grid on,',...'case 3,',...'delete(h),',...'y=sin(x)+cos(x);,',...'h=plot(x,y);,',...'grid on,',...'case 4,',...'delete(h),',...'y=exp(-sin(x));,',...'h=plot(x,y);,',...'grid on,',...'end'];set(hm,'callback',huidiao)set(gca,'position',[0.2 0.2 0.6 0.6]) title('弹出式菜单的使用')实例37：滑标的使用h0=figure('toolbar','none',...'position',[200 150 450 250],... 'name','实例37');[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;h0=mesh(x,y,z);h1=axes('position',...[0.2 0.2 0.5 0.5],...'visible','off');htext=uicontrol(gcf,...'units','points',...'position',[20 30 45 15],...'string','brightness',...'style','text');hslider=uicontrol(gcf,...'units','points',...'position',[10 10 300 15],...'min',-1,...'max',1,...'style','slider',...'callback',...'brighten(get(hslider,''value''))');实例38：多选菜单h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例38');[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;h0=mesh(x,y,z);hlist=uicontrol(gcf,'style','listbox',...'string','default|spring|summer|autumn|winter',... 'max',5,...'min',1,...'position',[20 20 80 100],...'callback',[...'k=get(hlist,''value'');,',...'switch k,',...'case 1,',...'colormap default,',...'case 2,',...'colormap spring,',...'case 3,',...'colormap summer,',...'case 4,',...'colormap autumn,',...'case 5,',...'colormap winter,',...'end']);实例39：菜单控制的使用h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例39');x=0:0.5:2*pi;y=cos(x);h=plot(x,y);grid onset(gcf,'toolbar','none')hm=uimenu('label','example');huidiao1=[...'set(hm_gridon,''checked'',''on''),',...'set(hm_gridoff,''checked'',''off''),',...'grid on'];huidiao2=[...'set(hm_gridoff,''checked'',''on''),',...'set(hm_gridon,''checked'',''off''),',...'grid off'];hm_gridon=uimenu(hm,'label','grid on',... 'checked','on',...'callback',huidiao1);hm_gridoff=uimenu(hm,'label','grid off',... 'checked','off',...'callback',huidiao2);实例40：UIMENU菜单的应用h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例40');h1=uimenu(gcf,'label','函数');h11=uimenu(h1,'label','轮廓图',...'callback',[...'set(h31,''checked'',''on''),',...'set(h32,''checked'',''off''),',...'[x,y,z]=peaks;,',...'contour3(x,y,z,30)']);h12=uimenu(h1,'label','高斯分布',...。
matlab编程常用指令
matlab编程常用指令
Matlab是一种功能强大的数值计算软件，许多科学家和工程师都喜欢使用它进行数据分析和可视化。
以下是一些Matlab编程常用指令：
1. clear - 清除变量
2. clc - 清除命令行窗口
3. close all - 关闭所有图形窗口
4. format - 控制命令窗口的显示格式
5. help - 查看函数的帮助文档
6. load - 加载数据文件
7. save - 保存变量到文件中
8. plot - 绘制二维图像
9. semilogx - 绘制x轴为对数坐标的二维图像
10. semilogy - 绘制y轴为对数坐标的二维图像
11. loglog - 绘制双对数坐标的二维图像
12. histogram - 绘制直方图
13. scatter - 绘制散点图
14. subplot - 在一个图像窗口中创建多个子图
15. xlabel - 坐标轴标签
16. ylabel - 坐标轴标签
17. title - 图像标题
18. legend - 图例
19. find - 找到满足条件的数组元素的索引
20. max - 求得矩阵的最大值
21. min - 求得矩阵的最小值
22. mean - 求得矩阵的平均值
23. median - 求得矩阵的中位数
24. std - 求得矩阵的标准差
25. sort - 对矩阵进行排序
这些是Matlab编程中常用的指令，当然还有更多的指令可以用来操作数据、绘图、数学计算等。
有了这些指令，我们可以更方便、更快捷地进行数据分析和可视化。
matlab实用程序百例（6）
matlab实用程序百例（6）可设置函数曲线光源的用户界面h0=figure('toolbar','none',...'position',[198****8468],...'name','实例53');h1=axes('parent',h0,...'position',[0.15 0.5 0.7 0.5],...'visible','off');[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;fh=surf(x,y,z);shading interpview([-60 30])camlight leftlightk=light('position',[0 -2 1]);button1=uicontrol('parent',h0,...'style','pushbutton',...'string','设置光线',...'position',[80 60 70 30],...'callback',[...'an1=inputdlg(''光线来源的X轴坐标'');,',... 'k1=str2num(an1{1});,',...'an2=inputdlg(''光线来源的Y轴坐标'');,',... 'k2=str2num(an2{1});,',...'an3=inputdlg(''光线来源的Z轴坐标'');,',... 'k3=str2num(an3{1});,',...'set(lightk,''position'',[k1 k2 k3]);,',...'set(edit1,''string'',[''['',num2str(k1),'''',num2str(k2),'' '',num2str(k3),'']'']);']);button2=uicontrol('parent',h0,...'style','pushbutton',...'string','关闭',...'position',[250 60 70 30],...'callback','close');edit1=uicontrol('parent',h0,...'style','edit',...'max',2,...'min',0,...'fontsize',15,...'backgroundcolor',[1 1 1],...'string','[0 -2 1]',...'position',[80 110 220 30]);text1=uicontrol('parent',h0,...'style','text',...'backgroundcolor',[0.75 0.75 0.75],...'fontsize',15,...'string','光线来源坐标',...'position',[80 140 220 30]);添加效果h0=figure('toolbar','none',...'position',[200 50 300 350],...'name','实例54');h1=axes('parent',h0,...'position',[0.2 0.4 0.6 0.6],...'visible','off');ezsurf('sin(sqrt(x.^2+y.^2))/sqrt(x.^2+y.^2)',[-6*pi,6*pi]) b1=uicontrol('parent',h0,...'units','points',...'tag','b1',...'style','pushbutton',...'string','设置',...'backgroundcolor',[0.75 0.75 0.75],... 'position',[40 50 50 20],...'callback',[...'view(0,75);,',...'shading interp;,',...'lightangle(-45,30);,',...'k=findobj(gca,''type'',''surface'');,'... 'set(k,''facelighting'',''phong'');,',... 'set(k,''ambientstrength'',0.3);,',...'set(k,''diffusestrength'',0.8);,',...'set(k,''specularstrength'',0.9);,',...'set(k,''specularexponent'',25);,',... 'set(k,''backfacelighting'',''unlit'')']); b2=uicontrol('parent',h0,...'units','points',...'tag','b2',...'style','pushbutton',...'string','关闭',...'backgroundcolor',[0.75 0.75 0.75],... 'position',[120 50 50 20],...'callback','close');查询日期h0=figure('toolbar','none',...'position',[198****8468],...'name','实例55');h1=axes('parent',h0,...'position',[0.15 0.5 0.7 0.5],...'visible','off');huidiao=[...'yearnum=str2num(get(edit1,''string''));,',...'monthnum=str2num(get(edit2,''string''));,',...'daynum=str2num(get(edit3,''string''));,',...'monthday=[0 31 28 31 30 31 30 31 31 30 31 30 31];,',... 'dyear=yearnum-2000;,',...'beishu=fix(dyear/4);,',...'yushu=rem(yearnum,4);,',...'if yushu==0,',...'monthday(3)=29;,',...'end,',...'mday=0;,',...'for i=1:monthnum,',...'mday=monthday(i)+mday;,',...'end,',...'yearday=mday+daynum-1;,',...'noweek=fix(yearday/7);,',...'set(edit5,''string'',[''第'',num2str(noweek),''周'']);,',...'if dyear>0,',...'if yushu==0,',...'beishu=beishu-1;,',...'end,',...'dday=yearday+365*dyear+beishu+1;,',...'end,',...'if dyear<=0,',...'dday=365*dyear+yearday+beishu;,',...'end,',...'mweek=rem(dday,7)+7;,',...'set(edit4,''string'',''Sunday'');,',...'end,',...'if mweek==9,',...'set(edit4,''string'',''Monday'');,',...'end,',...'if mweek==10,',...'set(edit4,''string'',''Tuesday'');,',...'end,',...'if mweek==11,',...'set(edit4,''string'',''Wednesday'');,',... 'end,',...'if mweek==12,',...'set(edit4,''string'',''Thursday'');,',... 'end,',...'if mweek==13,',...'set(edit4,''string'',''Friday'');,',...'end,',...'if mweek==7,',...'set(edit4,''string'',''Saturday'');,',... 'end,',...'if mweek==6,',...'set(edit4,''string'',''Friday'');,',...'end,',...'if mweek==5,',...'set(edit4,''string'',''Thursday'');,',... 'end,',...'if mweek==4,',...'set(edit4,''string'',''Wednesday'');,',... 'end,',...'set(edit4,''string'',''Tuesday'');,',... 'end,',...'if mweek==2,',...'set(edit4,''string'',''Monday'');,',... 'end,',...'if mweek==1,',...'set(edit4,''string'',''Sunday'');,',... 'end'];edit1=uicontrol('parent',h0,...'style','edit',...'horizontalalignment','right',...'position',[40 300 50 20]);text1=uicontrol('parent',h0,...'style','text',...'string','年',...'horizontalalignment','left',...'position',[90 300 50 20]);edit2=uicontrol('parent',h0,...'style','edit',...'horizontalalignment','right',...'position',[160 300 50 20]);text2=uicontrol('parent',h0,...'style','text',...'string','月',...'horizontalalignment','left',...'position',[210 300 50 20]);edit3=uicontrol('parent',h0,...'style','edit',...'horizontalalignment','right',...'position',[280 300 50 20]);text3=uicontrol('parent',h0,... 'style','text',...'string','日',...'horizontalalignment','left',...'position',[330 300 50 20]);edit4=uicontrol('parent',h0,... 'style','edit',...'horizontalalignment','left',...'position',[210 200 120 20]); text4=uicontrol('parent',h0,... 'style','text',...'string','查找的日期为',...'horizontalalignment','right',... 'position',[110 200 100 20]); edit5=uicontrol('parent',h0,... 'style','edit',...'horizontalalignment','left',...'position',[210 100 120 20]); text1=uicontrol('parent',h0,... 'style','text',...'string','该日处于',...'horizontalalignment','left',...'position',[160 100 50 20]); button1=uicontrol('parent',h0,... 'style','pushbutton',...'position',[80 40 80 30],...'string','开始',...'callback',huidiao);button2=uicontrol('parent',h0,...'style','pushbutton',...'position',[220 40 80 30],...'string','关闭',...'callback','close');图形效果（1）h0=figure('toolbar','none',...'position',[198****0468],...'name','实例56');h1=axes('parent',h0,...'position',[0.3 0.45 0.5 0.5],...'visible','off');l1=uimenu(gcf,'label','Draw figure',... 'tag','l1');huidiao=[...'if get(r1,''value'')==1,',...'shading faceted,',...'end,',...'if get(r2,''value'')==1,',...'shading flat,',...'end,',...'if get(r3,''value'')==1,',...'shading interp,',...'end,',...'k=get(p1,''value'');,',...'switch k,',...'case 1,',...'colormap(''cool''),',...'case 2,',...'colormap(''spring''),',...'case 3,',...'colormap(''summer''),',...'case 4,',...'colormap(''autumn''),',...'case 5,',...'colormap(''winter''),',...'end'];l11=uimenu('parent',l1,...'label','Surface',...'tag','l11',...'callback',[...'[x,y]=meshgrid(-8:0.5:8);,',... 'r=sqrt(x.^2+y.^2)+eps;,',... 'z=sin(r)./r;,',...'surf(x,y,z),',... huidiao]);l12=uimenu('parent',l1,...'label','Mesh',...'tag','l12',...'callback',[...'mesh(peaks),',... huidiao]);l13=uimenu('parent',l1,...'label','Membrane',...'tag','l13',...'callback',[...'mesh(membrane),',... huidiao]);f1=uicontrol('parent',h0,... 'units','points',...'listboxtop',0,...'position',[12 6 100 101],...'style','frame',...'tag','f1');r1=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 58.5 72.75 16.5],...'string','shading faceted',...'style','radiobutton',...'tag','r1',...'value',1,...'callback',[...'shading faceted,',...'set(r1,''value'',1);,',...'set(r2,''value'',0);,',...'set(r3,''value'',0);']);r2=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 35.25 78.75 18.75],...'string','shading flat',...'style','radiobutton',...'tag','r2',...'value',0,...'callback',[...'shading flat,',...'set(r2,''value'',1);,',...'set(r1,''value'',0);,',...'set(r3,''value'',0);']);r3=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 12.75 71.25 18.75],...'string','shading interp',...'style','radiobutton',...'tag','r3',...'value',0,...'callback',[...'shading interp,',...'set(r3,''value'',1);,',...'set(r1,''value'',0);,',...'set(r2,''value'',0);']);t1=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'fontsize',12,...'listboxtop',0,...'position',[14.25 75.75 90.75 22.5],...'string','平滑处理',...'style','text',...'tag','t1');t2=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'fontsize',12,...'listboxtop',0,...'position',[117 69 72.75 17.5],...'string','设置色调',...'style','text',...'tag','t2');p1=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],...'listboxtop',0,...'position',[116.25 39 72.75 20.25],...'string','Cool|Spring|Summer|Autumn|Winter',... 'style','popupmenu',...'tag','p1',...'value',1,...'callback',[...'k=get(p1,''value'');,',...'switch k,',...'case 1,',...'colormap(''cool''),',...'case 2,',...'colormap(''spring''),',...'case 3,',...'colormap(''summer''),',...'case 4,',...'colormap(''autumn''),',...'case 5,',...'colormap(''winter''),',...'end']);b1=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],...'listboxtop',0,...'position',[12 243 72.75 30.75],...'string','关闭',...'tag','b1',...'callback','close');b2=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[216.75 67.5 83.25 18.75],... 'string','Colorbar',...'tag','b2',...'callback','colorbar');图形效果h0=figure('toolbar','none',...'position',[168 94.5 315 289.5],...'name','实例57');h1=axes('parent',h0,...'position',[0.4 0.4 0.5 0.5],...'visible','off');f1=uicontrol('parent',h0,...'style','frame',...'position',[15 10 80 70],...'string','dull',...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'tag','r1',...'value',1,...'callback',[...'set(r1,''value'',1);,',...'set(r2,''value'',0);,',...'set(r3,''value'',0);,',...'material dull']);r1=uicontrol('parent',h0,...'style','radiobutton',...'position',[19.5 58.5 72.75 16.5],...'string','dull',...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'tag','r1',...'value',1,...'callback',[...'set(r1,''value'',1);,',...'set(r2,''value'',0);,',...'set(r3,''value'',0);,',...'material dull']);r2=uicontrol('parent',h0,...'style','radiobutton',...'position',[19.5 35.25 72.75 16.5],...'string','metal',...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'tag','r2',...'value',0,...'callback',[...'set(r2,''value'',1);,',...'set(r1,''value'',0);,',...'set(r3,''value'',0);,',...'material metal']);r3=uicontrol('parent',h0,...'style','radiobutton',...'position',[19.5 12.75 72.75 16.5],...'string','shiny',...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'tag','r3',...'value',0,...'callback',[...'set(r3,''value'',1);,',...'set(r1,''value'',0);,',...'set(r2,''value'',0);,',...'material shiny']);u1=uimenu('parent',h0,...'label','绘图',...'backgroundcolor',[0.753 0.753 0.753],... 'tag','u1',...'callback',[...'[x,y]=meshgrid(-8:0.5:8);,',...'r=sqrt(x.^2+y.^2)+eps;,',...'z=sin(r)./r;,',...'surf(x,y,z),',...'shading interp']);b1=uicontrol('parent',h0,...'style','pushbutton',...'position',[19.5 150 60 20],...'string','light',...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'tag','b1',...'callback','camlight headlight');b2=uicontrol('parent',h0,...'style','pushbutton',...'position',[19.5 100 60 20],...'string','关闭',...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'tag','b2',...'callback','close');可控制小球运动速度的用户界面h0=figure('toolbar','none',...'position',[198****0468],...'name','实例58');h1=axes('parent',h0,...'position',[0.25 0.45 0.5 0.5],...'visible','off');t=0:0.1:4*pi;x=sin(t);y=cos(t);plot(x,y)axis equalaxis offh=line('color',[1 0 0],...'linestyle','.',...'xdata',0,...'ydata',1,...'markersize',20,...'erasemode','xor');n=length(t);i=1;speed=0.01;k=0;b1huidiao=[...'k=0;,',...'while 1,',...'set(h,''xdata'',x(i),''ydata'',y(i));,',... 'drawnow,',...'pause(speed),',...'i=i+1;,',...'if i>n,',...'i=1;,',...'end,',...'if k==1,',...'break,',...'end,',...'end'];b1=uicontrol('parent',h0,...'units','points',...'tag','b1',...'style','pushbutton',...'string','开始',...'backgroundcolor',[0.75 0.75 0.75],... 'position',[30 80 50 20],...'callback',b1huidiao);b2=uicontrol('parent',h0,...'units','points',...'style','pushbutton',...'string','停止',...'backgroundcolor',[0.75 0.75 0.75],... 'position',[100 80 50 20],...'callback','k=1;');b3=uicontrol('parent',h0,...'units','points',...'tag','b3',...'style','pushbutton',...'string','关闭',...'backgroundcolor',[0.75 0.75 0.75],... 'position',[170 80 50 20],...'callback',[...'k=1;,',...'close']);s1=uicontrol('parent',h0,...'units','points',...'tag','s1',...'style','slider',...'value',50*speed,...'max',1,...'min',0,...'backgroundcolor',[0.75 0.75 0.75],... 'position',[30 115 190 20],...'callback',[...'m=get(s1,''value'');,',...'speed=m/50;']);t1=uicontrol('parent',h0,...'units','points',...'style','text',...'fontsize',15,...'string','小球运动速度',...'backgroundcolor',[0.75 0.75 0.75],...'position',[30 135 190 20]);设置坐标轴纵横轴比h0=figure('name','实例59');h1=axes('parent',h0,...'position',[0.3 0.45 0.5 0.5],...'visible','off');u1=uimenu('parent',h0,...'label','绘图',...'backgroundcolor',[0.753 0.753 0.753],... 'tag','u1',...'callback',[...'[x,y]=meshgrid(-8:0.5:8);,',...'r=sqrt(x.^2+y.^2)+eps;,',...'z=sin(r)./r;,',...'mesh(x,y,z),',...'shading interp,',...'axis normal']);f1=uicontrol('parent',h0,...'units','points',...'listboxtop',0,...'position',[12 6 100 150],...'style','frame',...'tag','f1');t1=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 130 72.75 16.5],...'string','坐标纵横比',...'style','text',...'tag','t1');r1=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 110 72.75 16.5],...'string','axis equal',...'style','radiobutton',...'tag','r1',...'value',1,...'callback',[...'set(r1,''value'',1);,',...'set(r2,''value'',0);,',...'set(r3,''value'',0);,',...'set(r4,''value'',0);,',...'set(r5,''value'',0);,',...'axis equal']);r2=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 85 72.75 16.5],...'string','axis square',...'style','radiobutton',...'tag','r2',...'value',0,...'callback',[...'set(r2,''value'',1);,',...'set(r1,''value'',0);,',...'set(r3,''value'',0);,',...'set(r4,''value'',0);,',...'set(r5,''value'',0);,',...'axis square']);r3=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 60 72.75 16.5],...'string','axis image',...'style','radiobutton',...'tag','r3',...'value',0,...'callback',[...'set(r3,''value'',1);,',...'set(r2,''value'',0);,',...'set(r1,''value'',0);,',...'set(r4,''value'',0);,',...'set(r5,''value'',0);,',...'axis image']);r4=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 35 72.75 16.5],...'string','axie vis3d',...'style','radiobutton',...'tag','r4',...'value',0,...'callback',[...'set(r4,''value'',1);,',...'set(r2,''value'',0);,',...'set(r3,''value'',0);,',...'set(r1,''value'',0);,',...'set(r5,''value'',0);,',...'axis vis3d']);r5=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[19.5 10 72.75 16.5],...'string','axis auto',...'style','radiobutton',...'tag','r5',...'value',0,...'callback',[...'set(r5,''value'',1);,',...'set(r2,''value'',0);,',...'set(r3,''value'',0);,',...'set(r4,''value'',0);,',...'set(r1,''value'',0);,',...'axis auto']);b1=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[12 243 72.75 30.75],...'string','关闭',...'tag','b1',...'callback','close');b2=uicontrol('parent',h0,...'units','points',...'backgroundcolor',[0.753 0.753 0.753],... 'listboxtop',0,...'position',[216.75 67.5 83.25 18.75],... 'string','Colorbar',...'tag','b2',...'callback','colorbar');动态文本显示h0=figure('toolbar','none',...'position',[198****0468],...'name','实例60');h1=axes('parent',h0,...'position',[0.25 0.45 0.5 0.5],...'visible','off');str1='当前阻尼比＝';z=0.52;t=0:0.1:10;y=step(1,[1 2*z 1],t);hline=plot(t,y);grid onr1=uicontrol('parent',h0,...'units','points',...'tag','r1',...'style','radio',...'string','grid on',...'position',[30 120 60 20],...'backgroundcolor',[0.75 0.75 0.75],... 'value',1,...'callback',[...'grid on,',...'set(r1,''value'',1);,',...'set(r2,''value'',0)']);r2=uicontrol('parent',h0,...'units','points',...'tag','r2',...'style','radio',...'string','grid on',...'position',[30 95 60 20],...'backgroundcolor',[0.75 0.75 0.75],... 'value',0,...'callback',[...'grid off,',...'set(r2,''value'',1);,',...'set(r1,''value'',0)']);s1=uicontrol('parent',h0,...'units','points',...'tag','s1',...'style','slider',...'value',z,...'position',[100 95 150 20],...'backgroundcolor',[0.75 0.75 0.75],... 'max',1,...'min',0,...'callback',[...'z=get(s1,''value'');,',...。
matlab程序大全
频率特性类题目1 一个系统的开环传递函数为 ，试绘制其当K=5、30时系统的开环频率特性Nyquist 图，并判断系统的稳定性。
w=linspace(0.5,5,1000)*pi; sys1=zpk([],[0 -10 -2],100) sys2=zpk([],[0 -10 -2],600) figure(1)nyquist(sys1,w)title('system nyquist charts with k=5') figure(2)nyquist(sys2,w)title('system nyquist charts with k=30')由图可知K=5时，开环Nyquist 曲线没有包围（-1，j0）点，所以系统稳定。
K=30时，开环Nyquist 曲线包围（-1，j0）点，所以系统不稳定。
2系统开环传递函数为 ，建立其零极点增益模型，然后分别绘制当K=5、K=30时系统的开环频率特性Bode 图，并判断系统的稳定性。
sys1=zpk([],[0 -10 -2],100) sys2=zpk([],[0 -10 -2],600) figure(1)[Gm1,Pm1,Wcg1,Wcp1]=margin(sys1) bode(sys1)title('system bode charts with k=5'),grid figure(2)[Gm2,Pm2,Wcg2,Wcp2]=margin(sys2) bode(sys2)title('system bode charts with k=30'),grid 因为K=5时，Wcg>Wcp,所以系统稳定。
K=10时，Wcg<Wcp,系统不稳定。
3 系统开环传递函数为 ，计算K=5和K=30时系统的幅值裕度与相位裕度。
sys1=zpk([],[0 -10 -2],100) sys2=zpk([],[0 -10 -2],600)[Gm1,Pm1,Wcg1,Wcp1]=margin(sys1) [Gm2,Pm2,Wcg2,Wcp2]=margin(sys2)4 已知某系统的闭环传递函数()s Φ如下，试用roots （）命令来判断系统的稳定性。
(word完整版)Matlab常用命令集合
Matlab常用命令集合1.常用命令—-〉管理命令和函数2.常用命令--〉管理变量和工作区(输入输出、内存管理等)3.常用命令—->管理命令控制窗口(command窗口）4.常用命令——〉文件和工作环境5.常用命令-—>启动和退出Matlab6.语言结构和调试命令--〉程序设计7.语言结构和调试命令--〉流程控制8.语言结构和调试命令—->交互输入9.语言结构和调试命令—-〉面向对象编程10.语言结构和调试命令-—〉程序调试11.语言结构和调试命令—->lasterr,lastwarn12.创建图形用户界面——〉对话框13.矩阵函数－数值线性代数-—〉矩阵分析14.矩阵函数－数值线性代数——>线性方程15.矩阵函数－数值线性代数-—〉特征值和奇异值16.矩阵函数－数值线性代数-—〉矩阵函数17.矩阵函数－数值线性代数--〉低级函数18.创建图形用户界面—-〉用户界面对象19.创建图形用户界面--〉其他函数20.创建图形用户界面-—〉其他函数21.图象可视化函数--〉三维绘图函数22.图象可视化函数—->绘制标注和网络23.图象可视化函数-—>体数据可视化24.图象可视化函数-—>表面、网格和轮廓绘制25.图象可视化函数-—〉域生成26.图象可视化函数--〉专门图形绘制27.图象可视化函数——〉视觉控制28.图象可视化函数-->颜色操作29.图象可视化函数—->打印函数30.图象可视化函数-->图形图象处理31.双重函数和非线性数值方法32.多项式和插值函数--〉多项式33.多项式和插值函数—-〉数据插值34.多项式命令补遗35.数组函数--〉单元数组函数36.数组函数——>多维数组函数37.位操作函数、结构和对象函数38.底层文件输入输出函数39.底层文件输入输出函数—-〉特殊函数40.字符串函数--〉常用函数41.字符串函数—->字符串操作42.字符串函数—->字符串和数值的转换43.字符串函数-->基转换44.声音处理函数-->常用声音函数45.声音处理函数-—>特殊声音函数46.声音处理函数-—〉WAV声音文件47.稀疏矩阵函数—->基本稀疏矩阵48.稀疏矩阵函数—-〉满阵和稀疏矩阵的转换49.稀疏矩阵函数—-〉稀疏矩阵的非零元素操作50.稀疏矩阵函数—->稀疏矩阵的可视化51.稀疏矩阵函数——〉排序算法52.稀疏矩阵函数—->范数、条件数和秩53.稀疏矩阵函数-—>线性方程的稀疏系统54.稀疏矩阵函数——〉线性方程的稀疏系统55.稀疏矩阵函数-->杂项函数56.数据分析和傅立叶变换—-〉基本运算57.数据分析和傅立叶变换--〉有限差分58.数据分析和傅立叶变换-—>相关59.数据分析和傅立叶变换——>滤波和卷积60.数据分析和傅立叶变换—-〉傅立叶变换61.数据分析和傅立叶变换-->向量函数62.数学函数和坐标变换—->基本数学函数63.数学函数和坐标变换--〉特殊函数64.数学函数和坐标变换-—〉坐标变换65.矩阵和矩阵操作基础--〉矩阵和数组基础66.矩阵和矩阵操作基础——〉特殊变量和常量67.矩阵和矩阵操作基础-—〉时间和日期68.矩阵和矩阵操作基础—->矩阵操作69.矩阵和矩阵操作基础-->特殊函数矩阵70.运算符和逻辑函数—->算数运算符71.运算符和逻辑函数—->关系运算符72.运算符和逻辑函数—-〉特殊函数73.运算符和逻辑函数——〉逻辑函数74.simulink的命令集75.rtw相关的命令1．常用命令——〉管理命令和函数addpath 添加目录到MATLAB搜索路径doc 在Web浏览器上现实HTML文档help 显示Matlab命令和M文件的在线帮助helpwinhelpdesk helplookfor 在基于Matlab搜索路径的所有M文件中搜索关键字partialpath 部分路径名path 所有关于路径名的处理pathtool 一个不错的窗口路径处理界面rmpath 删除搜索路径中指定目录type 显示指定文件的内容ver 版本信息version 版本号web 打开web页what 列出当前目录吓所有的M文件Mat文件和Mex文件whatsnew 显示readme文件which 显示文件位置（返回）2．常用命令-->管理变量和工作区（输入输出、内存管理等)clear 从内存中删除disp 显示文本或数组内容length 数组长度(最长维数）load 重新载入变量（从磁盘上)mlock 锁定文件,防止文件被错误删除munlock 解锁文件openvar 在数组编辑器中打开变量pack 整理内存空间save 保存变量到文件8＊)size 数组维数who 列出内存变量whos 列出内存变量,同时显示变量维数workspace 显示工作空间窗口（返回)3．常用命令—-〉管理命令控制窗口(command窗口）clc 清空命令窗口echo 禁止或允许显示执行过程format 设置输出显示格式home 光标移动到命令窗口左上角more 设置命令窗口页输出格式（返回)4．常用命令--〉文件和工作环境cd 改变工作目录copyfile 复制文件delete 删除文件和图形对象diary 把命令窗口的人机交互保存到文件dir 显示目录edit 编辑文本文件fileparts 返回文件的各个部分fullfile 使用指定部分建立文件全名inmem 返回内存(伪代码区)的matlab函数名ls 在unix系统中列出目录(win中亦可）matlabroot 根目录mkdir 新建目录open 打开文件pwd 显示当前目录tempdir 返回系统临时目录的名字tempname 随机给出一个临时字符串（可用作文件名）! 直接调用操作系统command命令(返回）5．常用命令-—〉启动和退出Matlabmatlabrc Matlab的启动M文件exit 退出Matlabquit 退出Matlabstartup 运行matlab启动文件（返回）6．语言结构和调试命令-—〉程序设计builtin 从可重载方法中调用内置函数eval 执行包含可执行表达式的字符串evalc 计算并返回表达式的值evalin 执行某个工作空间中的包含表达式的字符串feval 执行函数（从函数名或函数句柄)function 函数头global 定义全局变量nargchk 检查输入参数数目persistent 定义常量script 作为脚本的M文件（返回）7．语言结构和调试命令——>流程控制break 停止执行循环case case语句，switch语句一部分switch switch_exprcase case_exprstatmentscaseotherwiseendcatch try—catch语句一部分，捕捉程序else if条件语句一部分elseif if条件语句一部分end 终止for、while、switch、try和if语句error 显示错误信息for 循环语句一（确定次数）if if条件语句other switch语句一部分return 返回到调用函数switch 开关语句try try程序块warning 类似于disp，但可被禁止while 循环语句(次数不确定)（返回）8．语言结构和调试命令-->交互输入input 交互输入keyboard 在M文件中遇到keyboard时将在命令窗口产生交互直到输入return命令menu 为用户输入产生一个选择菜单pause 暂停（返回）9．语言结构和调试命令—->面向对象编程class 创建一个对象或者返回一个对象类double 转换为双精度inferiorto 亚类关系inline 创建一个内联函数int8，int16,int32 转换到符号整数isa 检查是否为所给类的对象loadobj load函数用户定义扩展saveobj save函数用户定义扩展single 转换为单精度superiorto 超类关系uint8,uint16，uint32 转换到无符号整数10语言结构和调试命令-—〉程序调试dbclear 断点清除dbcont 重新开始执行dbdown 改变当前工作空间dbmex 调试Mex文件dbquit 退出调试模式dbstack 显示函数调用堆栈dbstatus 列出所有断点dbstep 从断点处执行dbstop 设置断点dbtype 列出带行号的M文件内容dbup 改变当前工作空间（返回）11．语言结构和调试命令-—〉lasterr,lastwarnlasterr 返回Matlab中产生的最后一个异常信息。
matlab代码大全
matlab代码大全MATLAB主要命令汇总MATLAB函数参考附录1.1 管理用命令函数名功能描述函数名功能描述addpath 增加一条搜索路径 rmpath 删除一条搜索路径demo 运行Matlab演示程序 type 列出.M文件doc 装入超文本文档 version 显示Matlab的版本号help 启动联机帮助 what 列出当前目录下的有关文件lasterr 显示最后一条信息 whatsnew 显示Matlab的新特性lookfor 搜索关键词的帮助 which 造出函数与文件所在的目录path 设置或查询Matlab路径附录1.2管理变量与工作空间用命令函数名功能描述函数名功能描述clear 删除存中的变量与函数 pack 整理工作空间存disp 显示矩阵与文本 save 将工作空间中的变量存盘length 查询向量的维数 size 查询矩阵的维数load 从文件中装入数据 who,whos 列出工作空间中的变量名附录1.3文件与操作系统处理命令函数名功能描述函数名功能描述cd 改变当前工作目录 edit 编辑.M文件delete 删除文件 matlabroot 获得Matlab的安装根目录diary 将Matlab运行命令存盘 tempdir 获得系统的缓存目录dir 列出当前目录的容 tempname 获得一个缓存(temp)文件! 执行操作系统命令附录1.4窗口控制命令函数名功能描述函数名功能描述echo 显示文件中的Matlab中的命令 more 控制命令窗口的输出页面format 设置输出格式附录1.5启动与退出命令函数名功能描述函数名功能描述matlabrc 启动主程序 quit 退出Matlab环境startup Matlab自启动程序附录2 运算符号与特殊字符附录2.1运算符号与特殊字符函数名功能描述函数名功能描述+ 加 ... 续行标志- 减 , 分行符(该行结果不显示)* 矩阵乘 ; 分行符(该行结果显示).* 向量乘 % 注释标志^ 矩阵乘方 ! 操作系统命令提示符.^ 向量乘方矩阵转置kron 矩阵kron积 . 向量转置\ 矩阵左除 = 赋值运算/ 矩阵右除 == 关系运算之相等.\ 向量左除 ~= 关系运算之不等./ 向量右除 < 关系运算之小于: 向量生成或子阵提取 <= 关系运算之小于等于() 下标运算或参数定义 > 关系运算之大于[] 矩阵生成 >= 关系运算之大于等于{} & 逻辑运算之与. 结构字段获取符 | 逻辑运算之或. 点乘运算,常与其他运算符联合使用(如.\) ~ 逻辑运算之非xor 逻辑运算之异成附录2.2逻辑函数函数名功能描述函数名功能描述all 测试向量中所用元素是否为真is*(一类函数) 检测向量状态.其中*表示一个确定的函数(isinf)any 测试向量中是否有真元素*isa 检测对象是否为某一个类的对象exist 检验变量或文件是否定义 logical 将数字量转化为逻辑量find 查找非零元素的下标附录3 语言结构与调试附录3.1编程语言函数名功能描述函数名功能描述builtin 执行Matlab建的函数 global 定义全局变量eval 执行Matlab语句构成的字符串 nargchk 函数输入输出参数个数检验feval 执行字符串指定的文件 script Matlab语句及文件信息function Matlab函数定义关键词附录3.2控制流程函数名功能描述函数名功能描述break 中断循环执行的语句 if 条件转移语句case 与switch结合实现多路转移otherwise 多路转移中的缺省执行部分else 与if一起使用的转移语句 return 返回调用函数elseif 与if一起使用的转移语句 switch 与case结合实现多路转移end 结束控制语句块 warning 显示警告信息error 显示错误信息 while 循环语句for 循环语句附录3.3交互输入函数名功能描述函数名功能描述input 请求输入 menu 菜单生成keyboard 启动键盘管理 pause 暂停执行附录3.4面向对象编程函数名功能描述函数名功能描述class 生成对象 isa 判断对象是否属于某一类double 转换成双精度型 superiorto 建立类的层次关系inferiorto 建立类的层次关系unit8 转换成8字节的无符号整数inline 建立一个嵌对象附录3.5调试函数名功能描述函数名功能描述dbclear 清除调试断点 dbstatus 列出所有断点情况dbcont 调试继续执行 dbstep 单步执行dbdown 改变局部工作空间存 dbstop 设置调试断点dbmex 启动对Mex文件的调试sbtype 列出带命令行标号的.M 文件dbquit 退出调试模式 dbup 改变局部工作空间容dbstack 列出函数调用关系附录4 基本矩阵与矩阵处理附录4.1基本矩阵函数名功能描述函数名功能描述eye 产生单位阵 rand 产生随机分布矩阵linspace 构造线性分布的向量 randn 产生正态分布矩阵logspace 构造等对数分布的向量 zeros 产生零矩阵ones 产生元素全部为1的矩阵 : 产生向量附录4.2特殊向量与常量函数名功能描述函数名功能描述ans 缺省的计算结果变量 non 非数值常量常由0/0或Inf/Inf获得computer 运行Matlab的机器类型 nargin 函数中参数输入个数eps 精度容许误差(无穷小) nargout 函数中输出变量个数flops 浮点运算计数 pi 圆周率i 复数单元 realmax 最大浮点数值inf 无穷大 realmin 最小浮点数值inputname 输入参数名 varargin 函数中输入的可选参数j 复数单元 varargout 函数中输出的可选参数附录4.3时间与日期函数名功能描述函数名功能描述calender 日历 eomday 计算月末clock 时钟 etime 所用时间函数cputime 所用的CPU时间 now 当前日期与时间date 日期 tic 启动秒表计时器datenum 日期(数字串格式) toc 读取秒表计时器datestr 日期(字符串格式) weekday 星期函数datevoc 日期(年月日分立格式)附录4.4矩阵处理函数名功能描述函数名功能描述cat 向量连接 reshape 改变矩阵行列个数diag 建立对角矩阵或获取对角向量 rot90 将矩阵旋转90度fliplr 按左右方向翻转矩阵元素 tril 取矩阵的下三角部分flipud 按上下方向翻转矩阵元素 triu 取矩阵的上三角部分repmat 复制并排列矩阵函数附录5 特殊矩阵函数名功能描述函数名功能描述compan 生成伴随矩阵 invhilb 生成逆hilbert矩阵gallery 生成一些小的测试矩阵 magic 生成magic矩阵hadamard 生成hadamard矩阵 pascal 生成pascal矩阵hankel 生成hankel矩阵 toeplitz 生成toeplitz矩阵hilb 生成hilbert矩阵 wilkinson 生成wilkinson特征值测试矩阵附录6 数学函数附录6.1三角函数函数名功能描述函数名功能描述sin/asin 正弦/反正弦函数 sec/asec 正割/反正割函数sinh/asinh 双曲正弦/反双曲正弦函数 sech/asech 双曲正割/反双曲正割函数cos/acos 余弦/反余弦函数 csc/acsc 余割/反余割函数cosh/acosh 双曲余弦/反双曲余弦函数csch/acsch 双曲余割/反双曲余割函数tan/atan 正切/反正切函数 cot/acot 余切/反余切函数tanh/atanh 双曲正切/反双曲正切函数coth/acoth 双曲余切/反双曲余切函数atan2 四个象限反正切函数附录6.2指数函数函数名功能描述函数名功能描述exp 指数函数 log10 常用对数函数log 自然对数函数 sqrt 平方根函数附录6.3复数函数函数名功能描述函数名功能描述abs 绝对值函数 imag 求虚部函数angle 角相位函数 real 部函数conj 共轭复数函数附录6.4数值处理函数名功能描述函数名功能描述fix 沿零方向取整 round 舍入取整floor 沿-∞方向取整 rem 求除法的余数ceil 沿+∞方向取整 sign 符号函数附录6.5其他特殊数学函数函数名功能描述函数名功能描述airy airy函数 erfcx 比例互补误差函数besselh bessel函数(hankel函数) erfinv 逆误差函数bessili 改进的第一类bessel函数 expint 指数积分函数besselk 改进的第二类bessel函数 gamma gamma函数besselj 第一类bessel函数 gammainc 非完全gamma函数bessely 第二类bessel函数 gammaln gamma对数函数beta beta函数 gcd 最大公约数betainc 非完全的beta函数 lcm 最小公倍数betaln beta对数函数 log2 分割浮点数elipj Jacobi椭圆函数 legendre legendre伴随函数ellipke 完全椭圆积分 pow2 基2标量浮点数erf 误差函数 rat 有理逼近erfc 互补误差函数 rats 有理输出附录7 坐标转换函数名功能描述函数名功能描述cart2pol 笛卡儿坐标到极坐标转换 pol2cart 极坐标到笛卡儿坐标转换cart2sph 笛卡儿坐标到球面坐标转换 sph2cart 球面坐标到笛卡儿坐标转换附录8 矩阵函数附录8.1矩阵分析函数名功能描述函数名功能描述cond 求矩阵的条件数 rcond LINPACK倒数条件估计det 求矩阵的行列式 rref 矩阵的行阶梯型实现norm 求矩阵的数 rrefmovie 消元法解方程演示null 右零空间 subspace 子空间orth 正交空间 trace 矩阵的迹rank 求矩阵的秩附录8.2线性方程函数名功能描述函数名功能描述/,\ 线性方程求解 nnls 非零最小二乘chol Cholesky分解 pinv 求伪逆矩阵inv 矩阵求逆 qr 矩阵的QR分解lscov 最小二乘方差 qrdelete QR分解中删除一行lu 矩阵的LU三角分解 qrinsert QR分解中插入一行附录8.3特征值与奇异值函数名功能描述函数名功能描述banlance 改进特征值精度的均衡变换 qz QZ算法求矩阵特征值cdf2rdf 复块对角阵到实块对角阵转换 rdf2cdf 实块对角阵到复块对角阵转换eig 求矩阵的特征值和特征向量 schur Schur分解hess 求Hessenberg矩阵 svd 奇异值分解poly 求矩阵的特征多项式附录8.4矩阵函数函数名功能描述函数名功能描述expm 矩阵指数函数 logm 矩阵对数函数funm 矩阵任意函数 sqrtm 矩阵平方根附录9 数据分析与Fourier变换函数附录9.1基本运算函数名功能描述函数名功能描述cumprod 向量累积 prod 对向量中各元素求积cumsum 向量累加 sort 对向量中各元素排序max 求向量中最大元素 sortrows 对矩阵中各行排序min 求向量中最小元素 std 求向量中各元素标准差mean 求向量中各元素均值 sum 对向量中各元素求和median 求向量中中间元素 trapz 梯形法求数值积分附录9.2微分计算函数名功能描述函数名功能描述del2 离散Laplace变换 gradient 梯度计算diff 差分于近视微分附录9.3滤波与卷积函数名功能描述函数名功能描述Conv 卷给与多项式乘法 filter 一维数字滤波conv2 二维卷积 filter2 二维数字滤波Deconv 因式分解与多项式乘法附录9.4方差处理函数名功能描述函数名功能描述corrcoef 相关系数计算 cov 协方差计算附录9.5Fourier变换函数名功能描述函数名功能描述abs 绝对值函数 fftshift fft与fft2输出重排angle 相角函数 ifft 离散Fourier逆变换cplxpair 依共轭复数对重新排序 ifft2 二维离散Fourier逆变换fft离散Fourier变换 unwrap 相角矫正fft2 二维离散Fourier变换附录10 多项式处理函数附录10.1多项式处理函数名功能描述函数名功能描述conv 卷机与多项式乘法 polyfit 数据的多项式拟合deconv 因式分解与多项式乘法 polyval 多项式求值poly 求矩阵的特征多项式 polyvalm 多项式矩阵求值polyder 多项式求导 residue 部分分式展开polyeig 多项式特征值 roots 求多项式的根附录10.2数据插值函数名功能描述函数名功能描述griddata 数据网络的插值生成interpft 一维插值（FFT方法）interp1 一维插值（查表） interpn 多维插值（查表）。
matlab常用函数与常用指令大全[整理版]
![matlab常用函数与常用指令大全[整理版]](https://img.taocdn.com/s3/m/ff905ddadb38376baf1ffc4ffe4733687f21fc56.png)
matlab常用函数与常用指令大全matlab, 函数, 指令, 大全matlab常用函数与常用指令大全matlab常用函数- -1、特殊变量与常数ans 计算结果的变量名computer 确定运行的计算机eps 浮点相对精度Inf 无穷大I 虚数单位inputname 输入参数名NaN 非数nargin 输入参数个数nargout 输出参数的数目pi 圆周率nargoutchk 有效的输出参数数目realmax 最大正浮点数realmin 最小正浮点数varargin 实际输入的参量varargout 实际返回的参量操作符与特殊字符+ 加- 减* 矩阵乘法.* 数组乘（对应元素相乘）^ 矩阵幂 .^ 数组幂（各个元素求幂）\ 左除或反斜杠/ 右除或斜面杠./ 数组除（对应元素除）kron Kronecker张量积: 冒号() 圆括[] 方括 . 小数点.. 父目录 ... 继续, 逗号（分割多条命令）; 分号（禁止结果显示）% 注释! 感叹号' 转置或引用= 赋值== 相等<> 不等于& 逻辑与| 逻辑或~ 逻辑非xor 逻辑异或2、基本数学函数abs 绝对值和复数模长acos,acodh 反余弦，反双曲余弦acot,acoth 反余切，反双曲余切acsc,acsch 反余割，反双曲余割angle 相角asec,asech 反正割，反双曲正割secant 正切asin,asinh 反正弦，反双曲正弦atan,atanh 反正切，双曲正切tangent 正切atan2 四象限反正切ceil 向着无穷大舍入complex 建立一个复数conj 复数配对cos,cosh 余弦，双曲余弦csc,csch 余切，双曲余切cot,coth 余切，双曲余切exp 指数fix 朝0方向取整floor 朝负无穷取整*** 最大公因数imag 复数值的虚部lcm 最小公倍数log 自然对数log2 以2为底的对数log10 常用对数mod 有符号的求余nchoosek 二项式系数和全部组合数real 复数的实部rem 相除后求余round 取整为最近的整数sec,sech 正割，双曲正割sign 符号数sin,sinh 正弦，双曲正弦sqrt 平方根tan,tanh 正切，双曲正切3、基本矩阵和矩阵操作blkding 从输入参量建立块对角矩阵eye 单位矩阵linespace 产生线性间隔的向量logspace 产生对数间隔的向量numel 元素个数ones 产生全为1的数组rand 均匀颁随机数和数组randn 正态分布随机数和数组zeros 建立一个全0矩阵colon) 等间隔向量cat 连接数组diag 对角矩阵和矩阵对角线fliplr 从左自右翻转矩阵flipud 从上到下翻转矩阵repmat 复制一个数组reshape 改造矩阵roy90 矩阵翻转90度tril 矩阵的下三角triu 矩阵的上三角dot 向量点集cross 向量叉集ismember 检测一个集合的元素intersect 向量的交集setxor 向量异或集setdiff 向是的差集union 向量的并集数值分析和傅立叶变换cumprod 累积cumsum 累加cumtrapz 累计梯形法计算数值微分factor 质因子inpolygon 删除多边形区域内的点max 最大值mean 数组的均值mediam 中值min 最小值perms 所有可能的转换polyarea 多边形区域primes 生成质数列表prod 数组元素的乘积rectint 矩形交集区域sort 按升序排列矩阵元素sortrows 按升序排列行std 标准偏差sum 求和trapz 梯形数值积分var 方差del2 离散拉普拉斯diff 差值和微分估计gradient 数值梯度cov 协方差矩阵corrcoef 相关系数conv2 二维卷积conv 卷积和多项式乘法filter IIR或FIR滤波器deconv 反卷积和多项式除法filter2 二维数字滤波器cplxpair 将复数值分类为共轭对fft 一维的快速傅立叶变换fft2 二维快速傅立叶变换fftshift 将FFT的DC分量移到频谱中心ifft 一维快速反傅立叶变换ifft2 二维傅立叶反变换ifftn 多维快速傅立叶变换ifftshift 反FFT偏移nextpow2 最靠近的2的幂次unwrap 校正相位角多项式与插值conv 卷积和多项式乘法roots 多项式的根poly 具有设定根的多项式polyder 多项式微分polyeig 多项式的特征根polyfit 多项式拟合polyint 解析多项式积分polyval 多项式求值polyvalm 矩阵变量多项式求值residue 部分分式展开interp1 一维插值interp2 二维插值interp3 三维插值interpft 使用FFT的一维插值interpn 多维插值meshgrid 为3维点生成x和y的网格ndgrid 生成多维函数和插值的数组pchip 分段3次Hermite插值多项式ppval 分段多项式的值spline 3次样条数据插值绘图函数bar 竖直条图barh 水平条图hist 直方图histc 直方图计数hold 保持当前图形loglog x,y对数坐标图pie 饼状图plot 绘二维图polar 极坐标图semilogy y轴对数坐标图semilogx x轴对数坐标subplot 绘制子图bar3 数值3D竖条图bar3h 水平3D条形图comet3 3D慧星图cylinder 圆柱体fill3 填充的3D多边形plot3 3维空间绘图quiver3 3D震动（速度）图slice 体积薄片图sphere 球stem3 绘制离散表面数据wate***ll 绘制瀑布trisurf 三角表面clabel 增加轮廓标签到等高线图中datetick 数据格式标记grid 加网格线gtext 用鼠标将文本放在2D图中legend 图注plotyy 左右边都绘Y轴title 标题xlabel X轴标签ylabel Y轴标签zlabel Z轴标签contour 等高线图contourc 等高线计算contourf 填充的等高线图hidden 网格线消影meshc 连接网格/等高线mesh 具有参考轴的3D网格peaks 具有两个变量的采样函数surf 3D阴影表面图su***ce 建立表面低层对象surfc 海浪和等高线的结合surfl 具有光照的3D阴影表面trimesh 三角网格图1 常用指令(General Purpose Commands) 1.1 通用信息查询(General information) demo 演示程序help 在线帮助指令helpbrowser 超文本文档帮助信息helpdesk 超文本文档帮助信息helpwin 打开在线帮助窗info MA TLAB 和MathWorks 公司的信息subscribe MA TLAB 用户注册ver MA TLAB 和TOOLBOX 的版本信息version MA TLAB 版本whatsnew 显示版本新特征1.2 工作空间管理(Managing the workspace)clear 从内存中清除变量和函数exit 关闭MA TLABload 从磁盘中调入数据变量pack 合并工作内存中的碎块quit 退出MA TLABsave 把内存变量存入磁盘who 列出工作内存中的变量名whos 列出工作内存中的变量细节workspace 工作内存浏览器1.3 管理指令和函数(Managing commands and functions) edit 矩阵编辑器edit 打开M 文件inmem 查看内存中的P 码文件mex 创建MEX 文件open 打开文件pcode 生成P 码文件type 显示文件内容what 列出当前目录上的M、MA T、MEX 文件which 确定指定函数和文件的位置1.4 搜索路径的管理(Managing the seach patli) addpath 添加搜索路径rmpath 从搜索路径中删除目录path 控制MA TLAB 的搜索路径pathtool 修改搜索路径1.5 指令窗控制(Controlling the command window) beep 产生beep 声echo 显示命令文件指令的切换开关diary 储存MA TLAB 指令窗操作内容format 设置数据输出格式more 命令窗口分页输出的控制开关1.6 操作系统指令(Operating system commands) cd 改变当前工作目录computer 计算机类型copyfile 文件拷贝delete 删除文件dir 列出的文件dos 执行dos 指令并返还结果getenv 给出环境值ispc MA TLAB 为PC(Windows)版本则为真isunix MA TLAB 为Unix 版本则为真mkdir 创建目录pwd 改变当前工作目录unix 执行unix 指令并返还结果vms 执行vms dcl 指令并返还结果web 打开web 浏览器! 执行外部应用程序2 运算符和特殊算符(Operators and special characters)2.1 算术运算符(Arithmetic operators)+ 加- 减* 矩阵乘.* 数组乘^ 矩阵乘方.^ 数组乘方\ 反斜杠或左除/ 斜杠或右除./或.\ 数组除张量积[注]本表第三栏括号中的字符供在线救助时help 指令引述用2.2 关系运算符(Relational operators)= = 等号~= 不等号< 小于> 大于<= 小于或等于>= 大于或等于2.3 逻辑操作(Logical operators) & 逻辑与| 逻辑或~ 逻辑非xor 异或any 有非零元则为真all 所有元素均非零则为真2.4 特殊算符(Special characters) ：冒号( ) 圆括号[ ] 方括号{ } 花括号@ 创建函数句柄. 小数点. 构架域的关节点.. 父目录续行号, 逗号; 分号% 注释号! 调用操作系统命令= 赋值符号ˊ引号ˊ复数转置号.ˊ转置号[,] 水平串接[;] 垂直串接( ),{ },. 下标赋值( ),{ },. 下标标识subsindex 下标标识3 编程语言结构(Programming language constructs) 3.1 控制语句(Control flow)break 终止最内循环case 同switch 一起使用catch 同try 一起使用continue 将控制转交给外层的for 或while 循环else 同if 一起使用elseif 同if 一起使用end 结束for,while,if 语句for 按规定次数重复执行语句if 条件执行语句otherwise 可同switch 一起使用return 返回switch 多个条件分支try try-cathch 结构while 不确定次数重复执行语句3.2 计算运行(Evaluation and execution)assignin 跨空间赋值builtin 执行内建的函数eval 字符串宏指令evalc 执行MA TLAB 字符串evalin 跨空间计算串表达式的值feval 函数宏指令run 执行脚本文件3.3 脚本文件、函数及变量(Scripts,function,and variables) exist 检查变量或函数是否被定义function 函数文件头global 定义全局变量isglobal 若是全局变量则为真iskeyword 若是关键字则为真mfilename 正在执行的M 文件的名字persistent 定义永久变量script MA TLAB 命令文件3.4 宗量处理(Augument handling) inputname 实际调用变量名nargchk 输入变量个数检查nargin 函数输入宗量的个数nargout 函数输出宗量的个数nargoutchk 输出变量个数检查varagin 输入宗量varagout 输出宗量3.5 信息显示(Message display)disp 显示矩阵和文字内容display 显示矩阵和文字内容的重载函数error 显示错误信息fprintf 把格式化数据写到文件或屏幕lasterr 最后一个错误信息lastwarn 最后一个警告信息sprintf 按格式把数字转换为串warning 显示警告信息3.6 交互式输入(Interactive input)input 提示键盘输入keyboard 激活键盘做为命令文件pause 暂停uicontrol 创建用户界面控制uimenu 创建用户界面菜单4 基本矩阵函数和操作(Elementary matrices and matrix manipulation)4.1 基本矩阵(Elementary matrices)eye 单位阵linspace 线性等分向量logspace 对数等分向量meshgrid 用于三维曲面的分格线坐标ones 全1 矩阵rand 均匀分布随机阵randn 正态分布随机阵repmat 铺放模块数组zeros 全零矩阵: 矩阵的援引和重排4.2 矩阵基本信息(Basic array information)disp 显示矩阵和文字内容isempty 若是空矩阵则为真isequal 若对应元素相等则为1islogical 尤其是逻辑数则为真isnumeric 若是数值则为真length 确定向量的长度logical 将数值转化为逻辑值ndims 数组A 的维数size 确定矩阵的维数4.3 矩阵操作(Matrix manipulateion)blkdiag 块对角阵串接diag 创建对角阵，抽取对角向量end 数组的长度，即最大下标find 找出非零元素1 的下标fliplr 矩阵的左右翻转flipud 矩阵的上下翻转flipdim 交换对称位置上的元素ind2sub 据单下标换算出全下标reshape 矩阵变维rot90 矩阵逆时针90°旋转sub2idn 据全下标换算出单下标tril 抽取下三角阵triu 抽取上三角阵4.4 特殊变量和常数(Special variables and constants) ans 最新表达式的运算结果eps 浮点相对误差i,j 虚数单位inf 或Inf 无穷大isfinite 若是有限数则为真isinf 若是无穷大则为真isnan 若为非数则为真NaN 或nan 非数pi 3.1415926535897?.realmax 最大浮点数realmin 最小正浮点数why 一般问题的简明答案4.5 特殊矩阵(Specialized matrices) compan 伴随矩阵gallery 一些小测试矩阵hadamard Hadamard 矩阵hankel Hankel 矩阵hilb Hilbert 矩阵invhilb 逆Hilbert 矩阵magic 魔方阵pascal Pascal 矩阵rosser 典型对称特征值实验问题toeplitz T oeplitz 矩阵vander V andermonde 矩阵wilkinson Wilkinson’s 对称特征值实验矩阵5 基本数学函数(Elementary math functions)5.1 三角函数(T rigonometric) acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割asec 反正割asech 反双曲正割asin 反正弦asinh 反双曲正弦atan 反正切atanh 反双曲正切atan2 四象限反正切cos 余弦cosh 双曲余弦cot 余切coth 双曲余切csc 余割csch 双曲余割sec 正割sech 双曲正割sin 正弦sinh 双曲正弦tan 正切tanh 双曲正切5.2 指数函数(Exponential)exp 指数log 自然对数log10 常用对数log2 以2 为底的对数nestpow2 最近邻的2 的幂pow2 2 的幂sqrt 平方根5.3 复数函数(Complex)abs 绝对值angle 相角complex 将实部和虚部构成复数conj 复数共轭cplxpair 复数阵成共轭对形式排列imag 复数虚部isreal 若是实数矩阵则为真real 复数实部unwrap 相位角360°线调整5.4 圆整和求余函数(Rounding and remainder)ceil 朝正无穷大方向取整fix 朝零方向取整floor 朝负无穷大方向取整mod 模数求余rem 求余数round 四舍五入取整sign 符号函数6 特殊函数(Specialized math functions)cart2pol 直角坐标变为柱(或极)坐标cart2sph 直角坐标变为球坐标cross 向量叉积dot 向量内积isprime 若是质数则为真pol2cart 柱(或极)坐标变为直角坐标sph2cart 球坐标变为直角坐标7 矩阵函数和数值线性代数(Matrix functions－numerical linear algebra) 7.1 矩阵分析(Matrix analysis)det 行列式的值norm 矩阵或向量范数normest 估计2 范数null 零空间orth 值空间rank 秩rref 转换为行阶梯形trace 迹subspace 子空间的角度7.2 线性方程(Linear equations)chol Cholesky 分解cholinc 不完全Cholesky 分解cond 矩阵条件数condest 估计1-范数条件数inv 矩阵的逆lu LU 分解luinc 不完全LU 分解lscov 已知协方差的最小二乘积nnls 非负二乘解pinv 伪逆qr QR 分解rcond LINPACK 逆条件数\、/ 解线性方程7.3 特性值与奇异值(Eigenvalues and singular values) condeig 矩阵各特征值的条件数eig 矩阵特征值和特征向量eigs 多个特征值gsvd 归一化奇异值分解hess Hessenberg 矩阵poly 特征多项式polyeig 多项式特征值问题qz 广义特征值schur Schur 分解svd 奇异值分解svds 多个奇异值7.4 矩阵函数(Matrix functions)expm 矩阵指数expm1 矩阵指数的Pade 逼近expm2 用泰勒级数求矩阵指数expm3 通过特征值和特征向量求矩阵指数funm 计算一般矩阵函数logm 矩阵对数sqrtm 矩阵平方根7.5 因式分解(Factorization utility)cdf2rdf 复数对角型转换到实块对角型balance 改善特征值精度的平衡刻度rsf2csf 实块对角型转换到复数对角型8 数据分析和傅里叶变换(Date analysis and Fourier transforms)8.1 基本运算(Basic operations) cumprod 元素累计积cumsum 元素累计和cumtrapz 累计积分hist 统计频数直方图histc 直方图统计max 最大值mean 平均值median 中值min 最小值prod 元素积sort 由小到大排序sortrows 由小到大按行排序std 标准差sum 元素和trapz 梯形数值积分var 求方差8.2 有限差分(Finite differentces) del2 五点离散Laplaciandiff 差分和近似微分gradient 梯度8.3 相关(Correlation)corrcoef 相关系数cov 协方差矩阵subspace 子空间之间的角度8.4 滤波和卷积(Filtering and convoluteion) conv 卷积和多项式相乘conv2 二维卷积convn N 维卷积detrend 去除线性分量deconv 解卷和多项式相除filter 一维数字滤波器fliter2 二维数字滤波器8.5 傅里叶变换(Fourier transforms)fft 快速离散傅里叶变换fft2 二维离散傅里叶变换fftn N 维离散傅里叶变换fftshift 重排fft 和fft2 的输出ifft 离散傅里叶反变换ifft2 二维离散傅城叶反变换ifftn N 维离散傅里叶反变换ifftshift 反fftshift9 音频支持(Audio support)9.1 音频硬件驱动(Audio hardware drivers)sound 播放向量soundsc 自动标刻并播放waveplay 利用系统音频输出设配播放waverecor 利用系统音频输入设配录音9.2 音频文件输入输出(Audio file import and export) auread 读取音频文件(.au)auwrite 创建音频文件(.au)wavread 读取音频文件(.wav)wavwrite 创建音频文件(.wav)9.3 工具(Utilities)lin2mu 将线性信号转换为μ 一律编码的信号mu2lin 将μ 一律编码信号转换为线性信号10 插补多项式函数(Interpolation and polynomials) 10.1 数据插补(Data Interpolation)griddata 分格点数据griddata3 三维分格点数据griddatan 多维分格点数据interpft 利用FFT 方法一维插补interp1 一维插补interp1q 快速一维插补interp2 二维插补interp3 三维插补intern N 维插补pchip hermite 插补10.2 样条插补(Spline Interpolation)ppval 计算分段多项式spline 三次样条插补10.3 多项式(Polynomials)conv 多项式相乘deconv 多项式相除poly 由根创建多项式polyder 多项式微分polyfit 多项式拟合polyint 积分多项式分析polyval 求多项式的值polyvalm 求矩阵多项式的值residue 求部分分式表达roots 求多项式的根1#caidaoyiba11 数值泛函函数和ODE 解算器(Function functions and ODE solvers) 11.1 优化和寻根(Optimization and root finding)fminbnd 非线性函数在某区间中极小值fminsearch 单纯形法求多元函数极值点指令fzero 单变量函数的零点11.2 优化选项处理(Optimization Option handling)optimget 从OPTIONS 构架中取得优化参数optimset 创建或修改OPTIONS 构架11.3 数值积分(Numerical intergration)dblquad 二重(闭型)数值积分指令quad 低阶法数值积分quadl 高阶法数值积分11.4 绘图(Plotting)ezcontour 画等位线ezcontourf 画填色等位线ezmesh 绘制网格图ezmeshc 绘制含等高线的网格图ezplot 绘制曲线ezplot3 绘制3 维曲线ezpolar 采用极坐标绘图ezsurf 画曲面图ezsurfc 画带等位线的曲面图fplot 画函数曲线图11.5 内联函数对象(Inline function object) argnames 给出函数的输入宗量char 创建字符传输组或者将其他类型变量转化为字符串数组formula 函数公式inline 创建内联函数11.6 差微分函数解算器(Differential equation solvers) ode113 变阶法解方程ode15s 变阶法解刚性方程ode23 低阶法解微分方程ode23s 低阶法解刚性微分方程ode23t 解适度刚性微分方程odet23tb 低阶法解刚性微分方程ode45 高阶法解微分方程12 二维图形函数(Two dimensional graphs)12.1 基本平面图形(Elementary X－Y graphs) loglog 双对数刻度曲线plot 直角坐标下线性刻度曲线plotyy 双纵坐标图polar 极坐标曲线图semilogx X 轴半对数刻度曲线semilogy Y 轴半对数刻度曲线12.2 轴控制(Axis control)axes 创建轴axis 轴的刻度和表现box 坐标形式在封闭式和开启词式之间切换grid 画坐标网格线hold 图形的保持subplot 创建子图zoom 二维图形的变焦放大12.3 图形注释(Graph annotation)gtext 用鼠标在图上标注文字legend 图例说明plotedit 图形编辑工具text 在图上标注文字texlabel 将字符串转换为T ex 格式title 图形标题xlabel X 轴名标注ylabel Y 轴名标注12.4 硬拷贝(Hardcopy and printing)orient 设置走纸方向print 打印图形或把图存入文件printopt 打印机设置13 三维图形函数(Three dimensional graphs) 13.1 基本三维图形(Elementary 3－D plots) fill3 三维曲面多边形填色mesh 三维网线图plot3 三维直角坐标曲线图surf 三维表面图13.2 色彩控制(Color control)alpha 透明色控制brighten 控制色彩的明暗caxis (伪)颜色轴刻度colordef 用色风格colormap 设置色图graymon 设置缺省图形窗口为单色显示屏hidden 消隐shading 图形渲染模式whitebg 设置图形窗口为白底13.3 光照模式(Lighting)diffuse 漫反射表面系数light 灯光控制lighting 设置照明模式material 使用预定义反射模式specular 漫反射surfnorm 表面图的法线surfl 带光照的三维表面图13.4 色图(Color maps)autumn 红、黄浓淡色bone 蓝色调灰度图colorcube 三浓淡多彩交错色cool 青和品红浓淡色图copper 线性变化纯铜色调图flag 红－白－蓝黑交错色图gray 线性灰度hot 黑－红－黄－白交错色图hsv 饱和色彩图jet 变异HSV 色图lines 采用plot 绘线色pink 淡粉红色图prism 光谱色图spring 青、黄浓淡色summer 绿、黄浓淡色vga 16 色white 全白色winter 蓝、绿浓淡色13.5 轴的控制(Axis control)axes 创建轴axis 轴的刻度和表现box 坐标形式在封闭式和开启式之间切换daspect 轴的DataAspectRatio 属性grid 画坐标网格线hold 图形的保持pbaspect 画坐标框的PlotBoxAspectRatio 属性subplot 创建子图xlim X 轴范围ylim Y 轴范围zlim Z 轴范围zoom 二维图形的变焦放大13.6 视角控制(V iewpoint control)rotate3d 旋动三维图形view 设定3－D 图形观测点viewmtx 观测点转换矩阵13.7 图形注释(Graph annotation)colorbar 显示色条gtext 用鼠标在图上标注文字plotedit 图形编辑工具text 在图上标注文字title 图形标题xlabel X 轴名标注ylabel Y 轴名标注zlabel Z 轴名标注13.8 硬拷贝(Hardcopy and printing)orient 设置走纸方向print 打印图形或把图存入文件printopt 打印机设置verml 将图形保存为VRML2.0 文件14 特殊图形(Specialized graphs)14.1 特殊平面图形(Specialized 2－D graphs) area 面域图bar 直方图barh 水平直方图comet 彗星状轨迹图compass 从原点出发的复数向量图errorbar 误差棒棒图ezplot 画二维曲线ezpolar 画极坐标曲线feather 从X 轴出发的复数向量图fill 多边填色图fplot 函数曲线图hist 统计频数直方图pareto Pareto 图pie 饼形统计图plotmatrix 散点图阵列scatter 散点图stairs 阶梯形曲线图stem 火柴杆图14.2 等高线及二维半图形(Contour and 2－1/2D graphs) clabel 给等高线加标注contour 等高线图contourf 等高线图contour3 三维等高线ezcontour 画等位线ezcontourf 画填色等位线pcolor 用颜色反映数据的伪色图voronoi V oronoi 图14.3 特殊三维图形（Specialized 3－D graphs）bar3 三维直方图bar3h 三维水平直方图comet3 三维彗星动态轨迹线图ezgraph3 通用指令ezmesh 画网线图ezmeshc 画等位线的网线图ezplot3 画三维曲线ezsurf 画曲面图ezsurfc 画带等位线的曲面图meshc 带等高线的三维网线图meshz 带零基准面的三维网线图pie3 三维饼图ribbon 以三维形式绘制二维曲线scatter3 三维散点图stem3 三维离散杆图surfc 带等高线的三维表面图trimesh 三角剖分网线图trisurf 三角剖分曲面图waterfall 瀑布水线图14.4 内剖及向量视图(V olume and vector visualization) coneplot 锥体图contourslice 切片等位线图quiver 矢量场图quiver3 三维方向箭头图slice 切片图14.5 图像显示及文件处理(Image display and file I/O) brighten 控制色彩的明暗colorbar 色彩条状图colormap 设置色图contrast 提高图像对比度的灰色图gray 线性灰度image 显示图像imagesc 显示亮度图像imfinfo 获取图像文件的特征数据imread 从文件读取图像的数据阵（和伴随色图））imwrite 把强度图像或真彩图像写入文件14.6 影片和动画(Movies and animation)capture 当前图的屏捕捉frame2im 将影片动画转换为编址图像getframe 获得影片动画图像的帧im2frame 将编址图像转换为影片动画movie 播放影片动画moviein 影片动画内存初始化rotate 旋转指令14.7 颜色相关函数(Color related function)spinmap 颜色周期性变化操纵14.8 三维模型函数(Solid modeling)cylinder 圆柱面patch 创建块sphere 球面Surf2patch 将曲面数据转换为块数据15 句柄图形(Handle Graphics)15.1 图形窗的产生和控制(Figure window creation and control) clf 清除当前图close 关闭图形figure 打开或创建图形窗口gcf 获得当前图的柄openfig 打开图形refresh 刷新图形shg 显示图形窗15.2 轴的产生和控制(Axis creation and control) axes 在任意位置创建轴axis 轴的控制box 坐标形式在封闭式和开启式之间切换caxis 控制色轴的刻度cla 清除当前轴gca 获得当前轴的柄hold 图形的保持ishold 若图形处保持状态则为真subplot 创建子图15.3 句柄图形对象(Handle Graphics objects) axex 在任意位置创建轴figure 创建图形窗口image 创建图像light 创建光line 创建线patch 创建块rectangle 创建方surface 创建面text 创建图形中文本uicontextmenu 创建现场菜单对象uicontrol 用户使用界面控制uimenu 用户使用菜单控制15.4 句柄图形处理(Handle Graphics operations) copyobj 拷贝图形对象及其子对象delete 删除对象及文件drawnow 屏幕刷新findobj 用规定的特性找寻对象gcbf “正执行回调操作”的图形的柄gcbo “正执行回调操作”的控件图柄指令gco 获得当前对象的柄get 获得对象特性getappdat 获得应用程序定义数据isappdata 检验是否应用程序定义数据reset 重设对象特性rmappdata 删除应用程序定义数据set 建立对象特性setappdata 建立应用程序定义数据15.5 工具函数(Utilities)closereq 关闭图形窗请求函数ishandle 若是图柄代号侧为真newplot 下一个新图16 图形用户界面工具(Graphical user interface tools) align 对齐用户控件和轴cbedit 编辑回调函数ginput 从鼠标得到图形点坐标guide 设计GUImenu 创建菜单menuedit 菜单编辑propedit 属性编辑uicontrol 创建用户界面控制uimenu 创建用户界面菜单2#caidaoyiba17 字符串(Character string)17.1 通用字符串函数(General)blanks 空格符号cellstr 通过字符串数组构建字符串的元胞数组char 创建字符传输组或者将其他类型变量转化为字符串数组deblank 删除最后的空格double 把字符串变成ASCII 码值eval 执行串形式的MATLAB 表达式17.2 字符串查询(String tests)iscellstr 若是字符串组成的元胞数组则为真ischar 若是字符串则为真isletter 串中是字母则为真isspace 串中是空格则为真isstr 若是字符串则为真17.3 字符串操作(String operations)base2dec X－进制串转换为十进制整数bin2dec 二进制串转换为十进制整数dec2base 十进制整数转换为X 进制串dec2bin 十进制整数转换为二进制串dec2hex 十进制整数转换为16 进制串findstr 在一个串中寻找一个子串hex2dec 16－进制串转换为十进制整数hex2num 16－进制串转换为浮点数int2str 将整数转换为字符串lower 把字符串变成小写mat2str 将数组转换为字符串num2str 把数值转换为字符串strcat 把多个串连接成长串strcmp 比较字符串strcmpi 比较字符串（忽略大小写）strings MA TLAB 中的字符串strjust 字符串的对齐方式strmatch 逐行搜索串strnomp 比较字符串的前N 个字符strncmpi 比较字符串的前N 个字符（忽略大小写）strrep 用另一个串代替一个串中的子串strtok 删除串中的指定子串strvcat 创建字符串数组str2mat 将字符串转换为含有空格的数组str2num 将字符串转换为数值upper 把字符串变成大写18 文件输入/输出(File input/output)clc 清除指令窗口disp 显示矩阵和文字内容fprintf 把格式化数据写到文件或屏幕home 光标返回行首input 提示键盘输入load 从磁盘中调入数据变量pause 暂停sprintf 写格式数据到串sscanf 在格式控制下读串19 时间和日期(T ime and dates)clock 时钟cputme MA TLAB 战用CPU 时间date 日期etime 用CLOCK 计算的时间now 当前时钟和日期pause 暂停tic 秒表启动toc 秒表终止和显示20 数据类型(Data types and structures)20.1 数据类型(Data types)cell 创建元胞变量char 创建字符传输组或者将其他类型变量转化为字符串数组double 转化为16 位相对精度的浮点数值对象function handle 函数句柄inline 创建内联函数JavaArray 构建Java 数组JavaMethod 调用某个Java 方法JavaObject 调用Java 对象的构造函数single 转变为单精度数值sparse 创建稀疏矩阵struct 创建构架变量uint8(unit16、unit32) 转换为8(16、32)位无符号整型数int8(nit16、nit32) 转换为8(16、32)位符号整型数20.2 多维数组函数(Multi－dimensional array functions) cat 把若干数组串接成高维数组ndims 数组A 的维数ndgrid 为N－D 函数和插补创建数组ipermute 广义反转置permute 广义非共轭转置shiftdim 维数转换squeeze 使数组降维20.3 元胞数组函数(Cell array functions)cell 创建元胞变量celldisp 显示元胞数组内容cellfun 元胞数组函数cellplot 图示元胞数组的内容cell2struct 把元胞数组转换为构架数组deal 把输入分配给输出is cell 若是元胞则为真num2 cell 把数值数组转换为元胞数组struct2 cell 把构架数组转换为元胞数组20.4 构架函数(Structure functions)fieldnames 获取构架的域名getfield 获取域的内容isfield 若为给定构架的域名则为真isstruct 若是构架则为真rmfield 删除构架的域setfield 指定构架域的内容struct 创建构架变量20.5 函数句柄函数(Function handle functions)@ 创建函数句柄functions 列举函数句柄对应的函数func2str 将函数句柄数组转换为字符串str2func 将字符串转换为函数句柄20.6 面向对象编程(Object oriented programming functions) dlass 查明变量的类型isa 若是指定的数据类型则为真inferiorto 级别较低isjava 若是java 对象则为真isobject 若是对象则为真methods 显示类的方法名substruct 创建构架总量superiorto 级别较高21 示例(E xamples and demonstrations)demo 演示程序flow 无限大水体中水下射流速度数据intro 幻灯演示指令peaks 产生peaks 图形数据22 符号工具包(Symbolic Math T oolbox) 22.1 微积分(Calculus)diff 求导数limit 求极限int 计算积分jacobian Jacobian 矩阵symsum 符号序列的求和trylor T rylor 级数22.2 线性代数(Linear Algebra)det 行列式的值diag 创建对角阵，抽取对角向量eig 矩阵特征值和特征向量expm 矩阵指数inv 矩阵的逆jordan Jordan 分解null 零空间poly 特征多项式rank 秩rref 转换为行阶梯形svd 奇异值分解。
matlab语言常用算法程序集
matlab语言常用算法程序集MATLAB是一种常用的数值计算和科学计算语言，常用于工程、科学、统计分析和数据可视化等领域。
其功能强大且易于学习和使用，同时还提供了丰富的算法和函数库，用于实现各种常用算法。
以下是一些MATLAB常用算法程序集。
1.线性代数算法：MATLAB提供了一系列强大的线性代数函数，用于矩阵运算和线性方程求解。
包括矩阵乘法、求逆矩阵、求解线性方程组、特征值和特征向量计算等。
2.优化算法：MATLAB提供了多种优化算法，用于寻找函数的最小值或最大值。
包括最优化问题的求解、约束优化、非线性优化、整数规划等。
3.插值与拟合算法：MATLAB提供了多种插值和拟合函数，用于通过已知数据点求未知点的值或者近似函数。
常用的算法包括线性插值、多项式拟合、样条插值、最小二乘拟合等。
4.数值积分与微分算法：MATLAB提供了多种数值积分和微分函数，用于计算函数的定积分或求解微分方程。
包括常用的数值积分法（如梯形法、辛普森法）、微分方程求解（如欧拉法、四阶龙格库塔法）等。
5.图像处理算法：MATLAB提供了图像处理工具箱，包括图像读取、显示、滤波、变换等功能。
常用的算法包括图像平滑、边缘检测、图像分割、形态学运算等。
6.信号处理算法：MATLAB提供了信号处理工具箱，用于处理和分析信号数据。
包括滤波、傅里叶变换、频谱分析、信号降噪等。
7.统计分析算法：MATLAB提供了丰富的统计分析工具，用于数据分析和建模。
包括描述性统计、假设检验、ANOVA分析、回归分析、时间序列分析等。
8.机器学习算法：MATLAB提供了机器学习工具箱，用于实现各种机器学习算法。
包括监督学习算法（如线性回归、逻辑回归、支持向量机）、无监督学习算法（如聚类分析、主成分分析）等。
9.计算几何算法：MATLAB提供了计算几何工具箱，用于解决各种计算几何问题。
包括点线面的计算、多边形和曲线的操作、三维几何变换等。
10.数值分析算法：MATLAB提供了多种数值分析函数，用于数值计算和近似求解。
MATLAB如何使用-教程-初步入门大全资料
运算 数学表达式
加 a+b
减
a-b
乘 a×b
除 a÷ b
幂 a^b
MATLAB运算符
+ *
/(右除)或\(左除)
^
MATLAB表达式
a+b a-b a*b a/b或b\a a^b
示例
1+2 5-3 2*3
6/2或2\6 2^3
指出：右除相当于通常的除法。
22
七、MATLAB的变量与函数
1、变量 变量就是在程序的运行过程中，其数值可以变化的量
MATLAB是交互式的语言，输入命令即给出运算结 果。而命令窗口则是MATLAB的主要交互窗口，用 于输入和编辑命令行等信息，显示结果（图形除 外）。
当命令窗口中出现提示符“>>”时，表示MATLAB已 经准备好，可以输入命令、变量或运行函数。提示 符总是位于行首。
在每个指令行输入后要按回车键，才能使指令被 MATLAB执行。
28
矩阵的创建（续）
1、直接输入法－在命令窗口按规则输入方式创建矩阵
例1.在命令窗口创建简单的数值矩阵。
>>A=[1 3 2;3 1 0;2 1 5] 回车后在命令窗口显示如下结果
A=
132
310
215 例2.在命令窗口创建带运算表达式的矩阵，不显示结果。
>>y=[sin(pi/3),cos(pi/6);log(20),exp(2)]; 输入“y”回车，在命令窗口显示出来。
（3）在MATLAB安装目录\MATLAB6p5中双击 MATLAB快捷方式。
（4）在MATLAB安装目录\MATLAB6p5\bin\win32 中双击MATLAB.exe图标。
MATLAB操作命令大全
MATLAB操作命令大全1.基本操作- help：查看函数的帮助文档。
- save：将变量保存到文件中。
- load：从文件中加载变量。
- clear：清除当前工作空间中的变量。
- who：列出当前工作空间中的变量。
- whos：显示当前工作空间中变量的详细信息。
- quit：退出MATLAB。
2.变量操作-=：赋值操作，将值赋给变量。
- disp：显示变量的值。
- length：返回数组的长度。
- size：返回数组的大小。
- max：返回数组的最大值。
- min：返回数组的最小值。
- sum：返回数组元素的和。
3.数学操作-+：加法操作，将两个数值相加。
--：减法操作，将两个数值相减。
-*：乘法操作，将两个数值相乘。
-/：除法操作，将两个数值相除。
-^：指数操作，将一个数值提高到指定次幂。
- sqrt：返回一个数值的平方根。
- abs：返回一个数值的绝对值。
4.矩阵操作- eye：创建一个单位矩阵。
- zeros：创建一个全0矩阵。
- ones：创建一个全1矩阵。
- rand：创建一个0到1之间的随机矩阵。
- diag：返回对角线元素。
- inv：返回矩阵的逆矩阵。
- det：返回矩阵的行列式。
5.图形操作- plot：绘制二维线图。
- scatter：绘制散点图。
- bar：绘制柱状图。
- hist：绘制直方图。
- surf：绘制三维曲面图。
- contour：绘制等高线图。
- imagesc：绘制矩阵的颜色图。
6.控制流程操作- if：用于条件判断。
- for：用于循环操作。
- while：用于循环操作。
- switch：用于多条件判断。
- break：跳出循环。
- continue：跳过当前循环，并继续执行下一次循环。
7.文件操作- fopen：打开文件。
- fclose：关闭文件。
- fprintf：将数据写入文件。
- fscanf：从文件中读取数据。
- fseek：设置文件指针的位置。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
IIR滤波程序function y=iirfilter(b,a,x)nf=length(a);c=zeros(nf-1,nf-1);for i=nf-1:-1:1for j=1:ic(i,j)=b(i-j+1);d(i,j)=-a(i-j+1);endendd(1,1)=0;y(1)=b(1)*x(1);for i=2:nf-1ss=0;for j=2:iss=ss+x(j)*c(i,j)+y(j-1)*d(i,j-1);endy(i)=ss+x(1)*c(i,1);enda(1)=0;for i=nf:length(x)sum=0;for j=2:nfsum=sum+b(j)*x(i-j+1)-a(j)*y(i-j+1);endy(i)=sum+b(1)*x(i);end读写二进制文件写文件file_out = fopen('tdip.b', 'wb');fwrite(file_out,y,'double');%int8,int16,int24,double,float读文件[FileName,PathName] = uigetfile('*.*','Select the M-file');%选择需要处理的数据文件%% % 打开方式为二进制打开，其实'r'就行，matlab是默认二进制形式打开文件的file_id = fopen(FileName, 'rb');%设置输出文件名字jump_distance = 0; % 这里是跳过的字节数fseek(file_id, jump_distance, 'bof');i=0;week=0;week0=headt(7);while ~feof(file_id)% while(i<3)[data, ele_count] = fread(file_id, [4,fs/10] , 'int32' ) ;End批量处理文件编写程序时经常会遇到处理大量数据文件的情况，数据文件的命名比较相似，以一系列编号区分，如“a1.mat，a2.mat，... ，a100.mat " (假设其中的数据名称也为a1...a100) 。
为了可以批量处理这些文件，可以采用如下的方法：for i=1:1:100s=strcat('a', int2Str(i));load(strcat(s,'.mat'));x=eval(s); %将数据a1赋值给x，便于后继统一处理% ... %统一处理程序end写txt文件fid=fopen(‘1.txt’,’w t’);fprintf(fid,‘这是一个测试文件：\n’)fprintf(fid,’%f’,a)m=[A nskip nfft fxb real(pxb) imag(pxb)];fid = fopen('pxb.txt', 'at');fprintf(fid,'振幅开始时间fft点数激励频率pxb:\n');fprintf(fid,' %d %d %d %6.8f %6.8f %6.8f\n\n', m);% fprintf(fid,'l1 l2 l3 l4:\n')% fprintf(fid, ' %d %d %d %d %d %d %d %d\n ', l');fclose(fid);%a=’00:00:00’;% fprintf(fid,'%s ',a);用cell生成自动变量名并幅值clcclearfor i=1993:1996matName{i-1992}=strcat('Y',num2str(i)); % str=[num2str(i-1992) '.csv'];% matName{i-1992}=xlsread(str); matName{i-1992}=rand(1,i-1992);% eval([matName{i-1992} ' = rand(1,i);']); endmatName{1}(1)关于errorbar误差线：体现数据点置信水平和偏差MATLAB中errorbar的使用example#上下误差相等X = 0:pi/10:pi;Y = sin(X);E = 0:0.05:0.5;figureerrorbar(X,Y,E,'r');hold onplot(X,Y,'-mo',...'LineWidth',2,...'MarkerEdgeColor','k',...'MarkerFaceColor',[.49 1 .63],... 'MarkerSize',12);xlabel('x');ylabel('sin(x)');#上下误差不相等X = 0:pi/10:pi;Y = sin(X);L = 0:0.05:0.5;U=0.5:-0.05:0;figureerrorbar(X,Y,L,U,'r');hold onplot(X,Y,'-mo',...'LineWidth',2,...'MarkerEdgeColor','k',...'MarkerFaceColor',[.49 1 .63],... 'MarkerSize',12);xlabel('x');ylabel('sin(x)');误差棒端点拉长的问题（/blog/?p=45）在Matlab中有个画误差棒的函数errorbar，用在对数坐标下时，则会出现误差棒端点拉长的问题，严重影响画面（如下图左半轴所示）。
这个问题开始是在Matlab6.5中发现的，2006b中虽对此函数有所更新但并未解决这个问题。
幸运的是，这个函数没有嵌入内核，可以在命令窗口键入type errorbar ，直接察看其源代码。
找到定义误差棒端点长度的地方tee = (max(x(:))-min(x(:)))/100;% make tee .02 x-distance for error bars取误差棒端点大小为误差棒长度的2%xl = x – tee; 计算左端点xr = x + tee; 右端点找对地方，改起来就很简单了，在对数坐标中，点A变成了logA，要取到同一点，将A变成10^A 即可，后两行改为xl = x/10^tee;xr = x*10^tee;改名保存到工作目录，直接调用，效果如下图右半轴。
如果在这里加上个判断，就可以直接去目录MATLAB\R2006b\toolbox\matlab\specgraph下替换原始文件了:)时间做坐标轴clearclcx=load('a.txt');n=length(x);for i=1:nwhile(x(i,1)<10^5)x(i,1)=x(i,1)*10;endfor j=2:5while(x(i,j)<10^6)x(i,j)=x(i,j)*10;endendendx=x/10^6;ts=datenum('2013-03-25 19:30:06');% ¿ªÊ¼Ê±¼ätf=datenum('2013-03-26 08:06:06');% ½áÊøʱ¼ät=linspace(ts,tf,length(x));% 21minΪ×ܹ²µÄʱ¼äcolor=['r','g','m','k'];figure(1)plot(t,x(:,1),'LineWidth',2)axis([ts tf 0 4.5])set(gca,'Xtick',[ts:(tf-ts)/6:tf])datetick('x','yyyy/mm/DD HH:MM:SS','keepticks')title('1-5','Color','k','FontSize',16);ylabel('mV','FontName','Times New Roman','FontSize',20)for i=2:5hold onplot(t,x(:,i),color(i-1),'LineWidth',2)endlegend('SB','SC','SD','C1','C2')plot详细semilogx(p_result(1:f_num-1,3),p_result(1:f_num-1,1),'-bo',...'LineWidth',2.2,...'MarkerEdgeColor','k',...'MarkerFaceColor',[0 0 1],...'MarkerSize',6);% set(gca,'xlim',[1 floor(length(p))/2]);% set(gca,'ylim',[-10 35]);hold onsemilogx(p_result(1:f_num-1,3),p_result(1:f_num-1,2),'-bo',...'LineWidth',2.2,...'MarkerEdgeColor','k',...'MarkerFaceColor',[0 0 1],...'MarkerSize',6);添加当前文件夹及其子文件夹到路径currentFolder = pwd;addpath(genpath(currentFolder))在MATLAB中得到系统当前日期、时间也是经常用到的内容，由以下函数实现。