Теми рефератів
> Реферати > Курсові роботи > Звіти з практики > Курсові проекти > Питання та відповіді > Ессе > Доклади > Учбові матеріали > Контрольні роботи > Методички > Лекції > Твори > Підручники > Статті Контакти
Реферати, твори, дипломи, практика » Курсовые проекты » Чисельне рішення звичайних диференціальних рівнянь

Реферат Чисельне рішення звичайних диференціальних рівнянь





й довідник.- СПб .: Питер, 2002. - 528 с.



Додаток 1


Протокол програми (в М-файлі)

format long

% метод трапеції=1=1=1=3=100=abs (xn-x0)/n=int (x0, a, b); i=1: n - 1

x (i)=x0 + i * h;

s=s + 2 * int (x (i), a, b); ( метод трапецій )=(s + int (xn, a, b)) * h/2

disp ( аналітичні рішення )=intt (xn, a, b) -intt (x0, a, b)

% метод Симсона

j=0;=int (x0, a, b); j lt; n - 1=j + 1;=s + 4 * int (x0 + j * h, a, b); j == n - 1

continue=j + 1;=s + 2 * int (x0 + j * h, a, b); ( метод Симсона )=h * (s + int (xn, a , b))/3

t=[0.06943184 0.33000984 0.66999052 0.93056816];=[0.17392742 0.32607258 0.32607258 0.17392742];=0; i=1: 4

s=s + w (i) * int (x0 + (xn-x0) * t (i), a, b);

end ( квадратурні формули Гаусса )

s=(xn-x0) * s

tch=[0.083751 0.312730 0.5 0.68727 0.916249];=0; i=1: 5

s=s + 0.2 * int (x0 + (xn-x0) * tch (i), a, b); ( метод Чебишева )=(xn-x0) * s

function y=int (x, a, b)=1/((a * x + b) ^ 2); y=intt (x, a, b)=- 1/(a * (a * x + b));



Додаток 2


Протокол програми (в М-файлі) long=1

% кількість кроків=1000;

% обнуленіемассів=zeros (1, n);

y=zeros (1, n);

% н.у

x (1)=0; (1)=0.1;

% крок рішення

h=(xend-x (1))/n; i=2: n

k1=odddde (x (i - 1), y (i - 1));

k2=odddde (x (i - 1) + h/2, y (i - 1) + h * k1/2);

k3=odddde (x (i - 1) + h/2, y (i - 1) + h * k2/2);

k4=odddde (x (i - 1) + h, y (i - 1) + h * k3);

y (i)=y (i - 1) + h * (k1 + 2 * k2 + 2 * k3 + k4)/6;

x (i)=x (i - 1) + h;

% plot (x, y)=0.1;

[xm, ym]=ode45 (odddde, [0: 0.001: 1], y0); (x, y, xm, ym)



Додаток 3


Протокол програми (в М-файлі)

% Введення рядками=[2 - 4 3 1;

- +1 5 - 7 - 3;

10 - 2 4 4;

- 1 1 - 1 - 1];=[7;- 24; 34;- 6];=a b

% Введення стовпцями=[2 - 1 10 - 1;

- 4 +5 - 2 січня;

3 - 7 +4 - 1;

1 - 3 +4 - 1];=[7 - 24 34 - 6];=b/a=b * a ^ - 1=b * inv (a)



Додаток 4


Протокол програми (в М-файлі)

% початкове значення кореня=1;

% лічильник інтерація

n=0; abs (ur (x)/dur (x)) gt; 1e - 7

x=x-ur (x)/dur (x);

n=n + 1;

if n gt;=1000

disp ( рішення немає )

break

end(x)(n)y=ur(x)=x^4-5*x^3+8*x^2-5*x+1;y=dur(x)long=1;=7;(1)=x-h;(1)=ur(x(1));i=1:n;(i+1)=x(1)+i*h;(i+1)=ur(x(i+1));(x)(y)i=1:n- 1 (i)=y (i + 1) -y (i); i=1: n - 2y (i)=dy (i + 1) -dy (i); i=1: n - 3y (i )=d2y (i + 1) -d2y (i); i=1: n - 4y (i)=d3y (i + 1) -d3y (i);=(dy (1) + 0.5 * d2y (1)-d3y (1)/6 + d4y (1)/12)/h;


Назад | сторінка 9 з 9





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

  • Реферат на тему: Метод Ньютона (метод дотичних). Рішення систем нелінійних алгебраїчних рів ...
  • Реферат на тему: Рішення систем нелінійніх рівнянь. Метод ітерацій. Метод Ньютона-Канторов ...
  • Реферат на тему: Створення програми "Підрахунок частоти введених слів у текстовому файл ...
  • Реферат на тему: Прямі методи рішення лінійних систем. Метод квадратного кореня
  • Реферат на тему: Додаток, що реалізує метод заміни найменш значущих бітів для файлів з розши ...