Реферати українською » Информатика, программирование » Визначення мольной теплоємності методом інтерполяції


Реферат Визначення мольной теплоємності методом інтерполяції

Федеральне агентство за освітою

Федеральне державне освітнє установа

Вищої професійної освіти

">СИБИРСКИЙФЕДЕРАЛЬНЫЙ УНІВЕРСИТЕТ"

Інститут Космічних та інформаційних технологій

Кафедра системи штучного інтелекту

>КУРСОВАЯ РОБОТА

Тема:ОПРЕДЕЛЕНИЕМОЛЬНОЙТЕПЛОЕМКОСТИ МЕТОДОМИНТЕРПОЛЯЦИИ

Красноярськ, 2009


Зміст

1. Мета і завдання курсової роботи

2. Теоретичні основи курсової роботи

3. Масив вихідних даних

4. Математичні моделі, застосовувані для розрахунків

5. Результати підрахунків, котрі відрекомендовуються як таблиць і графіків

6. Текст програми

Висновок

Список літератури


1. Мета і завдання курсової роботи

Мета курсової роботи: закріплення навичок роботи із мовою високого рівня Сі, вміння писати цією мовою програми рішення технічних завдань (визначеннямольной теплоємності кисню, з допомогою методу інтерполяції).

Завдання: визначення приблизних значень теплоємності при високих температурах від 0 0 З до 1500 0із крокомDt=10 0 З, методами інтерполяції, що дозволяє дізнатися приблизні значення функції в проміжних точках.

2. Теоретичні основи курсової роботи

Результати експериментів найчастіше є таблицю наступного виду:

X

де Х - це то, можливо, наприклад, час, аf(X) швидкість чи, як у нашому прикладі Х – це температура, аf(X) це теплоємність.

З цієї таблиці, наприклад, відомі значення функціїf(X) в точках x0 і x1, але нічого не знаємо про значення, наприклад, у точці , проте, існують методи, дозволяють дізнатися приблизні значення функції в проміжних точках. До таких методам ставляться методи інтерполяції.

Визначення 1:Интерполяцией називається пошук наближеною функціїF(X), такий щоF(xі)=>f(xі), деi=0,1…n, af(xі) відомі значення функціїF(X) наотрезке[x0, xn]. Крапки, у якихF(xі)=>f(xі) називаються вузлами інтерполяції.

Визначення 2: Якщо знайденаинтерполяционная функціяF(X) для відрізка [x0, xn] має область визначення поза цього відрізка, тоді називатиметься екстраполяцією функціїf(x).

Однією з методів інтерполяції є методИнтерполяции статечниммногочленом

Шукатимемоинтерполяционную функціюF(X) як багаточлена ступеня n:

(*)

>МногочленPn(x) має n+1 коефіцієнт, отже, n+1 умова, накладене на багаточлен однозначно визначить його коефіцієнти, які можна отримані їх умови:

чи

>Разрешив неї щодо aі (>i=0,1…,n), одержимо аналітичне вираз дляполинома (*).


3. Масив вихідних даних

>Опитним шляхом знайдено дані істинноїмольной теплоємності киснюm порівн при постійному тискуP=const, за нормальної температуриt=0 0 З,t=500 0 З, іt=1000 0 З, представлені таблицею 1.

Таблиця 1.

№ варіанта

1 29.2741 33.5488 35.9144
2 29.2801 33.5501 35.9201
3 29.2729 33.5493 35.9167
4 29.30 33.5479 35.9251
5 29.2752 33.5485 35.9109
6 29.2748 33.5397 35.8999
7 29.2752 33.5501 35.9123
8 29.2744 33.5486 35.9128
9 29.2699 33.5484 35.9251
10 29.2742 33.5481 35.9109
11 29.2753 33.5399 35.9201
12 29.2748 33.5501 35.9167
13 29.2801 33.5493 35.9144
14 29.2729 33.5479 35.9201
15 29.2744 33.5485 35.9123
16 29.2699 33.5493 35.9128
17 29.2742 33.5479 35.9251
18 29.2753 33.5485 35.9109
19 29.2748 33.5397 35.9128
20 29.2752 33.5501 35.9251
21 29.2744 33.5486 35.9201
22 29.2741 33.5484 35.9167
23 29.2801 33.5481 35.9144
24 29.2753 33.5486 35.9201

>мольний теплоємність інтерполяція програма

У нашому випадку розглядаються дані варіанта №5.



№ варіанта

5 29.2752 33.5485 35.9109

4. Математичні моделі, застосовувані для розрахунків

>Интерполяционний багаточленm порівн=>f(t0), матиме такий вигляд:

,

з урахуванням нього складається система лінійних рівнянь, дозволивши яку щодо коефіцієнтів a, b,d, одержимоинтерполяционную функцію.Составим тих данихинтерполяционние рівняння:

1. 

2.   

3.   

4.   

 

 

5.

>y=29,2752+0,0104575*t-0,0000038218*t2


5. Результати підрахунків

>t°,C

>m порівн

>t°,C

