MATLAB命令汇总.doc

已学MATLAB命令汇总一、MATLAB初步1 基本操作 显示变量名 who 显示变量信息 whos清除变量 clear 保存变量 save装载变量 load 关闭图形 close清除图形 clf 设置显示格式 format帮助 help 注释 % 续行 演示 demo2 基本常数： pi I j inf eps NaN exp(1) 3 算术运算： + - * / sqrt .* ./ . 4 内部函数：sin(x) tan(x) asin(x) atan(x) abs(x) round(x) floor(x) ceil(x) log(x) sign(x) length(v) size(m) y, k=sort(x)二、 矩阵与数组1矩阵与数组生成 1 2 3;4 5 6 linspace(a,b,n) eye(n) zeros(m,n) ones(m,n) rand(m,n)magic(n) hilb(n) diag(A) diag ( a1, a2, . . . an )2矩阵与数组运算 A+B A-B A*B Ab b/A Ak A.*B A./B A.2 det(A) dot(a,b)cross(a,b) rref(A) trace(A) inv(A) rank(A) eig(A) null(A) orth(A) 注意：1A， sin(A)3矩阵与数组操作 A(i, j) =a B=A(1,3, 2,1) A(2, :)=1 2 3 A. A find(A0)A(A1) B=(A=1) flipud(A) tril(A,-k) triu(A, k)三、函数作图1平面曲线 plot(x,y) fplot(fun,a,b)ezplot(f, a,b) ezplot(x,y, a,b)2空间曲线 plot3(x,y,z) ezplot3(x,y,z,a,b)3空间曲面 mesh(x,y,z) meshz(x,y,z) meshc(x,y,z) surf(x,y,z) ezmesh(z,a,b,c,d) ezmesh(x,y,z,a,b,c,d)ezsurf(z,a,b,c,d) ezsurf(x,y,z,a,b,c,d)4动画 m(k)=getframe; movie(m, n);5画面安排 hold on; hold off; sunbplot; figure; 6旋转曲面 x,y,z=cylinder(r, n); 7球面 x,y,z=sphere(n)8椭球面 x,y,z=ellipsoid(x0, y0, z0, a, b, c, n)四、编程初步1 M文件：(1) 脚本文件（命令汇集，不调用参数）； (2) 函数文件（调用参数） 子函数 全局变量( global: k c ) 2 函数定义与调用(1) 函数文件： function z=fun(x,y)z=x.3+y.3-3*x.*y;(2) 内嵌函数： fun=inline( x.3+y.3-3*x.*y,x,y);(3) 匿名函数： fun=(x) x.3+y.3-3*x.*y;3三种结构(1) 顺序结构 (2) 循环结构 for end; while expression;(3) 分支结构 if else end; switchcase1case2otherwiseend五、符号运算符号变量 syms a b c d x 符号表达式 y = a*x2+b*x+c求反函数 g=finverse(f, x)函数复合 compose(f,g,x,y) 函数展开 expand(f)因式分解 factor(f)化简 simplify(f); simple(f) 符号矩阵 A=a+b c-d; x y 逆矩阵 inv(A) 行列式 det(A) 特征值 eig(A)变量代换 subs(y, x, a); subs(fun, x, y, 2, 5 )符号转数值 double(x)转字符串 char(x)类型转化 str2num; num2str 字符串求值 eval(s)函数求值 feval(f)计算精度 vap(f, n);优美格式 pretty(f)极限 limit(f, x, a); 左极限 limit(f, x, a,left); 右极限 limit(f, x, a,right)求和 s=symsym(f, n, n1, n2)微分 s=diff(f, var, n)不定积分 s=int(f, x);定积分 s=int(f, x, a, b)重积分 s=int (int (f, y, c, d), x, a, b) 泰勒级数 s=tailor(f, n, x0, x)解代数方程 x=solve(fun, x)解微分方程 t,y=dsolve(odefun, y0) 六、矩阵代数特定多项式 p=poly(x); 其中x=x1, x2, , xn多项式赋值 polyval(p,a)多项式分解 f = roots(p)多项式加减 p1+p2；p1-p2 多项式相乘 conv(p1,p2) 多项式相除 q,r=deconv(p1,p2) 多项式求导 polyder(p,k) 向量的模 norm(x) 矩阵最简型 rref(A) Ax=0的基础解系 null(A); 或有理基础解系 null(A,r)Ax=b的特解 x=AB （或最小二乘解） A的正交规范基 orth(A) A的列向量正交规范化 特征多项式 poly(A) 特征向量特征值 u, v=eig(A)Jordan标准形 jordan(A)七、方程求根多项式方程求根 x = roots(p)一般方程求根 x = fsolve(fun, x0); x = fzero(fun, x0); 迭代法求根： f(x) = 0 x= g(x) 一般迭代： xk+1= g(xk) （收敛条件：| g(x0) | 1e-6 c=(a+b)/2; if f(a)*f(c)0 a=c;else a=c;b=cende=e/2; k=k+1;end x=(a+b)/2; x, k八、数值微积分 差商 diff(y) ./ diff(x) (偏)导数 fx=gradient(f,x); fx,fy=gradient(f,x,y)向量场 quiver(x,y, fx,fy)定积分 trapz(y) quad(fun,a,b) quadl(fun,a,b)重积分 dblquad(fun,a,b,c,d) triplequad(fun,a,b,c,d,e,f)九、微分方程数值解 t,x=ode45(odefun,t1,t2,x0); 或ode23(odefun,t1,t2,x0);欧拉方法： y(k+1)=y(k)+h*f(x(k),y(k);改进欧拉方法： z=y(k)+h*f(x(k),y(k); y(k+1)=y(k)+h*( f ( x(k), y(k) )+f ( x(k+1), z ) )/2;龙格库达方法： z1=f(x(k), y(k); z2=f(x(k)+h/2, y(k)+h*z1/2); z3=f(x(k)+h/2, y(k)+h*z2/2); z4=f(x(k)+h, y(k)+h*z3); y(k+1)= y(k)+h*(z1+2*z2+2*z3+z4)/6;十、插值与拟合 一元插值 yi=interp1(x,y,ximethod) method： linear, cubic, spline 样条插值 pp=spline(x,y); yi=ppval(pp,xi)画样条图 fnplt(pp)二元插值 zi=interp2(x,y,z,xi,yi,method) method： linear, cubic, spline 杂乱数据插值 zi=griddata(x,y,z,xi,yi,method) method： 同上多项式拟合 p=polyfit(x, y, k); yi=polyval(p, xi)最小二乘拟合 lsqcurvefit(fun, c0, x, y); y = fun(c, x) 十一、优化问题有约束一元函数最小值 x,fval = fminbnd( fun, x1, x2 )无约束多元函数最小值 x,fval = fminsearch(fun, x0)线性规划 x, fval = linprog( f, A, b, Aeq, beq, lb, ub, options)二次规划 x, fval = quadprog( H, f, A, b, Aeq, beq, lb, ub, x0 )非线性规划 x,fval = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)极小极大 x,fval = fminimax(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)十二、统计分析 最大、最小 m1=max(data) m2=min(data) 和、积 s=sum(data) s=prod(data) 累积和、积 s=cumsum(data) s=cumprod(data) 样本均值 m=mean(data) 修正均值 m=trimmean(data, percent)样本中值 m=median(data)样本标准差 d=std(data)样本方差 d2=var(data)数据排序 y=sort(data); y, k=sort(data) 直方图 hist(data, k)随机矩阵 rand(m,n); 随机置换 randperm(n)正态检验函数 normplot(data)正态分布密度 p=normpdf(x,