Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теплоэнергетические процессы и установки. Общая энергетика

.pdf
Скачиваний:
0
Добавлен:
29.11.2025
Размер:
4.17 Mб
Скачать

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

f (x)

 

1

 

exp x m 2

/ 2 n2 ,

(1.6)

 

 

 

n

 

2

 

 

 

 

 

где f(x) – плотность вероятности измеряемых значений величины; m математическое ожидание, являющееся наиболее вероятным значением измеряемой величины х и представляющее собой

среднее арифметическое значение:

 

n

 

/ n ;

m

x

 

 

i

 

i 1

 

 

n среднее квадратическое отклонение измеряемой величины;

 

 

 

n

x

m 2

n

 

i

.

 

 

 

n 1

 

 

 

i 1

 

 

На рис. 1.1 представлена нормальная кривая распределения случайной величины, где по оси абсцисс отложены результаты измерений, а по оси ординат – плотность вероятности их появления. Площадь под кривой, соответствующая какому-либо интервалу по оси абсцисс, представляет собой вероятность Р попадания случайного результата измерения в этот интервал. После интегрирования (1.6) найдем

 

 

 

 

x

x m 2

 

P

 

2

 

2 n2 dx .

 

 

 

e

(1.7)

n

 

 

 

 

2 0

 

 

11

Рис. 1.1. Нормальное распределение погрешностей

Если взять интервал по оси абсцисс в долях от среднего квадратического отклонения, т. е. ввести новую переменную (х m)/ n = t, то выражение (1.7) примет вид

 

 

 

 

t

P

 

2

 

e t 2 / 2dt 2Ф t ,

n

 

 

 

 

2 0

где Ф(t) – так называемая функция нормального распределения, для которой составлены таблицы.

Таким образом, на участке т п, т + п размещаются 0,6228 всех произведенных измерений. В более широких границах, например т 2 п, т + 2 п размещаются уже 0,9546 всех измерений, а из «трехсигмовых» границ выходят лишь 0,0028 измерений.

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

12

Конечная цель анализа выполненных измерений состоит в определении погрешности среднего арифметического значения

 

n

 

/ n.

x

x

 

 

i

 

i 1

 

 

Согласно теории погрешностей оценкой точности измерения среднего арифметического значения, принимаемого за истинное значение измеряемой величины, является среднее квадратическое отклонение

 

 

1

n

m n /

n

xi x .

 

 

 

 

n(n 1) i 1

Таким образом, среднее квадратическое отклонение оценки сред-

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

Так как нормальный закон справедлив при бесконечном числе измерений (практически при n > 200), то для оценки доверительного интервала пользуются распределением Стьюдента, учитывающим влияние конечного числа измерений на величину доверительного интервала (при n распределение Стьюдента сходится с нормальным).

Границы доверительного интервала для заданного значения доверительной вероятности при ограниченном числе наблюдений записываются в виде

x t m x x t m ,

(1.8)

где t m – коэффициент Стьюдента, значения которого в зависи-

мости от числа измерений (п – 1) и доверительной вероятности приведены в табл.1.1.

13

Таблица 1.1 Значения коэффициента Стьюдента в зависимости от числа

измерений n – 1 и доверительной вероятности

п – 1

 

t для значений

 

0,683

0,9

0,95

0,997

 

1

1 833

6,314

12,710

234,8

2

1 283

2920

4,303

18,72

3

1 197

2,353

3,182

9,005

4

1 142

2,132

2,776

6,485

5

1,110

2,015

2,571

5,404

6

1,089

1,943

2,447

4,819

7

1,075

1,845

2,365

4,455

8

1,066

1,859

2,306

4,209

9

1,058

1,833

2,262

4,032

10

1,052

1,812

2,228

3,898

12

1,042

1,782

2,179

3,711

14

1,036

1,761

2,145

3,586

16

1,031

1,746

2,120

3,496

18

1,027

1,734

2,101

3,430

20

1,024

1,725

2,086

3,378

30

1,016

1,697

2,042

3,230

Следует иметь в виду, что, как было упомянуто выше, измерения, содержащие грубые промахи, должны быть исключены, так как не заслуживают доверия. На практике отбрасывают результаты из-

мерений, погрешности которых превышают значения (4–3) m. Случайная составляющая погрешностей косвенных измерений

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

14

 

n

f

2

 

y

 

 

,

 

 

xi

m

 

i 1

i

 

где mi – средние квадратические отклонения результатов измере-

ний величин xi .

Границы доверительного интервала при заданной доверительной вероятности определяются по формуле (1.8). Таким образом, обозначая для заданной доверительной вероятности погрешности из-

мерений величин xi через i , получим

 

 

n

f

2

 

y

 

 

(1.9)

 

 

xi

i .

 

i 1

 

 

Погрешности в (1.9) выражаются в тех же единицах, что и величина у. Если непосредственно измеряемые величины хi являются разно-

родными, то пользуются относительными погрешностями.

1.7. Оценка погрешностей измерительных систем при технических измерениях

Для технических измерений используют измерительные системы, как правило, состоящие из первичного преобразователя (датчика), промежуточного преобразователя (или линии связи) и измерительного прибора (или аналого-цифрового преобразователя перед входом в ПЭВМ). Для оценки погрешностей измерительных систем при технических измерениях используют два метода. В первом методе, аналогичном описанному в п. 1.5, погрешности измерительной системы оценивают пределами допустимых основной и дополнительной погрешностей измерений. Однако нельзя ожидать, что описанный в п. 1.5 неблагоприятный случай будет часто встречаться. Поэтому погрешности измерительной системы следует находить из суммы квадратов пределов допустимых погрешностей составляющих элементов:

