胡良剑 matlab 数学实验常用指令整理

matlab常用指令MATLAB是一款非常实用的科学计算软件，在使用过程中，一些常用的指令是非常必要的。

在本篇文章中，我们将会介绍MATLAB常用指令，以使你更加熟练掌握MATLAB的使用。

一、基本数学运算+ 加- 减* 乘/ 除^ 幂（指数）sqrt 平方根exp 取指数log 取自然对数log10 取以10为底的对数sin 正弦cos 余弦tan 正切asin 反正弦acos 反余弦atan 反正切abs 绝对值rem 模运算fix 向零取整floor 向负无穷取整ceil 向正无穷取整round 四舍五入mod 取摸余数二、变量与矩阵1、赋值：通过等号将数值赋给变量，如：a=3;b=2.1;c=2+3i;2、数列：建立一个等差数组，例如：d=1:10; %1到10的等差数列e=linspace(0,2*pi,100); %0到2*pi之间的100个等间距点 a=[1 2 3;4 5 6;7 8 9];b=zeros(2,3);c=ones(3,2);d=rand(3,3);e=eye(4);4、矩阵元素操作：通过下标访问矩阵中的元素，例如：a(1,2) %输出a矩阵第一行第二列的元素b(2,3)=7 %将b矩阵第二行第三列的元素赋为75、矩阵运算：矩阵加减乘除，如：a+b %对应元素相加a-b %对应元素相减a*b %矩阵乘法a/b %矩阵除法a' %矩阵转置6、矩阵函数：除了使用基本操作外，还能使用各种矩阵相关函数完成矩阵计算，例如：inv(a) %矩阵求逆det(a) %矩阵求行列式trace(a) %矩阵求迹eig(a) %求特征值rank(a) %矩阵的秩size(a) %返回矩阵的大小max(a) %求矩阵元素最大值min(a) %求矩阵元素最小值sum(a) %求矩阵元素的和prod(a) %求矩阵所有元素的乘积mean(a) %求矩阵元素的平均值三、绘图1、二维绘图：绘制二维函数的曲线、散点图等，例如：x=linspace(-3,3,100); %生成-3到3之间的100个等间距点y=sin(x);plot(x,y); %绘制正弦函数曲线plot(x,y,'r--'); %绘制红色的正弦函数曲线，形状为虚线xlabel('x values');ylabel('y values');title('sine function');grid on;四、数据处理1、数据导入：在MATLAB中，可以通过各种方式将数据导入，如：a=load('filename.txt'); %从文件中载入数据b=xlsread('filename.xls'); %从Excel文件中载入数据五、编程1、条件语句：通过条件语句实现程序的分支结构，例如：if(a<0)disp('a is negative');elseif(a==0)disp('a is zero');elsedisp('a is positive');endfor i=1:10disp(i);end3、函数：在MATLAB中，可以自定义函数，函数调用格式为：function [out1,out2,...]=function_name(in1,in2,...)%函数说明%计算过程end4、脚本：在MATLAB中，脚本是一些命令或函数的集合，可以将脚本保存到文件中执行，例如：%脚本说明a=1;b=2;c=a+b;disp(c);以上便是MATLAB一些常用指令的详细介绍。

