胡良剑 matlab 数学实验常用指令整理
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一些常用指令的详细介绍。
MATLAB常用指令(全)(print)
persistent 定义永久变量 script MATLAB 命令文件 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 提示键盘输入
singular values) condeig 矩阵各特征值的条件数 eig 矩阵特征值和特征向量 eigs 多个特征值 gsvd 归一化奇异值分解 hess Hessenberg 矩阵 poly 特征多项式
ispc MATLAB 为 PC(Windows)版本则为 ( ),{ },. 下标赋值
helpbrowser 超文本文档帮助信息
真
( ),{ },. 下标标识
helpdesk 超文本文档帮助信息
isunix MATLAB 为 Unix 版本则为真
subsindex 下标标识
helpwin 打开在线帮助窗
cplxpair 复数阵成共轭对形式排列 imag 复数虚部 isreal 若是实数矩阵则为真 real 复数实部 unwrap 相位角 360°线调整 5.4 圆整和求余函数(Rounding and
胡良剑 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常用指令
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 提供了广泛的功能和指令,可以满足各种数学计算、数据分析和可视化需求。
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))矩阵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)的导数。
(完整版)MATLAB数学实验第二版答案(胡良剑)
数学实验答案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操作命令大全
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命令一、常用对象操作:除了一般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指令整理精品资料
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命令大全
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卡水表的管段明设计外,其余均暗设。
胡良剑 matlab 数学实验常用指令整理
第七章
1、符号对象的定义 P117 s=sym(str) 将数值或字符串 str 转化为符号对象 s, 数值为有理表示; s= sym(num,’d’) 将数值表达式转化为符号表达式,数值用十进制表示; syms var1 var2 定义 var1, var2, 为符号变量 subs(s, old, new) 将符号表达式 s 中的符变量 old 用 new 代替 2、计算精度和数据类型转换 P118 digits(n) 将数值计算精度设为 n 位; x=vpa(s) 求 s 的数值结果 x=vpa(s,n) 采用 n 位计算精度求 s 的数值结果 double(s) 将符号对象转化为双精度数值 char(s) 将符号对象转化为字符串 数据类型转换图指令详解 P118 3、符号函数计算 P120 factor(expr) 对 expr 作因式分解 expand(expr) 将 expr 展开 collect(expr,v) 将 expr 按变量 v 合并同类项 simple(expr) 将 expr 化简 g=finverse(f,v) 求函数 f(v)的反函数 g(v) fg=compose(f, g) 求函数 f(v)和 g(v)的复合函数 f(g(v)) [n,d]=numden(expr) 分式通分,n 返回分子,d 返回分母 funtool 函数分析图形界面 4、极限和级数及微分 P122 limit(s,x,a) 返回符号表达式 s 当 x->a 时的极限 limit(s,x,a,’right’) 返回 s 当 x->a 时的右极限 limit(s,x,a,’left’) 返回 s 当 x->a 时的左极限 symsum(s,n,a,b) 返回符号表达式 s 表示的通项当自变量 n 由 a 到 b 的和 diff(s,x) 返回符号表达式 s 对 x 的导函数 diff(s,x,n) 返回 s 对 x 的 n 阶导函数。 taylor(s,n,a,x) 返回符号表达式 s 在 a 点 aylor 展开到 n-1 次式.自变量为 x taylortool Taylor 分析图形界面 jacobian(f,x) 返回向量函数 f 的 Jacobian 矩阵 5、积分 P124 int(s) 符号表达式 s 的不定积分.
matlab常用函数与常用指令大全[整理版]
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数学实验第二版答案(胡良剑)精编版
数学实验答案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指令大全
分享我的分享当前分享返回分享首页»分享 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函数指令大全
包括各种函数、命令含义、还有simulink模块的,希望有帮助A aabs 绝对值, 模acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割all 所有元素均非零则为真alpha 透明控制angle 相角ans 最新表达式的运算结果any 有非零元则为真area 面域图asec 反正割asech 反双曲正割asin 反正弦asinh 反双曲正弦atan 反正切atan2 四象限反正切atanh 反双曲正切autumn 红、黄浓淡色axis 轴的刻度和表现B bbar 直方图binocdf 二项分布概率binopdf 二项分布累积概率binornd 产生二项分布随机数组blanks 空格符号bode 给出系统的对数频率曲线bone 蓝色调浓淡色阵box 坐标封闭开关break 终止最内循环brighten 控制色彩的明暗butter ButterWorth低通滤波器C ccaxis (伪)颜色轴刻度cd 设置当前工作目录cdf2rdf 复数对角型转换到实块对角型ceil 朝正无穷大方向取整cell 创建单元数组char 创建字符串数组或者将其他类型变量转化为字符串数组charfcn Maple函数Children 图形对象的子对象clabel 等高线标注class 判别数据类别clc 清除指令窗中显示内容clear 从内存中清除变量和函数clf 清除当前图形窗图形close 关闭图形窗collect 合并同类项Color 图形对象色彩属性colorbar 显示色条colorcube 三浓淡多彩交错色colordef 定义图形窗色彩colormap 设置色图comet 彗星状轨迹图comet3 三维彗星动态轨迹线图compass 射线图;主用于方向和速度cond 矩阵条件数conj 复数共轭continue 将控制转交给外层的for或while循环contour 等高线图contourf 填色等高线图conv 卷积和多项式相乘cool 青和品红浓淡色图copper 线性变化纯铜色调图corrcoef 相关系数cos 余弦cosh 双曲余弦cot 余切coth 双曲余切cov 协方差矩阵csc 余割csch 双曲余割cumsum 元素累计和cumtrapz 梯形法累计积分D ddblquad 二重(闭型)数值积分指令deconv 解卷和多项式相除del2 计算曲率demos 演示函数det 行列式的值diag 创建对角阵,抽取对角向量diff 求导数,差分和近似微分digits 控制符号数值的有效数字位数dir 列出目录清单dirac 单位冲激函数disp 显示矩阵和文字内容disttool 概率分布计算交互界面doc 列出指定工具包中所有函数名docsearch 进行多词条检索double 把符号常数转化为16位相对精度的浮点数值对象drawnow 刷新屏幕dsolve 求解符号常微分方程E eedit 矩阵编辑器,打开M文件Ei maple 指数积分eig 矩阵特征值和特征向量end 数组的最大下标,结束for,while,if 语句eps 浮点相对误差EraseMode 图形对象属性error 显示错误信息exit 关闭MATLABexp 指数expand 对指定项展开expm 矩阵指数eye 单位阵ezcontour 画等位线ezcontourf 画填色等位线ezmesh 画网线图ezmeshc 画带等位线的网线图ezplot 绘制符号表达式的二维图形ezplot3 画三维曲线ezpolar 画极坐标曲线ezsurf 画曲面图ezsurfc 画带等位线的曲面图F ffactor 进行因式或因子分解false 按指定大小创建全0逻辑数组feather 从X轴出发的复数向量图,羽毛图feval 函数宏指令fill 多边形填色图find 寻找非数单下标标识findsym 确认表达式中自由符号变量fix 朝零方向取整flag 红-白-蓝-黑交错色图fliplr 矩阵的左右翻转flipud 矩阵的上下翻转floor 朝负无穷大方向取整fminbnd 非线性函数在某区间中极小值fminsearch 单纯形法求多元函数极值点指令for (end) 按规定次数重复执行语句format 设置数据输出格式fourier Fourier变换fsolve 解非线性方程组的最简单格式function 函数文件头functions 观察函数句柄内涵function handle 函数句柄funfun 数值泛函函数和ODE解算器funm 计算一般矩阵函数fzero 单变量函数的零点G ggallery 产生测试矩阵gca 获得当前轴的柄gcf 获得当前图的柄general 通用指令get 获得图柄getframe 获得影片动画图象的帧ginput 用鼠标在图上获取数据global 定义全局变量gradient 梯度gray 线性灰度grid on/off 画坐标网格线H hheaviside 单位阶跃函数help 在线帮助指令helpbrowser 超文本文档帮助信息helpdesk 超文本文档帮助信息helpwin 打开在线帮助窗hidden 网线图消隐开关hist 统计频数直方图histfit 带拟曲线的统计频数直方图hold on/off 图形的保持horner 转换成嵌套形式hot 黑-红-黄-白交错色图hsv 饱和色彩图I ii, j 虚数单位if end 条件执行语句if-else-end 程序分支控制ifourier Fourier反变换ilaplace Laplace反变换imag 复数虚部image 图像impulse 给出系统的冲激响应ind2sub 据单下标换算出全下标inf 或Inf 无穷大inline 创建内联函数input 提示键盘输入int 计算积分int2str 整数转换为字符串inv 矩阵的逆invhilb Hilbert矩阵法求逆阵isa 判断指定变量类别ischar 若是字符串则为真isempty 若是空矩阵则为真isfinite 若是有限数则为真isglobal 若是全局变量则为真ishandle 是否图柄isinf 若是无穷大则为真isletter 串中是字母则为真islogical 若是逻辑数则为真isnan 若为非数则为真isnumeric 若是数值则为真isolate maple的特殊指令isprime 是否质数isreal 若是实数矩阵则为真isspace 串中是空格则为真iztrans Z反变换J jjacobian Jacobian 矩阵jet 变异HSV色图jordan Jordan分解K kkeyboard 键盘获得控制权L llaplace Laplace变换legend 形成图例说明length 确定数组长度light 灯光控制lighting 设置照明模式limit 求极限line 创建线对象LineStyle 图形线对象属性-线型LineWidth 图形属性-线宽linmod2 从SIMULINK模型得到系统的状态方程linspace 线性等分向量load 从磁盘调入数据变量Location 图形对象属性-位置log 自然对数log10 常用对数log2 以2为底的对数logical 将数值转化为逻辑值logspace 对数等分向量lookfor 关键词检索M mmagic 魔方阵maple 进入MAPLE工作空间计算Marker 图形对象属性-点形状MarkerEdge-Color 图形对象属性-点边界色彩MarkerFace-Color 图形对象属性-点域色彩MaekerSize 图形对象属性-点大小material 对象材质max 最大值md 创建目录mean 平均值mesh 三维网线图meshgrid 用于三维曲面的分格线坐标mfun 对MAPLE中若干经典特殊函数的数值计算mfunlist MAPLE经典特殊函数列表mhelp 查阅Maple中的库函数及其调用方法min 最小值minreal 状态方程最小实现mkdir 创建目录mod 模数求余more 命令窗口分页输出的控制开关movie 播放影片动画moviein 影片动画内存初始化mtaylor Taylor级数展开N nNaN或nan 非数nargin 函数输入量的个数nargout 函数输出量的个数ndims 数组的维数norm 矩阵或向量范数normcdf 正态分布累计概率normpdf 服从N分布的随机变量取值x的概率密度normrnd 产生服从N分布的随机数组notebook 创建或打开M-book文件null 零空间num2str 把数值转换为字符串numden 提取公因式O oode45 高阶法解微分方程ones 全1 数组optimsetorth 值空间P ppack 合并工作内存中的碎块pascal Pascal 矩阵path 控制MATLAB的搜索路径pathtool 修改搜索路径pause 暂停pcolor 用颜色反映数据的伪色图peaks 产生peaks图形数据pi 3.1415926535897….pie 饼形统计图pink 淡粉红色图阵plot 直角坐标下线性刻度曲线plot3 三维直角坐标曲线图plotyy 双纵坐标图polar 极坐标曲线图poly 特征多项式,由根创建多项式poly2sym 将多项式转换为符号多项式polyfit 多项式拟合polyval 求多项式的值polyvalm 求矩阵多项式的值pow2 2的幂pretty 习惯方式显示prism 光谱色图阵prod 元素积Q qquad 低阶法数值积分quadl 高阶法数值积分quit 退出MATLABquiver 二维箭头图;主用于场强、流向R rrand 均匀分布随机数组randn 正态分布随机数组random 产生各种分布随机数组randsrc 产生均布数组rank 秩real 复数实部realmax 最大浮点数realmin 最小正浮点数rem 求余数repmat 铺放模块数组reshape 矩阵变维residue 求部分分式表达return 返回roots 求多项式的根rose 频数扇形图;主用于统计rot90 矩阵逆时针旋转90度rotate 旋转指令round 四舍五入取整rref 转换为行阶梯形S ssave 把内存变量存入磁盘sec 正割sech 双曲正割set 设置图形对象属性shading 图形渲染模式shg 显示图形窗sign 函数符号,符号函数simple 运用各种指令化简符号表达式simplify 恒等式简化simulink 打开SIMULINK集成环境sin 正弦sinh 双曲正弦size 确定数组大小slice 切片图solve 求解代数方程组sphere 产生球面数据spinmap 颜色周期性变化操纵spring 青、黄浓淡色sqrt 平方根square 轴属性为方型ss 产生状态方程LTI对象stairs 阶梯形曲线图std 标准差stem 杆图stem3 三维离散杆图str2func 创建函数句柄(punct)strcmp 比较字符串String 图形对象属性-字符串subexpr 运用符号变量置换子表达式subplot 创建子图subs 通用置换指令sum 元素和summer 绿、黄浓淡色图阵surf 三维表面图surfc 带等高线的三维表面图switch-case 多个条件分支sym 产生符号对象syms 定义基本符号对象symsum 符号序列的求和T ttan 正切tanh 双曲正切taylor Taylor级数text 图形上文字标注tf 产生传递函数LTI对象tfdata 从对象中提取传递函数分子分母多项式系数tic 秒表起动title 图形名toc 秒表终止和显示trace 迹trapz 梯形数值积分true 按指定大小创建全1逻辑数组triplequad 三重(闭型)数值积分指令type 显示文件内容V vvar 求方差version MATLAB 版本view 设定3-D图形观测点vpa 给出数值型符号结果W wwhat 列出当前目录上的M、MAT、MEX文件which 确定指定文件所在的目录while end 不确定次数重复执行语句whitebg 图形底色控制who 列出工作内存中的变量名whos 列出工作内存中的变量细节winter 蓝、绿浓淡色X xxlabel X轴名标注xor 异或Y yYcolor 图形对象属性-纵轴颜色ylabel Y轴名标注Z zzeros 全零矩阵zlabel Z轴名标注zoom 二维图形的变焦放大ztrans Z变换Simulink模块Add 求和模块Breaker 开关Current Measurement 电流测量器Dc Voltage Source 直流电压源Discrete Filter 离散滤波器模块Gain 增益模块In1 输入端口模块Integrator 连续函数积分Out1 输出端口模块Parallel RLC Branch RLC并联支路PowerguiProduct 乘法器Random SourceScope 示波模块Series RLC Branch RLC串联支路Simulink SIMULINK基本库Sine Wave 正弦波输出Step 阶跃输出Transfer Fcn 传递函数模块Voltage Measurement 电压测量器总黄酮生物总黄酮是指黄酮类化合物,是一大类天然产物,广泛存在于植物界,是许多中草药的有效成分。
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:设置文件指针的位置。
matlab数学实验胡良剑第五章
第五章%Exercise 1x=[0 4 10 12 15 22 28 34 40];y=[0 1 3 6 8 9 5 3 0];trapz(x,y)%Exercise 2x=[0 4 10 12 15 22 28 34 40];y=[0 1 3 6 8 9 5 3 0];diff(y)./diff(x)%Exercise 3xa=-1:0.1:1;ya=0:0.1:2;[x,y]=meshgrid(xa,ya);z=x.*exp(-x.^2 -y.^3);[px,py] = gradient(z,xa,ya);px%Exercise 4t=0:0.01:1.5;x=log(cos(t));y=cos(t)-t.*sin(t);dydx=gradient(y,x)plot(x,dydx) %dydx函数图,作图观察x=-1时,dydx的值约0.9 %以下是更精确的编程计算方法[x_1,id]=min(abs(x-(-1)));%找最接近x=-1的点,id为这个点的下标dydx(id)%Exercise 5(1)fun=@(x)1/sqrt(2*pi)*exp(-x.^2/2);quadl(fun,0,1)或用trapzx=linspace(0,1,100);y=1/sqrt(2*pi)*exp(-x.^2/2);trapz(x,y)%Exercise 5(2)fun=inline('exp(2*x).*cos(x).^3');quadl(fun,0,2*pi)或用trapzx=linspace(0,2*pi,100);y=exp(2*x).*cos(x).^3;trapz(x,y)%Exercise 5(3)fun=@(x)x.*log(x.^4).*asin(1./x.^2);quadl(fun,1,3)或用trapzx=1:0.01:3;y=feval(fun,x);trapz(x,y)%Exercise 5(4)fun=@(x)sin(x)./x;quadl(fun,1e-10,1) %注意由于下限为0,被积函数没有意义,用很小的1e-10代替%Exercise 5(5)fun=inline('x.^(-x)','x');quadl(fun,1e-10,1) %注意由于下限为0,被积函数没有意义,用很小的1e-10代替%Exercise 5(6)fun=inline('sqrt(1+r.^2.*sin(th))','r','th');dblquad(fun,0,1,0,2*pi)%Exercise 5(7)%先在Editer窗口建立90页函数dblquad2,再在Command窗口clear;fun=@(x,y)1+x+y.^2;clo=@(x)-sqrt(2*x-x.^2);dhi=@(x)sqrt(2*x-x.^2);dblquad2(fun,0,2,clo,dhi,100)%Exercise 6t=linspace(0,2*pi,100);x=2*cos(t);y=3*sin(t);dx=gradient(x,t);dy=gradient(y,t);f=sqrt(dx.^2+dy.^2);trapz(t,f)%Exercise 6另一解法%先写参数方程x=2*cos(t);y=3*sin(t);%计算x'(t)=-2*sin(t),y'(t)=3*cos(t)%4*sin(t)^2+9*cos(t)^2=4+5*cos(t)^2fun=@(t)sqrt(4+5*cos(t).^2);quadl(fun,0,2*pi)%Exercise 7%先算出z的梯度dz/dx=(1-2*x^2)*exp(-x^2-y^2),dz/dy=(1-2*y^2)*exp(-x^2-y^2);%根据曲面面积公式fun=@(x,y)sqrt(1+((1-2*x.^2).*exp(-x.^2-y.^2)).^2+((1-2*y.^2)*exp(-x.^2-y.^2)).^2);dblquad(fun,-1,1,0,2)%或者用下列纯粹离散化解法xa=linspace(-1,1);ya=linspace(0,2);[x,y]=meshgrid(xa,ya);z=x.*exp(-x.^2-y.^2);[zx,zy]=gradient(z,xa,ya);f=sqrt(1+zx.^2+zy.^2);s=0;for i=2:length(xa)for j=2:length(ya)s=s+(xa(i)-xa(i-1))*(ya(j)-ya(j-1))*(f(i,j)+f(i-1,j)+f(i,j-1)+f(i-1,j-1))/4;%每个近似长方体高用四顶点平均值endends%Exercise 8funl=inline('-(-x).^0.2.*cos(x)');funr=inline('x.^0.2.*cos(x)');quadl(funl,-1,0)+quadl(funr,0,1)%Exercise 9 (以I32为例)fun=@(x)abs(sin(x));h=0.1;x=0:h:32*pi;y=feval(fun,x);t1=trapz(x,y)h=pi;x=0:h:32*pi;y=feval(fun,x);t2=trapz(x,y)%步长与周期一致,结果失真q1=quad(fun,0,32*pi)q2=quadl(fun,0,32*pi)%Exercise 10(1)先在Editer窗口建立88页函数deriv,再在Command窗口fun=inline('x.^2.*sin(x.^2+3*x-4)','x');deriv(fun,[1.3 1.5],0.1,1e-3) %取0.1为初始步长%注:书后习题答案错,1.3处导数应为2.4177,1.5处导数应为-11.3330%Exercise 10(2)%先在程序编辑器,写下列函数,保存为ex5_10_2ffunction d=ex5_10_2f(fname,a,h0,e)h=h0;d=(fname(a+h)-2*fname(a)+fname(a-h))/(h*h);d0=d+2*e;while abs(d-d0)>ed0=d;h0=h;h=h0/2;d=(fname(a+h)-2*fname(a)+fname(a-h))/(h*h);end%再在指令窗口执行fun=@(x)x.^2*sin(x.^2-x-2);d=ex5_10_2f(fun,1.4,0.1,1e-3)%Exercise 11%提示:f上升时,f'>0;f下降时,f'<0; f极值,f'=0.%Exercise 12在程序编辑器,写下列函数,保存为ex5_12function I=ex5_12(fname,a,b,n)h=(b-a)/n;x=a:h:b;f=fname(x);I=f(1)+f(n+1);for i=2:nif i-2*floor(i/2)==0I=I+4*f(i);elseI=I+2*f(i);endendI=h/3*I;%再在指令窗口执行ex5_12(inline('1/sqrt(2*pi)*exp(-x.^2/2)'),0,1,50) %注:原题n=5改为偶数n=50%更加符合Matlab风格的编程ex5_12function I=ex5_12f(fname,a,b,n)x=linspace(a,b,n+1);f=fname(x);I=(b-a)/n/3*(f(1)+f(n+1)+2*sum(f(3:2:n))+4*sum(f(2:2:n)));%Exercise 13fun=inline('5400*v./(8.276*v.^2+2000)','v');quadl(fun,15,30)%Exercise 14重心不超过凳边沿。
MATLAB常用命令检索总结计划大全要点总结计划
MATLAB 命令MATLABCOMMANDSANDFUNCTIONSLISTAAABS绝对值、模、字符的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 控制轴刻度微风格的高层指令BBBAR 二维直方图BAR3 三维直方图BAR3H 三维水平直方图BARH 二维水平直方图BASE2DECX 进制变换为十进制BIN2DEC 二进制变换为十进制BLANKS 创立空格串BONE 蓝色彩黑白色图阵BOX 框状坐标轴BREAKWHILE 或FOR环中止指令BRIGHTEN 亮度控制CCCAPTURE (3版从前)捕捉目前图形CART2POL 直角坐标变成极或柱坐标CART2SPH 直角坐标变成球坐标CAT串接成高维数组CAXIS 色标尺刻度CD指定目前目录CDEDIT 启动用户菜单、控件回调函数设计工具CDF2RDF 复数特点值对角阵转为实数块对角阵CEIL 向正无量取整CELL 创立元胞数组CELL2STRUCT 元胞数组变换为构架数组CELLDISP 显示元胞数组内容CELLPLOT 元胞数组内部结构图示CHAR 把数值、符号、内联类变换为字符对象CHI2CDF 散布累计概率函数CHI2INV 散布逆累计概率函数CHI2PDF 散布概率密度函数CHI2RND 散布随机数发生器CHOLCHOLESKY 分解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 创立圆柱DDDBLQUAD二重数值积分DEAL 分派宗量DEBLANK 删去串尾部的空格符DEC2BASE 十进制变换为X进制DEC2BIN DEC2HEX DECONV 十进制变换为二进制十进制变换为十六进制多项式除、解卷DELAUNAYDELAUNAYDEL2 失散LAPLACIAN 三角剖分差分DEMOMATLAB演示DET队列式DIAG 矩阵对角元素提取、创立对角阵DIARYMATLAB 指令窗文本内容记录DIFF数值差分、符号微分DIGITS 符号计算中设置符号数值的精度DIR 目录列表DISP显示数组DISPLAY 显示对象内容的重载函数DLINMOD 失散系统的线性化模型DMPERM 矩阵DULMAGE-MENDELSOHN分解DOS履行DOS指令并返回结果DOUBLE把其余种类对象变换为双精度数值DRAWNOW 更新事件队列逼迫MATLAB刷新屏幕DSOLVE 符号计算解微分方程EEECHOM文件被履行指令的显示EDIT 启动M文件编写器EIG求特点值和特点向量EIGS求指定的几个特点值END控制流FOR等结构体的结尾元素下标EPS浮点相对精度ERROR显示犯错信息并中止履行ERRORTRAP 错误发生后程序能否持续履行的控制ERF偏差函数ERFC偏差补函数ERFCX 刻度偏差补函数ERFINV 逆偏差函数ERRORBAR 带偏差限的曲线图ETREEPLOT 画消去树EVAL 串演算指令EVALIN 跨空间串演算指令EXIST 检查变量或函数能否已定义EXIT 退出MATLAB环境EXP指数函数EXPAND 符号计算中的睁开操作EXPINT 指数积分函数EXPM 常用矩阵指数函数EXPM1PADE法求矩阵指数EXPM2TAYLOR 法求矩阵指数EXPM3 特点值分解法求矩阵指数EYE单位阵EZCONTOUR 画等位线的简捷指令EZCONTOURF 画填色等位线的简捷指令EZGRAPH3 画表面图的通用简捷指令EZMESH 画网线图的简捷指令EZMESHC 画带等位线的网线图的简捷指令EZPLOT 画二维曲线的简捷指令EZPLOT3 画三维曲线的简捷指令EZPOLAR 画极坐标图的简捷指令EZSURF画表面图的简捷指令EZSURFC 画带等位线的表面图的简捷指令FACTOR符号计算的因式分解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浮点运算次数FLOWMATLAB 供给的演示数据FMIN 求单变量非线性函数极小值点(旧版)FMINBND 求单变量非线性函数极小值点FMINS 纯真形法求多变量函数极小值点(旧版)FMINUNC 拟牛顿法求多变量函数极小值点FMINSEARCH 纯真形法求多变量函数极小值点FNDER 对样条函数求导FNINT 利用样条函数求积分FNVAL 计算样条函数区间内随意一点的值FNPLT绘制样条函数图形FOPEN翻开外面文件FOR组成FOR环用FORMAT 设置输出格式FOURIERFOURIER 变换FPLOT返函绘图指令FPRINTF 设置显示格式FREAD 从文件读二进制数据FSOLVE 求多元函数的零点FULL 把稀少矩阵变换为非稀少阵FUNM 计算一般矩阵函数FUNTOOL 函数计算器图形用户界面FZERO求单变量非线性函数的零点GGGAMMA 函数GAMMAINC 不完整函数GAMMALN 函数的对数GCA获取目前轴句柄GCBO获取正履行"回调"的对象句柄GCF获取目前图对象句柄GCO获取目前对象句柄GEOMEAN 几何均匀值GET获知对象属性GETFIELD 获知构架数组的域GETFRAME 获取电影的帧画面GINPUT 从图形窗获取数据GLOBAL 定义全局变量GPLOT 依图论法例绘图GRADIENT 近似梯度GRAY 黑白灰度GRID 画分格线GRIDDATA 规则化数据和曲面拟合GTEXT 由鼠标搁置说明文字GUIDE 启动图形用户界面交互设计工具HHHARMMEAN 调解均匀值HELP在线帮助HELPWIN 交互式在线帮助HELPDESK 翻开超文本形式用户指南HEX2DEC 十六进制变换为十进制HEX2NUM 十六进制变换为浮点数HIDDEN 透视和消隐开关HILBHILBERT 矩阵HIST频数计算或频数直方图HISTC 端点定位频数直方图HISTFIT 带正态拟合的频数直方图HOLD 目前图上重画的切换开关ILAPLACELAPLACEHORNER 分解成嵌套形式HOT黑红黄白色图HSV饱和色图IIIF-ELSE-ELSEIF 条件分支结构IFFT失散FOURIER反变换IFFT2 二维失散FOURIER反变换IFFTN 高维失散FOURIER反变换IFFTSHIFT 直流重量对中的谱的反操作IFOURIERFOURIER 反变换I,J缺省的"虚单元"变量反变换IMAG 复数虚部IMAGE 显示图象IMAGESC 显示亮度图象IMFINFO 获取图形文件信息IMREAD 从文件读取图象IMWRITE 把图象写成文件IND2SUB 单下标转变成多下标INF 无量大INFOMATHWORKS 企业网点地点INLINE 结构内联函数对象INMEM 列出内存中的函数名INPUT 提示用户输入INPUTNAME 输入宗量名INT 符号积分INT2STR 把整数数组变换为串数组INTERP1 一维插值INTERP2 二维插值INTERP3 三维插值INTERPNN维插值INTERPFT 利用FFT插值INTROMATLAB自带的入门指引INV 求矩阵逆INVHILBHILBERT 矩阵的正确逆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反变换JJ,KKJACOBIAN 符号计算中求JACOBIAN 矩阵JET蓝头红尾饱和色JORDAN 符号计算中获取JORDAN标准型KEYBOARD 键盘获取控制权KRONKRONECKER 乘法例则产生的数组LLLAPLACELAPLACE 变换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 解非线性最小二乘问题LULU分解MMMAD 均匀绝对值偏差MAGIC 魔方阵MAPLE&NB,SP; 运作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级数睁开NNNDIMS 求数组维数NAN 非数(预约义)变量NARGCHK 输入宗量数考证NARGIN 函数输入宗量数NARGOUT 函数输出宗量数NDGRID产生高维格点矩阵NEWPLOT准备新的缺省图、轴NEXTPOW2取最靠近的较大2次幂NNZ矩阵的非零元素总数NONZEROS矩阵的非零元素NORM矩阵或向量范数NORMCDF正态散布累计概率密度函数NORMEST预计矩阵2范数NORMINV正态散布逆累计概率密度函数NORMPDF正态散布概率密度函数NORMRND正态随机数发生器NOTEB OOK启动MATLAB和WORD的集成环境NULL零空间NUM2STR把非整数数组变换为串NUM DEN获取最小公分母和相应的分子表达式NZMAX指定寄存非零元素所需内存OOODE1非STIFF微分方程变步长解算器ODE15SSTIFF 微分方程变步长解算器ODE23T 适量STIFF微分方程解算器ODE23TBSTIFF 微分方程解算器ODE45非STIFF微分方程变步长解算器ODEFILEODE 文件模板ODEGET 获知ODE选项设置参数ODEPHAS2ODE 输出函数的二维相平面图ODEPHAS3ODE 输出函数的三维相空间图ODEPLOTODE 输出函数的时间轨迹图ODEPRINT 在MATLAB指令窗显示结果ODESET创立或改写ODE选项构架参数值ONES全1数组OPTIMSET 创立或改写优化泛函指令的选项参数值ORIENT 设定图形的排放方式ORTH值空间正交化PACK采集MATLAB 内存碎块扩大内存PAGEDLG 调出图形排版对话框PATCH创立块对象PATH设置MATLAB 搜寻路径的指令PATHTOOL 搜寻路径管理器PAUSE暂停PCODE创立预解译P码文件PCOLOR 伪彩图PEAKSMATLAB 供给的典型三维曲面PERMUTE 广义转置PI(预约义变量)圆周率PIE二维饼图PIE3三维饼图PINK 粉红色图矩阵PINV 伪逆PLOT 平面线图PLOT3 三维线图PLOTMATRIX 矩阵的散点图PLOTYY 双纵坐标图POISSINV 泊松散布逆累计概率散布函数POISSRND 泊松散布随机数发生器POL2CART 极或柱坐标变成直角坐标POLAR 极坐标图POLY矩阵的特点多项式、根集对应的多项式POLY2STR 以习惯方式显示多项式POLY2SYM 双精度多项式系数转变成向量符号多项式POLYDER 多项式导数POLYFIT 数据的多项式拟合POLYVAL 计算多项式的值POLYVALM 计算矩阵多项式POW22的幂PPVAL 计算分段多项式PRETTY 以习惯方式显示符号表达式PRINT 打印图形或SIMULINK 模型PRINTSYS 以习惯方式显示有理分式PRISM 光谱色图矩阵PROCREAD 向MAPLE输送计算程序PROFILE 函数文件性能评估器PROPEDIT 图形对象属性编写器PWD显示目前工作目录QUAD 低阶法计算数值积分QUAD8 高阶法计算数值积分(QUADL) QUIT 推出MATLAB 环境QUIVER 二维方向箭头图QUIVER3 三维方向箭头图RRRAND 产生均匀散布随机数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实数块对角阵转为复数特点值对角阵RSUMSRIEMANN 和SSSAVE把内存变量保存为文件SCATTER 散点图SCATTER3 三维散点图SEC正割SECH双曲正割SEMILOGXX 轴对数刻度坐标图SEMILOGYY 轴对数刻度坐标图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 创立多个符号对象TTTAN 正切TANH 双曲正切TAYLORTOOL 进行TAYLOR迫近剖析的交互界面TEXT 文字说明TF创立传达函数对象TIC 启动计时器TITLE 图名TOC封闭计时器TRAPZ 梯形法数值积分TREELAYOUT 睁开树、林TREEPLOT 画树图TRIL 下三角阵TRIM 求系统均衡点TRIMESH 不规则格点网线图TRISURF 不规则格点表面图TRIU 上三角阵TRY-CATCH 控制流中的TRY-CATCH 结构TYPE 显示M文件UUUICONTEXTMENU 创立现场菜单UICONTROL 创立用户控件UIMENU 创立用户菜单UNMKPP 逐段多项式数据的反清晰化UNWRAP 自然态相角UPPER变换为大写字母VVVAR方差VARARGIN 变长度输入宗量VARARGOUT 变长度输出宗量VECTORIZE 使串表达式或内联函数适于数组运算VER版本信息的获取VIEW 三维图形的视角控制VORONOIVORONOI 多边形VPA随意精度(符号类)数值WWWARNING 显示警示信息WHAT 列出目前目录上的文件WHATSNEW 显示MATLAB 中README 文件的内容WHICH 确立函数、文件的地点WHILE 控制流中的WHILE 环结构WHITE 全白色图矩阵WHITEBG 指定轴的背风景WHO 列出内存中的变量名WHOS列出内存中变量的详尽信息WINTER 蓝绿调冬色图WORKSPACE 启动内存阅读器XX,YY,ZZXLABELX 轴名XOR 或非逻辑YESINPUT 智能输入指令YLABELY 轴名ZEROS全零数组ZLABELZ 轴名ZOOM 图形的变焦放大和减小ZTRANS 符号计算Z变换。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第四章
1、多项式 P63 y=polyval (p, x) 求得多项式 p 在 x 处的值 y,x 可以是一个或多个点 p3=conv(p1,p2) 返回多项式 p1 和 p2 的乘积 [p3,r]=deconv(p1,p2) p3 返回多项式 p1 除以 p2 的商, r 返回余项 x=roots(p) 求得多项式 p 的所有复根. p=polyfit(x,y,k)用 k 次多项式拟合向量数据(x, y), 返回多项式的降幂系数 2、函数零点 P65 ( 1) 、x=fzero(Fun, x0) 一元函数求实根:返回一元 函数 Fun 的一个零点, 其 中 Fun 为函数 (三种形式) x0 为标量时, 返回函数在 x0 附近的零点; x0 为区间[a, b]时, 返回在[a,b]中的零点 (2) 、[x,f,h]=fsolve(Fun, x0), 方程(组)数值求解 P x: 返回一元或多元函数 Fun 在 x0 附近的一个零点,其中 x0 为迭代初值向量; f:返回 Fun 在 x 的函数值, 应该接近 0;
5、矩形区域重积分 P87 z=dblquad(Fun,a,b,c,d) 求得二元函数 Fun(x,y) 的重积分, a、b 为变量 x 的下、上限;c、d 为变量 x 的下、上限. z=triplequad(Fun,a,b,c,d,e,f) 求 得 三 元 函 数 un(x,y,z) 的 重 积分 , 格 式 类 似 dblquad。
第五章
1、数值差分 P84 n 维向量 x=(x1, x2, , xn) 的差分定义为 n-1 维向量 x = (x2-x1, x3-x2, , xnxn-1)。 diff(x) 如果 x 是向量,返回向量 x 的差分,如果 x 是矩阵,则按各列作差分。 diff(x, k) k 阶差分,即差分 k 次 例:x=[1 2 3 4 5]; diff(x), diff(x,2) A=[1 2;3 4; 5 6],diff(A), diff(A,2) 2、数值导数和梯度 P84 q=polyder(p) 求得由向量 p 表示的多项式导函数的向量表示 q. 例:p=[1 2 3 4];q=polyder(p) Fx=gradient(F,x) 返回向量 F 表示的一元函数沿 x 方向的导函数 F ’(x).其中 x 是与 F 同维数的向量. [Fx,Fy]=gradient(F,x,y) 返回矩阵 F 表示的二元函数的数值梯度 (F ’x,F’y), 当 F 为 m×n 矩阵时, x,y 分别为 n 维和 m 维的向量. 3、梯形积分法 z=trapz(x,y) 返回积分的近似值,其中 x 表示积分区间的离散化向量; y 是与 x 同维数的向量,表示被积函数 。 4、高精度数值积分 P86 quad(Fun,a, b) 自适应步长 Simpson 积分法求得 Fun 在区间[a, b]上的定积分 z=quadl(Fun,a,b) 高精度 Lobatto 积分法.格式同 quad. 注意:trapz,quad,quadl 均不能用于求解反常积分以及一些奇异积分。
h: 返回值如果大于 0,说明计算结果可靠,否则计算结果不可靠。 3、函数极值 P68 min(y) 返回向量 y 的最小值 max(y) 返回向量 y 的最大值 [x,f]=fminbnd(fun,a,b) x 返回一元函数 y=f(x)在区间[a,b]内的局部极小值点, f 返回局部极小值,fun 为函数句柄或匿名函数或 inline 函数。 [x,f]=fminsearch(fun,x0) x 返回多元函数 y=f(x)在初始值 x0 附近的局部极小值 点,f 返回局部极小值. x, x0 均为向量。 注意:在使用 fsolve 和 fminserrch 等指令时,多变量必须合写成一个向量变量, 如用 x(1),x(2) 4、最小二乘法 P69 c= lsqnonlin (Fun,c0) 使用迭代法搜索最优参数 c. 其中 Fun 是以参数 c(可以是向量)为自变量的函数, c0 为参数 c 的近似初值(与 c 同维向量), 表示误差向量 y-f(c,x)(x, y 为数据向量), 作为迭代初值。 c=lsqcurvefit(Fun2,c0, x, y) 从外部输入数据, 这里 Fun2 为两变量 c 和 x 的函数 f(c, x)。
第一章
1、数组表示方法 P9 a=[1,2,3;4,5,6;7,8,9] 或者 a=[1 2 3 456 7 8 9]; 2、等差数列构造一维数组 P9 b=0:3:10(0 3 6 9) b=0:10(0-10,增量默认为 1) b=linspace(0,10,4)(将区间[0,10]等分为 4-1=3 份) 3、数组的提取和拼接 P10 b([1,end]) 提取数组 b 的首尾元素 d=a(2,1:3) 提取 a 的第二行,1-3 列,可简写为 d=a(2,:) [s,t]=max(b) 》s=10 t=4,s 为最大值,t 为最大值的编址 4、数组运算 P12 数组加减 A+B,A-B 数组乘数组 A.*B 数组除法 左除 A.\B,右除 B./A 数与数组的混合运算 k+A,k-A,k*A,A*k,A.^k,k.^A,k./A 注意:点运算只有点乘、点乘方、点除三个,表示元素之间的运算。 左除 A.\B 是 B 除以 A 中对应元素,右除则反之 5、数学函数 P14 exp 指数函数 log 自然对数 ln log10 以 10 为底的对数
ห้องสมุดไป่ตู้
第七章
1、符号对象的定义 P117 s=sym(str) 将数值或字符串 str 转化为符号对象 s, 数值为有理表示; s= sym(num,’d’) 将数值表达式转化为符号表达式,数值用十进制表示; syms var1 var2 定义 var1, var2, 为符号变量 subs(s, old, new) 将符号表达式 s 中的符变量 old 用 new 代替 2、计算精度和数据类型转换 P118 digits(n) 将数值计算精度设为 n 位; x=vpa(s) 求 s 的数值结果 x=vpa(s,n) 采用 n 位计算精度求 s 的数值结果 double(s) 将符号对象转化为双精度数值 char(s) 将符号对象转化为字符串 数据类型转换图指令详解 P118 3、符号函数计算 P120 factor(expr) 对 expr 作因式分解 expand(expr) 将 expr 展开 collect(expr,v) 将 expr 按变量 v 合并同类项 simple(expr) 将 expr 化简 g=finverse(f,v) 求函数 f(v)的反函数 g(v) fg=compose(f, g) 求函数 f(v)和 g(v)的复合函数 f(g(v)) [n,d]=numden(expr) 分式通分,n 返回分子,d 返回分母 funtool 函数分析图形界面 4、极限和级数及微分 P122 limit(s,x,a) 返回符号表达式 s 当 x->a 时的极限 limit(s,x,a,’right’) 返回 s 当 x->a 时的右极限 limit(s,x,a,’left’) 返回 s 当 x->a 时的左极限 symsum(s,n,a,b) 返回符号表达式 s 表示的通项当自变量 n 由 a 到 b 的和 diff(s,x) 返回符号表达式 s 对 x 的导函数 diff(s,x,n) 返回 s 对 x 的 n 阶导函数。 taylor(s,n,a,x) 返回符号表达式 s 在 a 点 aylor 展开到 n-1 次式.自变量为 x taylortool Taylor 分析图形界面 jacobian(f,x) 返回向量函数 f 的 Jacobian 矩阵 5、积分 P124 int(s) 符号表达式 s 的不定积分.
第三章
1、特殊矩阵生成 P46 zeros(m,n) m 行 n 列的零矩阵; ones(m,n) m 行 n 列的元素全为 1 的阵;
eye(n) n 阶单位矩阵; rand(m,n) m 行 n 列[0,1]上均匀分布随机数矩阵 2、矩阵分析 P48 rank(A) 秩 det(A) 行列式; inv(A) 逆矩阵; null(A) Ax=0 的基础解系; 例:A=[1 2 3; 4 5 6; 7 8 9], rank(A), det(A), inv(A), null(A) 3、特征值与标准形 P49 eig(A) 方阵 A 的特征值 [V, D]=eig(A)返回方阵 A 的特征值和特征向量。 其中 D 为的特征值构成的对角阵, 每个特征值对应的 V 的列为属于该特征值的一个特征向量。 [V, J]=jordan(A) 返回 A 的相似变换矩阵和约当标准形 例:A=[1 2; 3 4]; t= eig(A), [V, D]=eig(A) 注意:矩阵是一个二维数组,它的加减数乘运算与数组运算是一致的,但有两点 需要注意:1、对于乘法、乘方和除法等三种运算,矩阵运算与数组运算符及含 义都不同。矩阵运算按线性变换定义,符号中无点;数组运算按对应元素运算定 义,使用点运算。2、数与矩阵加减、矩阵除法在数学上没有一一,MATLAB 中 定义了这两类运算,数与矩阵加减与数组相同,A\B 和 B/A 分别表示 AX=B 和 XA=B 的解。 4、判断一个矩阵是否可以相似对角化,求出其特征向量,若线性无关即可。
inline 内嵌函数 feval 函数求值 使用格式 fun=inline(‘expr’,arg1,arg2,……) fun=@(arg1,arg2,……)expr,这里 expr 为函数表达式,arg1,arg2 为自变量名字符 串。 如,计算 s
n
n 1
100
1
2
>>fname=inline(‘sum(1./(1:m).^2)’,’m’) >>feval(fname,100) 二、作图 1、曲线图 P31 plot(x,y) 以数据(x(i), y(i))为节点的折线图, 其中 x, y 为同长度的向量 fplot(‘fun’,[a,b]) 函数 fun 在区间[a,b] 上的函数图 plot3(x,y,z) 空间曲线图, 其中 x, y, z 为同长度的向量 图形导出: file→Save; Edit →Copy figure 折线图:plot([1 2 3 4],[2 1 8 3]) 曲线图: fplot(‘sin(x)’,[0,2*pi]) 曲线图: x=0:pi/50:2*pi; y=sin(x); plot(x,y) 2、曲面图 P33 [x,y]=meshgrid(xa,ya)当 xa, ya 分别为 m 维和 n 维行向量, 得到 x 和 y 均为 n 行 m 列矩阵。meshgrid 常用于生成 X-Y 平面上的网格数据。 mesh(x,y,z) 绘制网面图,是最基本的曲面图形命令,其中 x,y,z 是同阶矩阵,表示曲 面三维数据。 surf(x,y,z) 绘制曲面图,与 mesh 用法类似 contour(x,y,z) 绘制等高线图,与 mesh 用法类似。 contour3(x,y,z) 绘制三维等高线图,与 mesh 用法类似。 3、一些图例和说明 P34 title(‘字符串’) 图像标题说明; xlabel,ylabel,zlabel 坐标轴 x,y,z 说明; hold on/hold off 保留/释放现有图形 axis([a,b,c,d]) 确定坐标轴范围 a<x<b,c<y<d axis([a,b,c,d,e,f]) 定制 3 维坐标轴范围 figure\close 开\关一个新图形窗口 subplot(m,n,k) 将图形窗口分为 m*n 个子图,指向第 k 幅图 legend(str1,str2, ...) 图例