<<
>>

Практическое занятие

30 мин Упражнение 18.1. Простые вычисления с использованием программы Mathcad

Задача.

Найти ребро куба, равновеликого шару, площадь поверхности которого равна площади боковой поверхности прямого кругового конуса, у которого высота вдвое меньше, чем длина образующей. Объем этого конуса равен 1.

Анализ. Основные геометрические формулы, используемые при расчете.

1 1

Объем конуса — V=>-rrm\.

Площадь боковой поверхности конуса — S = ш1.

Соотношение в конусе между радиусом основания, высотой и длиной образующей — г[1] [2] [3] + К1 = Р.

Площадь поверхности шара — V=4nR2.

4 3 ч

Объем шара — V=^nr . Объем куба — V= a .

(Определить) на панели инструментов Calculator (Счет). Итак, надо ввести V: 1 В документе появится полноценный оператор присваивания.

Е:=1

4. Путем несложных преобразований получим, что радиус основания конуса можно вычислить по формуле

Вводить эту формулу следует слева направо. Порядок ввода этой формулы следующий. Сначала введите знак корня произвольной степени: кнопка Nth Root (Корень данной степени) на панели инструментов Calculator (Счет) или комбинация клавиш CTRL+V Щелкните на черном квадратике, стоящем на месте показателя степени, и введите цифру 3. Щелкните на квадратике, замещающем подкоренное выражение, нажмите клавиши V*.

Введите знак квадратного корня: кнопка Square Root (Квадратный корень) на панели инструментов Calculator (Счет) или клавиша\ — и цифру 3. Прежде чем вводить знаменатель, дважды нажмите клавишу ПРОБЕЛ. Обратите внимание на синий уголок, который указывает на текущее выражение. Предполагается, что знак операции связывает выбранное выражение со следующим. В данном случае это безразлично, но в целом этот прием позволяет вводить сложные формулы, избегая ручного ввода дополнительных скобок. Нажмите клавишу /. Чтобы ввести число я, можно воспользоваться комбинацией клавиш CTRL+SHIFT+P или соответствующей кнопкой на панели инструментов Calculator (Счет).

5. Введите формулы для вычисления длины образующей и площади боковой поверхности конуса:

1= -ОД ; S = п-г-1

Указание знака умножения между переменными обязательно, так как иначе программа МаЛеаЛ сочтет, что указана одна переменная с именем из нескольких букв.

рументов Calculator (Счет). После формулы появится знак равенства и вычисленный результат. а= 0.7102

ШИ Вычислять можно как действительные, так и комплексные выражения. Обозначение мнимой единицы (i) следует вводить непосредственно после числового коэффициента, который нельзя опускать, даже если он равен единице.

9. Вернитесь к самому первому выражению и отредактируйте его. Вместо значения 1 присвойте переменной значение 8. Сразу же перейдите к последней введенной формуле и обратите внимание, что результат расчета сразу же стал отражать новые начальные данные.

ЩМы познакомились с методикой простейших вычислений в программе MathCad. Описанная техника позволяет использовать эту программу как «интеллектуальный калькулятор» для автоматического расчета по известным формулам.

Особенностью программы Mathcad является возможность практически мгновенного перерасчета с другими начальными данными.

Упражнение 18.2. Физические вычисления с использованием единиц измерения

30 мин

Постановка задачи. Теплоизолированный космический аппарат, находящийся на орбите Земли, имеет на борту приборы с электрической мощностью, которая может изменяться в ходе работы от А) = 75 Вт (дежурный режим) до = 200 Вт (сеанс связи). С целью обеспечения предсказуемого теплового режима в теплоизоляции сделано отверстие площадью 5j ,на которое попадает поток солнечной энергии W— 1400 Вт/м2. Полученная энергия излучается аппаратом через это и дополнительное отверстие в теплоизоляции с площадью S2 в режиме «черного тела». Каковы должны быть площади отверстий, если допустимый диапазон температур для оборудования, расположенного в аппарате, составляет 20-30°С? Анализ задачи. Минимальная температура аппаратуры соответствует режиму минимального тепловыделения. В этом случае поступающая мощность где 7) — минимальная допустимая температура в градусах Кельвина. В условиях теплового баланса эти мощности должны быть равны.

