- •Федеральное агентство по образованию
- •2404000000-35 Удк 681.142:519.6
- •Оглавление
- •1. Основы работы в Mathcad 10
- •2. Роль численных методов 36
- •3. Методы аппроксимации и интерполирования 42
- •4. Лабораторная работа № 1. Интерполирование степенными многочленами 44
- •5. Лабораторная работа № 2. Параметрическая идентификация математических моделей методами аппроксимации 65
- •6. Решение систем линейных алгебраических уравнений 93
- •7. Лабораторная работа № 3. Решение систем линейных уравнений приближенными методами 109
- •8. Лабораторная работа № 4. Решение нелинейных уравнений приближенными методами 142
- •9. Лабораторная работа № 5. Решение систем нелинейных уравнений приближенными методами 174
- •10. Лабораторная работа № 6. Решение обыкновенных дифференциальных уравнений первого порядка численными методами 197
- •11. Лабораторная работа № 7. Численное интегрирование 212
- •12. Лабораторная работа № 8. Моделирование реактора идеального вытеснения для многостадийной химической реакции с линейной кинетикой 226
- •13. Лабораторная работа № 9. Расчет моделей процессов диффузии и теплопроводности с помощью явной разностной схемы 249
- •Предисловие
- •1. Основы работы в Mathcad
- •1.1. Панели инструментов
- •1.2. Ввод и вывод данных
- •1.3. Осуществление несложных вычислений
- •1.4. Построение и настройка графиков
- •1.5. Программирование в Mathcad
- •1.5.1. Программирование без программирования
- •1.5.2. Язык программирования Mathcad
- •1.5.3. Создание программы (Add Line)
- •1.5.4. Редактирование программы
- •1.5.5. Локальное присваивание ()
- •1.5.6. Условные операторы (if, otherwise)
- •1.5.7. Операторы цикла (for, while, break, continue)
- •1.5.8. Возврат значения (return)
- •1.5.9. Перехват ошибок (on error)
- •1.5.10. Примеры программирования
- •2. Роль численных методов
- •2.1. Этапы решения задачи на компьютере
- •2.2. Математические модели
- •2.3. Численные методы
- •3. Методы аппроксимации и интерполирования
- •4. Лабораторная работа № 1. Интерполирование степенными многочленами
- •4.1. Постановка задачи
- •4.2. Порядок выполнения работы
- •4.3. Краткие теоретические сведения
- •4.3.1. Метод неопределенных коэффициентов
- •4.3.2. Интерполяционный многочлен Лагранжа
- •4.3.3. Интерполяционные формулы Ньютона для равностоящих узлов
- •4.4. Примеры выполнения
- •4.4.1. Интерполирование степенными многочленами с использованием метода неопределенных коэффициентов
- •4.4.2. Интерполирование степенными многочленами с использованием второй интерполяционной формулы Ньютона
- •4.5. Требования к отчету
- •4.6. Контрольные вопросы и задания
- •4.7. Задания
- •5. Лабораторная работа № 2. Параметрическая идентификация математических моделей методами аппроксимации
- •5.1. Постановка задачи
- •5.2. Порядок выполнения работы
- •5.3. Краткие теоретические сведения
- •5.3.1. Метод выбранных точек
- •5.3.2. Метод средних
- •5.3.3. Метод наименьших квадратов
- •5.4. Примеры выполнения
- •5.4.1. Аппроксимация с использованием метода выбранных точек
- •5.4.2. Аппроксимация с использованием метода средних
- •5.4.3. Аппроксимация с использованием метода наименьших квадратов
- •5.4.4. Сравнительный анализ методов аппроксимации
- •5.5. Требования к отчету
- •5.6. Контрольные вопросы и задания
- •5.7. Задания
- •6. Решение систем линейных алгебраических уравнений
- •6.1. Общие положения
- •6.2. Точные методы решения систем линейных уравнений
- •6.2.1. Метод Крамера
- •6.2.2. Метод Гаусса
- •6.2.3. Метод обращения матриц
- •7. Лабораторная работа № 3. Решение систем линейных уравнений приближенными методами
- •7.1. Постановка задачи
- •7.2. Порядок выполнения работы
- •7.3. Краткие теоретические сведения
- •7.3.1. Математическое описание реактора идеального смешения непрерывного действия
- •7.3.2. Математическое описание кинетических закономерностей химических превращений
- •7.3.4. Приближенные методы решения систем линейных уравнений
- •7.3.4.1. Метод простых итераций
- •7.3.4.2. Метод Зейделя
- •7.4. Примеры выполнения
- •7.4.1. Пример выполнения задания точным методом
- •7.4.2. Пример выполнения задания методом итераций и методом Зейделя
- •7.5. Требования к отчету
- •7.6. Контрольные вопросы и задания
- •7.7. Задания
- •8. Лабораторная работа № 4. Решение нелинейных уравнений приближенными методами
- •8.1. Постановка задачи
- •8.2. Порядок выполнения работы
- •8.3. Краткие теоретические сведения
- •8.3.1. Этапы решения нелинейного уравнения
- •8.3.4. Метод деления отрезка пополам (вилки, дихотомии)
- •8.3.5. Метод Ньютона (метод касательных)
- •8.3.6. Метод простых итераций
- •8.4. Пример выполнения задания методом итераций
- •8.5. Требования к отчету
- •8.6. Контрольные вопросы и задания
- •8.7. Задания
- •9. Лабораторная работа № 5. Решение систем нелинейных уравнений приближенными методами
- •9.1. Постановка задачи
- •9.2. Порядок выполнения работы
- •9.3. Краткие теоретические сведения
- •9.3.1. Метод Ньютона
- •9.3.2. Метод итераций
- •9.4. Примеры выполнения
- •9.4.1. Метод Ньютона
- •9.4.2. Метод итераций
- •9.5. Требования к отчету
- •9.6. Контрольные вопросы и задания
- •9.7. Задания
- •10. Лабораторная работа № 6. Решение обыкновенных дифференциальных уравнений первого порядка численными методами
- •10. 1. Постановка задачи
- •10.2. Порядок выполнения работы
- •10.3. Краткие теоретические сведения
- •10.3.1. Метод Эйлера
- •10.3.2. Модифицированный метод Эйлера
- •10.3.3. Метод Эйлера-Коши
- •10.3.4. Метод Рунге-Кутта 4-го порядка
- •10.4. Примеры выполнения
- •10.4.1. Реализация метода Эйлера в математическом редактореMathcad
- •10.4.2. Решение обыкновенного дифференциального уравнения 1-го порядка с помощью функции rkfixed
- •10.5. Требования к отчету
- •10.6. Контрольные вопросы и задания
- •10.7. Задания
- •11. Лабораторная работа № 7. Численное интегрирование
- •11. 1. Постановка задачи
- •11.2. Порядок выполнения работы
- •11.3. Краткие теоретические сведения
- •11.3. 1. Метод прямоугольников
- •11.3.2. Метод трапеций
- •11.4. Пример выполнения
- •11.5. Требования к отчету
- •11.6. Контрольные задания
- •1 Таблица 101.7. Задания
- •12. Лабораторная работа № 8. Моделирование реактора идеального вытеснения для многостадийной химической реакции с линейной кинетикой
- •12. 1. Постановка задачи
- •12.2. Порядок выполнения работы
- •12.3. Краткие теоретические сведения
- •12.3.1. Математическая модель реактора идеального вытеснения
- •12.3.2. Численное решение систем дифференциальных уравнений
- •12.4. Пример выполнения
- •12.5. Проверка расчета с помощью функции rkfixed
- •12.6. Требования к отчету
- •12.7. Контрольные вопросы и задания
- •1 Таблица 112.8. Задания
- •13. Лабораторная работа № 9. Расчет моделей процессов диффузии и теплопроводности с помощью явной разностной схемы
- •13.3.2. Решение уравнений в частных производных
- •13.3.3.Метод сетки
- •13.3.4. Явная разностная схема
- •13.3.5. Условия устойчивости явной разностной схемы
- •13.4. Пример выполнения
- •13.5. Требования к отчету
- •13.6. Контрольные вопросы и задания
- •13.7. Задания
- •Библиографический список
- •Использование
8.7. Задания
Варианты 1.1.-1.5
Тело брошено под углом к горизонту с начальной скоростью0из положения с координатойy0(рис. 70). Изменение высоты телаyпри падении описывается законом:
, (68)
дальность падения тела sпвычисляется по формуле:
, (69)
скорость в конце полета можно рассчитать по формуле:
. (70)
В(68)-(70)g=9,81 – ускорение свободного падения, м/с2;t – время падения, с;y0– начальная координата, м;0– начальная скорость, м/с;– угол к горизонту, под которым брошено тело, …;sп– дальность падения тела, м.
Рис. 70. Изменение высоты тела при падении (к вариантам 1.1.-1.5)
1.1.Используя (68), найти времяt, через которое тело упадет на землю. Исходные данные:=30,0=10 м/с,y0=5 м, 0 t2 с.
1.2.Используя (68), найти угол к горизонту, под которым было брошено тело, если известно, что0=10 м/с,y0=8 м,t=2 с, 0 рад.
1.3.Используя (69), найти угол к горизонту, под которым было брошено тело, если известно, что0=10 м/с,Sп=5 м,t=2 с, 0 рад.
1.4. По формуле (70) найти начальную скорость брошенного тела0, если известно, что=/3,п=10 м/с,t=1c, 0020 м/с.
1.5.По формуле (70) найти угол к горизонту, под которым было брошено тело, если известно, что0=10 м/с,п=14 м/с,t=1,6c, 0 рад.
Варианты 1.6.-1.8
Брусок массыmнаходится на наклонной плоскости с углом наклона(рис. 71). На него воздействует силаF. Коэффициент трения бруска о наклонную плоскостьk. Брусок находится в состоянии покоя, если выполняется следующее условие:
. (71)
В (71) g=9,81 – ускорение свободного падения, м/с2;– угол наклона, …;k– коэффициент трения бруска о наклонную плоскость;m– масса бруска, кг;F– сила воздействия, Н.
Рис. 71. Силы, действующие на брусок (к вариантам 1.6.-1.8)
1.6. Рассчитать силу воздействияF, при которой брусок остается в покое с учетом, что угол наклона плоскости=30. Исходные данные:k=0,2,m=1 кг,g=9,81 м/с2, 0 F20 Н.
1.7. Рассчитать коэффициент тренияkбруска о наклонную плоскость, при котором брусок остается в покое с учетом, что=/6,m=1 кг,g=9,81 м/с2,F=3,3 Н, 0 k2.
1.8. Рассчитать массу брускаm, при котором брусок остается в покое, если известно, что=/6 ,k=0,2,g=9,81 м/с2,F=3,32 Н, 0m2 кг.
Варианты 1.9.-1.11
Тело совершает гармонические затухающие колебания. Амплитуда колебаний А. Смещение тела описывается законом:
, (72)
,,, (73)
где А0–начальная амплитуда колебаний, м;0– начальная фаза колебаний, …;T– период колебаний, с;– коэффициент затухания;r– коэффициент трения;m– масса тела, кг;зат– циклическая частота затухающих колебаний, с-1; ,0– собственная частота свободных незатухающих колебаний, с-1;t– время, с;x– величина смещения точки, м.
1.9.Определить момент времениtп, когда тело в первый раз с момента начала колебаний вернется в положение равновесия. Исходные данные:T=15 с,А0=10 м,0=30, r=0,2,m=5 кг, 0 t10с,=3,14 рад.
1.10. Определить коэффициент затухания ,если известно, что с момента начала колебаний до момента, когда тело в первый раз вернулось в положение равновесия, прошло 6,7 с. Исходные данные:T=15 с,А0=10 м,0=/6,=3,14 рад, 00,3.
1.11. Рассчитать массу телаm,если известно, что с момента начала колебаний до момента, когда точка в первый раз вернулась в положение равновесия прошло 6,67 с. Исходные данные:T=15 с,А0=10 м,0=/6, r=0,2,=3,14 рад, 2 m6 кг.
Вариант 2.1
Точка совершает гармонические колебания с периодом T. Амплитуда колебанийА. Смещение точки описывается законом:
, (74)
где 0– начальная фаза колебаний, …;t– время, с,T– период колебаний, с;А– амплитуда колебаний, м;x– величина смещения точки, м. Определить момент времениtп, когда точка в первый раз с момента начала колебаний вернется в положение равновесия. Исходные данные:T=2 с,А=0,1 м,0=60, 0 t10с,=3,14 рад.
Варианты 2.2.-2.4
Тело, расположенное на горизонтальной плоскости и присоединенное к пружине (рис. 72), совершает гармонические колебания. Трение о поверхность равно нулю. Смещение тела описывается законом:
,, (75)
гдеА– амплитуда колебаний, м;0– начальная фаза колебаний, …;– частота колебаний,с-1;k– жесткость пружины;m – масса тела, кг;t– время, с;x– величина смещения точки, м.
Рис. 72. Колебания тела (к вариантам 2.2.-2.4)
2.2.Определить момент времениtп, когда тело в первый раз с момента начала колебаний вернется в положение равновесия (точка О). Исходные данные:А=0,5 м,0=0,k=0,09 Н/м2,m=1 кг, 0 t 15с,=3,14 рад.
2.3. Рассчитать жесткость пружиныk, если известно, что с момента начала колебаний до момента, когда тело в первый раз вернулось в положение равновесия, прошло 5 с. Исходные данные:А=1,5 м,0=60,m=0,85 кг,=3,14 рад, 0k0,2.
2.4. Рассчитать массу тела, если известно, что с момента начала колебаний до момента, когда тело в первый раз вернулась в положение равновесия, прошло 15 с. Исходные данные:А=0,3 м,0=30,k=0,075 Н/м2,=3,14 рад, 0m5 кг.
Вариант 2.5
При сложении двух одинаково направленных гармонических колебаний с одинаковой частотой
,(76)
результирующее колебание описывается законом:
,, (77)
, (78)
, (79)
где А1, А2, А– амплитуды составляющих и результирующего колебаний, м;1, 2, – начальные фазы составляющих и результирующего колебаний, , T – период колебаний, с; – частота колебаний, с-1; t – время, с; x – величина смещения точки, м.
Определить момент времени tп, когда тело в первый раз с момента начала колебаний вернется в положение равновесия. Исходные данные:А1=3 м,А2=1 м,1=30,2=60,T=10 с, 0t3с,=3,14 рад.
Варианты 2.6.-2.9
Изменение заряда на обкладках конденсатора, входящего в состав колебательного контура (рис. 73), подчиняется закону:
, (80)
,, (81)
гдеq0– амплитудное значение заряда в момент времениt=0, Кл;0 – начальная фаза колебаний, …; – коэффициент затухания,зат – циклическая частота свободных затухающих колебаний, с -1;R– электрическое сопротивление контура, Ом;L– индуктивность контура, Гн;С– емкость конденсатора, Ф;t– время, с.
Рис. 73. Колебательный контур (к вариантам 2.6.-2.9)
2.6.Определить момент времениtп, когда заряд на обкладках конденсатора равен нулю. Исходные данные:R=1 Ом,L=0,1 Гн,С=4010-6Ф,0=60,q0=510-6Кл, 0 t0,007с,=3,14 рад.
2.7.Рассчитать индуктивность колебательного контураL, если известно, что полный разряд конденсатора происходит за время, равное 0,05 с. Исходные данные:R=0,5 Ом, ,С=2010-6Ф,0=60,q0=1,510-6Кл,=3,14 рад, 0,3L0,5 Гн.
2.8. Рассчитать сопротивление колебательного контураR, если известно, что полный разряд конденсатора происходит за время равное 0,0042 с. Исходные данные:L=0,25 Гн,С=0,410–6 Ф,0=60,q0=2,5 Кл,=3,14 рад, 750R2000 Ом.
2.9. Рассчитать емкость колебательного контураC, если известно, что полный разряд конденсатора происходит за время, равное 0,5 с. Исходные данные:R=0,5 Ом,L=0,25 Гн,0=30,q0=0,510-6Кл,=3,14 рад, 3010-6С3110-6Ф.
Варианты 2.10.-2.12
На спокойную воду спущен корабль с нулевой вертикальной скоростью. Если пренебречь силами, обусловленными вязкостью воды, свободные вертикальные колебания корабля в спокойной воде описываются законом:
, (82)
где M– масса корабля,– плотность воды,S– площадь горизонтальной проекции корабля,g– ускорение свободного падения,t – время.
2.10.Определить момент времениtп, когда корабль достигает своего первоначального положения. Исходные данные:M=500 т,=1 т/м3,S=250 м2,g=9,81 м/с2, 0 t1,5с,=3,14 рад.
2.11. Определить массу корабляM, если известно, что корабль достигает своего первоначального положения за 10 с. Исходные данные:=1,03 т/м3,S=500 м2,g=9,81 м/с2,=3,14 рад, 400M500 т.
2.12. Определить площадь горизонтальной проекции корабляS, если известно, что корабль достигает своего первоначального положения за 20 с. Исходные данные:=1 т/м3,M=500 т,g=9,81 м/с2,=3,14 рад, 200S300 м2.
Варианты 3.1.-3.4
Тело массы m находится на наклонной плоскости, составляющей угол с вертикалью(рис. 74). К телу прикреплена пружина, жесткость которойС. Пружина параллельна наклонной плоскости. В начальный момент времени телу сообщается скорость0, направленная вниз по наклонной плоскости. Если начало координат взять в положении статического равновесия, уравнение движения тела запишется:
,, (83)
гдеg– ускорение свободного падения, м/с2;t– время, с;m – масса тела, кг;С– жесткость пружины;– угол плоскости с вертикалью, …;0– начальная скорость тела, м/с;x– величина смещения точки, м.
Рис. 74. Тело на наклонной плоскости (к вариантам 3.1.-3.4)
3.1.Найти моменты времениt, когда тело возвращается в начальное положение. Исходные данные:С=0,09 Н/м2,m=1,5 кг,0=0,3 м/с,g=9,81 м/с2,=45, 0 t18с,=3,14 рад.
3.2. Рассчитать жесткость пружиныС, если известно, что тело возвращается в начальное состояние за 6,37 с. Исходные данные:m=1,5 кг,0=0,3 м/с,g=9,81 м/с2,=45,=3,14 рад, 0,02С0,5 Н/м2.
3.3. Рассчитать массу телаm, если известно, что тело возвращается в начальное состояние за 6,37 с. Исходные данные:С=0,09 Н/м20=0,3 м/с,g=9,81 м/с2,=45,=3,14 рад, 0,5 m 3 кг.
3.4. Рассчитать угол , если известно, что тело возвращается в начальное состояние за 6,37 с. Исходные данные: m=1,5 кг,С=0,09 Н/м2,0=0,3 м/с,g=9,81 м/с2,=3,14 рад, 0,2рад.
Вариант 3.5
На гладкой поверхности, наклоненной к горизонту под углом , находится прикрепленный к пружине груз (рис. 75). Статическое удлинение пружины равноf. В начальный момент пружина была растянута из ненапряженного состояния на длину, равную 3 f, и груз отпущен без начальной скорости. Колебания груза в этом случае могут быть описаны законом:
, (84)
гдеg– ускорение свободного падения, м/с2;t– время, с;f – статическое удлинение пружины, м;– угол наклона плоскости к горизонту, …;x– величина смещения точки, м.
Рис. 75. Прикрепленный к пружине груз (к варианту 3.5)
Найти момент времени tп, когда груз возвращается в начальное положение. Исходные данные:f=0,5 м,g=9,81 м/с2,=60, 0 tс,=3,14 рад.
Вариант 3.6
На два вращающиеся в противоположные стороны цилиндрические шкива одинакового радиуса свободно наложен однородный стержень (рис. 76). Центры шкивов находятся на горизонтальной прямой О1О2. РасстояниеО1О2=2l. Стержень приводится в движение силами трения, развивающимися в точках касания его со шкивами. Эти силы пропорциональны давлению стержня на шкив, причем коэффициент пропорциональности (коэффициент трения) равенf. При сдвиге стержня из положения симметрии наx0при начальной скорости0=0 уравнение движения стержня:
, (85)
гдеg– ускорение свободного падения, м/с2;t– время, с;x– величина смещения точки, м;f– коэффициент пропорциональности;x0 – начальный сдвиг стержня из положения симметрии, м.
Рис. 76. Вращающие шкивы (к варианту 3.6)
Найти моменты времени t, когда стержень пройдет положение симметрии. Исходные данные:f=0,252,g=9,81 м/с2,l=0,25 м, x0=0,1 м, 0 tс,=3,14 рад.
Варианты 3.7-3.8
Имеется шнур (рис. 77), закрепленный с одного конца и приводимый в колебания с другого (точка О). Колебания произвольной точкиМ, отстоящей на расстояниихот незакрепленного конца шнура длинойl, описывается уравнением плоской стоячей волны:
, (86)
,, (87)
где А– амплитуда колебаний, м;– дополнительное отставание по фазе, которое может возникать при отражении,;– частота колебаний, с-1;T– период колебаний, с;– скорость распространения волны, м/с;t– время, с;l– длина шнура, м.
Рис. 77. Колебания точки шнура (к вариантам 3.7.-3.8)
3.7.Определить момент времениtп, когда точкаМзанимает начальное положение. Исходные данные:А=0,3 м,=30,=2100 м/с,l=5 м,х=2 м, T=2 с, 0t0,1с,=3,14 рад.
3.8. Рассчитать скорость распространения волны,если известно, что точкаМ заняла первоначальное положение через 0,086 с. Исходные данные:А=0,3 м,=30,l=5 м,х=2 м, T=2 с,=3,14 рад, 10003200 м/с.
Варианты 3.9. - 3.11
Брусок массы mнаходится на наклонной плоскости с углом наклона(рис. 78). На него воздействует силаF. Коэффициент трения бруска о наклонную плоскостьk. Брусок находится в состоянии равномерного прямолинейного движения, если выполняется следующее условие:
. (88)
В (88) g=9,81 – ускорение свободного падения, м/с2;– угол наклона, …;k– коэффициент трения бруска о наклонную плоскость;m– массы бруска, кг;F– сила воздействия, Н.
3.9.Определить какова должна быть сила воздействияF, чтобы брусок оставался в состоянии равномерного прямолинейного движения при угле наклона плоскости=30. Исходные данные:k=0,3,m=2 кг,g=9,81 м/с2, 0 F30 Н.
Рис. 78. Силы, действующие на брусок (к вариантам 3.9.-3.11)
3.10. Рассчитать коэффициент тренияkбруска о наклонную плоскость, при котором брусок остается в состоянии равномерного прямолинейного движения при угле наклона плоскости=30с учетом, чтоm=2 кг,g=9,81 м/с2,F=20,82 Н, 0,01 k 1,5.
3.11. Рассчитать массу брускаm, при котором брусок остается в состоянии равномерного прямолинейного движения, если известно, что=30,k=0,3,g=9,81 м/с2,F=20,82 Н, 0,01 m 4 кг.