matlab与scilab的比较2.docVIP

matlab 使用老出错，决定试试scilab. 几个例子来学习 scilab例1： 画等值线的命令//定义一个向量x=[ 1 2 3 ; 1 3 5 ; 8 2 8];//和matlab 相同的命令a=1:3;b=1:3;nz=2;//画等值线图，nz 表示等值线的取值点数。contour2d(a,b,x,nz)上面的命令可以存在一个文件里，例如 a.sci, 然后在scilab的命令行里可以用exec(a.sci)来执行，就会看到画的图如下 例2： matlab的reshape命令，在scilab下面是//定义一个向量--x=[1 2 3 4 5 6 6 7 7 8]x = 1. 2. 3. 4. 5. 6. 6. 7. 7. 8.//向量的维数--size(x)ans = 1. 10.//定义个矩阵, 相当于matlab的 y=reshape(x,2,5)--y=matrix(x,2,5)y = 1. 3. 5. 6. 7. 2. 4. 6. 7. 8.例3：文件，输入输出// 新建一个文件，向文件中输入随机数u=file(open,aa.txt,unknown)for k=1:4// 生成 1-by-4 的随机矩阵a=rand(1,4)write(u,a)end//rewind 使得文件的指针回到文件头。file(rewind,u)//读出2-by-4的矩阵x=read(u,2,4)//关闭文件file(close,u)例4：和waterfall 和像的画图函数x 如果是一个矩阵，则mesh(x)会的得到一个图像，如下图 例5： 把生成的图像转化成eps格式保存//设置标题xtitle(Nk=1024, k=431:592)//存成eps格式的文件，文件名将会是f1024zoom.epsxs2eps(0,f1024zoom) 例6: 打开和写入文件可以使用 C 程序的格式p=mopen(s2xx1024.sci,r);Nx=81;Nk=1024;d=512; d1=Nk/2-d+1; d2=Nk/2+d;//文件保存的是Nx 行 Nk 列的矩阵//下面的f读入所有的数据, %lg 表示读入 double 精度f=mfscanf(-1,p,%lg);// 把上面读入的一维的向量转成矩阵，要注意行列参数的顺序， 转成 Nk 行 Nx 列的矩阵f=matrix(f,Nk,Nx);//再次转置才得到我们想要的 Nx 行 Nk 列的矩阵f=f;mclose(p);p1=mopen(k1.txt,w);k1=1:Nk;k=(k1-Nk/2-0.5)/Nk;kk=k(d1:d2);b = kk;pos = 1;for posy=f(pos,1:Nk);yy=y(d1:d2);mv = mean(yy);b = [b yy];mfprintf(p1,%16.15e %16.15e\n,b)file(close,p1);例7: 找一个好的spline 曲线来拟合数据// this is an artifical example where the datas xd and yd// are build from a perturbed sin functiona = 0; b = 2*%pi;sigma = 0.1; // standard deviation of the gaussian noisem = 200; // number of experimental points//下面的函数给出[a,b]上200个等分点，包括两个端点。xd = linspace(a,b,m);yd = sin(xd) + grand(xd,nor,0,sigma);//我们要构造一个spline曲线，使得曲线在点//x(i) , i=1,2,...，6的函数值等于 y(i),//曲线在点 x(i) 的导数值等于 d(i)n = 6; // number of breakpointsx = linspace(a,b,n);// compute the spline//lsq_splin 函数就是计算这样的spline 曲线，记住 spline 曲线和// {y(i), d(i)， = 1,2,。。。，6} 有这一一对应的关系，lsq_splin//计算出来的一条是原来给出