最近我想实现的巷道评估ISO分类。 请参阅本文中使用车辆的加速度测量来估算路面不平度 。 但是,我碰到了大约绘制ISO分类图的一些问题。 的代码如下所示:
%% generate ISO Clasification Curve (m^2/(rad/m))
AngSpaFre = 10^(-3):0.01:10^2; % spatial frequency (rad/m)
ParamTable = [0 1*10^(-6) 2*(10^(-6)); % ISO classification parameters
2*(10^(-6)) 4*(10^(-6)) 8*(10^(-6)); % (lower bound, geometric average, upper bound)
8*(10^(-6)) 16*(10^(-6)) 32*(10^(-6));
32*(10^(-6)) 64*(10^(-6)) 128*(10^(-6));
128*(10^(-6)) 256*(10^(-6)) 512*(10^(-6));
512*(10^(-6)) 1024*(10^(-6)) 2048*(10^(-6));
2048*(10^(-6)) 4096*(10^(-6)) 8192*(10^(-6));
8192*(10^(-6)) 16384*(10^(-6)) 32768*(10^(-6))];
class_num = 8;
len = length(AngSpaFre);
iso_table = zeros(class_num, len);
for i = 1:class_num
geo_mean = ParamTable(i, 2);
geo_upbound = ParamTable(i, 3);
for j = 1:len
if(AngSpaFre(j) <= 1)
iso_table(i,j) = geo_mean*(AngSpaFre(j)^(-2));
else
iso_table(i,j) = geo_mean*(AngSpaFre(j)^(-1.5));
end
end
end
figure,
ht = loglog(AngSpaFre, iso_table(1,:), AngSpaFre, iso_table(2,:), AngSpaFre, iso_table(3,:), AngSpaFre, iso_table(4,:), AngSpaFre, iso_table(5,:), AngSpaFre, iso_table(6,:), AngSpaFre, iso_table(7,:), AngSpaFre, iso_table(8,:));
hY = get(gca,'ylim');
hold on;
loglog([1 1], hY);
xlabel('Spatial Frequency \Omega (rad/m)');
ylabel('PSD (m^2/(rad/m))');
title('ISO Classification (log-log scale)');
而目前的结果显示为图波纹管:
然后2个问题来自于我的实现,并希望任何一个可以给我一些解决方案或简单的例子来使问题清晰。
Q1:由于分类曲线应该认识到道路的质量水平。 所以我想添加的每个线以下的标签,以显示像下面的图中的水平。 我怎么能认识到这一点在MATLAB?
Q2:另外,我想在“双对数”情节的右侧添加其他y轴没有像下图数刻度值。 如何实现这一目标? 我指的是官方网站的例子很多,但大多数集中在“阴谋” /“plotyy”的例子。 :(