Пакет Mathcad
 Скачать 96.5 Kb.
|
|
МО и ПО РФ НГТУ Реферат по информатике на тему «Пакет MathCAD» факультет: АВТ группа: А - 514 студент: Коваленко С.А. преподаватель: Мамонова В.Г. Новосибирск - 1997 Содержание :
2.2 Документ 1______________________________________________5 2.3 Графические возможности _______________________________6
4. MathCAD или программирование на языках высокого уровня _18 1. Введение Математические и научно - технические расчеты являются важной сферой применения персональных компьютеров . Часто они выполняются с помощью программ , написанных на языке высокого уровня, например Бейсике или Паскале. Сегодня эту работу нередко выполняет обычный пользователь ПК. Для этого он вынужден изучать языки программирования и многочисленные, подчас весьма тонкие капризные численные методы математических расчетов. Нередко при этом из под руки способного физика, химика или инженера выходят далёкие от совершенства программы. Это не вполне нормальное положение может изменить к лучшему применение интегрированных программных систем автоматизации математических расчетов (Eureka, MathCAD, MatLab и др.). Здесь рассматриваются возможности и эволюция одной из таких систем - MathCAD. Фирма MathSoft Inc.(США) выпустила первую версию системы в 1986 г. Главная отличительная особенность системы MathCAD заключается в её входном языке, который максимально приближён к естественному математическому языку, используемому как в трактатах по математике, так и вообще в научной литературе. В ходе работы с системой пользователь готовит так называемые документы. Они одновременно включают описания алгоритмов вычислений, программы управляющие работой систем, и результат вычислений. По внешнему виду тексты мало напоминают обычной программы . 2. Возможности системы MathCAD объединяет в себе простой текстовый редактор, математический интерпретатор и графический процессор. Система ориентирована на IBM - совместимые компьютеры. Упомянутые выше документы содержат текстовые, формульные и графические блоки. На экране дисплея они занимают прямоугольные области, границы которых обычно не видны (но при введении в блок курсора противолежащие углы прямоугольных областей отмечаются прямоугольниками). Блоки выполняются слева направо и сверху вниз. Подготовка и исполнение документов MathCAD может осуществляется с помощью: главного меню и спускающихся подменю (для их появления необходимо нажать клавишу Из режимов работы необходимо особо отметить режимы auto и manual. Режим auto обеспечивает автоматические вычисления сразу после загрузки документа по мере его прокрутки (скроллинга) на экране дисплея. В этом режиме скроллинг ощутимо замедлен, особенно при использовании системы на ПК класса IBM PC XT без математического сопроцессора. Режим manual (ручной) позволяет осуществлять быстрый скроллинг без выполнения документа. Для проведения вычислений от начала документа и до конца видимой на экране дисплея его части нужно нажать клавишу Текстовые блоки являются не более чем комментариями. Их назначение - пояснить сухое и лаконичное математическое описание, представленное на входном языке системы. Текстовые блоки могут быть полноформатными (на всю длину строки) и в виде прямоугольников ограниченных размеров. Если ввести знак «кавычки», то на экране дисплея появится пара кавычек, между которыми вводить и редактировать текст в обычном порядке. Текстовый редактор системы не обладает всеми возможностями специализированных редакторов текста, однако позволяет корректировать тексты, выравнивать их по краю, перемещать текстовые блоки в любое место документа и т.д. Весьма удобны средства редактирования документов, позволяющие, в частности, стирать указанный курсором блок (клавиша   Математический интерпретатор системы - наиболее интересная её часть. Математические формулы, подлежащие интерпретации, записываются в общепринятом виде. Например, вычисление квадратного корня из двух в системе MathCAD задаётся как û2 =, а не в виде PRINT SQR (2) , как это делается, скажем, на Бейсике. Для ввода формул используются шаблоны, вводимые определёнными комбинациями клавиш. Имеется возможность изменения формата представления чисел, например числа знаков после разделительной точки, погрешности вычислений и обозначения мнимой единицы (i на j и наоборот) при операциях с комплексными числами.Чтобы вывести на экран дисплея графический блок, необходимо установить курсор на место левого верхнего будущего графика и ввести знак @. На экране дисплея появится прямоугольник - шаблон будущего графика. Начиная с версии 2.0, масштаб можно и не указывать - он вычисляется автоматически. Возможности системы поясняют следующие конкретные примеры.  Название операции Задание операции Численный пример 2 Возведение в степень X×Y 3 = 9 Вычисление факториала X! 4! = 24 Вычисление квадратного корня \X Ö9 = 3 Вычисление абсолютного значения ½X ½-5½ = 5 Сложение X+Y 2 + 3 = 5 Умножение X*Y 2×3 = 6 12 Деление X/Y `3 = 4 Присваивание значений переменной X:Y X : = 8 Вывод значения переменной X = X = 8 Задание циклической переменной i : = N1..N2 i : = 1,...,5 Суммирование членов ряда i \ X åX = 2i i i 3 Перемножение членов ряда i # X Õ X = 3.84×10 i i 1 Вычисление определённого интеграла x&f(x) ó Ö2xdx = 0.93 õ 0 Задание функции пользователя f(X) : ... f(x): = sin x Инициализация переменной х : ... x : = 1, f (x) = 0.841 df(x) = 0.54 Вычисление производной x ? f (x) ` dx Задание и проверка неравенства X > Y 5 > 4 = 1, 4 > 5 = 0 Вычисление специальных функций Jn(x) = Jn(1,.5) = 0. 242 Бесселя и интеграла вероятности J1(x) = J1(.5) = 0.242 erf(x) = erf(1) = 0.843 Документ 1. Примеры выполнения математических операций . Документ 1 иллюстрирует задание и исполнение в системе MathCAD ряда математических действий. Среди них вычисление алгебраических, тригонометрических и гиперболических функций, сумм и произведений рядов, определённого интеграла и производной. В MathCAD предусмотрены средства для решения нелинейных уравнений, не имеющих аналитических решений. Так , функция root (f(x,y,z,),x) ищет значение переменной x, при котором f(x,y,z) = 0. Более сложные вычисления (решение систем нелинейных уравнений, минимизация функций нескольких переменных и др.) обеспечиваются организацией вычислительного блока, открываемого словом Given. Документ 2 показывает решение актуальной для садовода или любителя баньки задачи: как, сгибая железный лист, получить ящик заданного объёма. Оказывается, есть три решения. Ящик может быть неглубоким, но с дном большой площади или глубоким, но с дном малой площади. Третье решение физически нереально. Во второй части документа показано решение задачи о конструировании ящика максимального объёма, ещё более актуальной при нашем дефиците строительных материалов. Специалистов в электротехнике и радиотехнике наверняка привлечёт способность системы MathCAD выполнять все предусмотренные в ней вычисления как с действительными, так и с комплексными числами. В документе 3 приведены примеры операций с комплексными числами, начиная от простых и кончая сложными. К последним относится вычисление комплексного кругового интеграла, в ходе которого выполняется численное интегрирование и дифференцирование с комплексными аргументами. Начиная с версии 2.0 в MathCAD введён функционально полный набор векторных и матричных операций. Это существенно облегчает решение задач линейной алгебры. В качестве примера в документе 3 даётся решение системы линейных уравнений с комплексными коэффициентами, в ходе которого производится обращение комплексной матрицы. К таким уравнениям приводит анализ электрических и электронных цепей на переменном токе. Весьма привлекательны средства линейной и сплайн-интерполяции и экстраполяции данных. Линейная интерполяция графически означает просто соединение узловых точек графика отрезками прямых. В отличии от неё сплайн-интерполяция напоминает соединение этих точек с помощью гибкой линейки. Строго математически это означает проведение через каждые три точки линии, описываемой кубическим полиномом. При этом во всех стыкуемых точках обеспечивается непрерывность как первой , так и второй производной каждого из полиномов. Сплайн-интерполяция - это мощное средство представления данных, заданных небольшим числом узловых точек. Документ 4 демонстрирует задание в виде векторов напряжений и токов N - образной вольт - амперной характеристики туннельного диода. Затем проводиться интерполяция-экстраполяция этой характеристики описанными в документе способами . Можно заметить ,что сплайн- интерполяция в данном случае предпочтительнее линейной . MathCAD имеет обширный набор статических операторов и функций ,обеспечивающих генерацию случайных чисел ,а также вычисление среднего ,дисперсии и вариации ,коэффициентов линейной регрессии, ряда специальных математических функций .Документ 5 организует генерацию 200 случайных чисел ,их представление на плоскости ,вычисление ряда статистических параметров и построение гистограммы распределения . Графические возможности системы. Мы уже отмечали графические возможности системы. MathCAD позволяет строить самые разнообразные графики: в декартовой и в полярной системе координат, с масштабной сеткой и без неё, с линейным и логарифмическим масштабом, с отметкой линий прямоугольниками, крестами, ромбами и т.д. Задание вида и размера графика осуществляется вводом соответствующего формата. Для задания формата можно ввести внутрь шаблона графика курсор и нажать клавишу logs = 0,0 subdivs = 1,1 size = 5,15 type = 1 Если параметры logs - нули, график строится с линейным масштабом, иначе- с логарифмическим (в этом случае параметры указывают число делений шкалы в пределах декады) . Параметры subdivs задают число делений шкалы, а параметры size - размеры графика, выраженные в знакоместах. Во всех этих случаях первый параметр относится к оси Y графика, второй - к оси Х. Параметр type описывает указание о типе графика в виде малой или большой латинской буквы. Например, указание L задаёт сплошной график, d «строит» точки в узлах и т.д. Возможна комбинация таких указаний.  Конструирование железного ящика заданного объёма VO : = 7.5 Заданный объём ящика W : = 4 Ширина листа L : = 8 Длина листа X : = 0,0.2 .. 5 Расстояние от линии отгиба листа V(X) : = (L - 2X)×(W - 2X)×X Объём ящика F(X) : = V(X) -VO Основное уравнение Графическое решение задачи   60F(X), 0   20 0 X 5 Решение задачи численным методом Инициализация Решение Комментарий X : = 0 root(V(X) - VO, X) = 0.297 Плоский ящик X : = 1 root(V(X) - VO, X) = 1.5 Глубокий ящик X : = 4 root(V(X) - VO, X) = 4.203 Решение физически нереально (X > W/2) Конструирование железного ящика максимального объёма X : = 1 Инициализация Given Начало блока решения V(X) : = (L - 2X)×(W - 2X)×X Основное уравнение V(X) » 100 Объём, заведомо превышающий требуемый X M : = minerr ( X ) Поиск оптимального значения Х X M = 0.848 Найденное оптимальное значение Х V(X M ) = 12.317 Максимально возможный объём ящика Документ 2. Задача о конструирование железного ящика заданного и максимального объёма  Задание мнимой единицы i : = Ö-1 Арифметические операции Z1 : = 2 + 3i Z2 : = 4 + 5i Z : = Z1 + Z2 Z = 6 + 8i Re ( Z ) = 6 Im ( Z ) = 8 sin ( Z1 ) = 9.154 - 4. 169i Вычисление комплексных корней квадратного уравнения х : = 0 + 3i ( Инициализация первого корня ) 2 root(x + 2x +15, x) = - 1 + 3.742i ( Первый корень ) x : = 0 - 3i ( Инициализация второго корня ) 2 root(x + 2x +15, x) = - 1 + 3.742i ( Второй корень ) Решение систем линейных уравнений с комплексными коэффициентами é10 + 200i 0 - 200i ù é5 + 0iù А : =ê ú B : = ê ú ë0 - 200i 0 + 170iû ë0 + 0iû -1 X : = A B (Решение с помощью матричных операторов) é0.037 + 0.131i ù Х : =ê ú (Вектор решения) ë0.044 + 0.154i û Вычисление комплексного кругового интеграла 1 f(x) : = ¾ z(t) : = cos(t) + sin(t) x -4 te : = 6.2832 TOL : = 10 (Погрешность) te ó éd ù ôf(z(t))ê ¾ z(t)ú dt = 6.283
0 Документ 3. Примеры операций с комплексными числами.  Задание ВАХ туннельного диода  0 0.2 50 ( Векторы исходных данных, .4 20 содержащие координаты U: = .6 I : = 3 семи узловых точек ВАХ ) .8 4 1.0 14 1.2 55 Линейная интерполяция ВАХ linterp(U, I, 0.15 ) = 37.5 (Примеры интерполяции ВАХ) linterp(U , I, 0.5) = 11.5 J(V) : = linterp(U, I, V) (Задание функции J(V) ВАХ ) V : = -0.05, - 0.025 .. 1.2  40 На графике ВАХ при линейнойинтерполяции отчётливо видны отрезки прямых, и кривая В АХ  J(V),0 неестественно-40 -0.05 V 1.2 Интерполяция кубическими сплайнами IS : = cspline(U ,I) (Векторы вторых производных) interp (IS, U, I, 0.15) = 49.493 (Примеры сплайн - интерполяции) interp (IS, U, I, 0.5) = 8.191 J(V) : = interp (IS, U, I, V ) (Задание функции J(V) ВАХ) V : = -0.05, - 0.025 .. 1.2 40 Кривая ВАХ при сплайн -интерполяции отличается плавностью и похожа на J(V),0 реальную кривую ВАХ-40 -0.05 V 1.2 Документ 4. Линейная и сплайн - интерполяция N - образной вольтамперной характеристики (ВАХ) туннельного диода .  Генерация 200 случайных чисел с равномерным распределениемi : = 1..200 x : = rnd ( 10 ) i Графическое представление случайных чисел  10 х График наглядно показывает i равномерность распределения случайных чисел 0 1 i 200 Вычисление основных статистических параметров массива х mean(x) = 4.619 var(x) = 8.869 max (x) = 9.95 min (x) = 0 stdev (x) = 2.978 Подготовка данных к построению гистограммы N : = 10 j : = 0..N k : = 0..N - 1 10 intervals : = 1 + j ¾ P : = hist (intervals, x ) j N Гистограмма распределения чисел в массиве х  30     P k0 0 intervals 10 k |