Потоковые ленты иллюстрируют направление потока, похожего на линии потоков, но могут также показать вращение вокруг оси потока путем скручивания соединительной линии ленточной формы. 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;