Режим максимального тепловыделения соответствует максимальной температуре аппаратуры. В этом случае W5) + А2 — стГ24 (5) + S2).

1. Запустите программу Mathcad.

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

W: =1400

М:= 75-watt N2 := 200-watt П := (20 + 273)-К Т2 := (30 + 273)-К

3. Обозначения физических единиц присоединяйте к соответствующим значениям через знак умножения. Если нужное обозначение неизвестно, используйте команду Insert • Unit (Вставка • Единица измерения).

Измеряемая величина выбирается в списке Dimension (Размерность), а нужная единица измерения — в списке Unit (Единица измерения).

4. Присвойте переменной а значение постоянной Стефана-Больцмана

.

м •tv

Чтобы ввести греческую букву, используйте панель инструментов Greek (Греческий алфавит) или введите соответствующую латинскую букву (в данном случае «5») и сразу же нажмите комбинацию клавиш CTRL+G. Так как специальной единицы для размерности этой константы не существует, ее следует составить из стандартных единиц с помощью умножения и деления.

5. Введите полученные в ходе анализа формулы для вычисления площадей отверстий, присвоив полученные значения переменным 51 и 52.

(N2-Ti4-Nl-T24) W- (N2 - N1)- c-(N2- Tl-m-T24)

^94 ••=-------------------------- . ссо.- , ,

И'-(724-Г14) gW-П/- П4)

6. Изменение значений параметров, заданных в условии задачи, приводит к автоматическому перерасчету формул. В частности, исследуйте, изменяя значение переменной W, как изменяются требования к такому методу терморегуляции при удалении аппарата от Солнца и приближении к нему

Вт Вт

W= 2700 ——;на орбите Марса W= 500 —п~). м*- м41

7. Обратите внимание, что результат содержит единицы измерения в соответствии с системой единиц СИ. Используемая система единиц отображается в диалоговом окне Insert Unit (Вставка единиц измерения).

8. Чтобы изменить используемую систему единиц, дайте команду Math • Options (Математика • Параметры) и в открывшемся диалоговом окне Math Options (Параметры расчета) выберите вкладку Unit System (Система единиц). Выберите систему CGS и посмотрите, как изменились результаты (они теперь выражаются в квадратных сантиметрах). Если, например, выбрать американскую систему единиц (U.S.), то результат будет выражен в квадратных футах.

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

Это позволяет немедленно получать результат в наиболее удобной форме.

Упражнение 18.3. Векторы и матрицы 30 мин

1

1 по нормированным собственным векторам мат- 1

Анализ. Первый этап решения задачи состоит в нахождении собственных значений и собственных векторов данной матрицы. Затем необходимо найти вектор Т, такой, что5Т=У где S — матрица, столбцы которой представляют собой собственные вектора матрицы М.

1. Запустите программу Mathcad.

2. Создайте матрицу М. Начните запись оператора присваивания, а для ввода правой части нажмите комбинацию клавиш CTRL+M, воспользуйтесь командой Insert • Matrix (Вставка • Матрица) или щелкните на кнопке Matrix or Vector (Матрица или вектор) на панели инструментов Matrix (Матрица).

3. В открывшемся диалоговом окне Insert Matrix (Вставка матрицы) укажите число строк и столбцов (по три) и щелкните на кнопке ОК

4. Введите значения элементов матрицы в отведенные места.

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

6. Собственные значения квадратной матрицы можно получить при помощи функции eigenvals. Результатом ее работы является вектор собственных значений, присвойте его переменной L.

7. Функция eigenvec позволяет получить собственный вектор, соответствующий данному собственному значению.

Ей нужны два параметра: матрица, для которой ищется собственный вектор, и собственное значение, которому он соответствует. Чтобы записать собственные вектора в качестве столбцов матрицы S, надо присвоить вычисленное значение столбцу матрицы. Столбцы матрицы в программе Mathcad выбираются специальным верхним индексом, заключенным в угловые скобки. Чтобы ввести номер столбца, нажмите комбинацию клавиш CTRL+6 или щелкните на кнопке Matrix Column (Столбец) на панели инструментов Matrix (Матрица), после чего введите номер нужного столбца матрицы. Будьте внимательны — столбцы и строки матрицы нумеруются начиная с нуля.