>m порівн

 0 29.2752  400 32.8467
 10 29.3794  410 32.9203
 20 29.4828  420 32.9932
 30 29.5855  430 33.0653
 40 29.6874  440 33.1366
 50 29.7885  450 33.2072
 60 29.8889  460 33.2770
 70 29.9885  470 33.3460
 80 30.0873  480 33.4143
 90 30.1854  490 33.4818
 100 30.2827  500 33.5485
 110 30.3793  510 33.6145
 120 30.4551  520 33.6797
 130 30.5701  530 33.7441
 140 30.6643  540 33.8078
 150 30.7578  550 33.8707
 160 30.8506  560 33.9329
 170 30.9425  570 33.9943
 180 31.0337  580 34.0549
 190 31.1242  590 33.1148
 200 31.2138  600 34.1739
 210 31.3027  610 34.2322
 220 31.3909  620 34.2897
 230 31.4783  630 34.3466
 240 31.5649  640 34.4026
 250 31.6507  650 34.4579
 260 31.7358  660 34.5124
 270 31.8201  670 34.5661
 280 31.9037  680 34.6191
 290 31.9865  690 34.6713
 300 32.0685  700 34.7228
 310 32.1497  710 34.7735
 320 32.2302  720 34.8234
 330 32.3100  730 34.8725
 340 32.3890  740 34.9209
>t°,C

>m порівн

>t°,C

>m порівн

750 34.9686  1150 36.2470
760 35.0154  1160 36.2633
770 35.0615  1170 36.2788
780 35.1069  1180 36.2936
790 35.1514  1190 36.3076
800 35.1952  1200 36.3208
 810 35.2383  1210 36.3333
 820 35.2806  1220 36.3450
 830 35.3221  1230 36.3559
 840 35.3628  1240 36.3661
 850 35.4028  1250 36.3755
 860 35.4420  1260 36.3842
 870 35.4805  1270 36.3920
 880 35.5185  1280 36.3992
 890 35.5551  1290 36.4055
 900 35.5913  1300 36.4111
 910 35.6267  1310 36.4159
 920 35.6613  1320 36.4200
 930 35.6952  1330 36.4233
 940 35.7283  1340 36.4258
 950 35.7607  1350 36.4276
 960 35.7922  1360 36.4286
 970 35.8230  1370 36.4288
 980 35.8531  1380 36.4283
 990 35.8824  1390 36.4270
 1000 35.9109  1400 36.4250
 1010 35.9387  1410 36.4222
 1020 35.9656  1420 36.4186
 1030 35.9919  1430 36.4142
 1040 36.0173  1440 36.4091
 1050 36.0420  1450 36.4032
 1060 36.0660  1460 36.3966
 1070 36.0891  1470 36.3892
 1080 36.1116  1480 36.3810
 1090 36.1332  1490 36.3721
 1100 36.1541  1500 36.3624
 1110 36.1742
 1120 36.1935
 1130 36.2121
 1140 36.2299

Графік:

6. Текст програми

#>include<stdio.h>

#>include<conio.h>

#>include<math.h>

>floatandrey (>floatc1,floatc2,floatm);

>voidmain()

{>clrscr();

>floatp1,p2,b,d;

>intt1=500,i;

>floatk1=29.2752,k2=33.5485,k3=35.9109;

>p1=(k2-k1)/t1;

>p2=(k3-k1)/(2*t1);

>d=-(p1-p2)/t1;

>b=p1-t1*d;

>printf ("nb=%f",b);

>printf ("nd=%f",d);

>andrey (>b,d,k1);}

>floatandrey (>floatc1,floatc2,floatm)

{>clrscr();

>floatt[1000];

>floaty[1000];

>inth=10,i;

>for (>t[0]=0,i=0;i<=150;i++)

{>t[i]=t[0]+i*h;

>y[i]=m+c1*t[i]+c2*t[i]*t[i];

>printf ("nt[%i]=%7.2fy[%i]=%7.2f",i,t[i],i,y[i]);}

>getch();}


Висновок

Дані істинноїмольной теплоємності киснюm порівн ,знайдені дослідним шляхом при постійному тискуP=const, за нормальної температуриt=0 0 З,t=500 0 З, іt=1000 0 З, збіглися зm порівн, які я знайшла з допомогою мови Сі. Отже, метод інтерполяції спрацював.


Список літератури:

1.Паппас Кріс Мюрей. Програмування мовоюС++:-К.: Видавнича групаBHV, 2000. -320с.

2.Крячков А.В.,Сухинина І.В.,ТомшинВ.К. Програмування на З повагою та З++.Практикум:Учеб. посібник для вузів/Крячков А.В.,Сухинина І.В.,ТомшинВ.К.: Під ред.Томшина – 2-ге вид.испр. – М.: Гаряча лінія – телеком. 2000 – 344 з.: мул.

3.Подбельский В.В., ФомінС.С. Програмування мовою Сі:Учеб. посібник – 2-гедоп. вид. – М.: Фінанси і статистика, 2000 – 600 з.: мул.

4.ГутерР.С.,ОвчинскийБ.В. Елементи чисельного аналізу та математичного опрацювання результатів досвіду. – 2-ге вид., перераб. – М.: Наука, 1970, 432 з.

5. ВолковЕ.А.Численние методи. – 2-ге вид.испр. – М.: Наука, 1987, 248 з.

6.Мудров А.Є.Численние методи для ПЕОМ мовами Бейсік, Фортран і Паскаль – Томськ: ">РАСКО", 1991, - 272 з.: мул.

7.Плис А.І.,Сливина Н.А. Лабораторний практикум з вищої математики.:Учеб.пособ. для втузів. . – 2-ге вид., перераб. ідоп. – М.:Висш. шк., 1994. – 416 з.


Схожі реферати:

Навігація