(一)P20.T1(1）指令行：[1 2;3 4]+10-2i(3) 指令行：[1 2;3 4].\[20 10;9 2] (5) 指令行：exp([1 2;3 4])(7) 指令行：prod([1 2;3 4])(9) 指令行：abs([1 2;3 4]-pi)(11) 指令行：find([10 20;30 40]>=[40,30;20,10]) (13) 指令行：all([1 2;3 4]>1)(15) 指令行：linspace(3,4,5)P20.T2(1) 指令行：clear;a=1,b=num2str(a),c=a>0,a==b,a==c,b==c(2) 指令行：clear;fun='abs(x)',x=-2,eval(fun),double(fun)P20.T3指令行：T='log(r)/(n*log(1+0.01*p))';r=2;p=0.5;n=12;eval(T)(二)P40.T2指令行：clear;s=0;n=0;s=s+log(1+n); while s<=100n=n+1;s=s+log(1+n); m=n;endm结果：m =37P40.T3指令行：clear; F(1)=1;F(2)=1;k=2;x=0; e=1e-8; a=(1+sqrt(5))/2;while abs(x-a)>ek=k+1;F(k)=F(k-1)+F(k-2);x=F(k)/F(k-1);enda,x,kP40.T5指令行：clear;x=0:24;y=[15 14 14 14 14 1516 18 20 22 23 2528 31 32 31 29 2725 24 22 20 18 17 16];plot(x,y);P40.T6(1)指令行：x=-2:0.1:2y=x.^2.*sin(x.^2.-x-2);plot(x,y,':ro');hold on;fplot(inline('x^2*sin(x^2-x-2)'),[-2,2 ]);hold off;(2)指令行：xa=-3:3;ya=-3:3;[x,y]=meshgrid(xa,ya); z=x.^2+y.^2;mesh(x,y,z)(5) 指令行：clear;close;t=0:0.1:2*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z);(7)指令行：x=0:0.1:2*pi;subplot(1,3,1);y1=sin(x);plot(x,y1);subplot(1,3,2);y2=sin(x).*sin(10*x);plot(x,y2);subplot(1,3,3);y3=-sin(x);plot(x,y3);P41.T7指令行：x=-2:0.05:2;y=1.1*(x>1.1)+x.*(x<=1.1).*(x>=-1.1)-1 .1*(x<-1.1);plot(x,y)（三）P59 T2指令行： A=[4 1 -1;3 2 -6;1 -5 3];B=[9;-2;1];A\B;[rank(a),rank([a,b])]由于方程有3个未知数，且)()(ArAr==3,所以方程解唯一。

matlab常用指令1. **基本数学运算**：+ `+`：加+ `-`：减+ `*`：乘+ `/`：除+ `^`：次方2. **矩阵和向量操作**：+ `eye(n)`：创建 n×n 单位矩阵+ `ones(n)`：创建 n×n 全 1 矩阵+ `zeros(n)`：创建 n×n 全 0 矩阵+ `A * B`：矩阵乘法+ `A\B`：矩阵求逆+ `size(A)`：获取矩阵 A 的大小+ `length(v)`：获取向量 v 的长度3. **绘图指令**：+ `plot(x, y)`：绘制 x 对 y 的曲线+ `subplot(m, n, p)`：创建 m×n 子图，并选择第 p 个子图 + `title('')`：设置图形标题+ `xlabel('')`：设置 x 轴标签+ `ylabel('')`：设置 y 轴标签+ `legend('')`：添加图例4. **数据分析**：+ `mean(x)`：计算向量 x 的平均值+ `std(x)`：计算向量 x 的标准差+ `sort(x)`：对向量 x 进行排序+ `max(x)`：找出向量 x 中的最大值+ `min(x)`：找出向量 x 中的最小值5. **文件操作**：+ `load('file.mat')`：加载 MAT 文件+ `save('file.mat', A)`：将矩阵 A 保存到 MAT 文件+ `fopen('file.txt', 'r')`：以只读方式打开文本文件+ `fscanf(fid, '%f', [n, m])`：从文件中读取数据并存储在矩阵中这只是 MATLAB 中一些常用指令的一小部分，MATLAB 提供了广泛的功能和指令，可以满足各种数学计算、数据分析和可视化需求。

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=1k=0;while b-a>epsx=(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)的导数。

数学实验答案Chapter 1Page20,ex1(5) 等于[exp(1),exp(2);exp(3),exp(4)](7) 3=1*3, 8=2*4(8) a为各列最小值，b为最小值所在的行号(10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture(11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10)(12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10)Page20, ex2(1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码Page20,ex3>> r=2;p=0.5;n=12;>> T=log(r)/n/log(1+0.01*p)Page20,ex4>> x=-2:0.05:2;f=x.^4-2.^x;>> [fmin,min_index]=min(f)最小值最小值点编址>> x(min_index)ans =0.6500 最小值点>> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点f1 =0.0328x1_index =24>> x(x1_index)ans =-0.8500>> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点f2 =0.0630x2_index =65>> x(x2_index)ans =1.2500>> z=magic(10)z =92 99 1 8 15 67 74 51 58 4098 80 7 14 16 73 55 57 64 414 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 3417 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6679 6 13 95 97 29 31 38 45 7210 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59>> sum(z)>> sum(diag(z))>> z(:,2)/sqrt(3)>> z(8,:)=z(8,:)+z(3,:)Chapter 2Page 45 ex1先在编辑器窗口写下列M函数，保存为eg2_1.m function [xbar,s]=ex2_1(x)n=length(x);xbar=sum(x)/n;s=sqrt((sum(x.^2)-n*xbar^2)/(n-1));例如>>x=[81 70 65 51 76 66 90 87 61 77];>>[xbar,s]=ex2_1(x)Page 45 ex2s=log(1);n=0;while s<=100n=n+1;s=s+log(1+n);endm=nPage 40 ex3clear;F(1)=1;F(2)=1;k=2;x=0;e=1e-8; a=(1+sqrt(5))/2;while abs(x-a)>ek=k+1;F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1);enda,x,k计算至k=21可满足精度clear;tic;s=0;for i=1:1000000s=s+sqrt(3)/2^i;ends,toctic;s=0;i=1;while i<=1000000s=s+sqrt(3)/2^i;i=i+1;ends,toctic;s=0;i=1:1000000;s=sqrt(3)*sum(1./2.^i);s,tocPage 45 ex5t=0:24;c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ...31 32 31 29 27 25 24 22 20 18 17 16];plot(t,c)Page 45 ex6(1)x=-2:0.1:2;y=x.^2.*sin(x.^2-x-2);plot(x,y)y=inline('x^2*sin(x^2-x-2)');fplot(y,[-2 2]) (2)参数方法t=linspace(0,2*pi,100);x=2*cos(t);y=3*sin(t); plot(x,y)(3)x=-3:0.1:3;y=x;[x,y]=meshgrid(x,y);z=x.^2+y.^2;surf(x,y,z)(4)x=-3:0.1:3;y=-3:0.1:13;[x,y]=meshgrid(x,y);z=x.^4+3*x.^2+y.^2-2*x-2*y-2*x.^2.*y+6;surf(x,y,z)(5)t=0:0.01:2*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z)(6)theta=linspace(0,2*pi,50);fai=linspace(0,pi/2,20); [theta,fai]=meshgrid(theta,fai);x=2*sin(fai).*cos(theta);y=2*sin(fai).*sin(theta);z=2*cos(fai);surf(x,y,z)(7)x=linspace(0,pi,100);y1=sin(x);y2=sin(x).*sin(10*x);y3=-sin(x);plot(x,y1,x,y2,x,y3)page45, ex7x=-1.5:0.05:1.5;y=1.1*(x>1.1)+x.*(x<=1.1).*(x>=-1.1)-1.1*(x<-1.1);plot(x,y)page45,ex9clear;close;x=-2:0.1:2;y=x;[x,y]=meshgrid(x,y);a=0.5457;b=0.7575;p=a*exp(-0.75*y.^2-3.75*x.^2-1.5*x).*(x+y>1);p=p+b*exp(-y.^2-6*x.^2).*(x+y>-1).*(x+y<=1);p=p+a*exp(-0.75*y.^2-3.75*x.^2+1.5*x).*(x+y<=-1);mesh(x,y,p)page45, ex10lookfor lyapunovhelp lyap>> A=[1 2 3;4 5 6;7 8 0];C=[2 -5 -22;-5 -24 -56;-22 -56 -16];>> X=lyap(A,C)X =1.0000 -1.0000 -0.0000-1.0000 2.0000 1.0000-0.0000 1.0000 7.0000Chapter 3Page65 Ex1>> a=[1,2,3];b=[2,4,3];a./b,a.\b,a/b,a\bans =0.5000 0.5000 1.0000ans =2 2 1ans =0.6552 一元方程组x[2,4,3]=[1,2,3]的近似解ans =0 0 00 0 00.6667 1.3333 1.0000矩阵方程[1,2,3][x11,x12,x13;x21,x22,x23;x31,x32,x33]=[2,4,3]的特解Page65 Ex 2(1)>> A=[4 1 -1;3 2 -6;1 -5 3];b=[9;-2;1];>> rank(A), rank([A,b]) [A,b]为增广矩阵ans =3ans =3 可见方程组唯一解>> x=A\bx =2.38301.48942.0213(2)>> A=[4 -3 3;3 2 -6;1 -5 3];b=[-1;-2;1];>> rank(A), rank([A,b])ans =3ans =3 可见方程组唯一解>> x=A\bx =-0.4706-0.2941(3)>> A=[4 1;3 2;1 -5];b=[1;1;1];>> rank(A), rank([A,b])ans =2ans =3 可见方程组无解>> x=A\bx =0.3311-0.1219 最小二乘近似解(4)>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1 2 3]';%注意b的写法>> rank(a),rank([a,b])ans =3ans =3 rank(a)==rank([a,b])<4说明有无穷多解>> a\bans =110 一个特解Page65 Ex3>> a=[2,1,-1,1;1,2,1,-1;1,1,2,1];b=[1,2,3]';>> x=null(a),x0=a\bx =-0.62550.6255-0.20850.4170x0 =11通解kx+x0Page65 Ex 4>> x0=[0.2 0.8]';a=[0.99 0.05;0.01 0.95];>> x1=a*x, x2=a^2*x, x10=a^10*x>> x=x0;for i=1:1000,x=a*x;end,xx =0.83330.1667>> x0=[0.8 0.2]';>> x=x0;for i=1:1000,x=a*x;end,xx =0.83330.1667>> [v,e]=eig(a)v =0.9806 -0.70710.1961 0.7071e =1.0000 00 0.9400>> v(:,1)./xans =1.17671.1767 成比例，说明x是最大特征值对应的特征向量Page65 Ex5用到公式(3.11)(3.12)>> B=[6,2,1;2.25,1,0.2;3,0.2,1.8];x=[25 5 20]'; >> C=B/diag(x)C =0.2400 0.4000 0.05000.0900 0.2000 0.01000.1200 0.0400 0.0900>> A=eye(3,3)-CA =0.7600 -0.4000 -0.0500-0.0900 0.8000 -0.0100-0.1200 -0.0400 0.9100>> D=[17 17 17]';x=A\Dx =37.569625.786224.7690Page65 Ex 6(1)>> a=[4 1 -1;3 2 -6;1 -5 3];det(a),inv(a),[v,d]=eig(a) ans =-94ans =0.2553 -0.0213 0.04260.1596 -0.1383 -0.22340.1809 -0.2234 -0.0532v =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170d =-3.0527 0 00 3.6760 00 0 8.3766(2)>> a=[1 1 -1;0 2 -1;-1 2 0];det(a),inv(a),[v,d]=eig(a) ans =1ans =2.0000 -2.0000 1.00001.0000 -1.0000 1.00002.0000 -3.0000 2.0000v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i-0.5773 0.5774 0.5774-0.5774 0.5773 - 0.0000i 0.5773 + 0.0000id =1.0000 0 00 1.0000 + 0.0000i 00 0 1.0000 - 0.0000i(3)>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]A =5 76 57 10 8 76 8 10 95 7 9 10>> det(A),inv(A), [v,d]=eig(A)ans =1ans =68.0000 -41.0000 -17.0000 10.0000-41.0000 25.0000 10.0000 -6.0000-17.0000 10.0000 5.0000 -3.000010.0000 -6.0000 -3.0000 2.0000v =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209d =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887(4)(以n=5为例)方法一（三个for）n=5;for i=1:n, a(i,i)=5;endfor i=1:(n-1),a(i,i+1)=6;endfor i=1:(n-1),a(i+1,i)=1;enda方法二（一个for）n=5;a=zeros(n,n);a(1,1:2)=[5 6];for i=2:(n-1),a(i,[i-1,i,i+1])=[1 5 6];enda(n,[n-1 n])=[1 5];a方法三（不用for）n=5;a=diag(5*ones(n,1));b=diag(6*ones(n-1,1));c=diag(ones(n-1,1));a=a+[zeros(n-1,1),b;zeros(1,n)]+[zeros(1,n);c,zeros(n-1,1)] 下列计算>> det(a)ans =665>> inv(a)ans =0.3173 -0.5865 1.0286 -1.6241 1.9489-0.0977 0.4887 -0.8571 1.3534 -1.62410.0286 -0.1429 0.5429 -0.8571 1.0286-0.0075 0.0376 -0.1429 0.4887 -0.58650.0015 -0.0075 0.0286 -0.0977 0.3173>> [v,d]=eig(a)v =-0.7843 -0.7843 -0.9237 0.9860 -0.92370.5546 -0.5546 -0.3771 -0.0000 0.3771-0.2614 -0.2614 0.0000 -0.1643 0.00000.0924 -0.0924 0.0628 -0.0000 -0.0628-0.0218 -0.0218 0.0257 0.0274 0.0257d =0.7574 0 0 0 00 9.2426 0 0 00 0 7.4495 0 00 0 0 5.0000 00 0 0 0 2.5505Page65 Ex 7(1)>> a=[4 1 -1;3 2 -6;1 -5 3];[v,d]=eig(a)v =0.0185 -0.9009 -0.3066-0.7693 -0.1240 -0.7248-0.6386 -0.4158 0.6170d =-3.0527 0 00 3.6760 00 0 8.3766>> det(v)ans =-0.9255 %v行列式正常, 特征向量线性相关，可对角化>> inv(v)*a*v 验算ans =-3.0527 0.0000 -0.00000.0000 3.6760 -0.0000-0.0000 -0.0000 8.3766>> [v2,d2]=jordan(a) 也可用jordanv2 =0.0798 0.0076 0.91270.1886 -0.3141 0.1256-0.1605 -0.2607 0.4213 特征向量不同d2 =8.3766 0 00 -3.0527 - 0.0000i 00 0 3.6760 + 0.0000i>> v2\a*v2ans =8.3766 0 0.00000.0000 -3.0527 0.00000.0000 0.0000 3.6760>> v(:,1)./v2(:,2) 对应相同特征值的特征向量成比例ans =2.44912.44912.4491(2)>> a=[1 1 -1;0 2 -1;-1 2 0];[v,d]=eig(a)v =-0.5773 0.5774 + 0.0000i 0.5774 - 0.0000i-0.5773 0.5774 0.5774-0.5774 0.5773 - 0.0000i 0.5773 + 0.0000id =1.0000 0 00 1.0000 + 0.0000i 00 0 1.0000 - 0.0000i>> det(v)ans =-5.0566e-028 -5.1918e-017i v的行列式接近0, 特征向量线性相关，不可对角化>> [v,d]=jordan(a)v =1 0 11 0 01 -1 0d =1 1 00 1 10 0 1 jordan标准形不是对角的，所以不可对角化(3)>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]A =5 76 57 10 8 76 8 10 95 7 9 10>> [v,d]=eig(A)v =0.8304 0.0933 0.3963 0.3803-0.5016 -0.3017 0.6149 0.5286-0.2086 0.7603 -0.2716 0.55200.1237 -0.5676 -0.6254 0.5209d =0.0102 0 0 00 0.8431 0 00 0 3.8581 00 0 0 30.2887>> inv(v)*A*vans =0.0102 0.0000 -0.0000 0.00000.0000 0.8431 -0.0000 -0.0000-0.0000 0.0000 3.8581 -0.0000-0.0000 -0.0000 0 30.2887本题用jordan不行, 原因未知(4)参考6(4)和7(1)Page65 Exercise 8只有(3)对称, 且特征值全部大于零, 所以是正定矩阵. Page65 Exercise 9(1)>> a=[4 -3 1 3;2 -1 3 5;1 -1 -1 -1;3 -2 3 4;7 -6 -7 0]>> rank(a)ans =3>> rank(a(1:3,:))ans =2>> rank(a([1 2 4],:)) 1,2,4行为最大无关组ans =3>> b=a([1 2 4],:)';c=a([3 5],:)';>> b\c 线性表示的系数ans =0.5000 5.0000-0.5000 1.00000 -5.0000Page65 Exercise 10>> a=[1 -2 2;-2 -2 4;2 4 -2]>> [v,d]=eig(a)v =0.3333 0.9339 -0.12930.6667 -0.3304 -0.6681-0.6667 0.1365 -0.7327d =-7.0000 0 00 2.0000 00 0 2.0000>> v'*vans =1.0000 0.0000 0.00000.0000 1.0000 00.0000 0 1.0000 v确实是正交矩阵Page65 Exercise 11设经过6个电阻的电流分别为i1, ..., i6. 列方程组如下20-2i1=a; 5-3i2=c; a-3i3=c; a-4i4=b; c-5i5=b; b-3i6=0;i1=i3+i4;i5=i2+i3;i6=i4+i5;计算如下>> A=[1 0 0 2 0 0 0 0 0;0 0 1 0 3 0 0 0 0;1 0 -1 0 0 -3 0 0 0; 1 -1 0 0 0 0 -4 0 0;0 -1 1 0 0 0 0 -5 0;0 1 0 0 0 0 0 0 -3; 0 0 0 1 0 -1 -1 0 0;0 0 0 0 -1 -1 0 1 0;0 0 0 0 0 0 -1 -1 1];>>b=[20 5 0 0 0 0 0 0 0]'; A\bans =13.34536.44018.54203.3274-1.18071.60111.72630.42042.1467Page65 Exercise 12>> A=[1 2 3;4 5 6;7 8 0];>> left=sum(eig(A)), right=sum(trace(A))left =6.0000right =6>> left=prod(eig(A)), right=det(A) 原题有错, (-1)^n应删去left =27.0000right =27>> fA=(A-p(1)*eye(3,3))*(A-p(2)*eye(3,3))*(A-p(3)*eye(3,3)) fA =1.0e-012 *0.0853 0.1421 0.02840.1421 0.1421 0-0.0568 -0.1137 0.1705>> norm(fA) f(A)范数接近0ans =2.9536e-013Chapter 4Page84 Exercise 1(1)roots([1 1 1])(2)roots([3 0 -4 0 2 -1])(3)p=zeros(1,24);p([1 17 18 22])=[5 -6 8 -5];roots(p)(4)p1=[2 3];p2=conv(p1, p1);p3=conv(p1, p2);p3(end)=p3(end)-4; %原p3最后一个分量-4roots(p3)Page84 Exercise 2fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x');fzero(fun,2)Page84 Exercise 3fun=inline('x^4-2^x');fplot(fun,[-2 2]);grid on;fzero(fun,-1),fzero(fun,1),fminbnd(fun,0.5,1.5)Page84 Exercise 4fun=inline('x*sin(1/x)','x');fplot(fun, [-0.1 0.1]);x=zeros(1,10);for i=1:10, x(i)=fzero(fun,(i-0.5)*0.01);end;x=[x,-x]Page84 Exercise 5fun=inline('[9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);16*x(1)-x(1)^3-2*x(2)^ 2-16*x(3)^2]','x');[a,b,c]=fsolve(fun,[0 0 0])Page84 Exercise 6fun=@(x)[x(1)-0.7*sin(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))];[a,b,c]=fsolve(fun,[0.5 0.5])Page84 Exercise 7clear; close; t=0:pi/100:2*pi;x1=2+sqrt(5)*cos(t); y1=3-2*x1+sqrt(5)*sin(t);x2=3+sqrt(2)*cos(t); y2=6*sin(t);plot(x1,y1,x2,y2); grid on; 作图发现4个解的大致位置，然后分别求解y1=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.5,2])y2=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.8,-2])y3=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[3.5,-5])y4=fsolve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[4,-4])Page84 Exercise 8(1)clear;fun=inline('x.^2.*sin(x.^2-x-2)');fplot(fun,[-2 2]);grid on; 作图观察x(1)=-2;x(3)=fminbnd(fun,-1,-0.5);x(5)=fminbnd(fun,1,2);fun2=inline('-x.^2.*sin(x.^2-x-2)');x(2)=fminbnd(fun2,-2,-1);x(4)=fminbnd(fun2,-0.5,0.5);x(6)=2feval(fun,x)答案: 以上x(1)(3)(5)是局部极小，x(2)(4)(6)是局部极大，从最后一句知道x(1)全局最小，x(2)最大。

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命令一、常用对象操作：除了一般windows窗口的常用功能键外。