8. В правой части оператора присваивания надо указать собственное значение матрицы. Собственные значения являются элементами вектора L. Номер элемента указывается как нижний индекс. Для ввода нижнего индекса нажмите клавишу [ или воспользуйтесь кнопкой Subscript (Индекс) на панели инструментов Matrix (Матрица). Итоговый оператор для первого собственного вектора будет выглядеть следующим образом:

5 :=eigenvec(M,I0)

Аналогично задайте операторы для второго и третьего собственных значений.

9. Для нахождения коэффициентов при собственных векторах в разложении необходимо решить систему линейных уравнений. Ее удобно записать в матричной форме. Создайте вектор Г с тремя элементами. Величины этих элементов значения не имеют.

10. Запишите ключевое слово given.

11. Ниже запишите матричное уравнениеS-T= V. Знак логического равенства введите с помощью комбинации клавиш CTRL+=.

12. Найдите коэффициенты в разложении при помощи функции/ind.

0, 836

0, 334

0, 148

Мы научились производить операции с векторами и матрицами, использовать соответствующие функции, выделять столбцы матриц и отдельные элементы. Матричная запись часто позволяет представить задачу в более удобной форме.

Упражнение 18.4. Аналитические вычисления

45 мин

Задача 1. На приведенной схеме сопротивление ЕЕ является переменным. Определить, как меняется ток между точками А и В в зависимости от величины этого сопротивления.

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

д-/4-д-/3= О

Эту систему надо решить, не подставляя конкретных значений вместо параметров R, RR и Е,

1. Запустите программу Mathcad.

2. Введите ключевое слово given.

3. Введите уравнения системы, полученной в ходе анализа. Обозначьте неизвестные токи переменными /О, II, 12,13,14. Фиксированное сопротивление Я обозначьте переменной #0. Обратите внимание, что присваивать начальные значения токов или задавать значения переменных Я0, ЛЛ иЕне требуется.

4. Введите функцию find, перечислив в качестве параметров неизвестные /О, Л, 12,13,14. Затем введите оператор аналитического вычисления, который выглядит как стрелка, направленная вправо, и вводится комбинацией клавиш CTRL+. или кнопкой Evaluate Symbolically (Вычислить аналитически) на панели инструментов Evaluation (Вычисление).

5. Щелкните за пределами данного блока, и программа Ма^еаё произведет аналитическое решение системы уравнений.

^ (3-RR+R0)

(R0-(5-RR+3-RQ))

4---------- 5-------

(5-RR+3-R0))

(-RR+R0)

Е• -------- ‘---------------------

(R0-(3-RR+3-R0))

^ (RR+3-R0)

(R0-(5-RR+3-R0))

(RR+R0)

2-Е------------------------------

(ДО-(5-ДД+3*йО))

Полученный результат позволяет провести полный анализ схемы.

Задача 2. Найти все корни уравнения:

