Все входные и выходные данные в заданиях этой группы являются вещественными числами.
Begin1°. Дана сторона квадрата a. Найти его периметр P = 4·a.
Begin2°. Дана сторона квадрата a. Найти его площадь S = a2.
Begin3°. Даны стороны прямоугольника a и b. Найти его площадь S = a·b и периметр P = 2·(a + b).
Begin4°. Дан диаметр окружности d. Найти ее длину L = π·d. В качестве значения π использовать 3.14.
Begin5°. Дана длина ребра куба a. Найти объем куба V = a3 и площадь его поверхности S = 6·a2.
Begin6°. Даны длины ребер a, b, c прямоугольного параллелепипеда. Найти его объем V = a·b·c и площадь поверхности S = 2·(a·b + b·c + a·c).
Begin7°. Найти длину окружности L и площадь круга S заданного радиуса R:
L = 2·π·R, S = π·R2.
В качестве значения π использовать 3.14.
Begin8°. Даны два числа a и b. Найти их среднее арифметическое: (a + b)/2.
Begin9°. Даны два неотрицательных числа a и b. Найти их среднее геометрическое, т. е. квадратный корень из их произведения: (a·b)1/2.
Begin10°. Даны два ненулевых числа. Найти сумму, разность, произведение и частное их квадратов.
Begin11°. Даны два ненулевых числа. Найти сумму, разность, произведение и частное их модулей.
Begin12°. Даны катеты прямоугольного треугольника a и b. Найти его гипотенузу c и периметр P:
c = (a2 + b2)1/2, P = a + b + c.
Begin13°. Даны два круга с общим центром и радиусами R1 и R2 (R1 > R2). Найти площади этих кругов S1 и S2, а также площадь S3 кольца, внешний радиус которого равен R1, а внутренний радиус равен R2:
S1 = π·(R1)2, S2 = π·(R2)2, S3 = S1 − S2.
В качестве значения π использовать 3.14.
Begin14°. Дана длина L окружности. Найти ее радиус R и площадь S круга, ограниченного этой окружностью, учитывая, что L = 2·π·R, S = π·R2. В качестве значения π использовать 3.14.
Begin15°. Дана площадь S круга. Найти его диаметр D и длину L окружности, ограничивающей этот круг, учитывая, что L = π·D, S = π·D2/4. В качестве значения π использовать 3.14.
Begin16°. Найти расстояние между двумя точками с заданными координатами x1 и x2 на числовой оси: |x2 − x1|.
Begin17°. Даны три точки A, B, C на числовой оси. Найти длины отрезков AC и BC и их сумму.
Begin18°. Даны три точки A, B, C на числовой оси. Точка C расположена между точками A и B. Найти произведение длин отрезков AC и BC.
Begin19°. Даны координаты двух противоположных вершин прямоугольника: (x1, y1), (x2, y2). Стороны прямоугольника параллельны осям координат. Найти периметр и площадь данного прямоугольника.
Begin20°. Найти расстояние между двумя точками с заданными координатами (x1, y1) и (x2, y2) на плоскости. Расстояние вычисляется по формуле
((x2 − x1)2 + (y2 − y1)2)1/2.
Begin21°. Даны координаты трех вершин треугольника: (x1, y1), (x2, y2), (x3, y3). Найти его периметр и площадь, используя формулу для расстояния между двумя точками на плоскости (см. задание Begin20). Для нахождения площади треугольника со сторонами a, b, c использовать формулу Герона:
S = (p·(p − a)·(p − b)·(p − c))1/2,
где p = (a + b + c)/2 — полупериметр.
Begin22°. Поменять местами содержимое переменных A и B и вывести новые значения A и B.
Begin23°. Даны переменные A, B, C. Изменить их значения, переместив содержимое A в B, B — в C, C — в A, и вывести новые значения переменных A, B, C.
Begin24°. Даны переменные A, B, C. Изменить их значения, переместив содержимое A в C, C — в B, B — в A, и вывести новые значения переменных A, B, C.
Begin25°. Найти значение функции y = 3x6 − 6x2 − 7 при данном значении x.
Begin26°. Найти значение функции y = 4(x−3)6 − 7(x−3)3 + 2 при данном значении x.
Begin27°. Дано число A. Вычислить A8, используя вспомогательную переменную и три операции умножения. Для этого последовательно находить A2, A4, A8. Вывести все найденные степени числа A.
Begin28°. Дано число A. Вычислить A15, используя две вспомогательные переменные и пять операций умножения. Для этого последовательно находить A2, A3, A5, A10, A15. Вывести все найденные степени числа A.
Begin29°. Дано значение угла α в градусах (0 ≤ α < 360). Определить значение этого же угла в радианах, учитывая, что 180° = π радианов. В качестве значения π использовать 3.14.
Begin30°. Дано значение угла α в радианах (0 ≤ α < 2·π). Определить значение этого же угла в градусах, учитывая, что 180° = π радианов. В качестве значения π использовать 3.14.
∙1.10. Содержание отчета
1. название и цель работы
2. номер варианта для выполнения задания и условия своего варианта
3. блок-схемы решения задачи
4. тексты программ
5. полученные при расчетах численные результаты
6. письменные ответы на контрольные вопросы
1.11 Контрольные вопросы
1.11.1 Каков порядок создания программы в интегрированной среде ABC-Pascal?
1.11.2 Какие основные функции выполняет система ABC-Pascal (по главному меню)?
1.11.3 Какие операции позволяет выполнять текстовый редактор ABC-Паскаля при подготовке программы? Примеры операций редактирования?
1.11.4 Как выполнить компиляцию программы с сохранением исполняемого модуля?
1.11.5 Как сохранить программу под другим именем?
1.11.6 Как открыть окно?
1.11.7 Какими способами можно выйти из среды?
1.11.8 Как создать исполняемый модуль?
1.11.9 Как вызвать контекстную помощь?
1.11.10 Какова структура Паскаль-программы?
1.11.11 Какие операторы используются для ввода (вывода) значений переменных?
1. найти корреляционную функцию сигнала RX(τ), определить дисперсию сигнала σX2 и стандартное отклонение σX, построить график корреляционной функции и объяснить характер ее изменения при изменении параметров энергетического спектра сигнала;
2. найти вероятность того, что реализация случайного сигнала через промежуток времени t0 после прохождения через нуль примет значение, большее σX;
3. определить шаг дискретизации h, достаточный для восстановления дискретизированного сигнала X(ih) c помощью двух первых полиномов Лежандра (линейная аппроксимация) с собственной средней квадратической погрешностью, не превышающей 1% от стандартного отклонения сигнала X(t);
4. найти значение коэффициента корреляции между сечениями сигнала, разделенными шагом дискретизации.
Исходные данные (вариант 18):
Энергетический спектр сигнала Х(t) - .
Варианты задания №1
№№
№№
0,01 с
2,0 с
5 мс
10 мс
2,5 с
0,35 c
70 мс
0,7 mc
3 с
1,2 mc
0,15 с
10 mc
7 c
20 mc
0,5 mc
3,5 c
6 c
7 mc
0,4 mc
23
8 mc
0,2 c
0,1 c
1 c
7 mc
2 c
0,5 c
0,35 c
0,1 c
0,24 c
4,5 mc
20 c
0,08 c
14 mc
6 mc
Решение:
1. Нахождение корреляционной функции сигнала RX(τ), определение дисперсии сигнала σX2 и стандартного отклонения σX, построение графика корреляционной функции и объяснение характера ее изменения при изменении параметров энергетического спектра сигнала.
Случайной называется такая функция, которая может принять тот или иной конкретный вид, неизвестно заранее какой именно. Этот вид называется реализацией. В данном случае случайная функция зависит от времени.
Под стационарным понимается такой случайный процесс, n – мерная плотность распределения вероятности которого не меняется при любом сдвиге всей группы точек вдоль оси времени.
Стационарный случайный процесс называется эргодическим, если любая его реализация имеет одни и те же статистические свойства. Следовательно, для эргодических случайных процессов любая статистическая характеристика, получаемая усреднением по множеству реализаций, может быть получена также усреднением за достаточно большой промежуток времени по одной реализации случайного процесса. Т.о. свойство эргодичности позволяет рассматривать не множество, а всего одну реализацию.
Случайная величина X(t) подчинена нормальному закону распределения, если плотность вероятности ее имеет вид
, (1.1)
где а – математическое ожидание X(t)
σ – стандартное отклонение X(t), σ>0.
Если X(t) центрированная величина, то а=0 и плотность распределения вероятности ее будет иметь вид
(1.2)
Функция распределения X(t) имеет следующий вид
(1.3)
Одномерная плотность вероятности и связанные с ней числовые характеристики позволяют получить важную информацию о свойствах случайного процесса. Однако для решения многих задач таких сведений не достаточно, так как они дают вероятностное представление о случайном процессе X(t) только в отдельные моменты времени и не определяют его изменение во времени. Для описания временных характеристик случайного процесса необходимо использовать корреляционную функцию и спектральных характеристики случайного процесса.
Начальный корреляционный момент двух значений случайной функции X(t) и X(t1), взятых в моменты времени t и t1 носит название корреляционной функции. Т.е. корреляцононная функция показывает степень сходства между сигналом и его сдвинутой копией.
В случае стационарности процесса корреляционная функция не будет зависеть от текущего значения времени t, и будет определяться только временным сдвигом . С учетом эргодичности стационарного процесса корреляционной функцией можно назвать среднее по времени от произведения X(t) и X(t+τ)
(1.4)
Для стационарного процесса корреляционная функция определяет зависимость случайной величины x в последующий момент времени (t+ τ) от предшествующего значения в момент t.
Корреляционная функция стационарного процесса является четной функцией. При τ=0 значение корреляционной функции является ее наибольшим значением. Значение корреляционной функции чаще всего будет тем меньше, чем больше промежутки времени τ , так как связь между далеко отстоящими друг от друга значениями х будет обычно слабее. Чем менее инерционен (более подвижен) объект наблюдения, тем быстрее убывает R(τ) c увеличение τ.
Задан энергетический спектр сигнала Х(t) (рис.1) - Он также носит название спектральной плотности сигнала и является четной функцией своего аргумента. Обладает важным свойством: интегрирование спектральной плотности по всем частотам от дает средний квадрат исходной функции времени Х(t). По своему физическому смыслу спектральная плотность есть величина, которая пропорциональна средней мощности процесса в интервале частот от ω до ω+dω.
Рисунок 1 Энергетический спектр сигнала X(t)
Спектральная плотность и корреляционная функция случайных процессов представляют собой взаимные преобразования Фурье, т.е. они связаны интегральными зависимостями. Т.к. спектральная плотность и корреляционная функция представляют собой четные вещественные функции, то их можно представить следующими соотношениями
(1.5)
Результаты интегрирования для наиболее распространенных случаев приведены в специальных таблицах, из которой имеем
Корреляционная функция (в соответствии с таблицей для определенного вида энергетического спектра)
Таким образом, график корреляционной функции имеет следующий вид
Рисунок 2 График корреляционной функции
Рассмотрим характер изменения корреляционной функции сигнала при изменении параметров энергетического спектра, а именно, коэффициента, определяющего ширину полосы частот, и стандартного отклонения по сравнению с полученной функцией. Наглядный результат представлении ниже на рисунках 3 и 4.
- полученная функция
- функции с измененным стандартным отклонением
- функции с измененным коэффициентом, определяющим ширину полосы частот
Рисунок 3 Изменение стандартного отклонения
Рисунок 4 Изменение коэффициента, определяющего ширину полосы частот
Связь между спектральной плотностью и видом функции времени Х(t) заключается том, что чем «уже» график спектральной плотности, т.е. чем меньшие частоты представлены в спектральной плотности, тем медленнее изменяется величина Х во времени. Связь между видом спектральной плотности и видом функции времени получается обратной по сравнению со связью между корреляционной функцией и самим процессом (чем менее инерционен (более подвижен) объект наблюдения, тем быстрее убывает R(τ) c увеличение τ). Отсюда вытекает, что более «узкому» графику спектральной плотности соответствует более «широкий» график корреляционной функции и наоборот.
2.Нахождение вероятности того, что реализация случайного сигнала через промежуток времени t0 после прохождения через нуль примет значение, большее σX .
Как уже было сказано ранее, случайной называется такая функция, которая может принять тот или иной конкретный вид, неизвестно заранее какой именно. Этот вид называется реализацией. В практике анализа случайных сигналов применение операции сдвига во времени позволяет получить два сечения случайного процесса в произвольные моменты времени, а именно, в нашем случае имеем два сечения случайного сигнала X(t) в моменты времени , где .
Совокупность этих двух сечений образуют двумерную случайную величину , которая описывается двумерной плотностью вероятности (совместной плотностью распределения значений сигнала X(t)) , которая имеет следующий вид
, (1.6)
Где ,
- значение нормированной корреляционной функции RX(τ) при τ=t0 .
Условная плотность распределения значений сигнала X(t) в момент при условии, что в сечении t1 сигнал имел значение , находится из следующего выражения
(1.7)
где - плотность распределения вероятности значений сигнала X(t) в первом сечении, находится в соответствии с нормальным законом распределения по формуле (1.2)
(1.8)
Найдем условную плотность распределения значений сигнала X(t) в момент . Для этого в выражение (1.7) подставим выражения (1.6) и (1.8)
(1.9)
Далее необходимо найти вероятность того, что реализация случайного сигнала через промежуток времени t0 после прохождения через нуль примет значение, большее σX. Для этого воспользуемся связью интегральной функции распределения F(x) (формула (1.3)) с функцией Лапласа Ф (функция верхнего предела, значения которой табулированы) для нормального распределения
(1.10)
Т.к. , то, пользуясь выражением (1.10) и найденной условной плотностью распределения (1.9), найдем искомую вероятность
Таким образом, вероятность того, что реализация случайного сигнала через промежуток времени t0 после прохождения через нуль примет значение, большее σX, равна
.
3. Определение шага дискретизации h, достаточного для восстановления дискретизированного сигнала X(ih) c помощью двух первых полиномов Лежандра (линейная аппроксимация) с собственной средней квадратической погрешностью, не превышающей 1% от стандартного отклонения сигнала X(t).
Наиболее известными специальными функциями, широко используемыми при разложении непрерывного сигнала в обобщенный ряд Фурье, являются полиномы Лежандра, которые имеют вид
По условию, для восстановления дискретизированного сигнала X(ih) необходимо воспользоваться двумя первыми полиномами Лежандра
(1.12)
Полномы Лежандра ортонормированны на интервале [-1;1] с весом ρ(t)=1. Модель экспериментального процесса X(t), построенная как ряд Фурье по полиномам Лежандра, имеет вид
(1.13)
Коэффициенты Фурье-Лежандра вычисляют по формуле:
(1.14)
Модель в виде ряда Фурье по ортогональным функциям для экспериментальной функции X(t) должна строиться на интервале ортонормированности базисных функций. Это значит, что если экспериментальная реализация X(t) определена на интервал [0,T], то чтобы построить ее модель в виде ряда Фурье по полиномам Лежандра нужно перейти от этого интервала к интервалу [-1;1] путем замены переменных:
(1.15)
Ряд Фурье – Лежандра примет вид:
(1.16)
Коэффициенты Фурье-Лежандра вычисляют по формуле:
(1.17)
Часто экспериментальная реализация представляет собой множество дискретных значений , снимаемых в течении определенного интервала времени Т, в точках , где i – номер отсчета, i=0,1,…N.
N – количество отсчетов или длина реализации, h – интервал между соседними отсчетами или шаг дискретизации, который вычисляется по формуле
(1.18)
где Т – длительность интервала наблюдения.
Для такого дискретизированного сигнала:
(1.19)
Тогда ряд Фурье – Лежандра также представляет собой набор дискретных значений, каждое из которых можно получить по формуле
(1.20)
Коэффициенты Фурье-Лежандра вычисляют по формуле:
(1.21)
Преимущество обобщенных рядов Фурье состоит не только в том, что они применяются при моделировании сложных непериодических процессов, но и в том, что для достижения достаточно хорошего приближения к исходной функции можно ограничиться малым числом слагаемых ряда m, то есть использовать полиномы низших порядков.
Точность моделирования определяется среднеквадратической погрешностью аппроксимации
(1.22)
Пользуясь выражениями (1.20), (1.21) и (1.12) запишем ряд Фурье - Лежандра для нашего дискретизированного сигнала X(ih)
За xi принимаем дискретизированный энергетический спектр сигнала
Тогда (1.23)
Выражение (1.23) определяется количеством отсчетов N, которое напрямую связано с шагом дискретизации h и длительностью интервала наблюдения соотношением (1.18) . По условию, необходимо найти такой шаг дискретизации h, при котором среднеквадратическая погрешность аппроксимации (соотношение (1.22)) не превышала бы .
Для этого смоделируем выражение (1.23) и (1.22) в программе MathCAD с возможностью варьирования вышеуказанными параметрами:
Рассматривая энергетический спектр сигнала (рис.1), можно увидеть, что оба конца графика уходят в бесконечность, следовательно, при расчете шага дискретизации необходимо было подобрать такою длительность интервала наблюдения, чтобы среднеквадратическая погрешность аппроксимации не превышала порядка заданного значения. С помощью построенной модели было установлено, что при T=14000 порядок погрешности соответствует заданному. Далее с помощью той же модели был установлен оптимальный шаг дискретизации, при котором среднеквадратическая погрешность аппроксимации не превышает заданного значения, а дальнейшее уменьшение шага дискретизации не приводит к значительному уменьшению погрешности.
Таким образом, полученный шаг дискретизации равен при средней квадратической погрешности аппроксимации .
4. Нахождение значения коэффициента корреляции между сечениями сигнала, разделенными шагом дискретизации.
Коэффициент корреляции – мера статистической связи между случайными величинами, а в нашем случае между сечениями сигнала, разделенными шагом дискретизации. В п.1 была найдена корреляционная функция случайного сигнала X(t):
, в которой .
Коэффициент корреляции или нормированная корреляционная функция находится по формуле
.
23
вдоль оси времени.
, (1.1)
(1.2)
(1.3)
. С учетом эргодичности стационарного процесса корреляционной функцией можно назвать среднее по времени от произведения X(t) и X(t+τ)
(1.4)
Он также носит название спектральной плотности сигнала и является четной функцией своего аргумента. Обладает важным свойством: интегрирование спектральной плотности по всем частотам от 
дает средний квадрат исходной функции времени Х(t). По своему физическому смыслу спектральная плотность есть величина, которая пропорциональна средней мощности процесса в интервале частот от ω до ω+dω.
(1.5)
- полученная функция
- функции с измененным стандартным отклонением
- функции с измененным коэффициентом, определяющим ширину полосы частот
, где 
.
, которая описывается двумерной плотностью вероятности (совместной плотностью распределения значений сигнала X(t)) 
, которая имеет следующий вид
, (1.6)
,
- значение нормированной корреляционной функции RX(τ) при τ=t0 
.
при условии, что в сечении t1 сигнал имел значение 
, находится из следующего выражения
(1.7)
- плотность распределения вероятности значений сигнала X(t) в первом сечении, находится в соответствии с нормальным законом распределения по формуле (1.2)
(1.8)
(1.9)
(1.10)
, то, пользуясь выражением (1.10) и найденной условной плотностью распределения (1.9), найдем искомую вероятность
.
(1.12)
(1.13)
(1.14)
(1.15)
(1.16)
(1.17)
, снимаемых в течении определенного интервала времени Т, в точках 
, где i – номер отсчета, i=0,1,…N.
(1.18)
(1.19)
(1.20)
(1.21)
(1.22)
(1.23)
.
при средней квадратической погрешности аппроксимации 
.
(1.24)