Function q(t:real):real;
BEGIN
q:=x[1]+x[2]*t+x[3]*t*t+x[4]*t*t*t;
END;
procedure simpson(a1,b1:real;n2:integer;var s:real);
var
hs:real;e,i:integer;
begin
hs:=(b1-a1)/(2*n2);
s:=q(a1)+q(b1);
e:=1;
for i:=1 to 2*n2-1 do
begin
s:=s+(3+e)*q(a1+i*hs);
e:=-e;
end;
s:=s*hs/3;
end;
BEGIN
Clrscr;
write('vedit kilkistto4ok tabli4noi fynksi N=');readln(n1);
{ctepin polimena} m1:=3;
writeln('vedit vixidni to4ki');
for i3:=1 to n1 do begin
read(x1[i3]); read(y1[i3]);
WRITELN;
end;
writeln('_________________________');
writeln('vihidni dani');
for i3:=1 to n1 do
begin
write(' ',x1[i3]:5:2); write(' ',y1[i3]:5:2);
WRITELN;
end;
c[1]:=0; b[1]:=0;
MNK (n1,x1,y1,a);
for i2:=1 to m1+1 do n:=n+1;
Gaus(n,a,b,x);
writeln('*** vektor rozvyazky ***');
writeln;
for i1:=1 to n do writeln('x[',i1,']=',x[i1]:7:3);
{kofisienti zapicyutsya a0+a1*x+a2*x^2+a3*x^3+...=const}
write(' q(t)'); write(' y');
writeln;
for i3:=1 to n1 do
begin
e1:=e1+Sqr(q(x1[i3])-y1[i3]);
write(' ',q(x1[i3]):7:3); write(' ',y1[i3]:7:3);
writeln;
end;
writeln('poxibka aproksimasyi',e1:7:5);
write('введіть b2=');read(b2);
a1:=0; b1:=b2/10;
eps:=0.001;
repeat
k:=2;
simpson(a1,b1,k,s);
repeat
s0:=s; k:=k*2;
simpson(a1,b1,k,s);
writeln('b1=',b1:5:2,' s=',s:5:2);
until (abs(s-s0)<=eps)or(k>500);
b1:=b1+b2/10;
until b1>b2;
END.
3.7 Відлагодження програми
Виправлення синтаксичних помилок на підставі повідомлень транслятора
Після написання програми транслятор знайшов ні помилку у 149-му рядку з повідомленням «Error 26: type mismatsh» I, j, k, n2: integer; Після заміни real на integer, програма була запущена транслятором на виконання.
Таблиця 3.3 – Результати контролю правильності розв’язання задачі
Вихідні дані | Результати розрахунків | ||||||
ППП MathCad 2000 | Програма на мові Turbo Pascal | ||||||
400 800 1200 1600 2000 2400 2800 3200 3600 4000 4400 4800 5200 5600 | 35,3 43,5 52,6 58,1 58,5 61,6 67,9 72,3 67,5 64,4 65,3 57,9 55,6 54,3 | 363 352 351 342 338 335 329 324 334 340 341 356 363 370 | 34.47-0.123x+0.028x2-0.0012x3 | 421.53-2.37x+0.452x2+0.072x3 | 31.71-0.12x+0.027x2-0.001x3 | 417.32-2.41x+0.457x2+0.07x3 |
В результаті виконаної роботи були розраховані основні параметри підшипника коробки передач автомобіля, за результатами розрахунків ми дізнались параметри підшипника. Основна мета було знайти динамічну завантаженість ми отримали 3.17 H
Склавши програму по розрахунку системи лінійних рівнянь методом Гауса розрахували коефіцієнти системи. Знайдені коефіціенти
х
=-3.9, х =0.63, х =-4.63, х =9.98, х =1.094,Таким чином у нас є 5 чисел які заміняють невідомі.
Розрахували швидкісну характеристику автомобіля і склали таблицю яку можна роздрукувати. З таблиці видно що з більшання обертів і зростає навантаження.
ПЕРЕЛІК ПОСИЛАНЬ
1. Смирнов Г.А. Теория движения колесных машин. – М.: Машиностроение, 1981.- 271с.
2. Литвинов А.С., Фаробин Я.Е. Автомобиль. Теория експлуатационных свойств.- Москва: Машиностроение, 1989.- 237с.
3. Дьконов В. MATHCAD 8/2000: специальный справочник – СПб: Издательство «Питер», 2000.- 592с.
4. Кошарний М.Ф. Основи механіки та енергетики автомобіля. – Київ: Вища школа, 1992. – 200с.
5. Turbo Pascal 7.0 for users. BHV, 1999. – 448p.
6. Т.А.Павловськая Pascal Издательство «Питер», 2003.-291с.
7. П.П.Овчинников, Ф.П.Яремчук, В.М.Михайленко -Вища математика Частина-1 : Київ «Техніка»2003 -597с.
8. Савуляк В.І. Навчальний посібник
9. Кишеніна Н.В., Кишені В.О., «основи інформаційний технологій та програмування « частина-1 Вінниця 2003.