(1 +y-y2fy= 2

Анализ. Это уравнение четвертого порядка. Легко подобрать один корень {у = 1). Остающееся уравнение третьего порядка не имеет рациональных корней, так что поиск других корней этого уравнения — дело непростое. Неясно даже, сколько еще действительных корней имеет данное уравнение. Результаты численного решения зависят от подбора начального приближения и поэтому не гарантируют отыскания всех корней уравнения. Мы же решим это уравнение аналитически.

6. Введите заданное уравнение. Чтобы раскрыть скобки, дайте команду Symbolics • Simplify (Аналитические вычисления • Упростить).

7. Выделите в полученном уравнении независимую переменную (в данном случае#) и дайте команду Symbolics • Variable • Solve (Аналитические вычисления • Переменная • Решить).

Программа Mathcad выдаст вектор, элементами которого являются корни данного уравнения.

8. Полученный результат содержит сложные комплексные выражения, и его невозможно применить с пользой (все еще непонятно, являются ли корни действительными или комплексными). Чтобы разделить действительную и мнимую части, выделите результат вычисления целиком и дайте команду Symbolics • Evaluate • Complex (Аналитические вычисления • Вычислить • В комплексном виде). Если программа Mathcad не справится с преобразованием всего набора корней целиком, выполните преобразование корней поочередно: дайте команду Symbolics • Evaluate • Complex (Аналитические вычисления • Вычислить • В комплексном виде), поочередно выделив каждый из корней, записанных в комплексном виде.

9. Теперь полученное выражение надо упростить. Выделив его целиком, дайте команду Symbolics • Simplify (Аналитические вычисления • Упростить). Выражение станет существенно проще, причем станет понятно, что все корни уравнения действительные (все мнимые компоненты сократятся).

10. Последний шаг — раскрытие скобок, в данном случае упрощение аргументов тригонометрических функций. Для этого примените команду Symbolics • Expand (Аналитические вычисления • Раскрыть). Полученная запись—наилучшее представление точного решения, которое можно получить с помощью программы Mathcad.

11. Чтобы получить результат в числовом виде, достаточно ввести в конце выражения (итогового или на любой из предыдущих стадий) команду вычисления (=).

ИИ Мы научились использовать программу Mathcad для выполнения аналитических вычислений. Это позволяет получать точные решения задач, содержащих переменные параметры, анализировать полученные результаты, а также получать полный набор решений для некоторыхтиповуравнений.

Упражнение 18.5. Анализ результатов испытаний

30 мин

Задача. К пружине последовательно подвешивали грузы массой 1, 2, 3, 20 кг.

В результате был получен список величин удлинения пружины (в миллиметрах). Определить основные статистические параметры полученного набора измерений. Рассчитать жесткость пружины и массу узла, использованного для крепления грузов к пружине, воспользовавшись методом наименьших квадратов.

Таблица измерений:

Анализ. Для решения этой задачи достаточно использовать стандартные средства статистических вычислений, имеющиеся в программе Mathcad, Теоретически, растяжение пружины определяется формулой к • х ={т + т0) ■ g. Если определить статистическими методами коэффициенты а и b уравнении х-а • т + Ь, то получим:

а и а

1. Запустите программу Mathcad.

2. Введите таблицу данных, предназначенных для статистического анализа, как матрицу с двумя столбцами, первый из которых содержит веса грузов, а второй — значения растяжения пружины.

3. Определите число точек в наборах данных с помощью функции rows.

п := rows(data) n = 20

4. Вычислите среднее растяжение пружины в ходе эксперимента с помощью функции mean.

Y= data^ mean(Y) = 31.695

1. Вычислите медиану значений растяжения пружины при помощи функции median.

median(Y) = 30.7

2. Вычислите среднеквадратичное отклонение и дисперсию величины растяжения пружины при помощи функции stdev.

stdev( У) = 17.39625 stdev(Y)2= 302.62947

3. Определите коэффициенты линейного уравнения являющегося наилучшим приближением для данных наборов данных. Функция slope позволяет вычислить коэффициент наклона прямой, а функция intercept — свободный член.

Х:= data

bQ := intercepted Y) bQ = 0.07421

Ьл := slope(d^) b{ = 3.0115

4. Определите жесткость пружины.

k= 3,254-103 (Н/м)

5. Определите массу узла крепления.

т = 24,64 (г)

6. Сохраните созданный документ для использования в следующем упражнении.

[У| Мы научились применять функции, используемые для статистического анализа данных. Программа Mathcad содержит и другие функции аналогичного назначения, которые можно использовать для интерполяции и экстраполяции данных, а также их сглаживания.

Упражнение 18.6. Построение графиков

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

Анализ. Для построения графика можно использовать функцию, заданную набором данных или формулой. Формулы для функций, полученных в результате проделанных расчетов, необходимо определить, прежде чем их можно будет использовать при построении графика.

1. Запустите программу Mathcad.

2. Загрузите документ, созданный в предыдущем упражнении.

3. Переместите точку ввода в нижнюю часть документа.

4. Запишите формулу функции г(х) для определения координат точек, лежащих на аппроксимирующей прямой. Коэффициенты соответствующего уравнения были получены в предыдущем упражнении.

r(x) := b0+ Ьх

5. Нажмите клавишу @, щелкните на кнопке X-Y Plot (Декартовы координаты) на панели инструментов Graph (График) или дайте команду Insert • Graph • X-Y Plot (Вставка • График • Декартовы координаты). В документе появится область для создания графика.

6. Вместо заполнителя в нижней части графика укажите в качестве независимой переменной первый столбец матрицы data (datai0> или X).

7. Вместо заполнителя слева от графика укажите, что по вертикальной оси должны откладываться значения из второго столбца матрицы data и определенная выше линейная функция г(Х). В качестве разделителя используется запятая. Диапазон значений для осей координат выбирается программой Mathcad автоматически.

8. Чтобы изменить вид автоматически построенного графика, дважды щелкните внутри него. Откроется диалоговое окно Formatting Currently Selected X-Y Plot (Форматирование графика в декартовых координатах). Первая запись в списке на вкладке Traces (Кривые) соответствует первой отображенной кривой. Для изменения записи используются поля под списком.\

Рис. 18.5. Задание способа отображениялиний графика

9. Под столбцом Legend Label (Подпись) введите название графика.

10. В раскрывающемся списке под столбцом Symbol (Маркер) выберите способ обозначения для отдельных точек.

11. Под столбцом Туре (Вид линии) укажите, что необходимо пометить отдельные точки (points), а не провести непрерывную линию.

12. Выберите в списке вторую кривую и настройте ее отображение по своему вкусу.

13. Установите флажок Hide Arguments (Скрыть параметры), чтобы не отображать названия осей.

14. Сбросьте флажок Hide Legend (Скрыть подписи), чтобы включить отображение под графиком заданных подписей кривых.

15. В поле Title (Заголовок) на вкладке Labels (Надписи) задайте название графика и включите режим его отображения: флажок Show Title (Показать заголовок).

16. Постройте график, на котором отображалась бы величина отклонения экспериментальных точек от линии приближения Y- (Ь0 + Ьх-Х). Отформатируйте его, используя те же средства, что и в предыдущем случае.

Я Заголовок и подписи, использующие русские буквы, могут отображаться неправильно. Коррекцию обеспечивает выбор шрифта, правильно воспроизводящего кириллицу. Дайте команду Format • Equation (Формат • Выражение), в раскрывающемся списке Style Name (Имя стиля) выберите пункт Variables (Переменные) и щелкните на кнопке Modify (Изменить). Для задания шрифта используйте Полесо списком Шрифт.

Упражнение 18.7. Построение трехмерных графиков

30 мин

Задача. Изобразить на графике приблизительную форму электронных облаков в атомах. Анализ. По современным представлениям, электронные уровни в атоме определяются четырьмя квантовыми числами. Форма электронного облака определяется двумя из этих чисел:

• число / определяет тип орбитали (значения 0-3 соответствуют я-, р-, й- и /орбиталям);

• число т определяет магнитный момент электрона и может изменяться в диапазоне от -/до /.

При т=0 форма электронного облака определяется на основе многочленов Лежандра первого рода:

В этом случае ДФ) - \

Параметрическое задание соответствующей x(0,0)=Y^)-sin ф'соэ 9 у(0,ф)=У(ф)-зт ф-sin 9

2(0,ф)=У(ф)-СО5 ф

Углы 0, ф изменяются в диапазоне от 0 до 2л.

1. Запустите программу МаШеай.

2. Определите переменную /, которая укажет тип орбитали.

/:= 3

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

Удобно определить границы диапазона в целых числах (через точку с запятой, на экране изображаются две точки), а затем произвести перемасштабирование на отрезок [0; 2я].

¿:= 0..100 7 := 0..100

. 6. В появившейся области графика вместо заполнителя укажите имена отображаемых матриц через запятую, заключив все их в скобки: (X0,Y0,Z0).

7. Чтобы изменить формат построенного графика, дважды щелкните на его области. Откроется диалоговое окно 3-D Plot Format (Формат трехмерного графика).

8. На вкладке General (Общие) установите флажок Equal Scales (Равный масштаб), чтобы обеспечить одинаковый масштаб по осям координат.

9. На вкладке Appearance (Оформление) установите переключатель Fill Surface (Заливка поверхности), чтобы обеспечить заливку построенной поверхности.

10. На вкладке Lighting (Подсветка) включите режим освещения поверхности. Установите флажок Enable Lighting (Включить подсветку), отключите все источники света, кроме первого.

11. На панели Light Location (Размещение источника) задайте координаты источника света. Используйте кнопку Применить, чтобы сразу видеть последствия сделанных настроек. По окончании настройки закройте диалоговое окно щелчком на кнопке ОК

12. Путем протягивания мыши в области графика измените направление осей координат, чтобы изображение было видно наиболее отчетливо.

13. Изменяя значение /, можно увидеть форму электронных облаков для разных орбиталей, в том числе и не встречающихся в природе.

Мы научились строить трехмерные графики с изображением объемных поверхностей, заданных параметрически. Это фактически означает умение изображать любые фигуры, которые могут потребоваться в ходе практической работы.

Упражнение 18.8. Решение дифферен циальных уравнений

Задача. Найти функцию у(х), удов- летворяющую дифференциальному с1у

уравнению------- + у = хбшх и имею-

dx

щую значение 0 при х - 0 .

Анализ. Это простое дифференциальное уравнение допускает точное

аналитическое решение. В данном Рис. 18.7. Трехмерное изображение электронной упражнении предполагается исполь- /-орбитали

зование стандартной функции программы Mathcad, осуществляющей численное решение данного уравнения. Результат вычислений можно после этого сравнить с точным решением.

1. Запустите программу Mathcad.

2. Задайте начальное значение функции как элемент вектора у, размерность которого соответствует числу решаемых уравнений (в данном случае единице): г/0 := 0.

3. Создайте функцию Т(х, у), которая вычисляет значение производной при заданных значениях независимой переменной и неизвестной функции:

Т(х,у) :=-у0+лт-со8(*)

4. Определите начальное (точка 0) и конечное значение отрезка интегрирования.

я:=0, Ъ:= 12-71

5. Укажите число шагов интегрирования.

К =20

6. Вычислите численное решение уравнения при помощи функции гкДхвё.

Z := ткйхеё(у,а,Ь,К,Т)

Результат вычислений — матрица Z с двумя столбцами, первый из которых содержит значения независимой переменной, а второй — соответствующие значения функции.

7. Постройте график полученного решения.

8. Определите аналитическое решение данного уравнения при тех же начальных условиях.

9. Нанесите аналитическую кривую на тот же график и сравните поведение численного и точного решения.

Рис. 18.8. Графики численного и точного решения дифференциального уравнения

10. Измените число шагов, на которые делится отрезок интегрирования, и исследуйте, как изменяется результат расчета при уменьшении и увеличении этого параметра.

¡Ц Мы научились численно решать дифференциальные уравнения первого порядка с помощью программы Ма^саФ Использованный метод без изменений переносится на системы, содержащие два или большее число дифференциальных уравнений. Увеличение величины шага интегрирования ускоряет получение результата, но снижает его точность. При слишком большой величине шага результат расчетов может вообще не соответствовать реальному решению.

<< | >>
Источник: Под ред. С. В. Симоновича. Информатика. Базовый курс. 2-е издание. 2004

Еще по теме Практическое занятие:

  1. Занятие любовью
  2. Занятие 6. Образ «Я»
  3. Занятие 14. Я — организатор
  4. Статья 203-2. Занятие игорным бизнесом
  5. 1- Занятия не должны мешать семье.
  6. Взгляд в будущее... О занятиях и карьере
  7. Место занятий - весь мир
  8. Схема занятия по психогимнастике
  9. Посетите занятия по импровизационной комедии
  10. Занятие 13. Я в ситуациях взаимодействия
  11. Организация отдыха, досуга и занятости несовершеннолетних
  12. Занятие 7. Самораскрытие. Принятие себя