15

 

n

f

2

 

y

 

 

(1.10)

 

 

xi

xi .

 

i 1

 

 

Второй метод (вероятностно-статистический) позволяет получить оценки, близкие к действительным. Для того, чтобы воспользоваться этим методом, необходимо знать статистические характеристики систематической и случайной составляющих погрешностей.

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

погрешностей xи.с будет равно сумме систематических составляющих погрешностей отдельных преобразователей:

n xи.с i .

i 1

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

n

и.с Di ,

i 1

где Di = i2 – дисперсии случайных составляющих средств измерений.

16

2. СТАТИСТИЧЕСКИЕ МЕТОДЫ ОБРАБОТКИ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ

2.1.Средние арифметические

исредние квадратические отклонения

Значение искомого параметра, вычисленное на основе ограниченного числа опытов, всегда будет содержать элемент случайности. Такое приближенное значение называется оценкой параметра. Например, как было показано в п. 1.6, оценкой для математического ожидания служит среднее арифметическое x наблюдавшихся значений случайной величины xi в n независимых опытах:

 

n

 

/ n .

x

x

 

 

i

 

i 1

 

 

Оценкой для среднего квадратического отклонения является ста-

тистическое среднее квадратическое отклонение s (стандартное отклонение), представляющее меру ширины гистограммы – кривой распределения, построенной по экспериментальным данным:

 

1

n

xi

x 2 .

s

 

 

 

n 1 i 1

 

 

При очень большом числе опытов x

и s будут с большой веро-

ятностью весьма близки к m и (оценки, обладающие такими свойствами, называются состоятельными). Кроме того, желательно, чтобы, пользуясь величинами x вместо m и s вместо , мы не делали систематических ошибок в сторону завышения или занижения (такие оценки называются несмещенными). Наконец, выбранные несмещенные оценки должны обладать по сравнению с другими оценками минимальным средним квадратическим отклонением. Оценки, обладающие таким свойством, называются эффективны-

ми. В связи с этим Гаусс предложил метод наименьших квадратов

(точнее, минимума суммы квадратических ошибок). Следуя Гауссу,

17

определим vi xi m и обозначим квадратными скобками сумму по i от 1 до п, например

n

n

2 .

v v

, vv v

i

i

 

i 1

i 1

 

Для того, чтобы оценка для m была эффективной, значениеvv xx 2m x nm2 должно быть минимальным. Минимум

этого выражения относительно m получаем путем приравнивания к нулю первой производной:

2 x 2nm 0 ,

откуда

m x1 x2 ... xn / n x .

Следовательно, смысл среднего арифметического заключается в том, что оно позволяет свести к минимуму стандартное отклонение.

2.2. Закон распределения ошибок Гаусса

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

Сначала решим эту задачу для функции одной переменной у = f (х). Для величины х измеряемые значения х1, х2,…, хп дают среднее арифметическое х = [х] /п. Для каждого xi вычисляется соответствующее значение функции уi = f (хi). Так как xi x vi , то, полагая vi достаточно малым, получим

yi

f (x vi ) f x vi df / dx

 

 

 

 

 

 

y vi df / dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y

x, y .

 

 

 

 

 

 

 

 

 

 

 

18

Далее получаем

yi yi y vi df / dx x, y ,

y y df / dx 2 vv .

Следовательно:

sy df / dx sx .

Вслучае большого числа переменных вычисления совершенно аналогичны. Пусть z = f (х, у,…) и наблюдаемыми величинами яв-

ляются xi, yi . Тогда

zi f xi , yi ,... f x vi , y ui ,...

f x, y,... vi df / dx x ,y ui df / dy x, y ...,

откуда следует, что

z zi z vi df / dx ui df / fy ...

Пренебрегая удвоенными произведениями и членами более высокого порядка, получим

z z df / dx 2 vv df / dy 2 uu ....

Откуда следует закон распространения ошибок:

sz fx2sx2 f y2sy2 ... ,

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

19

2.3.Взвешенное среднее

Вэкспериментальных исследованиях при определении среднего арифметического часто приходится иметь дело с результатами измерений различной точности. Для определения погрешностей изме-

рений в этом случае вводятся веса wi так, чтобы измерениям большей точности соответствовали большие веса. Тогда среднее арифметическое сформулируется в виде

 

n

 

 

 

 

 

wi xi

 

wx

 

wx ,

x

i 1

 

 

n

w

 

 

 

w0

 

wi

 

 

 

 

 

i 1

 

 

 

 

где значение w0 считается постоянным. Согласно закону распределения ошибок Гаусса, стандартное отклонение среднего арифметического

s2

 

1

w2s2

w2s2

... w2s2

.

 

x

 

w02

1 1

2 2

n n

 

 

 

 

 

 

 

Веса выбираются таким образом, чтобы sx2 было минимальным. Из последнего уравнения следует, что

sx2

1

w12s12 w22s22 ... w0 w1 ...wn 1 2 sn2 .

w02

 

 

Это выражение становится минимальным, если sx2 / wi 0 , т. е.

2w s2

2w s2

0

или w

w

1

s2

/ 1

s2

.

(2.1)

i i

n n

 

i

n

 

i

 

n

 

 

Следовательно, веса должны быть обратно пропорциональны квадратам стандартных отклонений.

20