
- •Введение
- •1. Запуск математического пакета maple, структура окна
- •2. Описание основных команд–процедур пакета maple
- •2.1. Арифметические операции, типы данных в Maple
- •2.2. Структура команд Maple
- •Команда (параметр_1, параметр_2, …, параметр_n);
- •2.3. Преобразование математических выражений
- •2.3.1. Выделение частей арифметических выражений
- •2.3.2. Тождественные преобразования арифметических выражений
- •Collect(выражение, имя переменной относительно которой приводятся подобные);
- •3. Функции maple, решение уравнений и неравенств
- •3.1. Способы задания функций и замена переменных
- •Piecewise(cond_1,f1,cond_2,f2, …).
- •3.2. Операции оценивания
- •3.3. Решение уравнений и их систем
- •3.4. Решение неравенств
- •4. Графические и геометрические возможности пакета maple
- •4.1. Графические возможности пакета Maple
- •4.1.1. Двумерная графика
- •Polarplot ([переменная, выражение,параметры переменной],опции);
- •Pointplot(data_list,опции);
- •Sphereplot ((выражение), параметры_, параметры_, опции);
- •Sphereplot ([r_выражение, _выражение, _выражение], параметр_1, параметр_2, опции);
- •4.1.2. Трехмерная графика
- •4.2. Анимация
- •4.3. Геометрические пакеты Maple
- •4.3.1. Стереометрия
- •4.3.2. Планиметрия
- •5. Дифференциальное и интегральное исчисление функции одной переменной
- •5.1. Вычисление пределов в Maple для некоторых математических операций существует по две команды: прямого и отложенного исполнения. Имена команд состоят из одинаковых букв за исключением первой:
- •5.2. Дифференцирование
- •5.2.1. Вычисление производных
- •5.2.2. Дифференциальный оператор
- •5.3. Исследование функций
- •5.3.1. Непрерывность функции и точки разрыва
- •5.3.2. Экстремумы. Наибольшее и наименьшее значение функции
- •5.3.3. Исследование функции по общей схеме
- •{Arctan(0)}
- •6. Интегрирование
- •6.1. Аналитическое и численное интегрирование
- •6.2. Определенные интегралы
- •Is assumed to be: RealRange(Open(-1),3)
- •6.3. Основные методы интегрирования
- •7. Дифференциальные уравнения
- •7.1. Аналитическое решение дифференциальных уравнений
- •7.1.1. Общее решение дифференциальных уравнений
- •7.1.2. Фундаментальная (базисная) система решений
- •7.1.3. Решение задачи Коши или краевой задачи
- •7.1.4. Системы дифференциальных уравнений
- •7.1.5. Приближенное решение дифференциальных уравнений с помощью степенных рядов
- •7.2. Численное решение дифференциальных уравнений
- •7.2.1. Формат команд dsolve и odeplot для нахождения численного решения дифференциальных уравнений в графическом виде
- •Italic, 12]):
- •Italic, 12]):
- •7.2.2. Пакет графического представления решений дифференциальных уравнений dEtools
- •7.2.3. Построение фазовых портретов систем дифференциальных уравнений
- •8. Функции многих переменных, векторный анализ, ряды, интегральные преобразования
- •8.1. Дифференциальное исчисление функций многих переменных
- •8.1.1. Частные производные
- •8.1.2. Локальные и условные экстремумы функций многих переменных
- •Maximize(f,{x1,…,xn},range);
- •Minimize(f,{x1,…,xn}, range);
- •8.2. Интегральное исчисление функций многих переменных
- •8.3. Векторный анализ
- •8.4. Ряды и произведения
- •8.4.1. Вычисление суммы ряда и произведений
- •8.4.2. Разложение функции в степенной ряд и ряд Тейлора
- •8.5. Интегральные преобразования
- •8.5.1. Преобразование Фурье
- •8.5.2. Преобразование Лапласа
- •9. Линейная алгебра
- •9.1. Векторная алгебра
- •Warning, the protected names norm and trace have been redefined and unprotected
- •9.2. Действия с матрицами
- •9.2.1. Определение матрицы
- •9.2.2. Арифметические операции с матрицами
- •9.2.3. Определители, миноры и алгебраические дополнения, ранг и след матрицы
- •9.2.4. Обратная и транспонированная матрицы
- •Evalm(1/a);
- •Inverse(a).
- •9.2.5. Определение типа матрицы
- •1*Sqrt(3)/2,-1/2]]);
- •9.2.6. Функции от матриц
- •9.3. Спектральный анализ матрицы
- •9.4. Системы линейных уравнений, матричные уравнения
- •For имя переменной цикла in выражение 1 do выражение 2 od;
- •10.3. Процедуры-функции
- •10.4. Процедуры
- •Writeto (“имя файла”)
- •Appendto (“имя файла”)
- •Print(список Maple-выражений, перечисляемых через запятую);
- •Lprint(список Maple-выражений, перечисляемых через запятую);
- •Readdata(“имя файла”, тип переменной(integer/float – последний тип устанавливается по умолчанию),счетчик чисел);
- •Readlib(readdata):
- •12. Использование математического пакета maple для научных исследований
- •12.1. Исследование влияния изменяемых параметров плоской помольной камеры мельницы противоточного действия на скорость энергоносителя
- •12.1.1. Постановка задачи
- •12.1.2. Решение задачи
- •12.2. Определение коэффициента полезного заполнения материалом деформируемых бочкообразных камер устройств для тонкого и свертонкого помола материалов
- •12.2.1. Расчетная схема
- •12.2.2. Решение задачи
- •12.3. Исследование влияния основных конструктивных и технологических параметров на режимы работы трубной мельницы
- •Formula(1.8)
- •Formula(1.9)
- •Formula(2.1)
- •Formula(2.2)
- •Formula(2.3)
- •Formula(3.1)
- •Formula(3.2)
- •Formula(3.3)
- •13. Формат наиболее используемых команд аналитического пакета maple
- •Заключение
- •Библиографический список
- •Оглавление
Polarplot ([переменная, выражение,параметры переменной],опции);
Пример.
Построить график функции sin(3*r) при r=0..7 (рис. 4.6).
[> restart; with(plots):
Warning, the name changecoords has been redefined
[>polarplot([r, sin(3*r),r=0..7],color=black,
thickness=3);
Рис. 4.6. График функции, построенный в полярной системе координат
В графическом пакете plot есть команда для построения графика по координатам точек, для этого нужно ввести координаты точек командой
data_list:=[[x1,y1],…,[xn,yn]];
а затем выполнить команду
Pointplot(data_list,опции);
Пример.
Построить график по координатам точек [2,4],[-4,4],[-2,3],[-1,4], [0,5] (рис. 4.7).
[> restart;
[> with(plots):
Warning, the name changecoords has been redefined
[> data_list:=[[2,4],[-4,4],[-2,3],[-1,4], [0,5]];
[>pointplot(data_list,style=line,color=black, thickness=3);
Рис. 4.7. График, построенный по координатам точек
Формат команды построения графиков в координат:
Sphereplot ((выражение), параметры_, параметры_, опции);
или
Sphereplot ([r_выражение, _выражение, _выражение], параметр_1, параметр_2, опции);
Пример.
Построить график функции в сферической системе координат (4/3)^theta*sin(phi), при theta=-1..2*Pi, phi=0..2*Pi (рис. 4.8).
[> restart;
[> with(plots):
Warning, the name changecoords has been redefined
[> sphereplot ((4/3)^theta*sin(phi),theta=-1..2*Pi, phi=0..2*Pi);
Рис. 4.8. График функции, построенный в сферической системе координат
4.1.2. Трехмерная графика
Команды трехмерной графики аналогичны командам двумерной графики и имеют окончания “3d”. Причем число параметров, как правило, больше на единицу, а точки определяются тремя координатами. Общий вид команды:
plot3d ({выражение_1, выражение_2,…}, переменная_1=а..b, переменная_2=c..d, опции);
Опции этой команды совпадают с параметрами построения двумерной графики. Кроме того, к специфическим опциям пакета plot3d следует отнести:
light=[x1,y1,c1,c2,c3] – задание подсветки поверхности, создаваемой источником света из точки со сферическими координатами x1, y1, цвет которой задается долями красного (с1), зеленого (с2), синего (с3) цветов из интервала [0,1];
style – задает стиль рисунка;
POINT – точки;
LINE – линии;
HIDDEN – сетка с удалением невидимых линий;
PATCH – заполнитель;
WIREАFRAME – сетка с выводом невидимых линий;
CONTOUR – линии уровня;
PATCHCONTOUR – заполнитель и линии уровня.
Параметр shading – задает опцию интенсивности заполнения, NONE – без раскраски.
Данная команда позволяет выводить на одном рисунке несколько поверхностей, задаваемые однотипными Maple выражениями, которые зависят от двух переменных переменная_1[а,b], переменная _2[c,d]. Вид рисунка можно менять при помощи опций.
Пример.
Построить поверхности функций xsin(2y)+ycos(3x), sqrt(x^2+y^2)-10, при x=-Pi..Pi, y=-Pi/2..Pi/2 (рис. 4.9).
[> restart;
[>plot3d({x*sin(2*y)+y*cos(3*x), sqrt(x^2+y^2)-10}, x=-Pi..Pi, y=-Pi/2..Pi/2, grid=[20,20],title="2 poverhnosti", axes=FRAMED, orientation=[20,60], color=x+y);
Рис. 4.9. Поверхности трехмерных функций
Если поверхность задана параметрически: x=x(u,v), y=y(u,v), z=z(u,v), построение графика задается командой
plot3d([x(u,v), y(u,v), z(u,v)],u=u1..u2, v=v1..v2);
В случае если поверхность задана неявно F(x,y,z)=0, то построение поверхности осуществляется с помощью следующей команды:
implicitplot3d(F(x,y,z)=0, x=x1..x2, y=y1..y2, z=z1..z2);
Кроме того, в пакете plot имеется команда spacecurve для построения пространственной кривой, заданной параметрически. Даная команда имеет следующую структуру:
spacecurve([x(t),y(t),z(t)],t=t1..t2);
Пример.
Построить поверхности функций, заданных параметрически x^2-y^2+z^2=4 при x=-10..10, y=-10..10, z=-10..10 (рис. 4.10, а).
[> restart;
[> with(plots):
Warning, the name changecoords has been redefined
[> implicitplot3d(x^2-y^2+z^2=4, x=-10..10, y=-10..10, z=-10..10, scaling=CONSTRAINED);
Если уменьшить границы просмотра графика ( с x=-10..10, y=-10..10, z=-10..10 до x=-2..2, y=-2..2, z=-2..2), то это приведет к его сильному искажению (рис. 4.10, б).
[> restart;
[> with(plots):
Warning, the name changecoords has been redefined
[> implicitplot3d(x^2-y^2+z^2=4, x=-2..2, y=-2..2, z=-2..2, scaling=CONSTRAINED);
а)
б)
Рис. 4.10. Поверхность функции, заданной параметрически
(а –x=-10..10, y=-10..10, z=-10..10;
б – x=-2..2, y=-2..2, z=-2..2)
Пример.
Построить трехмерное отображение изменения давления (Р=1,45,25,50,70,100,150,180 МПа) по уравнению регрессии, полученного экспериментально от трех переменных (х1, х2, х3).
Р=5.47-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+ +5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+ +5.02*x3^2.
Чтобы получить уравнения для построения трехмерного отображения нужно указанные в задании величины давления вычесть из первого члена уравнения, таким образом, получим восемь уравнений регрессии
[> restart;
[> with(plots):
Warning, the name changecoords has been redefined
[> implicitplot3d({4.47-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2,
-19.53-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+
5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2,
-39.53-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+
5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2,
-44.53-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+ 5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2,
-64.53-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+ 5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2,
-94.53-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+ 5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2,
-144.53-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+ 5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2,
-174.53-32.19*x1-22.09*x2-11.32*x3+0.625*x1*x2+ 5.63*x1*x3+0.625*x2*x3+38.26*x1^2+16.33*x2^2+5.02*x3^2 }, x1=-1..1, x2=-1..1 ,x3=-1..1, scaling=CONSTRAINED, orientation=[-56,49], axes=boxed);
На графике (рис. 4.11) отображено шесть поверхностей зависимости давления от трех факторов от трех переменных (х1, х2, х3) вместо восьми поверхностей, это объясняется тем, что давление больше 150 нецелесообразно использовать.
Рис. 4.11. Поверхности зависимости давления от трех факторов от трех переменных (х1, х2, х3).