Метод Адамса

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

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

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

Программа №1

Matlab M

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

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)?

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

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

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