Пакет математических вычислений MAPLE.



 

Меню системы. Справка. Аналитические вычисления. Способы определения функций. Визуализация данных средствами пакета MAPLE.

Дополнительные настройки графических процедур.

 

Цель работы:

Ознакомление с интерфейсом системы математических вычислений MAPLE. Освоение приемов работы в системе. Изучение основных структур и типов данных, обрабатываемых системой.

Ознакомление с графическими возможностями и графическими процедурами системы математических вычислений MAPLE. Освоение приемов работы в системе с графическими процедурами. Использование возможностей пакета для визуализации данных.

 

Порядок выполнения задания

 

Ознакомиться с интерфейсом системы математических вычислений MAPLE, с пунктами меню системы, с использованием справочной системы пакета и выполнением примеров справки.

 

1. Запустить пакет MAPLE, ознакомиться с разделами главного меню системы и пунктами этих разделов. Просмотреть всплывающие подсказки к кнопкам инструментальной панели. Переключиться в текстовый режим, набрать несколько строк текста, выполнить форматирование текста, просмотреть кнопки появившейся контекстной панели инструментов.

2. Перейти в режим вычислений. Набрать несколько аналитических выражений по правилам синтаксиса системы MAPLE и освоить работу с контекстным меню (правая кнопка на результате вычислений).

3. Вызвать справочную систему через раздел главного меню Help. Ознакомиться со структурой справочной системы, с путеводителем по справке. Освоить механизм выполнения примеров из справки (пункт меню Edit\ Copy Examples).

4. Перейти в первоначальное окно, подключить пакет расширений «student», вызвать справку по процедуре пакета «middlebox» и выполнить в новом окне пример из справки.

5. Сохранить вычисления первоначального окна в личной папке. Удалить вычисления с рабочего листа (пункт меню Edit \ Remove Output ), сохранить полученный файл в личной папке.

6. Выполнить вычисления в новом рабочем окне с данными разных типов.

 

Образец выполнения задания по пункту «аналитические вычисления»:

1. Найти неопределенный интеграл

> restart:with(student):Int((3*x+4)*exp(3*x),x);

> F(x):=value(%);

Проверка дифференцированием

> f(x):=diff(F(x),x);

2. Вычислить определенный интеграл

> Int((x^2+4*x+3)*cos(x),x=-1..0);

> value(%);

> evalf(%);

 Первообразная

> F(x):=int((x^2+4*x+3)*cos(x),x);

Интеграл=F(0)-F(-1)

> R:=subs(x=0,F(x));S:=subs(x=-1,F(x));

3. Вычислить площадь фигуры

 

> y1:=4-x^2;y2:=x^2-2*x;

> plot([y1,y2],x=-2..3);

Найдем абсциссы точек пересечения

> solve(y1=y2,x);

Следовательно, площадь вычисляется  как

> Int(y1-y2,x=-1..2);

> value(%);

4. Найти расстояние от точки до плоскости

> restart;

Составим уравнение плоскости

> A:=matrix(3,3,[[x+3,y+1,z-1],[-9+3,1+1,-2-1],[3+3,-5+1,4-1]]);

> with(linalg):det(A)=0;

Warning, the protected names norm and trace have been redefined and unprotected

 

Это уравнение плоскости с коэффициентом =0 при у

Расстояние от точки до плоскости вычислим по формуле

> A:=-6;B:=0;C:=12;x[0]:=-7;y[0]:=0;z[0]:=-1;

 

> d := abs(A*x[0]+B*y[0]*C*z[0]*D)/sqrt(A^2+B^2+C^2);

5. Найти точку пересечения прямой и плоскости

Положим

выразим x, y, z через t и подставим в уравнение плоскости

> P:=x-3*y+7*z-24=0;

> x:=1-t;y:=4*t-5;z:=2*t+1;

> P;

Отсюда t = 1; и координаты точки пересечения

x = 0,y = -1, z = 3

 

> t:=1;x;y;z;

Визуализация данных.

 

Ознакомиться с пакетом графики «plots» и внутренними процедурами пакета с помощью справочной системы. Выполнить построение двумерных графиков функций с установкой дополнительных свойств графики: вставка текста, вид линий, надписи и легенда рисунка; построение поточечных графиков; анимация.   

 

1. Вызвать справку по процедуре “plot” пакета, выполнить примеры справки. Ознакомиться с дополнительными опциями команды “plot”.

2. Построить график затухающего сигнала, использовав процедуры построения графиков из справочной системы.

3. Построить график ломаной линии, заданной координатами точек. Вставить в рисунок надписи, установить стиль линии, цвет, толщину линии, надписи на осях как в образце (файл: «lab8.mws»).

4. Построить анимацию затухающего сигнала с изменением частоты гармоники.

5. Построить графики функций, заданных аналитически. В случае разрывных функций подобрать параметры процедуры “plot” для правильного отображения графиков.

6. Рассмотреть примеры построения графиков в полярной системе координат.

7. Ознакомиться с примерами трехмерной графики и анимации, вызвав соответствующие разделы справки и выполнив примеры из справки.

 

Образец выполнения задания (файл « lab 8. mws »):

 

Первые шаги в графике. Наберем в командном режиме следующие команды:

> restart;with(plots);

Ознакомимся с названиями процедур пакета. Для вызова справки по любой процедуре необходимо выполнить щелчок по названию процедуры и нажать клавишу "F1".

Построим график функции, задав ее формулой.

>

> sin(x)/x;

> smartplot(sin(x)/x);

Правой кнопкой укажем на формулу и выберем действие построения графика "plot2d".

