Набор данных вектора MATLAB® под названием wind
представляет воздушные потоки по Северной Америке. Этот пример использует комбинацию методов:
Линии потоков, чтобы проследить скорость ветра
Плоскости разбиения, чтобы показать перекрестные виды в сечении данных
Контуры в плоскостях разбиения, чтобы улучшить видимость окраски плоскости разбиения
Загрузите данные и определите минимальные и максимальные значения, чтобы определить местоположение плоскостей разбиения и контурных графиков (load
min
Max
).
load wind xmin = min(x(:)); xmax = max(x(:)); ymax = max(y(:)); zmin = min(z(:));
Вычислите величину векторного поля (который представляет скорость ветра) сгенерировать скалярные данные для slice
команда. Создайте плоскости разбиения вдоль оси X в xmin
, 100, и
xmax
, вдоль оси Y в ymax
, и вдоль оси z в zmin
. Задайте интерполированную поверхность, окрашивающую, таким образом, срез, окрашивающий, указывает на скорость ветра, и не чертите ребра (sqrt
срез
'FaceColor'
'EdgeColor'
).
wind_speed = sqrt(u.^2 + v.^2 + w.^2); hsurfaces = slice(x,y,z,wind_speed,[xmin,100,xmax],ymax,zmin); set(hsurfaces,'FaceColor','interp','EdgeColor','none') colormap jet
Чертите светло-серые линии контура в плоскостях разбиения, чтобы помочь определить количество отображения цвета (contourslice
'EdgeColor'
'LineWidth'
).
hcont = ... contourslice(x,y,z,wind_speed,[xmin,100,xmax],ymax,zmin); set(hcont,'EdgeColor',[0.7 0.7 0.7],'LineWidth',0.5)
В этом примере все линии потоков запускаются в значении оси X 80 и порождают линейную оболочку столбцов 20 - 50 в направлении Y и от 0 до 15 в z-направлении. Сохраните указатели линий потоков и установите ширину линии и цвет (meshgrid
поток
'LineWidth'
, Color
).
[sx,sy,sz] = meshgrid(80,20:10:50,0:5:15); hlines = streamline(x,y,z,u,v,w,sx,sy,sz); set(hlines,'LineWidth',2,'Color','r')
Настройте представление, расширив ось z, чтобы облегчить читать график (view
daspect
ось
).
view(3) daspect([2,2,1]) axis tight
Смотрите coneplot
для примера тех же данных, отображенных на графике с конусами.