ВУЗ:
Составители:
Рубрика:
7) Построить пластину:
> g1: = line([0,0,0], [0,-2*r,0],
color = red, thickness = 2):
8) Получить картину обтекания:
> g2: = complexplot3d(Im(chi), zeta = -4-4*I..4 + 0*I,
style = contour, orientation = [-90,0],
axes = boxed,color = blue):
9) Отобразить совмещенные графические данные:
> display(g1,g2);
10) Выделить компоненту скорости вдоль оси x:
> ux: = Re(diff(chi,zeta));
:= ux +
1
2
1
2
ℜ − −
ζ
+ ζ
2
4
2
−
1
2
ζ
2 + ζ
2
4
−
ζ
2
+ ζ
2
4
2
2
10) Выделить компоненту скорости вдоль оси y:
> uy: = -Im(diff(chi,zeta));
:= uy −
1
2
ℑ − −
ζ
+ ζ
2
4
2
−
1
2
ζ
2 + ζ
2
4
−
ζ
2
+ ζ
2
4
2
2
11) Оформить выражение для расчета компонент скорости потока в точках с координатами (x, y) в виде
функций.
> uxf: = (x,y)->subs(zeta = x + I*y,ux);
:= uxf
→ (),
x
y
()subs ,
=
ζ
+
x
y
I
ux
> uyf: = (x,y)->subs(zeta = x + I*y,uy);
:= uyf
→ (),
x
y
()subs ,
=
ζ
+
x
y
I
uy
12) Рассчитать значения проекций скорости на координатные оси в конкретной точке.
> evalf(uxf(-r,-r));
0.5688644811
> evalf(uyf(-r,-r));
у
х
