scilab中文帮助手册.doc

（1） 程序设计programming：abort: 中断计算说明：abort命令中断当前计算并给出Scilab的提示符(“-”)。如果用户使用pause命令进入高层界面，可使用abort命令回到初始层界面，出现Scilab的提示符(“-”)。Ans: 计算结果说明：如果表达式的值没有被赋与任何变量，那么Scilab自动产生一个变量Ans，用来存储当前表达式的值。backslash ()：矩阵左除调用格式：如x=Ab说明：backslash ()代表矩阵左除，x=Ab是方程A*x=b的解。如果A是一个方阵或者可逆矩阵（nonsingular或者invertible），那么x=Ab是方程A*x=b的唯一确定的解，且等于x=inv(A)*b.（通俗的说，就是方程数目等于未知数个数，所以解唯一。如A=1 2 1;3 2 2;2 3 1,b=8;13;11）如果A不是一个方阵，那么x是最小二乘解，此时，该解x使norm(A*x-b)最小。（通俗的说，就是方程数目不等于未知数个数。如果方程数目大于未知数个数，本来是没有准确解的，但可以根据最小二乘法原理拟合出最佳解。如A=rand(3,2);b=1;1;1。此时，由于A不是一个方阵，无法使用inv(A)命令，因此x=Ab等效于x=pinv(A)*b. 如果方程数目小于未知数个数，应该有无穷个解，用x=Ab命令得出的解只是其中的一个解，而用pinv(A)*b命令得出的解则比较粗糙，没有用x=Ab命令得出的解精确。）如果A是一个列满秩矩阵（矩阵的秩等于列数），那么此最小二乘解x=Ab是唯一能使norm(A*x-b)最小的解。如果A不是一个列满秩矩阵，那么此最小二乘解x=Ab不能使norm(A*x-b)最小，使norm(A*x-b)最小的解应该是x=pinv(A)*b）。inv(A)要求A是方阵，pinv(A)则不要求A是方阵bool2s ：将布尔矩阵转化为0、1矩阵调用格式：bool2s(x)参数：x是一个布尔向量或者布尔矩阵或者常数矩阵。说明：如果x是一个布尔矩阵，bool2s(x)返回一个0、1矩阵，1代表真，0代表假。如果x是一个普通矩阵，bool2s(x)也返回一个0、1矩阵，1代表非0的数字。boolean：它表示Scilab对象或者布尔型变量和操作符 与或非& | 说明：布尔型变量为%T（代表真）%F (代表假)。这些变量可以用来定义矩阵。布尔型型矩阵可以像普通矩阵操作，如元素的抽取、插入与连接，但不能进行如( +, *, -, , )等运算，布尔型型矩阵有三个特殊的操作符号：与或非& | brackets括号：包括左右圆括号()，左右方括号.方括号调用格式：a11,a12,.;a21,a22,.;.圆括号调用格式: s1,s2,.=func(.) 参数a11,a12,. : 任何合适的矩阵 (实数, 多项式, 有理数等);参数s1,s2,. : 任何变量名;说明： 用于向量定义或者矩阵连接或者用在调用函数的左边存储变量，而（）用于函数调用。 在用于矩阵连接时，空格、逗号用来隔开每“行”中的元素，分号、回车则用来隔开各“行”。 用在调用函数的左边存储变量时，必须用逗号隔开个变量。break -结束循环说明：在for 或者 while 循环中，break命令强制结束循环试试看k=0; while 1=1, k=k+1; if k 100 then break, end; endcallFortran or C 用户程序调入调用格式有两种长型 y1,.,yk=call(ident,x1,px1,tx1,.,xn,pxn,txn,out,ny1,my1,py1,ty1,.,nyl,myl,pyl,tyl) 短型 y1,.,yk=call(ident,x1.,xn) 参数说明：ident：字符串，代表所要调用的子程序的名字。Xi：实数矩阵或者字符串，代表输入变量pxi,pyi：整数，代表各输入变量在调入ident程序时的相对位置txi,tyi：字符”d”r” “c” “i”，代表输入变量的数据类型，d代表双精度型，r代表实型（浮点型），c代表字符型，I代表整型。Out：关键词，用来格开输入变量于输出变量。当有Out时,表明是长型调入,而没有Out是短型调入.ny1, my1是输出变量的大小功能说明：call命令用于在scilab于调入Fortran or C 用户程序，所要调入的程序应该预先与scilab链接，链接的方法如下：1.软链接:在scilab进程中使用命令link链接2.硬链接:先用scilab下已有的程序(default/Ex-fort.f)作为调入程序,再在文件(default/Flist)中增加一个入口,然后再使用命令make bin/scilex与scilab连接.call语句有两种调用形式,长型与短型.短型调用语法简单,能够很快地编译代码,但必须写一个小小的C或者F接口程序.长型调用则不需要,但语法复杂,代码解释速度也比较慢.例子:C语言/简单的C语言 f1=#include void fooc(c,a,b,m,n) double a,*b,c; int *m,*n; int i; for ( i =0 ; i k则为空集). j:d:k生成一个向量j, j+d, ., j+m*d.该操作符可以用来挑选向量或者矩阵中的某一行,某一列或者某一个元素A(:) 是指A中的元素按从左到右从上到下的顺序排成的一列.A(:,j) 是指A中的第j列元素. A(j:k) 是指A(j),A(j+1),.,A(k) A(:,j:k) 是指A(:,j),A(:,j+1),.,A(:,k) A(:)=w 是指用w将矩阵A填满.comma - (,) 列,指令与参数分隔符说明:该操作符用来分隔函数中的参数,或者用来分隔单行矩阵中的列,空格也可以用来分隔单行矩阵中的列,但用逗号更加合适.另外也可以用来分隔指令,分号;也可以用来分隔指令,但不会将结果显示在屏幕上.例子:a=1,2,3;4,5,6;a=1,b=1;c=2comments -注释注释行以/开始,它不会为scilab所编译continue -控制语句中的关键词,使程序转向下一个循环.说明:在for和while循环中,命令continue使程序转向下一个循环,而忽略该次循环中continue语句和end语句之间的指令.例子1:for k=1:10,K=k;if k2&kj+1&k0,n就是捕获的错误类型号；如果n0,n就是所被清除的错误类型号；如果n0,表示所有类型的错误都会被清除（默认）；error_table - 错误信息表说明：这一页列出预先定义的错误信息表，以及它们对应的序号。其中一些错误信息是Scilab软件自己的剖析器错误(parser errors)或者特殊内建错误(specific builtin errors),其他的错误信息比较通用，可以Scilab函数中使用。带*号的是语句error(n,pos)中句柄有关的。（表略）error - error messages调用格式：error(message ,n) error(n) error(n,pos)参数：message：一个字符串，要显示的错误信息。n：整数。错误信息的序号。pos : 整数。错误信息的参数。说明：error函数允许发布一个错误信息来处理错误，默认情况下error函数停止执行当前命令返回提示符界面。这个默认设置可以用errcatch或execstr(.,errcatch)函数更改。error(message)显示message所代表的错误信息，与此错误信息对应的序号是10000error(message,n)显示message所代表的错误信息，与此错误信息对应的序号是n，n应该大与10000error(n) 显示已预定义好了的错误信息，该错误信息的序号是n有一些预定义的错误信息要求有参数，在这种情况下，必须用error(n,pos)中的pos来指定参数值，在其他情况下，pos参数可以不写。etime - 耗时调用格式：e = etime(t2,t1)参数：t2：是一个6-10个值的向量t1：是一个6-10个值的向量e : t2与t1之间的秒数说明：当t1和t2都是含有10个值的向量时，他们的格式必须与 getdate函数返回的格式一样。在这种情况下， 它们的第三个，第四个，第五个值都被忽略。当t1和t2都是含有6个值的向量时，他们的格式必须Year Month Day Hour Minute Second格式一样，其中秒可以精确到毫秒。t1和t2必须大小相同，t1和t2也可以是矩阵，但是该矩阵每一行都必须符合以上所述格式。例子：t1=2004 06 10 17 00 12.345t2=2004 06 10 17 01 13.965E1=etime(t2,t1)t1=2004 06 24 162 5 10 17 00 12 345t2=2004 06 24 162 5 10 17 01 13 965E2=etime(t2,t1)注释：当当t1和t2都是含有10个值的向量时，如2004 06 24 162 5 10 17 00 12 345，2004表示年份，06表示月份， 24表示周数， 162表示天数， 5表示本周中的第几天（以周日为第一天），10表示本月中的第几日， 17表示时， 00表示分， 12表示秒的整数部分，345表示秒的小数部分。evstr - 计算表达式的值调用格式：H=evstr(Z)H,ierr=evstr(Z)参数Z : 字符串矩阵M或者列表list(M,Subexp)M : 字符串矩阵Subexp : 字符串向量H : 矩阵 ierr : 整数,误差指示说明：该函数的运算结果返回字符串矩阵M中表达的的计算结果。矩阵的每一个元素必须定义一个有效的Scilab表达式。如果对M中表达式的计算有错，那么，对于返回一个值的表达式H=evstr(M)来说，它对错误不予理睬(raises the error as usual)；对于返回二个值的表达式H,ierr=evstr(M)来说，虽然显示出错，但会返回一个出错序号给ierr。如果z是一个列表，subexp是一个字符串向量，用来定义子字符串，对这些子字符串的计算是在对M的计算之前。这些子字符串必须要在M中以%(k)的形式提及，其中k是字符串向量subexp中子表达式的序号。evstr(a=1)是无效语句，请用execstr代替。例子：a=1; b=2; Z=a,b;H=evstr(Z) /对字符串矩阵中的元素a,b求值，并将结果传给矩阵H，结果：H=1 2a=1; b=2; Z=list(%(1),%(1)-%(2),a+1,b+1);/其中M%(1),%(1)-%(2) ;Subexp=a+1,b+1;/%(1)就表示a+1,%(1)-%(2)就表示(a+1)-(b+1)evstr(Z)/计算a+1与(a+1)-(b+1)结果：H=2 -1execstr - 执行scilab字符串形式的代码调用格式：execstr(instr) ierr=execstr(instr,errcatch ,msg) 参数：instr : 字符串向量, 需要执行的Scilab指令. ierr : 整数, 0或者错误号. msg : 值为m或n的字符串.默认值为n.说明：该命令执行由参数instr所给的指令。如果errcatch没有出现，出错处理照常进行。如果errcatch出现了，并且在执行指令的过程中产生了错误，那么execstr函数本身不发布出错信息，但它会在出错点上中断执行instr指令，并启用错误号为ierr的出错信息，在这种情况下，出错信息的显示由msg选项控制。m: 表示错误信息会显示而且会被记录。n: 不显示出错信息,但错误信息也会记录(参考lasterror命令)，这个是默认选项.ierr=execstr(instr,errcatch)可以处理语法错误。 例子：execstr(a=1) / 执行a=1execstr(1+1) / 该语句其实没有用途execstr(if %t then; a=1; b=a+1; else b=0 end)execstr(a=zzzzzzz,errcatch) /不显示出错信息execstr(a=zzzzzzz,errcatch,m) /显示出错信息/句法错误举例execstr(a=1?02,errcatch) /不显示出错信息lasterror(%t) /查看出错信息execstr(a=1 2 3),errcatch) /不显示出错信息lasterror(%t) /查看出错信息exec - 执行脚本文件调用格式：exec(path ,mode) exec(fun ,mode) ierr=exec(path,errcatch ,mode) ierr=exec(fun,errcatch ,mode) 参数：path : 字符串, 脚本文件的路径 mode : 一个整数, 执行模式（如下）fun : scilab函数 ierr : 整数, 0 或者出错号说明：exec(path ,mode)连续地执行由路径path给定的scilab脚本文件，,mode是可选项，有以下几种不同的模式：0 : 默认选项-1 : 不显示任何信息1 : 显示出每一行命令2 : 显示出提示符3 : 既显示出每一行命令，又显示出提示符（312）4 : 在显示出提示符之前暂停. 当按回车之后又继续执行命令7 : 既又暂停，又有提示符，还会显示出每一行命令，这对演示很有用途（7124）exec(fun ,mode)命令以脚本形式执行函数fun，没有输入变量，没有输出变量，也没有特殊的变量环境。这种形式很有效，因为脚本代码可以被预编译（参考getf, comp）。这种脚本运算方法允许在库中以函数的形式存储脚本。如果在执行的时候遇到了错误，并且errcatch标志出现了，exec不会发布出错信息，但它会中断执行指令，并启用错误号为ierr的出错信息，如果errcatch标志没有出现，标准的出错处理照常进行。例子：mputl(a=1;b=2,c:/exam)/mputl把a=1;b=2转化为ascii码，并存储为c:/exam.exec(c:/exam)/执行该脚本文件whos -name a /查看以a为首字母的变量，deff(y=foo(x),a=x+1;y=a2)clear a bfoo(1)/ a是一个在函数foo的环境中创建的变量.当foo返回一个值时，它就会被破坏，对比前后两个whos-name a 所产生的结果就可以知道whos -name a x=1exec(foo)/a和y在当前环境又产生了whos -name a /对比前两个whos -name a 所产生的结果exists - 检查变量是否存在调用格式exists(name ,where) 参数：name : 字符串 where : 可选择的字符串，默认为all说明：exists(name)如果又名为name的变量则返回1，否则返回0注意：exists函数的返回值可能与环境有关exists(name,local)如果在当前函数的环境下存在名为name的变量则返回1，其它情况下返回0例子：deff(foo(x),.disp(exists(a12),exists(a12,local) disp(exists(x),exists(x,local)foo(1)a12=;foo(1)exit - 结束当前Scilab进程说明：结束当前Scilab进程。external - Scilab对象，外部函数与对象说明：特殊命令下外部函数与程序。一个外部对象（external）就是一个用作高级命令（some high-level primitives）参数的函数或者程序，例如ode , optim , schur等命令外部对象（函数或者程序）的调用次序是由高级命令设置相关参数来决定。例如：比如说外部函数costfunc就是高级命令optim的一个参数，它的调用次序一定要如此：f,g,ind=costfunc(x,ind)，然后高级最优化命令optim按以下方式调用：optim(costfunc,.)这里costfunc计算函数分f（x）的值和它在x点处的梯度（ind是一个整数，它的用处是使结果更加精确化）如果外部对象还需要其它值，这些变量可以在它的环境下定义，当然，他们也可以放在一个列表中。比如说，外部函数f,g,ind=costfunc(x,ind,a,b,c)。该函数对optim是有效的，不过得保证这个外部对象是一个列表list(costfunc,a,b,c)，并且对optim调用格式如下：optim(list(costfunc,a1,b1,c1),.一个外部对象也可能是Fortran或者C程序，这对于提高运算速度有帮助。这些程序的名字以字符串的形式传给高级命令。程序的调用次序也是规定的，这方面的例子在routines/default 目录之下。外部Fortran或者C程序也能够被动态的调用。extraction - 矩阵或列表元素的提取调用格式x(i,j) x(i) .=l(i) .=l(k1).(kn)(i) or .=l(list(k1,.,kn,i) l(k1).(kn)(i,j) or l(list(k1,.,kn,list(i,j) 参数：x : 各种类型的矩阵 l : 列表变量i,j : 索引k1,.kn : 索引说明：矩阵：i和j可以实数标量，向量，或者正元素矩阵* r=x(i,j)命令建立一个矩阵r，其中元素r(l,k)=x(int(i(l),int(j(k)，其中int(i(l)是取整命令，即取i(l)的整数部分。l的取值从1到向量i的总个数size(i,*)，k的取值从1到向量j的总个数size(j,*)，注意i和j的最大值必须小于或者等于矩阵x的行数（size(x,1)）与列数（size(x,2)）. 例子：x=1.23 4.6 6.7 9.5;6.8 6.21 5.14 1.47;.32 5.1 4.9 8.32;5.6 .47 1.56 7.56 i=1:3 j=2:4r=x(i,j)* 如果x是一个1x1矩阵，那么r=x(i)命令建立这样一个矩阵r，其中i必须是一个元素为1的矩阵，而r(l,k)=x(int(i(l),int(i(k)，l的取值从1到矩阵i的行数（size(i,1) ），k的取值从1到矩阵i的列数（size(x,2)）例子：x=5.6i=1 1 1;1 1 1 x(i)* 如果x是一个行向量，那么r=x(i)命令建立这样一个向量r，其中i也是一个向量（如i=1:3），r(l)=x(int(i(l) ，l的取值从1到向量i的总个数size(i,*)， i的最大值必须小于或者等于向量x的总个数size(x,*) . x=1.23 4.6 6.7 9.5 6.8 6.21 5.14 1.47 i=1:3r=x(i)* 如果x是一个多行多列的矩阵，那么r=x(i)命令建立这样一个列向量r，i也是一个向量（如i=1:7），r(l) (l的取值从1到向量i的总个数size(i,*)，至于如何赋值，请参见例子i的最大值必须小于或者等于向量x的总个数size(x,*) .x=1.23 4.6 6.7 9.5;6.8 6.21 5.14 1.47;.32 5.1 4.9 8.32;5.6 .47 1.56 7.56 i=1:7r=x(i)*（: ）冒号代表所有的元素* r=x(i,:) 建立一个矩阵 r ，其中r(l,k)=x(int(i(l),k)，l的取值从1到向量i的总个数size(i,*)，k 的取值从1到矩阵x的总列数size(x,2)x=1.23 4.6 6.7 9.5;6.8 6.21 5.14 1.47;.32 5.1 4.9 8.32;5.6 .47 1.56 7.56 i=1:3 r=x(i,:)* r=x(:,j) 建立一个矩阵 r，其中 r(l,k)=x(l,int(j(k)，l的取值从1到矩阵x的总行数size(x,1)，k的取值从1到向量j的总个数size(j,*) . x=1.23 4.6 6.7 9.5;6.8 6.21 5.14 1.47;.32 5.1 4.9 8.32;5.6 .47 1.56 7.56 j=1:3 r=x(:,j)*r=x(:)建立一个列向量，这个列向量是按从左到右的顺序把每一列串连起来组成的，它等效于matrix(x,size(x,*),1) .-布尔型向量. 如果索引变量(i or j)是布尔型向量，那么scilab执行find(i)或者find(j)命令，找出其中的为真的序号，且看以下例子：x=1.23 4.6 6.7 9.5;6.8 6.21 5.14 1.47;.32 5.1 4.9 8.32;5.6 .47 1.56 7.56 i=%t %t;%f %f;%t %t;%f %f;%t %t;%f %f find(i)x(i)-多项式. 如果索引变量(i or j)是多项式向量或者隐形多项式向量，那么scilab执行horner(i,m)或者horner(j,n)命令，其中m和n与x的维数有关，为了增加可读性，我们建议使用$符号。例子/ MATRIX CASEa=1 2 3;4 5 6a(1,2)a(1 1,2)a(:,1)a(:,3:-1:1)a(1)a(6)a(:)a(%t %f %f %t)a(%t %f,2 3)a(1:2,$-1)a($:-1:1,2)a($)/x=testx(1 1;1 1;1 1)/b=1/%s,(%s+1)/(%s-1)b(1,1)b(1,$)b(2) / the numerator/ LIST OR TLIST CASEl=list(1,qwerw,%s)l(1)a,b=l(3 2)l($)x=tlist(l(2:3) /form a tlist with the last 2 components of l/dts=list(1,tlist(x;a;b,10,2 3);dts(2)(a)dts(2)(b)(1,2)a,b=dts(2)(a,b)feval - 多次求值调用格式z=feval(x,y,f) z=feval(x,f)参数x,y : 两个向量f : 汗水或者字符串(对Fortran或者C调用)一个向量型参量或者两个向量型参量的函数多次求值z=feval(x,f)返回向量z,z由z(i)=f(x(i) 决定z=feval(x,y,f)返回向量z,其中z(i,j)=f(x(i),y(j)f是一个外部对象(函数或者程序),它接受一个或者两个实型参数,f返回的是实数或者复数.在Fortran调用的情况下,函数f必须在子程序ffeval.f中定义(在文件夹SCIDIR/routines/default中)例子:deff(z=f(x,y),z=x2+y2);feval(1:10,1:5,f)deff(z=f(x,y),z=x+%i*y);feval(1:10,1:5,f)feval(1:10,1:5,parab) feval(1:10,parab)find - 查找布尔型向量或者矩阵的元素索引调用格式ii=find(x ,nmax) i1,i2,.=find(x ,nmax) 参数x : 可能是一个布尔型向量, 布尔型矩阵, 布尔型超级矩阵, 标准矩阵或者超级矩阵.nmax : 返回一个最多的索引个数. 默认值为-1, 代表搜索所有元素.这个选项对于避免搜索所有索引很有效. ii, i1, i2, . : 整型索引向量或者空矩阵.说明:如果x 是一个布尔型矩阵ii=find(x)返回值为真的元素索引,如果没有元素值为真,则返回一个空矩阵.i1,i2,.=find(x)返回索引向量i1(行)与 i2(列),如果没有符合条件的元素,则返回一个空矩阵.如果x 是一个标准矩阵或者超级矩阵,find(x)则等效于find(x0)find()返回 例子:ExamplesA=rand(1,20);w=find(A100)B=rand(1,20);w=find(B0.9) H(i(1),j(1),k(1)format - 数字打印与显示格式调用格式：format(type,long) format() 参数：type : 字符串long : 整数 (最大显示位数，默认为10) 说明：参数type设置当前打印格式，它是以下参数之一v: 变量格式 (默认) e: 科学计数格式long：定义显示的最大显示位数，默认为10.format()返回一个包含当前格式信息的向量，第一个分量是格式的种类（0表示v，1表示e），第二个分量是显示位数。例子：Examplesx=rand(1,5);format(v,10);xformat(20);xformat(e,10);xformat(20);xx=100 %eps;format(e,10);xformat(v,10);xformat(v)fort - Fortran or C语言拥护程序调用 （暂缺）for - 循环语句的关键词说明：for用来定义循环语句.它的语法是: for variable=expression ,instruction, . ,instruction,endfor variable=expression do instruction, ,instruction,end如果expression是矩阵或者一个行向量，那么variable依次取矩阵每一列的值有用的例子：for variable=n1:step:n2, .,end如果expression是一个列表变量，那么variable依次取列表中的元素。例子：n=5;for i = 1:n, for j = 1:n, a(i,j) = 1/(i+j-1);end;endfor j = 2:n-1, a(j,j) = j; end; afor e=eye(3,3),e,end for v=a, write(6,v),end for j=1:n,v=a(:,j), write(6,v),end for l=list(1,2,example); l,endfunptr - coding of primitives ( wizard stuff ) （暂缺，高级用户使用）getdate - 获得日期与时间信息调用格式dt=getdate() x=getdate(s) dt=getdate(x) 参数dt : 包含10个整型元素的向量 x : 一个整数，表示自1970年1月1日到某一天的秒数说明：dt=getdate() 按以下格式返回当前的时间信息dt(1): 年份从 0000 到 9999. dt(2): 月份从 01 到 12. dt(3): 按ISO 8601标准计算的周数从01 到 53. dt(4): 按古罗马日历计算的天数字从001 到 366. dt(5): 周数，1代表星期天. dt(6): 日数从 01 到 31. dt(7): 时数从 00 到 23. dt(8): 分数从 00 and 59. dt(9): 秒数从 00 and 59. dt(10): 毫秒数 从 000 and 999. x=getdate(s) 返回一个整数，它表示自1970年1月1日00:00 UTC(Unix Time Convention)以来所经过的秒数.dt=getdate(x) 结合上面一句，运行一下这个命令就知道什么意思了：）例子：w=getdate()mprintf(Year:%d,Month:%d,Day:%d,w(1),w(2),w(6);x=getdate(s)getdate(x)getenv - 获得环境变量的值调用格式：env=getenv(str , rep ) 参数：str : 字符串，用来指定环境变量的值rep : 可选字符串. 当使用了这个选项，函数getenv没有找到str所指定环境变量， 则返回rep的值env : 字符串，包含了环境变量的值说明：返回已存在的环境变量的值例子：getenv(SCI)getenv(FOO,foo)getfield - 列表域的提取调用格式x,.=getfield(i,l)参数x : 任何类型的矩阵l : list, tlist 或者mlist 变量i : 域索引, 参考extraction获得更多细节. 说明：这个函数等效于x,.=l(i)，都是用来提取表中的域，两者唯一不同是前者可以用于mlist对象.getpid -获得Scilab进程号调用格式id=getpid() 说明返回一个表示Scilab进程号的整数例子：d=SD_+string(getpid()+_getversion - 获得Scilab版本信息调用格式 version=getversion() version,opts=getversi