MATLAB对数刻度绘图 （不等间距坐标绘图）
主要概括
这两天师兄算了一个风速，有了采样时间和风速大小，让我帮忙绘制风速的功率密度谱函数，到最后看到不等比的坐标轴怎么画，最后搜出了，写在这边以供查找。
原图
原图如下所示，主要是横纵坐标为10的对数：


方法
semilogx(x,y)  绘图，x轴对数刻度，y轴线性刻度，

semilogy(x,y)  绘图，x轴线性刻度，y轴对数刻度，

loglog(x,y)       绘图，x,y轴都为对数刻度
结果

本实验取材于中南大学《MATLAB与科学计算》

今天学习主要有两种内容，第一、对数坐标图，第二、统计图

一、对数坐标图
semilogx(x1,y1,选项1,x2,y2,选项2,...)
semilogy(x1,y1,选项1,x2,y2,选项2,....)
loglog(x1,y1,选项1,x2,y2,选项2,...)

例子：绘出1/x的函数图像



代码
x=0:0.1:10;
y=1./x;
subplot(2,2,1);
plot(x,y)
title('plot(x,y)');
subplot(2,2,2);
semilogx(x,y)
title('semilogx(x,y)');
grid on
subplot(2,2,3);
semilogy(x,y);
title('semilogy(x,y)');
grid on
subplot(2,2,4);
loglog(x,y)
title('loglog(x,y)');
grid on


总结：其实就是对数的转换，你想转换x那就semilogx，如果想转换y那就semilogy，如果想都转换一样，那就是loglog
二、统计图

1、条形图

格式
bar(y,style)

ps:style可以有两类，grouped簇状分组 stacked 为堆积分组
代码：
hold off
x=[2015,2016,2017];
y=[68,80,115,98,102;75,88,102,99,110;81,86,125,105,115];
bar(x,y)
title('Group')

图片



2、直方图

格式
hist(x,y)

代码
>> y=randn(500,1);
>> subplot(2,1,1);
>> hist(y);
>> title('高斯分布直方图');
>> subplot(2,1,2);
>> x=-3:0.2:3;
>> hist(y,x);
>> title('指定区间中心点的直方图')
>> 

图片



3、饼图

格式：
pie(x,y)

代码
score=[5,17,23,9,4];
ex=[0,0,0,0,1];
hold
pie(score,ex)
legend('优秀','良好','中等','及格','不及格','location','eastoutside')


图片

close all;
x = [256 512 1024 2048 4096 8192];
y1 = [0.017360, 0.043902, 0.185930, 0.787196, 3.266647, 45.045046]
y2 = [0.000510, 0.000625, 0.000573, 0.000602, 0.000635, 0.000714]
semilogy(x,y1,'b*-',x,y2,'ro-','linewidth',2);
hold on
% semilogy(x,y2,'r.-','linewidth',2);
grid on;
legend('SSDFT','New approcah')


title('Initial Step');           %设置标题行
xlabel('Image Dimension (pixels)');     %设置横坐标
ylabel('Time (s)');   %设置纵坐标

PROBLEM：
为什么loglog 、semilogy不能画出对数图呢？
KEY： hold on必须在绘图命令之后，而不能直接接在figure之后。否则在绘制特殊坐标的时候，仍然显示的是直角坐标。
错误代码：
xlab=0:10;
figure;hold on;grid on;
semilogy(xlab,10.^(xlab),'-s');
semilogy(xlab,10.^(xlab),'-s');
运行结果：

结果显然不是对数坐标。
解决过程：
调试。
在m文件中：先设置一个断点breakpoint，然后在命令窗中运行这个m文件，程序会运行到断点处，停住。

然后按F10，程序就执行一步，按一下走一步。

如果运行到某个函数，你想进入这个函数一步一步调试，那么按F11，就是单步调试，进入函数。当然这里不需要。可以直接看出来，hold on默认直角坐标并保持之。
正确代码：
xlab=0:10;
figure;
semilogy(xlab,10.^(xlab),'-s');
hold on;
semilogy(xlab,10.^(xlab),'-s');
grid on;