Форум программистов, компьютерный форум, киберфорум
Наши страницы
Matlab
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.61/18: Рейтинг темы: голосов - 18, средняя оценка - 4.61
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 292
1

Методы Эйлера и Рунге-Кутта

04.06.2013, 17:50. Просмотров 3266. Ответов 8
Метки faq+ (Все метки)

Подскажите, пожалуйста, где ошибка. Выдает следующие сообщения после запуска на выполнение:
Matlab M
1
2
3
4
5
6
7
8
9
10
11
Error using feval
Undefined function 'f4' for input arguments of type 'double'.
 
Error in Euler (line 30)
    Y1(j+1)=Y1(j)+h*(feval(f,T(j),Y(j))-b1*Y1(j)-c1*Y(j)); % решение первой системы
 
Error in reshenie (line 36)
        [Y]=Euler(f,a,b,a1,b1,c1,ya1,ya2,h,p) % вызов функции Euler
 
Error in resh_DU (line 15)
[Y]=reshenie('f4',a1,b1,a,b,c,ya1,ya2,h); % вызов функции reshenie
Программа должна была решать ДУ методами Эйлера и Рунге-Кутта 4-го порядка (выбирает пользователь). Помогите, чем можете. Очень надо до завтра.
0
Вложения
Тип файла: rar Euler_RK4.rar (3.6 Кб, 99 просмотров)
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
04.06.2013, 17:50
Ответы с готовыми решениями:

Решить задачу Коши двумя методами: Эйлера и Рунге-Кутта
решить задачу коши "Численное решение задачи Коши для дифференциальных...

Решение уравнения численно, методом Рунге-Кутта, или Эйлера
Подскажите пожалуйста, как решить задачу. Как я понял, нужно всего лишь...

Решение дифференциального уравнения методом Эйлера и методом Рунге-кутта 4 порядка
Помогите пожалуйста решить уравнение y''-4y'+5y=2x2ex , методом Эйлера и...

Решение дифференциального уравнения методом Эйлера и методом Рунге-кутта
Помогите пожалуйста решить уравнение у' = 1 + х sin y, y(π) = 2π , методом...

Метод Рунге-Кутта
Здравствуйте, помогите пожалуйста решить проблему. Решаю простой дифур...

8
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 292
05.06.2013, 07:36  [ТС] 2
Чем могут быть вызваны подобные ошибки?
0
Зосима
4929 / 3301 / 312
Регистрация: 02.04.2012
Сообщений: 6,207
Записей в блоге: 15
Завершенные тесты: 1
05.06.2013, 09:12 3
У тебя не определена функция f4, ее нет в архиве!
Кроме того (я написал простенькую f4 и запустил) не хватает функции toch resh, однако создать такую, с пробелом в имени невозможно!
Но кое-что можно придумать, если будем знать f4
1
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 292
05.06.2013, 12:36  [ТС] 4
Сделала вот такую функцию. Но все равно не работает(((
Matlab M
1
2
3
function f=f4(x)
f=(exp(3*x))*(13*x-7);
end
Добавлено через 4 минуты
Помогите, пожалуйста, чтобы прога хоть как-то работала. Лишь бы ошибок не было. У меня зачет через час(((((
0
Зосима
4929 / 3301 / 312
Регистрация: 02.04.2012
Сообщений: 6,207
Записей в блоге: 15
Завершенные тесты: 1
05.06.2013, 12:41 5
ща подшаманю, выложу
1
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 292
05.06.2013, 12:45  [ТС] 6
Спасибо большущее, жду
0
Зосима
4929 / 3301 / 312
Регистрация: 02.04.2012
Сообщений: 6,207
Записей в блоге: 15
Завершенные тесты: 1
05.06.2013, 13:01 7
Так-с первым делом поясню: функция toch_resh - это точное решение данного ДУ.
Я его считал отдельно так:
Matlab M
1
2
R = dsolve('D2y+4*y=exp(3*t)*(13*t-7)','y(0)=1','Dy(0)=0')
ft = simplify(R) % упрощаем
Получилось:
Matlab M
1
2
3
4
5
R =
 2*cos(2*t) + sin(2*t) - (cos(2*t)*exp(3*t)*(2*cos(2*t) - 2*sin(2*t) - 2*t*cos(2*t) + 3*t*sin(2*t)))/2 - (sin(2*t)*exp(3*t)*(2*cos(2*t) + 2*sin(2*t) - 3*t*cos(2*t) - 2*t*sin(2*t)))/2
 
ft =
2*cos(2*t) - exp(3*t) + sin(2*t) + t*exp(3*t)
Последнее выражение я и ввел в toch_resh:
Matlab M
1
2
3
function y = toch_resh(t)
y = 2*cos(2*t) - exp(3*t) + sin(2*t) + t.*exp(3*t);
end
Вообщем получаются такие рисуночки

Методы Эйлера и Рунге-Кутта
Методы Эйлера и Рунге-Кутта


Теперь держи архивчик с файлами
Euler_RK4.rar
2
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 292
05.06.2013, 13:05  [ТС] 8
Супер! Спасибо огромное! Всё поняла. Вы мой спаситель!!!
1
Зосима
4929 / 3301 / 312
Регистрация: 02.04.2012
Сообщений: 6,207
Записей в блоге: 15
Завершенные тесты: 1
05.06.2013, 13:45 9
И тебе спасибо за готовые программки решения
0
05.06.2013, 13:45
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.06.2013, 13:45

Метод Рунге-Кутта
Необходимо решить дифф. уравнение методом Рунге-Кутта 4-го порядка....

Метод Рунге-Кутта 4 порядка
Помогите исправить ошибки % моя функция для вычисления значений function res...

Диф.ур. метод Рунге-Кутта
Здравствуйте. Необходимо реализовать метод Рунге-Кутта 4 порядка и решить...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
9
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru