Метод Адамса

@Стелла21 · Регистрация: 07.06.2014

Студворк — интернет-сервис помощи студентам

Привет,помогите исправить ошибку в коде,я знаю точно что они в строках с out,но не могу разобраться вот код написан в матлабе

Программа №1

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
function [out] = Adams()
out = [];
xi0=0;
yi0=0.1; 
h=0.02;
N=10;
out=[out;xi0,yi0];
k1 = h*F(xi0, yi0);
k2 = h*F(xi0+h/2.0, yi0+k1/2.0);
k3 = h*F(xi0+h,yi0+2.0*k2-k1);
yi1=yi0+(k1+4.0*k2+k3)/6.0;
xi1= xi0+h;
 
k11 = h*F(xi1, yi1);
k22 = h*F(xi1+h/2.0, yi1+k11/2.0);
k33 = h*F(xi1+h,yi1+2*k22-k11);
yi2=yi1+(k11+4.0*k22+k33)/6.0;
xi2= xi1+h;
 
for xi0=3:N
out=[out,xi0,yi2];
xi3= xi2+h;
yi3= yi2 + h/12*(23*F(xi2,yi2)- 16*F(xi1,yi1)+5*F(xi0,yi0));
yi0=yi1;
yi1=yi2;
yi2=yi3;
xi0=xi1;
xi1=xi2;
xi2=xi3;
 
end
plot(out(:,1),out(:,2),'.--g')
hold on
l=dsolve('Dy = 30*y*(x-0.2)*(x-0.7)','y(0)=0.1','x');
ezplot(l,[0 1])
end
function [out] = F(x,y)
out =(30*y*(x-0.2)*(x-0.7));
end

@Dimas_rdn · 08.06.2014, 10:57

Не знаю, сработало без ошибок, ниже результат

@Стелла21 · 08.06.2014, 14:38 **[ТС]**

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

Программа №2

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
function [out] = Untitled()
function [out] = F(x,y)
out =(30*y*(x-0.2)*(x-0.7));
end
out = [];
 
xi0=0;
yi0=0.1; 
h=0.02;
N=10;
out=[out;xi0,yi0];
k1 = h*F(xi0, yi0);
k2 = h*F(xi0+h/2.0, yi0+k1/2.0);
k3 = h*F(xi0+h,yi0+2.0*k2-k1);
yi1=yi0+(k1+4.0*k2+k3)/6.0;
xi1= xi0+h;
 
k11 = h*F(xi1, yi1);
k22 = h*F(xi1+h/2.0, yi1+k11/2.0);
k33 = h*F(xi1+h,yi1+2*k22-k11);
yi2=yi1+(k11+4.0*k22+k33)/6.0;
xi2= xi1+h;
 
for i=3:N
 
xi3= xi2+h;
yi3= yi2 + h/12*(23*F(xi2,yi2)- 16*F(xi1,yi1)+5*F(xi0,yi0));
yi0=yi1;
yi1=yi2;
yi2=yi3;
xi0=xi1;
xi1=xi2;
xi2=xi3;
out=[xi0,yi0;xi2,yi1;xi2,yi2;];
 
end
plot(out(:,1),out(:,2),'.--g')
hold on
l=dsolve('Dy = 30*y*(x-0.2)*(x-0.7)','y(0)=0.1','x');
ezplot(l,[0 1])
end

спасибо зарание

Добавлено через 7 минут
препод сказал что: С ошибкой пока не разобрались, поэтому график приближенного решения неверный. 2 кода и оба не верны,в чем может быть ошибка в строках с out?

@S_el · 08.06.2014, 14:47

Стелла21, а препод не сказал какая конкретно ошибка?

@Стелла21 · 08.06.2014, 15:04 **[ТС]**

В коде метода Адамса есть ошибки. Проверьте, пожалуйста, строки с out-вот все что он сказал,это еще с 1 кода

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

Добавлено через 8 минут
и еще может быть в этом?
plot(out(:,1),out(:,2),'.--g')
я это где то в инете прочитала но так и не поняла толком,знаю что рисует,может не правильно out(:,1),out(:,2)?

Новые блоги и статьи Все статьи Все блоги /
Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .

@Dimas_rdn 41 / 41 / 9 Регистрация: 22.10.2012 Сообщений: 91
	08.06.2014, 10:57
	Не знаю, сработало без ошибок, ниже результат Миниатюры 1

@S_el 2444 / 1842 / 406 Регистрация: 15.12.2013 Сообщений: 8,243
	08.06.2014, 14:47
	Стелла21, а препод не сказал какая конкретно ошибка? 0

@Стелла21 1 / 1 / 0 Регистрация: 07.06.2014 Сообщений: 23
	08.06.2014, 15:04 [ТС]
	В коде метода Адамса есть ошибки. Проверьте, пожалуйста, строки с out-вот все что он сказал,это еще с 1 кода Добавлено через 1 минуту как я его поняла отправив ему 2 код я от какой то ошибки избавилась,осталась 1 ошибка из за который график ине верен я уже не знаю что делать помогите пожалуйста уже месяца 2 пытаюсь сдать Добавлено через 8 минут и еще может быть в этом? plot(out(:,1),out(:,2),'.--g') я это где то в инете прочитала но так и не поняла толком,знаю что рисует,может не правильно out(:,1),out(:,2)? 0