добавим график функции sqrt(abs(x)), пересадив формулу-отклик на график левой кнопкой

Ø sqrt(abs(x));

 

 

 

Найдем на рисунке точку пересения графиков функций и для сравнения вычислим ее

Ø fsolve(sqrt(abs(x))=sin(x)/x,x=2 );

 

 

Рисуем график функции (1-x)/(1-x)

> y:=(1-x)/(1+x);

> plot(y,x=-6..3,-5..5,discont=true,color=blue);

Справка по команде plot

> ?plot

Скопируем примеры командой меню Edit \Copy Examples и вставим в новый рабочий лист

Выполним все команды примеров Edit \ execute\ worksheet.

 

Пример на анимацию:

 

Включим панели символов View \ Palettes | Show all palettes

будем с панели символов снимать греческие буквы

> f:=exp(-alpha*t)*sin(omega*t);

Зафиксируем

> alpha:=2;

> animate(f,t=0..3,omega=1..6);

Спрвка по команде 3-х мерной графики

> ?plot3d

Скопируем примеры и посмотрим 3-х мерные графики и анимацию (справка по команде animate)

Строим график ломаной линии

> M:=[[0,0],[-3,3],[4,4],[6,0],[4,-3],[-4,-3],[0,0]];

> plot(M,linestyle);

Устанавливаем атрибуты точечного рисунка и захватываем его без прорисовки, присваивая объект переменной без вывода результата (двоеточие в конце оператора)

> Ris1:=plot(M,style=point, symbol = circle, symbolsize =20, color=blue):;

Аналогично рисуем ломаную с установкой стиля линии

> plot(M,color=brown,linestyle=3,thickness=3,labels=[t,W],labelfont=[SYMBOL,18],xtickmarks=3,labeldirections=[HORIZONTAL, VERTICAL],title=Ломаная, titlefont=[TIMES,BOLDITALIC,18]);

Захватим и этот рисунок и объединим на одном

> Ris2:=plot(M,color=brown,linestyle=3,thickness=3,labels=[t,W],labelfont=[SYMBOL,18],xtickmarks=3,labeldirections=[HORIZONTAL, VERTICAL],title=Ломаная, titlefont=[TIMES,BOLDITALIC,18]):

> display(Ris1,Ris2);

 

дополнительные надписи на рисунках

> textplot([1,2,`это ломаная линия`],align={ABOVE,RIGHT});

> Ris3:=textplot([1,2,`это ломаная линия`],align={ABOVE,RIGHT}):

Ø display(Ris1,Ris2,Ris3);

Ø

Совместный график нескольких функций

 

> F1:=sin(x)/x;F2:=x^2*exp(-x);

> plot([F1,F2],x=-1..4,color=[blue,red]);

Графики в полярных координатах

Архимедова спираль

> rho:=phi;

> plot(rho,phi=0..20*Pi,coords=polar,color=brown);

Кардиоида

> rho:=4*(1-cos(phi));

> plot(rho,phi=0..2*Pi,coords=polar);

Улитка Паскаля

> alpha:='alpha':#освобождаем фиксированную переменную

> r:=2*(1+cos(alpha));

> plot(r,alpha=0..10,coords=polar);

Розы

> r:=sin(3*phi);:

> plot(r,phi=0..2*Pi,coords=polar);

> r:=5/(1-cos(alpha));

> plot(r,alpha=0..2*Pi,coords=polar);

Лемниската в неявной форме

освободим назначение переменной y

> y:='y':

> Z:=(x^2+y^2)^2-2*(x^2-y^2);

> implicitplot(Z,x=-5..5,y=-5..5,numpoints=5000);

Линии второго порядка

> Z:=5*x^2+4*x*y+8*y^2-32*x-56*y+80=0;

> implicitplot(Z,x=-5..8,y=0..8,numpoints=5000);

>

Анимация в полярных координатах

> rho:=sin(t*phi);

> animate([rho,phi,phi=-2*Pi..2*Pi], t=1..5, coords=polar, frames=100,numpoints=100);

>

 

 

Вычисление пределов

> restart;with(student);# подключение пакета студенческих вычислений

> Limit((x^2-4*x+1)/(2*x+1),x=2);

> value(%);# вычислить значение предыдущей команды

> Limit((5*x^3-7*x)/(1-2*x^3),x=infinity);value(%);

> y:=1/(x-1)-2/(x^2-1);

рисуем график этой функции, учитывая разрыв

> plot(y,x=-2..2,-5..5,discont=true,color=green,thickness=3);

> Limit(y,x=1);value(%);

Дифференцирование

> f:=sqrt(1-x^2);

> y1:=diff(f,x);

> y2:=diff(f,x$2);

> simplify(%);

Примеры функций, для которых требуется построить графики, подобрав параметры процедуры «plot» в случае разрывной функции.

>

>

>

>

>

>

>

>

>

>

Исследование функций

> y:=x/(1+x^2);

> plot(y,x=-5..5);

Точки пересечения с осью Ох

> solve(y,x);

Производная

> y1:=diff(y,x);Y1:=simplify(%);

Точки предположительного экстремума как корни производной

> solve(Y1,x);

Значения функции в этих точках

> subs(x=-1,y);subs(x=1,y);

Вторая производная

> y2:=diff(y,x$2);Y2:=simplify(%);

Знаки второй производной в точках экстремума

> subs(x=-1,Y2);

минимум

> subs(x=1,Y2);

максимум

Точки перегиба

> solve(Y2,x);

 


Практическая работа 9


Дата добавления: 2018-10-27; просмотров: 528; Мы поможем в написании вашей работы!

Поделиться с друзьями:






Мы поможем в написании ваших работ!