1、!dir 可以查看当前工作目录的文件。

!dir& 可以在dos状态下查看。

2、who 可以查看当前工作空间变量名， whos 可以查看变量名细节。

3、功能键：功能键快捷键说明方向上键 Ctrl+P 返回前一行输入方向下键 Ctrl+N 返回下一行输入方向左键 Ctrl+B 光标向后移一个字符方向右键 Ctrl+F 光标向前移一个字符Ctrl+方向右键 Ctrl+R 光标向右移一个字符Ctrl+方向左键 Ctrl+L 光标向左移一个字符home Ctrl+A 光标移到行首End Ctrl+E 光标移到行尾Esc Ctrl+U 清除一行Del Ctrl+D 清除光标所在的字符Backspace Ctrl+H 删除光标前一个字符Ctrl+K 删除到行尾Ctrl+C 中断正在执行的命令4、clc可以命令窗口显示的内容，但并不清除工作空间。

二、函数及运算1、运算符：＋：加，－：减， *：乘， /：除， \：左除 ^：幂，‘：复数的共轭转置，（）：制定运算顺序。

2、常用函数表：sin( ) 正弦（变量为弧度）Cot( ) 余切（变量为弧度）sind( ) 正弦（变量为度数）Cotd( ) 余切（变量为度数）asin( ) 反正弦（返回弧度）acot( ) 反余切（返回弧度）Asind( ) 反正弦（返回度数）acotd( ) 反余切（返回度数）cos( ) 余弦（变量为弧度）exp( ) 指数cosd( ) 余弦（变量为度数）log( ) 对数acos( ) 余正弦（返回弧度）log10( ) 以10为底对数acosd( ) 余正弦（返回度数）sqrt( ) 开方tan( ) 正切（变量为弧度）realsqrt( ) 返回非负根tand( ) 正切（变量为度数）abs( ) 取绝对值atan( ) 反正切（返回弧度）angle( ) 返回复数的相位角atand( ) 反正切（返回度数）mod(x,y) 返回x/y的余数sum( ) 向量元素求和3、其余函数可以用help elfun和help specfun命令获得。

