我们在做数据绘图或函数图像时经常需要使用对数坐标系。尤其是数据的范围跨越很多个数量级时,通常的线性坐标系下无法表现出数据特征。
Scilab 中Plot函数无法画出对数坐标。需要使用 plot2d 函数。
plot2d 函数的基本用法如下:
plot2d([logflag,][x,],y[,style[,strf[,leg[,rect[,nax]]]]])
plot2d([logflag,][x,],y,<opt_args>)
下面是一个简单的例子:
~~~
iter = linspace(0,10,11);
err = 10.^(-iter);
plot2d("nl", iter, err, style=2);
set(gca(),"grid",[1 1]);
~~~
![](https://box.kancloud.cn/2016-01-24_56a4234437d63.jpg)
这个例子如果在普通的坐标系下看,是这个样子的:
~~~
plot(iter,err);
set(gca(),"grid",[1 1]);
~~~
![](https://box.kancloud.cn/2016-01-24_56a42344461ed.jpg)
由于数据很快就很接近0了,在图中很难看出后面的趋势。
下面来详细的讲解一下plot2d函数。
plot2d("nl", iter, err, style=2);
“nl” 表示,横坐标为正常的模式(normal),纵坐标为对数(log).
Style = 2 表示的是曲线的颜色。2 表示的是colormap 中的第二项,也就是蓝色。
如果是负数,则表示用不同的线型。如果既要设置曲线的颜色,又要设置线型,那么。。。暂时还没搞定。
下面再给一个例子,通过rect参数控制xy的范围:
~~~
x=[0:0.1:2*%pi]';
plot2d(x,[sin(x) sin(2*x) sin(3*x)],rect=[0,0,6,0.5]);
~~~
![](https://box.kancloud.cn/2016-01-24_56a4234454376.jpg)
有点跑题了,接着介绍对数坐标系绘图。下面再给一个例子:
~~~
ind = linspace(0,6,7);
iter = 10.^ind;
err1 = 10.^(-ind);
err2 = (10.^(-ind)).^2;
xset("font size", 4);
plot2d("ll", iter, err1, style=2);
plot2d("ll", iter, err2, style=3);
title("Loglog","fontsize", 4);
xlabel("Iterations","fontsize", 4);
ylabel("Error","fontsize", 4);
set(gca(),"grid",[5 5]);
legend(['error1';'error2'],"in_lower_left");
~~~
![](https://box.kancloud.cn/2016-01-24_56a42344647a8.jpg)
这个图是双对数坐标,同时还调整了图上的文字。需要注意的是:
xset("font size", 4);
语句一定要在
legend(['error1';'error2'],"in_lower_left");
语句之前调用,否则得到的图形的legend 的字号会有问题,下面是个实验,先执行如下语句:
~~~
ind = linspace(0,6,7);
iter = 10.^ind;
err1 = 10.^(-ind);
err2 = (10.^(-ind)).^2;
plot2d("ll", iter, err1, style=2);
plot2d("ll", iter, err2, style=3);
title("Loglog","fontsize", 4);
xlabel("Iterations","fontsize", 4);
ylabel("Error","fontsize", 4);
set(gca(),"grid",[5 5]);
legend(['error1';'error2'],"in_lower_left");
~~~
![](https://box.kancloud.cn/2016-01-24_56a4234476f85.jpg)
这个图形是和我们的预期一样的。标题和Label的字号变大了,刻度和Legend的字号还是原来的大小。
接着执行:
~~~
xset("font size", 4);
~~~
![](https://box.kancloud.cn/2016-01-24_56a42344844ca.jpg)
结果是刻度上的字号更新为正确的大小了,可legend 的字号没变。看来这个是 scilab 的一个bug。
所以我们需要先设置字号,然后调用legend 函数。
- 前言
- gnuplot 中的数据平滑
- 利用 Maxima 求解常微分方程
- Maxima 的绘图功能 1
- Maxima 的绘图功能 2
- Maxima 的绘图功能 3
- Maxima 矩阵及矢量运算 1
- Maxima 矩阵及矢量运算 3
- Maxima 矩阵及矢量运算 4
- maxima 学习笔记(coeff 函数的应用)
- maxima 代数表达式变换
- maxima 学习笔记(基本微积分)
- maxima 代数方程求解
- Maxima 中的复数运算
- Maxima 的三角函数化简功能
- Maxima 的基本微积分操作
- 利用 gnuplot 绘制时间序列图
- gnuplot 入门教程 1
- gnuplot 入门教程 2
- gnuplot 入门教程 3
- gnuplot 入门教程 4
- gnuplot 读取逗号分隔的数据文件
- scilab 读取处理 wav 文件
- scilab 读取处理 wav 文件 (2)
- Scilab 处理声音数据(补充)
- 利用 SCILAB 设计 FIR 滤波器(窗函数法)
- 利用 SCILAB 设计 FIR 滤波器(Minimax法)
- 利用 SCILAB 设计 iir 滤波器设计(模拟滤波器双线性变换法)
- Scilab 的绘图函数(1)
- Scilab 的绘图函数(2)
- Scilab 的绘图函数(3)
- Scilab 的绘图函数(4)