几个命令.doc

matlab文件读写1. textread最常用的方式 A,B,C,. = textread(filename,format)这样可以直接将文件的内容全部读出，按顺序依次放入相应变量中，得到向量结果如果文件太大，则一则读出时间太长，二则memory可能受不了，可以指定读出N个数据，N 0 y = x * sin(x)endy = x2 + exp(x) 但此例中end后面的语句y = x2 + exp(x)无论是否大于0，都是要运行的，不能满足题目要求，故应用下面的格式2。 格式2： if 逻辑表达式 程序语句组1else 程序语句组2end if语句计算所跟的逻辑表达式的值，若值为真，就做其下面的程序语句组1，然后跳到end后的程序继续做。否则，值为假就做else后程序语句组2。然后再接着做end后的程序。 例7 当从键盘输入自变量的值，由分段函数给出的值。 解 程序如下： x = input(x=) %屏幕提示x=，由键盘输入值赋给xif x 0 y = x3;else y = 5 * x2;endx,y2多条件语句 格式： if 逻辑表达式1程序语句组1elseif 逻辑表达式2程序语句组2elseif 逻辑表达式n程序语句组nelse程序语句组n+1end if语句组判断所跟的逻辑表达式1的值，若值为真，就做其下面的程序语句组1，然后跳到end后的程序继续做。否则，值为假就再判断elseif后的逻辑表达式2的值。若值为真，就做其下面的程序语句组2，否则继续做下面的程序语句。 例8 当从键盘输入自变量的值，由分段函数给出的值。 解 程序如下： x = input(x=) %屏幕提示x=，由键盘输入值赋给xif x 0 y = -1;elseif x = 0 y = 0;else y = 1;endx, y2.13循环语句1for-end循环 格式： for 循环变量 = 初值 : 步长 : 终值循环体语句组end 例9 求自然数前10项之和。 解 程序如下： sum=0; %变量sum表示和，初值为0for i=1:10 sum = sum + i;endsum 运行结果如下： sum = 552while-end循环 格式： while 逻辑表达式 循环体语句组 end例10 求自然数前n项和，项数n由键盘输入。 解 程序如下： n = input(n=) sum = 0; k = 1; %给和与循环变量赋初值while k = n sum = sum + k; k = k + 1;endsum 运行程序时输入n = 100，运行结果如下： n = 100sum = 50502.14关系与逻辑运算 在关系运算与逻辑运算中，非零数为真，零值为假。在输出时对于真值输出为1，对于假值输出为0。1关系运算（表2.3）表2.3 关系运算关系操作符对应函数说明=Eq(A, B)等于=ne(A, B)不等于gt(A, B)大于=ge(A, B)大于等于 例11 A=3 4 5 6 7 8; B=2 4 6 7 5 8; C=eq(A, B) 运行结果如下：C = 0 1 0 0 0 1 D=ne(A,B) 运行结果如下：D = 1 0 1 1 1 0 E=lt(A,B) 运行结果如下：E = 0 0 1 1 0 0 例12 A=4; B=2 3 4 5 6; C=eq(A,B)运行结果如下：C = 0 0 1 0 0 D=AB运行结果如下：D = 1 1 0 0 0 E=A A=1 2 0; 0 4 0; B=3 2 0; 0 0 0; C=and(A, B) 运行结果如下：C = 1 1 0 0 0 0 E=A|B运行结果如下：E = 1 1 0 0 1 0 D=A 运行结果如下：D = 0 0 1 1 0 1 w=xor(A, B)运行结果如下：w = 0 0 0 0 1 0 H=any(A)运行结果如下：H = 1 1 0 K=all(A)运行结果如下：K = 0 1 02.15 MATLAB帮助1help 命令（1）当不知系统有何帮助内容时，可直接输入help以寻求帮助: help（回车）（2）当想了解某一主题的内容时，如输入： help syntax （了解Matlab的语法规定）（3）当想了解某一具体的函数或命令的帮助信息时，如输入： help sin （了解函数sin的相关信息）将显示图2.11帮助信息：图2.11 函数sin帮助信息点击图2.11中的doc sin将显示函数sin更详细的使用帮助，如图2.12所示。图2.12 函数sin详细帮助信息2lookfor命令现需要完成某一具体操作，不知有何命令或函数可以完成，如输入： lookfor line （查找与直线、线性问题有关的函数）3模糊查询 输入命令的前几个字母，然后按Tab键，系统就会列出所有以这几个字母开头的命令。如输入si，按Tab键，系统将列出所有以si开头的函数，如图2.13。图2.13 模糊查询 求极限 求极限之前应先使用syms命令将函数自变量声明为符号变量 例1 定义x，y为符号变量。 syms x y; （2）MATLAB求极限的命令见表 。表 求极限命令数学运算MATLAB命令limit(f)limit(f, x, a)limit(f, x, a, left)limit(f, x, a, right) 说明：inf表示无穷大；NaN代表不定值，即0/0或/；pi表示。 例 求极限：；limit(x(1/2)-2(-1/x),x,0,right) ans = 0若不加入选项right，则返回NaN。求导数的命令 （1）diff(f) %关于符号变量求函数的一阶导数 （2）diff(f, x) %求函数对变量的一阶导数 （3）diff(f,n) %关于符号变量求函数的阶导数（4）diff(f,x,n) %关于变量对函数求阶导数求函数的导数。syms xy=10x+x10+log(x)diff(y) ans = 10x*log(10)+10*x9+1/x求偏导：syms x y;Z=y*log(x)+x*log(y)diff(diff(Z,x),y) 设，求。 解 程序如下： syms x;y=log(1+x);a=diff(y,x,2)a = -1/(1+x)2subs(a,x,1)（同下面等效）x=1; eval(a) %eval将符号表达式转换成数值表达式 ans = -0.2500 隐函数的导数设，求。 解 ，先求，再求，则。 程序如下： syms x ydf_dx=diff(exp(y)+x*y-exp(1),x);df_dy=diff(exp(y)+x*y-exp(1),y);-df_dx/df_dy 运行结果为： ans =-y/(exp(y)+x) 不定积分与定积分格式 int(f) %对关于符号变量求不定积分： int(f, v) %对关于变量求不定积分： int(f, a, b) %对关于符号变量从到求定积分： int(f, v, a, b) %对关于变量从到求定积分： 注：软件运行的不定积分的运行结果中省略了加任意常数，定积分中的用inf表示，用-inf表示。 计算下列不定积分： （1）； （2）； symx x %先定义符号变量（1）int(sin(x)10)（2）int(exp(-2*x)*sin(3*x) 例2 求定积分 （1）； （2）。int(sin(x)/x,x,0,1); vpa(sinint(1),4) syms x %先定义符号变量（1）int(abs(x-2),x,0,4)ans =4（2）int(exp(-x2),x,0,1)ans =1/2*erf(1)*pi(1/2) 其中erf是误差函数，它不是初等函数，可使用eva命令将其转化成数字表达式： eval(ans)ans = 0.7468 求微分方程通解的命令函数dsolve( )用于求解微分方程。表示：（为缺省的自变量），表示对的阶导数。 （1）单个方程 调用格式： y = dsolve(Dy = 1 + y2) %求一阶微分方程的通解 y = dsolve(Dy = 1 + y2, y(0) = 1) %求一阶微分方程带初始条件的特解y = dsolve(D2y = cos(2*x), x) %求二阶微分方程的通解 y = dsolve(D2y = cos(2*x), y(0) = 1, Dy(0) = 0, x) %求二阶微分方程带初始条件%的特解 （2）方程组命令格式：eqal 表示一个微分方程condition表示：初始条件u, v=dsolve(eqal1, eqal2 , condition1,condition2) %两个方程，两个输出 S = dsolve(Df = g, Dg = h, Dj = f) %三个方程，结构输出 输出成员：S.f，S.g，S.h 例 求微分方程的通解。 解 程序如下： y=dsolve(Dy+2*x*y=x*exp(-x2),x) 运行结果如下：y =(1/2*x2+C1)*exp(-x2) 例3求微分方程在初始条件下的特解。 解 程序如下： y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) 运行结果如下：y =(exp(x)+exp(1)/x 例4 求微分方程的通解。 解 程序如下： y=dsolve(