Курсовая работа: Визначення динамічних похибок вимірювань
Название: Визначення динамічних похибок вимірювань Раздел: Рефераты по коммуникации и связи Тип: курсовая работа |
Міністерство освіти і науки України Вінницький національний технічний університет Інститут автоматики, електроніки та комп’ютерних систем управління Кафедра МПА ВИЗНАЧЕННЯ ДИНАМІЧНИХ ПОХИБОК ВИМІРЮВАНЬ Пояснювальна записка з дисципліни “Основи теорії похибок та обробки результатів вимірювань” до курсової роботи за спеціальністю 8.091302 “Метрологія та вимірювальна техніка” 08-03.ОТПОРВ.007.00.000 ПЗ Вінниця ВНТУ 2007 Зміст Вступ 1. Характеристики точності та правильності вимірювань 2. Практична частина Висновки Перелік посилань Додатки Додаток А. Розв’язок диференційного рівняння в пакеті Maple 7 ВСТУП При проведенні вимірювань завжди виникає перехідний процес, при якому сигнал на виході засобу вимірювання суттєво змінюється в часі. Це пояснюється інерційними властивостями засобу вимірювання, які зумовлюють виникнення динамічної похибки – складової похибки вимірювання, що виникає додатково до статичної похибки при проведенні динамічних вимірювань. Динамічна похибка представляється через динамічні характеристики та визначається як миттєва різниця значення вихідного сигналу, розрахованого за вхідним сигналом, і миттєвого значення вихідного сигналу в даний момент часу. Для визначення динамічної похибки необхідно знайти повні динамічні характеристики засобу вимірювання, а саме диференціальне рівняння, перехідну, імпульсну, амплітудно-частотну та фазочастотну характеристики, оскільки саме повні динамічні характеристики однозначно визначають зміни вихідного сигналу засобу вимірювання при будь-яких змінах у часі вхідного сигналу або впливних величин. 1. ХАРАКТЕРИСТИКИ ТОЧНОСТІ ТА ПРАВИЛЬНОСТІ ВИМІРЮВАНЬ Точність вимірювань – це характеристика її якості, що відображає близькість результату вимірювань до істинного значення вимірюваної фізичної величини. Точність виражають оберненим значенням модуля відносної похибки вимірювання
Точність засобу вимірювальної техніки характеризується його класом точності – узагальненою характеристикою, що визначається границями допустимої основної та зведеної похибок, а також іншими характеристиками, що впливають на його точність, значення яких регламентується. Засобам вимірювань з двома чи більшою кількістю діапазонів вимірювань даної фізичної величини допускається присвоювати два і більше класів точності. Засобам вимірювання, що призначені для вимірювання двох і більше фізичних величин, також допускається присвоювати різні класи точності для кожної фізичної величини. Границі допустимих основної та додаткової похибок засобів вимірювання встановлюють у формі абсолютних, зведених або відносних значень, залежно від характеру їх зв’язку з інформативним параметром вхідного чи вихідного сигналів. Границі допустимої абсолютної основної похибки встановлюють за формулами
де Границі основної зведеної похибки встановлюються за формулою
де Значення, які стоять у дужках, для засобів вимірювань, які розробляються заново, не використовуються. Нормувальне значення Для засобів вимірювання фізичної величини, для яких прийнята шкала з умовним нулем, нормувальне значення дорівнює різниці границь вимірювань. Для вимірювальних приладів з суттєво нерівномірною шкалою нормувальне значення встановлюють рівним всій довжині шкали або її частині, що відповідає діапазону вимірювань. При цьому границі абсолютної похибки встановлюють, як і довжину шкали, в одиницях довжини. Границі допустимої відносної похибки встановлюються за формулою:
де Для цифрових приладів клас точності переважно відображений двома числами, що записані через косу риску
Для проміжних показів ( Відповідно до означення зведеної похибки за класом точності (
Границі допустимої додаткової похибки засобів вимірювання можна встановити в формі , що відрізняється від форми встановлення границь допустимої основної похибки. Їх встановлюють: У вигляді сталого значення для всього діапазону значень виливної величини або сталих значень для певних інтервалів цього діапазону. Зазначенням відношення границі допустимої додаткової похибки, що відповідає регламентованому інтервалу значень впливної величини, до ширини цього інтервалу. Наведенням граничної функції впливу як залежності границь допустимої додаткової похибки від впливних величин. Наведенням функціональної залежності границь допустимих відхилень від номінальної функції впливу. Правильність вимірювання – це характеристика його якості, що відображає близькість до нуля систематичної похибки в його результаті. Якщо значення систематичної складової
Поправлення – значення фізичної величини, яке додається до результату вимірювання, щоб виключити систематичну похибку. Поправлення підсумовують із значенням міри, показом вимірювального приладу тощо. Іноді замість поправлення значення фізичної величини користуються поправляльним множником
Сукупними характеристиками точності і правильності вимірювань є їх відтворюваність і збіжність. Відтворюваність вимірювань – це близькість результатів вимірювань сталої фізичної величини отриманих у різних умовах, різними методами, засобами вимірювань, експериментаторами незалежно від місця та часу їх здійснення, тобто нерівноточних вимірювань, в яких систематична складова їх похибки стає випадковою. Збіжністю вимірювань називають їх відтворюваність в однакових умовах (рівноточні вимірювання). Збіжність результатів вимірювань відображає близькість до нуля випадкової похибки вимірювань. Збіжність може бути оцінена кількісно дисперсією результатів вимірювання. Виділяють ще термін роздільної (подільної) здатності засобу вимірювальної техніки як кількості вірогідно розрізнюваних значень вимірюваної фізичної величини, які вписуються в зону їх невизначеності в процесі вимірювань. Для такого означення терміну роздільної здатності засобу вимірювальної техніки визначають також ступінь вірогідності, з якою встановлюється ширина зони невизначеності Рисунок 1.1 – Графічна інтерпретація роздільної здатності Роздільна здатність вимірювань є теоретичною характеристикою, що залежить від точності та діапазону вимірювань. 2. ПРАКТИЧНА ЧАСТИНА Динамічне рівняння пов’язує вихідну величину
В операторній формі
Диференціальне рівняння динамічної системи є вичерпною її характеристикою, але його коефіцієнти важко піддаються експериментальному визначенню. Тому як характеристики перетворення в часовій області використовуються імпульсна перехідна (вагова) Імпульсна функція
Перехідна функція
є відгуком динамічної системи на вхідну дію у вигляді одиничної функції
З характеристиками перетворення у часовій області однозначно пов’язані характеристики перетворення в частотній області, що є наслідком дуальності часу і частоти. Усталена реакція на синусоїдний вхідний сигнал у загальному випадку є складною функцією параметрів засобу вимірювальної техніки і описується відповідними амплітудно-частотною та фазочастотною характеристиками, які можуть бути одержані з диференціального рівняння в результаті нижчеподаних математичних дій. Застосувавши до диференціального рівняння при початкових нульових умовах перетворення Лапласа, одержимо передаточну функцію
де Заміна оператора Лапласа в передаточній функції на
Комплексна частотна характеристика є вихідною для визначення амплітудно-частотної
Згідно індивідуального завдання необхідно знайти розв’язок диференціального рівняння другого порядку
Підставимо (2.13) в (2.12) і отримаємо:
Розв’язком даного рівняння буде функція
графічне зображення якої подано на рисунку 2.1. Рисунок 2.1 – Графічне представлення розв’язку диференціального рівняння Для знаходження перехідної характеристики підставимо в (2.12) як вхідний сигнал
Отримаємо розв’язок:
Графічно перехідна характеристика зображена на рисунку 2.2. Рисунок 2.2 – Перехідна характеристика Для знаходження імпульсної характеристики підставимо в (2.12) як вхідний сигнал
Отримаємо розв’язок:
Графічно імпульсна характеристика зображена на рисунку 2.3. Рисунок 2.3 – Імпульсна характеристика Знайдемо передатну функцію заданого диференціального рівняння
Замінимо оператор Лапласа в передатній функції на
Виділимо дійсну та уявну частини в знаменнику:
Помножимо чисельник та знаменник дробу на вираз, комплексно спряжений до знаменника, для того, щоб позбутись ірраціональності в знаменнику. В результаті отримаємо
З даного виразу маємо дійсну
частини комплексної частотної характеристики. Знайдемо амплітудно-частотну характеристику як корінь із суми піднесених до квадрату дійсної та уявної частин комплексної частотної характеристики:
Замінимо
Графічно амплітудно-частотну характеристику наведено на рисунку 2.4. Рисунок 2.4 – Амплітудно-частотна характеристика Знайдемо фазочастотну характеристику як мінус арктангенс відношення уявної частини комплексної частотної характеристики до дійсної
Після заміни
Графік фазочастотної характеристики наведено на рисунку 2.5. Рисунок 2.5 – Фазочастотна характеристика Отже, в даному розділі було знайдено розв’язок диференціального рівняння другого порядку, отримано перехідну, імпульсну, амплітудно-частотну та фазочастотну характеристики. Всі розв’язки отримані за допомогою пакету прикладних програм Maple 7 і наведені в додатку А. ВИСНОВКИ В даній курсовій розглянуто питання визначення динамічних похибок вимірювання за допомогою динамічних характеристик засобу вимірювання. В першому розділі розглянуто характеристики точності та правильності вимірювань, дано інтерпретацію понять роздільної здатності та класу точності засобу вимірювання, наведено методики визначення класу точності для різних типів засобів вимірювання. В другому розділі було знайдено розв’язок диференціального рівняння другого порядку, що описує залежність вихідного сигналу засобу вимірювання від вхідного, отримано перехідну, імпульсну, амплітудно-частотну та фазочастотну характеристики, оскільки саме вони як повні динамічні характеристики дозволяють визначити динамічну похибку засобу вимірювання. Всі розв’язки отримані за допомогою пакету прикладних програм Maple 7 і наведені в додатку А. ПЕРЕЛІК ПОСИЛАНЬ 1. Поліщук Є.С., Дорожовець М.М., Яцук В.О., та ін. Метрологія та вимірювальна техніка: Підручник / Є.С.Поліщук, М.М.Дорожовець, В.О.Яцук, В.М.Ванько, Т.Г.Бойко; За ред. проф. Є.С.Поліщука. – Львів: Видавництво “Бескид Біт”, 2003. 2. ДСТУ 2681-94. Метрологія. Терміни та визначення. – К.: Держстандарт України, 1994. 3. Володарський Є.Т., Кухарчук В.В, Поджаренко В.О., Сердюк Г.Б. Метрологічне забезпечення вимірювань і контролю. Навчальний посібник. – Вінниця ВДТУ, 2001. 4. Кухарчук В.В., Кучерук В.Ю., Долгополов В.П., Грумінська Л.В. Метрологія та вимірювальна техніка. Навчальний посібник. – Вінниця: УНІВЕРСУМ-Вінниця, 2004. Додаток А Розв’язок диференційного рівняння в пакеті Maple 7 x(t):=cos(t); ode1:=1.1*diff(diff(y(t),t),t)+2.5*diff(y(t),t)+3.1*sin(2.1*t)= =0.5*x(t); dsolve({ode1,y(0)=0,D(y)(0)=0}); y1(t) := 34100/115861*sin(21/10*t)+775000/2433081*cos(21/10*t)- -55/746*sin(t)-125/746*cos(t)+103507151/432161530*exp(-25/11*t)- -41/105; with(plots):plot({cos(t),y1(t)},t=0..50); dsolve({1.1*diff(diff(y(t),t),t)+2.5*diff(y(t),t)+ +3.1*sin(2.1*t)=0.5*Heaviside(t),y(0)=0, D(y)(0)=0}); y2(t) := 34100/115861*sin(21/10*t)+775000/2433081*cos(21/10*t)+ +1/5*Heaviside(t)*t+11/125*exp(-25/11*t)*Heaviside(t)- -11/125*Heaviside(t)+157542/579305*exp(-25/11*t)-62/105; plot({y2(t)},t=0..20); dsolve({1.1*diff(diff(y(t),t),t)+2.5*diff(y(t),t)+ +3.1*sin(2.1*t)=0.5*Dirac(t),y(0)=0, D(y)(0)=0}); y3(t) := -1/5*exp(-25/11*t)*Heaviside(t)+1/5*Heaviside(t)+ +34100/115861*sin(21/10*t)+775000/2433081*cos(21/10*t)+ +273403/579305*exp(-25/11*t)-83/105; plot({y3(t)},t=0..20); S(f):=(4.41-(2*Pi*f)^2)/(2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+6.51- -5*sqrt(-1)*(2*Pi*f)^3+22.05*2*Pi*f*sqrt(-1)); S(f):=(4.41-(2*Pi*f)^2)*((2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+6.51) + sqrt(-1)*(5*(2*Pi*f)^3+22.05*2*Pi*f))/ ((2.2*(2*Pi*f)^4- -9.3*(2*Pi*f)^2+6.51)^2 -( sqrt(-1)*(5*(2*Pi*f)^3+ +22.05*2*Pi*f))^2); A(f):=(4.41-(2*Pi*f)^2)*(2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+ +6.51) / ((2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+6.51)^2 -( sqrt(- -1)*(5*(2*Pi*f)^3+22.05*2*Pi*f))^2); B(f):=(4.41-(2*Pi*f)^2)*(5*(2*Pi*f)^3+ +22.05*2*Pi*f) / ((2.2*(2*Pi*f)^4-9.3*(2*Pi*f)^2+6.51)^2 – -( sqrt(-1)*(5*(2*Pi*f)^3+22.05*2*Pi*f))^2); K(f):=sqrt( (A(f))^2+(B(f))^2 ); plot(K(f),f=0..1); Q(f):=-arctan((B(f)/A(f))); plot(Q(f),f=0..1); |