Matlab常用指令总结.docxVIP

概论format long显示更多位数format short显示少位数format hex将数字显示为十六进制浮点数，（3fb999999999999a，a-f代表十六进制数，前三字符为）double()将数字转化为双精度浮点数ezplot(f,0,4)绘制f函数[0,4]图像plot(x,y)绘制点或者连线zeros(n,1)产生n维0向量fibonacci(n)产生N个斐波那契数求解矩阵sum(A)对每一列的矩阵元素求和a’为将矩阵A转置sum(diag(A))矩阵Ａ主对角元素之和flipud(A)将A翻过来（上下）det(A)求A行列式的解inv(A)求A的逆矩阵norm(A)矩阵范数eig(A)特征值svd(A)奇异值a=A(:,[1 3 2 4]) 交换矩阵二、三列AX=B可以为X=A\BxA=B可以为X=B/Anorm(x,1)计算x的一阶范数即所有数的和norm(x,2)方均根norm(x,inf)取最大值插值v=polyinterp(x,y,u)（基于拉格朗日）可以计算xy组成的点阵插值结果，u为输出v对应的取得x 的值u=[,]完整次数插值symx=sym(‘x’)创建符号变量symx=polyinterp(x,y,symx)利用符号变量求解pretty(p)让p变得规整simplify(p)化简pv=piecelin(x,y,u)线性分段插值v=pchip(x,y,u)v=pchiptx(x,y,u)为分段三次埃米特插值及其简化版v=spline(x,y,u)v=splinetx(x,y,u)为三次样条插值及其简化版interpgui(x,y)可以绘制上述四种插值图像，直观方程求根【M=2a=1b=2k=0;while b-aepsx=(a+b)/2;if x^2 Mb=xelsea=xendk=k+1;end】牛顿法，求一次导法k=0;x=10;/初值xprev=12;/初值的初值while abs(x-xprev)eps*abs(x)xprev=x;x=0.5*(x+2/x);/x-f(x)/ f(x)的导数。k=k+1;endfzero(f(x),初值) /将二分法和割线法和IQI算法收敛速度结合fzerotx(f(x),[,])简化版本（fzerogui(f(x),[,])）fmintx(f(x),a,b) ab为区间，求区间最小值最小二乘法b1= polyval(a1,x)将x值带入方程a1c=polyfit(x,y,n) 将x,y带入并用N次多项式拟合，输出高次到第次。X = lsqcurvefit(FUN,X0,XDATA,YDATA)FUN:用于拟合的函数X0:迭代初始值XDATA:要拟合的所有点的X的坐标YDATA:要拟合的所有点的Y的坐标最小二乘法手算公式在doc 11/13f=inline(a(1)*x.^2+a(2)*x+a(3),a,x)建立多元方程的简便方法。ff=@(x,y) x^2+y^2;积分：y=quad(@aaa,0,1) 辛普森法则应用于0-1区间函数要加.y=quadtx(f,0,1) 简化版int(‘fun’,a,b) int(sym(x^2),0,1)直接出公式的，如果输入数则出结果，输入syms a b 则出公式diff(a*x^2,x)数值微分diff(x^2,2) 为二次导ode23 显示的龙格-库塔公式ode45 比较四阶和五阶的公式p9欧拉法：【其实就是循环语句构成X=(3:0.01:3.2);k=1:1:20;Y(k)=0;Y(1)=-1;for j=1:20Y(j+1)=Y(j)+0.01*(X(j)^2+Y(j)^2);endplot(X,Y);】[x,y]=ode23(ff,[x1,x2],y0)可以直接显示出x和y的值高阶常微分解见第七章p19顶部和函数my fun3边界初值问题见p21底部编程，绘图，设置坐标轴，做出自己的坐标轴1. axis([xminxmaxyminymax])设置当前图形的坐标范围，分别为x轴的最小、最大值，y轴的最小最大值2. V=axis返回包含当前坐标范围的一个行向量3. axis auto将坐标轴刻度恢复为自动的默认设置4. axis manual冻结坐标轴刻度，此时如果hold被设定为on，那么后边的图形将使用与前面相同的坐标轴刻度范围5. axis tight将坐标范围设定为被绘制的数据范围6. axis fill这是坐标范围和屏幕的高宽比，使得坐标轴可以包含整个绘制的区域。该选项只有在PlotBoxaApectRatio或DataAspectRatioMode被设置为‘manual’模式才有效7. ax