Потоковые ленты иллюстрируют направление потока, похожего на линии потоков, но могут также показать вращение вокруг оси потока путем скручивания соединительной линии ленточной формы. streamribbon
функция позволяет вам задать угол скручивания (в радианах) для каждой вершины в потоковых лентах.
Когда используется в сочетании с curl
функция, streamribbon
полезно для отображения вихревой скорости вращения векторного поля. Следующий пример иллюстрирует этот метод.
Загрузите и выберите необходимую область в wind
использование набора данных subvolume
. Графический вывод полного набора данных сначала может помочь вам выбрать необходимую область.
load wind lims = [100.64 116.67 17.25 28.75 -0.02 6.86]; [x,y,z,u,v,w] = subvolume(x,y,z,u,v,w,lims);
Вычислите вихревую скорость вращения и скорость ветра.
cav = curl(x,y,z,u,v,w); wind_speed = sqrt(u.^2 + v.^2 + w.^2);
Использование meshgrid
создать массивы начальных точек для потоковых лент. Смотрите Начальные точки Определения для Потока информация о Plotsfor об определении массивов начальных точек.
stream3
вычисляет вершины линии потоков с размером шага .5
.
streamribbon
масштабирует ширину ленты на коэффициент 2
улучшать видимость скручивания (который указывает на вихревую скорость вращения).
streamribbon
возвращает указатели объектов подложки, которые это создает, которые затем используются, чтобы выбрать цвет к красному (FaceColor
), цвет поверхностных ребер к светло-серому (EdgeColor
), и немного увеличьте яркость рассеянного света, отраженного, когда подсветка будет применена (AmbientStrength
).
[sx sy sz] = meshgrid(110,20:5:30,1:5); verts = stream3(x,y,z,u,v,w,sx,sy,sz,.5); h = streamribbon(verts,x,y,z,cav,wind_speed,2); set(h,'FaceColor','r',... 'EdgeColor',[.7 .7 .7],... 'AmbientStrength',.6)
volumebounds
команда обеспечивает удобный способ установить axis
и цветные пределы.
Добавьте a grid
и набор view
для 3-D (streamribbon
не изменяет текущее представление).
camlight
создает свет, расположенный справа от точки зрения и lighting
устанавливает метод подсветки на Гуро.
axis(volumebounds(x,y,z,wind_speed)) grid on view(3) camlight right;