MATLAB指令整理1 常用指令(General Purpose Commands)bar3 三維直方圖4.3 特殊三維圖形（Specialized 3－D graphs）bar3h 三維水準直方圖comet3 三維彗星動態軌跡線圖ezgraph3 通用指令voronoi Voronoi 圖ezcontourf 畫填色等位線contourf 等高線圖contour3 三維等高線ezcontour 畫等位線ezmesh 畫網線圖pcolor 用顏色反映資料的偽色圖ezplot3 畫三維曲線scatter3 三維散點圖ribbon 以三維形式繪製二維曲線stem3 三維離散杆圖surfc 帶等高線的三維表面圖trimesh 三角剖分網線圖pie3 立體圓形圖meshz 帶零基準面的三維網線圖contour 等高線圖ezsurf 畫曲面圖ezsurfc 畫帶等位線的曲面圖meshc 帶等高線的三維網線圖ezmeshc 畫等位線的網線圖4.2 等高線及二維半圖形(Contour and 2－1/2D graphs)comet 彗星狀軌跡圖barh 水準直方圖compass 從原點出發的複數向量圖errorbar 誤差棒棒圖ezplot 畫二維曲線bar 直方圖area 面域圖printopt 印表機設置verml 將圖形保存為VRML2.0 檔4 特殊圖形(Specialized graphs)4.1 特殊平面圖形(Specialized 2 D graphs)ezpolar 畫極座標曲線feather 從X 軸出發的複數向量圖scatter 散點圖stairs 階梯形曲線圖stem 火柴杆圖trisurf 三角剖分曲面圖plotmatrix 散點圖陣列pie 餅形統計圖fill 多邊填色圖fplot 函數曲線圖hist 統計頻數直方圖pareto Pareto 圖clabel 給等高線加標注4.4 內剖及向量視圖(Volume and vector visualization)5.1 圖形窗的產生和控制(Figure window creation and control)5 控制碼圖形(Handle Graphics)clf 清除當前圖close 關閉圖形打開或創建圖形視窗Surf2patch 將曲面資料轉換為塊資料sphere 球面4.7 顏色相關函數(Color related function)spinmap 顏色週期性變化操縱4.8 三維模型函數(Solid modeling)cylinder 圓柱面gcf 獲得當前圖的柄openfig 打開圖形caxis 控制色軸的刻度cla 清除當前軸gca 獲得當前軸的柄hold 圖形的保持box 座標形式在封閉式和開啟式之間切換axis 軸的控制refresh 刷新圖形shg 顯示圖形窗5.2 軸的產生和控制(Axis creation and control)axesrotate 旋轉指令moviein 影片動畫記憶體初始化brighten 控制色彩的明暗4.5 圖像顯示及檔處理(Image display and file I/O)colorbar 色彩條狀圖colormap 設置色圖contrast 提高圖像對比度的灰色圖slice 切片圖quiver3 三維方向箭頭圖print 列印圖形或把圖存入檔coneplot 錐體圖contourslice 切片等位線圖quiver 向量場圖gray 線性灰度顯示圖像frame2im 將影片動畫轉換為編址圖像getframe 獲得影片動畫圖像的幀im2frame 將編址圖像轉換為影片動畫movie 播放影片動畫capture 當前圖的屏捕捉4.6 影片和動畫(Movies and animation)imagesc 顯示亮度圖像imfinfo 獲取圖像檔的特徵資料imread 從檔讀取圖像的資料陣（和伴隨色圖））imwrite 把強度圖像或真彩圖像寫入檔waterfall 瀑布水線圖3.8 硬拷貝(Hardcopy and printing)legend 圖例說明座標形式在封閉式和開啟詞式之間切換texlabel 將字串轉換為Tex 格式3 三維圖形函數(Three dimensional graphs)3.1 基本三維圖形(Elementary 3－D plots)軸控制(Axis control)semilogy Y 軸半對數刻度曲線plot 直角坐標下線性刻度曲線plotyy 雙縱坐標圖polar 極座標曲線圖semilogx Xfill3 三維曲面多邊形填色mesh 三維網線圖graymon 設置缺省圖形視窗為單色顯示幕hidden 消隱shading 圖形渲染模式whitebg 設置圖形視窗為白底colordef 用色風格caxis (偽)顏色軸刻度plot3 三維直角坐標曲線圖surf 三維表面圖3.2 色彩控制(Color control)alpha 透明色控制loglog 雙對數刻度曲線基本平面圖形(Elementary X－Y graphs)繪製曲線繪製含等高線的網格圖繪製3 維曲線採用極座標繪圖畫函數曲線圖繪製網格圖1.4 繪圖(Plotting)1.3 數值積分(Numerical intergration)dblquad 二重(閉型)數值積分指令quad 低階法數值積分quadl 高階法數值積分1.5 內聯函數物件(Inline function object)argnames 給出函數的輸入宗量ode23t 解適度剛性微分方程odet23tb 低階法解剛性微分方程ode45 高階法解微分方程二維圖形函數(Two dimensional graphs)ode23s 低階法解剛性微分方程ode23 低階法解微分方程formula 函數公式1.6 差微分函數解算器(Differential equation solvers)ode113 變階法解方程ode15s 變階法解剛性方程3.3 光照模式(Lighting)diffuse 漫反射表面係數zlim Z 軸範圍ylim Yzoom 二維圖形的變焦放大3.6 視角控制(Viewpoint control)rotate3d 旋動三維圖形xlim Xpbaspect 畫座標框的PlotBoxAspectRatio 屬性創建軸軸的刻度和表現daspect 軸的DataAspectRatio 屬性grid 畫座標格線view 設定3圖形觀測點viewmtx 觀測點轉換矩陣xlabel Xylabel Yzlabel Z 軸名標注ishold 若圖形處保持狀態則為真title 圖形標題在圖上標注文字3.7 圖形注釋(Graph annotation)顯示色條gtext 用滑鼠在圖上標注文字plotedit 圖形編輯工具3.5 軸的控制(Axis control)winter 藍、綠濃淡色autumn 紅、黃濃淡色3.4 色圖(Color maps)bone 藍色調灰度圖colorcube 三濃淡多彩交錯色cool 青和品紅濃淡色圖surfl 帶光照的三維表面圖surfnorm 表面圖的法線燈光控制lighting 設置照明模式material 使用預定義反射模式specular 漫反射copper 線性變化純銅色調圖flag 紅－白－藍黑交錯色圖spring 青、黃濃淡色summer 綠、黃濃淡色vga 16 色white 全白色prism 光譜色圖pink 淡粉紅色圖hot 黑－紅－黃－白交錯色圖hsv 飽和色彩圖jet 變異HSV 色圖lines 採用plot 繪線色orient 設置送紙方向subplot 創建子圖0.4 構架函數(Structure functions)struct2 cell 把構架陣列轉換為元胞陣列fieldnames 獲取構架的功能變數名稱getfield 獲取域的內容isfield 若為給定構架的功能變數名稱則為真num2 cell 把數值陣列轉換為元胞陣列is cell 若是元胞則為真cellfun 元胞陣列函數cellplot 圖示元胞陣列的內容cell2struct 把元胞陣列轉換為構架陣列deal 把輸入分配給輸出isstruct 若是構架則為真rmfield 刪除構架的域0.6 面向物件編程(Object oriented programming functions)str2func 將字串轉換為函數控制碼dlass 查明變數的類型isa 若是指定的資料類型則為真inferiorto 級別較低func2str 將函數控制碼陣列轉換為字串functions 列舉函數控制碼對應的函數setfield 指定構架域的內容struct 創建構架變數0.5 函數控制碼函數(Function handle functions)@ 創建函數控制碼celldisp 顯示元胞陣列內容cell 創建元胞變數JavaArray 構建Java 陣列inline 創建內聯函數JavaMethod 調用某個Java 方法JavaObject 調用Java 物件的構造函數single 轉變為單精確度數值function handle 函數控制碼double 轉化為16 位相對精度的浮點數值物件0 資料類型(Data types and structures)0.1 資料類型(Data types)char 創建字元傳輸組或者將其他類型變數轉化為字串陣列sparse 創建稀疏矩陣uint8(unit16、unit32) 轉換為8(16、32)位元無符號整型數permute 廣義非共軛轉置shiftdim 維數轉換squeeze 使陣列降維0.3 元胞陣列函數(Cell array functions)ipermute 廣義反轉置ndgrid 為N－D 函數和插補創建陣列int8(nit16、nit32) 轉換為8(16、32)位元符號整型數0.2 多維陣列函數(Multi－dimensional array functions)cat 把若干陣列串接成高維陣列ndims 陣列A 的維數isjava 若是java 對象則為真isobject 若是對象則為真expand 對指定項展開collect 合併同類項factor 進行因式或因數分解horner 轉換成嵌套形式numden 提取公因式2.3 化簡(Simplification)triu 抽取上三角陣rank 秩rref 轉換為行階梯形svd 奇異值分解tril 抽取下三角陣simple 運用各種指令化簡符號運算式simplify 恒等式簡化fminunc 擬牛頓法求多元函數極值點fsolve 解非線性方程組lsqnonlin 解非線性最小二乘問題solve 求解方程組finverse 求反函數dsolve 求解符號常微分方程subexpr 運用符號變數置換子運算式subs 通用置換指令22.4 方程求解(Solution of Equation)compose 求復函數poly 特徵多項式null 零空間peaks 產生peaks 圖形資料intro 幻燈演示指令2 符號工具包(Symbolic Math Toolbox)2.1 微積分(Calculus)diff 求導數flow 無限大水體中水下射流速度資料demo 演示程式methods 顯示類的方法名substruct 創建構架總量superiorto 級別較高1 示例(E xamples and demonstrations)limit 求極限int 計算積分eig 矩陣特徵值和特徵向量expm 矩陣指數inv 矩陣的逆jordan Jordan 分解diag 創建對角陣，抽取對角向量det 行列式的值jacobian Jacobian 矩陣symsum 符號序列的求和trylor Trylor 級數2.2 線性代數(Linear Algebra)toc 碼錶終止和顯示tic 碼錶啟動6 圖形用戶介面工具(Graphical user interface tools)newplot 下一個新圖align 對齊用戶控制項和軸cbedit 編輯回調函數ginput 從滑鼠得到圖形點座標ishandle 若是圖柄代號側為真closereq 關閉圖形窗請求函數rmappdata 刪除應用程式定義資料set 建立物件特性setappdata 建立應用程式定義資料5.5 工具函數(Utilities)guide 設計GUImenu 創建菜單blanks 空白字元號cellstr 通過字串陣列構建字串的元胞陣列deblank 刪除最後的空格把字串變成ASCII 碼值7.1 通用字串函數(General)7 字串(Character string)menuedit 功能表編輯propedit 屬性編輯uicontrol 創建用戶介面控制uimenu 創建用戶介面功能表reset 重設物件特性isappdata 檢驗是否應用程式定義資料rectangle 創建方patch 創建塊surface 創建面text 創建圖形中文本uicontextmenu 創建現場功能表物件line 創建線light 創建光5.3 控制碼圖形物件(Handle Graphics objects)axex 在任意位置創建軸figure 創建圖形視窗image 創建圖像用戶使用介面控制用戶使用功能表控制gcbo “正執行回調操作”的控制項圖柄指令gco 獲得當前對象的柄get 獲得物件特性getappdat 獲得應用程式定義資料gcbf “正執行回調操作”的圖形的柄findobj 用規定的特性找尋對象5.4 控制碼圖形處理(Handle Graphics operations)copyobj 拷貝圖形物件及其子物件delete 刪除物件及檔drawnow 螢幕刷新eval 執行串形式的MATLAB 運算式7.2 字串查詢(String tests)8 檔輸入/輸出(File input/output)upper 把字串變成大寫clc 清除指令視窗disp 顯示矩陣和文字內容fprintf 把格式化資料寫到檔或螢幕str2num 將字串轉換為數值str2mat 將字串轉換為含有空格的陣列strncmpi 比較字串的前N 個字元（忽略大小寫）strrep 用另一個串代替一個串中的子串strtok 刪除串中的指定子串strvcat 創建字串陣列home 游標返回行首input 提示鍵盤輸入date 日期etime 用CLOCK 計算的時間now 當前時鐘和日期pause 暫停cputme MATLAB 佔用CPU 時間clock 時鐘load 從磁片中調入資料變數sprintf 寫格式資料到串sscanf 在格式控制下讀串19 時間和日期(Time and dates)strnomp 比較字串的前N 個字元strmatch 逐行搜索串bin2dec 二進位串轉換為十進位整數base2dec X－進制串轉換為十進位整數dec2base 十進位整數轉換為X 進制串dec2bin 十進位整數轉換為二進位串dec2hex 十進位整數轉換為16 進制串7.3 字串操作(String operations)isstr 若是字串則為真iscellstr 若是字串組成的元胞陣列則為真ischar 若是字串則為真isletter 串中是字母則為真isspace 串中是空格則為真findstr 在一個串中尋找一個子串hex2dec 16－進制串轉換為十進位整數strcmp 比較字串strcmpi 比較字串（忽略大小寫）strings MATLAB 中的字串strjust 字串的對齊方式strcat 把多個串連接成長串num2str 把數值轉換為字串hex2num 16－進制串轉換為浮點數int2str 將整數轉換為字串lower 把字串變成小寫mat2str 將陣列轉換為字串optimset 創建或修改OPTIONS 構架1.2 優化選項處理(Optimization Option handling)global 定義總體變數function 函數檔頭isglobal 若是總體變數則為真iskeyword 若是關鍵字則為真mfilename 正在執行的M 檔的名字exist 檢查變數或函數是否被定義腳本檔、函數及變數(Scripts,function,and variables)evalc 執行MATLAB 字串evalin 跨空間計算串運算式的值feval 函數巨集指令run 執行腳本檔persistent 定義永久變數script MATLAB 命令文件varagout 輸出宗量varagin 輸入宗量資訊顯示(Message display)display 顯示矩陣和文字內容的重載函數error 顯示錯誤資訊nargoutchk 輸出變數個數檢查nargout 函數輸出宗量的個數宗量處理(Augument handling)inputname 實際調用變數名nargchk 輸入變數個數檢查nargin 函數輸入宗量的個數字串巨集指令builtin 執行內建的函數break 終止最內迴圈控制語句(Control flow)case 同switch 一起使用catch try 一起使用continue 將控制轉交給外層的for 或while 迴圈編程語言結構(Programming language constructs)subsindex 下標標識[,] 水準串接[;] 垂直串接下標賦值( ),{ },. 下標標識else 同if 一起使用elseif 同if 一起使用try try-cathch 結構while 不確定次數重複執行語句計算運行(Evaluation and execution)assignin 跨空間賦值switch 多個條件分支return 返回結束for,while,if 語句for 按規定次數重複執行語句if 條件執行語句otherwise 可同switch 一起使用lasterr 最後一個錯誤資訊lastwarn 最後一個警告資訊reshape 矩陣變維ind2sub 據單下標換算出全下標rot90 矩陣逆時針90°旋轉sub2idn 據全下標換算出單下標特殊變數和常數(Special variables and constants)flipdim 交換對稱位置上的元素flipud 矩陣的上下翻轉blkdiag 塊對角陣串接end 陣列的長度，即最大下標find 找出非零元素1 的下標fliplr 矩陣的左右翻轉ans 最新運算式的運算結果eps 浮點相對誤差pi 3.1415926535897?.realmax 最大浮點數realmin 最小正浮點數why 一般問題的簡明答案NaN 或nan 非數isnan 若為非數則為真i,j 虛數單位inf 或Inf 無窮大isfinite 若是有限數則為真isinf 若是無窮大則為真矩陣操作(Matrix manipulateion)size 確定矩陣的維數linspace 線性等分向量eye 單位陣logspace 對數等分向量meshgrid 用於三維曲面的分格線座標ones 全1 矩陣基本矩陣(Elementary matrices)基本矩陣函數和操作(Elementary matrices and matrix manipulation) 按格式把數位轉換為串warning 顯示警告資訊3.6 互動式輸入(Interactive input)keyboard 啟動鍵盤做為命令檔rand 均勻分佈隨機陣randn 正態分佈隨機陣islogical 尤其是邏輯數則為真isnumeric 若是數值則為真length 確定向量的長度logical 將數值轉化為邏輯值isequal 若對應元素相等則為1isempty 若是空矩陣則為真repmat 鋪放模組陣列zeros 全零矩陣: 矩陣的援引和重排矩陣基本資訊(Basic array information).ˊ轉置號ˊ複數轉置號path 控制MATLAB 的搜索路徑rmpath 從搜索路徑中刪除目錄pathtool 修改搜索路徑指令窗控制(Controlling the command window)beep 產生beep 聲addpath 添加搜索路徑搜索路徑的管理(Managing the seach patli)pcode 生成P 碼檔type 顯示檔內容what 列出當前目錄上的M、MAT、MEX 文件which 確定指定函數和檔的位置echo 顯示命令檔指令的切換開關diary 儲存MATLAB 指令窗操作內容刪除檔dir 列出的文件dos 執行dos 指令並返還結果getenv 給出環境值copyfile 檔拷貝computer 電腦類型format 設置資料輸出格式more 命令視窗分頁輸出的控制開關作業系統指令(Operating system commands)cd 改變當前工作目錄open 打開文件mex 創建MEX 文件ver MATLAB 和TOOLBOX 的版本資訊subscribe MATLAB 用戶註冊version MATLAB 版本whatsnew 顯示版本新特徵1.2 工作空間管理(Managing the workspace)info MATLAB 和MathWorks 公司的資訊helpwin 打開線上幫助窗通用資訊查詢(General information)help 線上幫助指令helpbrowser 超文本文檔幫助資訊helpdesk 超文本文檔幫助資訊clear 從記憶體中清除變數和函數exit 關閉MATLAB管理指令和函數(Managing commands and functions)矩陣編輯器edit 打開M 文件inmem 查看記憶體中的P 碼檔workspace 工作記憶體流覽器whos 列出工作記憶體中的變數細節pack 合併工作記憶體中的碎塊quit 退出save 把記憶體變數存入磁片who 列出工作記憶體中的變數名ispc MATLAB 為PC(Windows)版本則為真isunix MATLAB 為Unix 版本則為真all 所有元素均非零則為真any 有非零元則為真特殊算符(Special characters)：冒號( ) 圓括號xor 異或~ 邏輯非>= 大於或等於邏輯操作(Logical operators)& 邏輯與| 邏輯或[ ] 方括號{ } 花括弧% 注釋號! 調用作業系統命令= 賦值符號引號; 分號, 逗號小數點. 構架域的關節點.. 父目錄? 續行號<= 小於或等於> 大於2.1 算術運算符(Arithmetic operators)2 運算符和特殊算符(Operators and special characters)+ 加- 減* 矩陣乘執行外部應用程式web 打開web 流覽器mkdir 創建目錄pwd 改變當前工作目錄unix unix 指令並返還結果vms 執行vms dcl 指令並返還結果.* 陣列乘^ 矩陣乘方關係運算符(Relational operators)= = 等號~= 不等號< 小於[注]本表第三欄括弧中的字元供線上救助時help 指令引述用張量積.^ 陣列乘方\ 反斜杠或左除/ 斜杠或右除./或.\ 陣列除4特殊矩陣(Specialized matrices)compan 伴隨矩陣相關(Correlation)gradient 梯度corrcoef 相關係數cov 協方差矩陣subspace 子空間之間的角度差分和近似微分del2 五點離散Laplaciansum 元素和trapz 梯形數值積分var 求方差有限差分(Finite differentces).4 濾波和卷積(Filtering and convoluteion)卷積和多項式相乘8.5 傅裏葉變換(Fourier transforms)fft 快速離散傅裏葉變換fft2 二維離散傅裏葉變換fftn N 維離散傅裏葉變換fliter2 二維數字濾波器filter 一維數字濾波器conv2 二維卷積convn N 維卷積detrend 去除線性分量解卷和多項式相除std 標準差sortrows 由小到大按行排序7.5 因式分解(Factorization utility)sqrtm 矩陣平方根cdf2rdf 複數對角型轉換到實塊對角型balance 改善特徵值精度的平衡刻度rsf2csf 實塊對角型轉換到複數對角型logm 矩陣對數funm 計算一般矩陣函數7矩陣函數(Matrix functions)expm1 矩陣指數的Pade 逼近expm2 用泰勒級數求矩陣指數expm3 通過特徵值和特徵向量求矩陣指數資料分析和傅裏葉變換(Date analysis and Fourier transforms)基本運算(Basic operations)median 中值min 最小值prod 元素積sort 由小到大排序mean 平均值max 最大值cumprod 元素累計積cumsum 元素累計和cumtrapz 累計積分histc 直方圖統計fftshift 重排fft 和fft2 的輸出ifft 離散傅裏葉反變換多項式(Polynomials)spline 三次樣條插補conv 多項式相乘deconv 多項式相除由根創建多項式ppval 計算分段多項式樣條插補(Spline Interpolation)interp2 二維插補interp3 三維插補intern N 維插補pchip hermite 插補polyder 多項式微分polyfit 多項式擬合1.1 優化和尋根(Optimization and root finding)fminbnd 非線性函數在某區間中極小值fminsearch 單純形法求多元函數極值點指令fzero 單變數函數的零點數值泛函函數和ODE 解算器(Function functions and ODE solvers) roots 求多項式的根polyint 積分多項式分析polyval 求多項式的值polyvalm 求矩陣多項式的值residue 求部分分式表達interp1q 快速一維插補interp1 一維插補waveplay 利用系統音頻輸出設配播放soundsc 自動標刻並播放waverecor 利用系統音頻輸入設配錄音.2 音頻檔輸入輸出(Audio file import and export)auread 讀取音頻檔(.au)sound 播放向量.1 音頻硬體驅動(Audio hardware drivers)ifft2 二維離散傅城葉反變換ifftn N 維離散傅裏葉反變換ifftshift 反fftshift音頻支援(Audio support)auwrite 創建音頻檔(.au)wavread 讀取音頻檔(.wav)griddata 分格點數據griddata3 三維分格點數據griddatan 多維分格點數據interpft 利用FFT 方法一維插補數據插補(Data Interpolation)插補多項式函數(Interpolation and polynomials)wavwrite 創建音頻檔(.wav)9.3 工具(Utilities)lin2mu 將線性信號轉換為μ 一律編碼的信號mu2lin 將μ 一律編碼信號轉換為線性信號svds 多個奇異值schur Schur 分解sec 正割csch 雙曲余割sech 雙曲正割sin 正弦sinh 雙曲正弦csc 余割coth 雙曲餘切atan2 四象限反正切cos 余弦cosh 雙曲余弦cot 餘切tan 正切tanh 雙曲正切pow2 2 的冪sqrt 平方根複數函數(Complex)abs 絕對值nestpow2 最近鄰的2 的冪log2 以2 為底的對數指數函數(Exponential)exp 指數log 自然對數log10 常用對數atanh 反雙曲正切atan 反正切rosser 典型對稱特徵值實驗問題pascal Pascal 矩陣toeplitz Toeplitz 矩陣vander Vandermonde 矩陣wilkinson Wilkinson’s 對稱特徵值實驗矩陣magic 魔方陣invhilb 逆Hilbert 矩陣gallery 一些小測試矩陣hadamard Hadamard 矩陣hankel Hankel 矩陣hilb Hilbert 矩陣基本數學函數(Elementary math functions)三角函數(Trigonometric)asec 反正割asech 反雙曲正割asin 反正弦asinh 反雙曲正弦acsch 反雙曲余割acsc 反余割acos 反余弦acosh 反雙曲余弦acot 反餘切acoth 反雙曲餘切angle 相角complex 將實部和虛部構成複數condest 估計1-範數條件數cond 矩陣條件數lu LU 分解luinc 不完全LU 分解lscov 已知協方差的最小二乘積cholinc 不完全Cholesky 分解chol Cholesky 分解normest 2 範數orth 值空間trace 跡子空間的角度線性方程(Linear equations)nnls 非負二乘解pinv 偽逆gsvd 歸一化奇異值分解hess Hessenberg 矩陣polyeig 多項式特徵值問題qz 廣義特徵值eigs 多個特徵值condeig 矩陣各特徵值的條件數qr QR 分解rcond LINPACK 逆條件數\、/ 解線性方程特性值與奇異值(Eigenvalues and singular values)norm 矩陣或向量範數矩陣分析(Matrix analysis)ceil 朝正無窮大方向取整5圓整和求餘函數(Rounding and remainder)fix 朝零方向取整floor 朝負無窮大方向取整mod 模數求餘unwrap 相位角360°線調整real 複數實部conj 複數共軛cplxpair 複數陣成共軛對形式排列imag 複數虛部isreal 若是實數矩陣則為真rem 求餘數round 四捨五入取整isprime 若是質數則為真pol2cart 柱(或極)座標變為直角坐標sph2cart 球座標變為直角坐標矩陣函數和數值線性代數(Matrix functions－numerical linear algebra) dot 向量內積cross 向量叉積sign 符號函數6特殊函數(Specialized math functions)cart2pol 直角坐標變為柱(或極)座標cart2sph 直角坐標變為球座標optimget 從OPTIONS 構架中取得優化參數。

MATLAB命令大全原文已完。

下文为附加文档，如不需要，下载后可以编辑删除，谢谢！施工组织设计本施工组织设计是本着“一流的质量、一流的工期、科学管理”来进行编制的。

编制时，我公司技术发展部、质检科以及项目部经过精心研究、合理组织、充分利用先进工艺，特制定本施工组织设计。

一、工程概况：西夏建材城生活区27#、30#住宅楼位于银川市新市区，橡胶厂对面。

本工程由宁夏燕宝房地产开发有限公司开发，银川市规划建筑设计院设计。

本工程耐火等级二级，屋面防水等级三级，地震防烈度为8度，设计使用年限50年。

本工程建筑面积:27#楼3824.75m2;30#楼3824.75 m2。

室内地坪±0.00以绝对标高1110.5 m为准，总长27#楼47.28m；30#楼47.28 m。

总宽27#楼14.26m；30#楼14.26 m。

设计室外地坪至檐口高度18.6 00m，呈长方形布置，东西向，三个单元。

本工程设计屋面为坡屋面防水采用防水涂料。

外墙水泥砂浆抹面，外刷浅灰色墙漆。

内墙面除卫生间200×300瓷砖，高到顶外，其余均水泥砂桨罩面，刮二遍腻子；楼梯间内墙采用50厚胶粉聚苯颗粒保温。

地面除卫生间200×200防滑地砖，楼梯间50厚细石砼1：1水泥砂浆压光外，其余均采用50厚豆石砼毛地面。

楼梯间单元门采用楼宇对讲门，卧室门、卫生间门采用木门，进户门采用保温防盗门。

本工程窗均采用塑钢单框双玻窗，开启窗均加纱扇。

本工程设计为节能型住宅，外墙均贴保温板。

本工程设计为砖混结构，共六层。

基础采用C30钢筋砼条形基础，上砌MU30毛石基础，砂浆采用M10水泥砂浆。

一、二、三、四层墙体采用M10混合砂浆砌筑MU15多孔砖；五层以上采用M7.5混合砂浆砌筑MU15多孔砖。

本工程结构中使用主要材料：钢材：I级钢，II级钢；砼：基础垫层C10，基础底板、地圈梁、基础构造柱均采用C30，其余均C20。

本工程设计给水管采用PPR塑料管，热熔连接；排水管采用UPVC硬聚氯乙烯管，粘接；给水管道安装除立管及安装IC卡水表的管段明设计外，其余均暗设。