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

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

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

Подскажите, пожалуйста, где ошибка. Выдает следующие сообщения после запуска на выполнение:
http://www.cyberforum.ru/matlab/thread888582.html
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
Я подобрал для вас темы с готовыми решениями и ответами на вопрос Методы Эйлера и Рунге-Кутта (Matlab):

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

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

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

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

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

8
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 290
05.06.2013, 07:36  [ТС] #2
Чем могут быть вызваны подобные ошибки?
0
Зосима
4915 / 3284 / 308
Регистрация: 02.04.2012
Сообщений: 6,191
Записей в блоге: 15
Завершенные тесты: 1
05.06.2013, 09:12 #3
У тебя не определена функция f4, ее нет в архиве!
Кроме того (я написал простенькую f4 и запустил) не хватает функции toch resh, однако создать такую, с пробелом в имени невозможно!
Но кое-что можно придумать, если будем знать f4
1
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 290
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
Зосима
4915 / 3284 / 308
Регистрация: 02.04.2012
Сообщений: 6,191
Записей в блоге: 15
Завершенные тесты: 1
05.06.2013, 12:41 #5
ща подшаманю, выложу
1
yuppi
11 / 3 / 1
Регистрация: 22.01.2012
Сообщений: 290
05.06.2013, 12:45  [ТС] #6
Спасибо большущее, жду
0
Зосима
4915 / 3284 / 308
Регистрация: 02.04.2012
Сообщений: 6,191
Записей в блоге: 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
Сообщений: 290
05.06.2013, 13:05  [ТС] #8
Супер! Спасибо огромное! Всё поняла. Вы мой спаситель!!!
1
Зосима
4915 / 3284 / 308
Регистрация: 02.04.2012
Сообщений: 6,191
Записей в блоге: 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 порядка
Помогите исправить ошибки % моя функция для вычисления значений function res...

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

Дифференциальное уравнение методом Рунге-Кутта.
нужно решить дфиеренциальное уравнение в matlab методом рунге-кутта. вот оно:...

Уравнение, метод Рунге-Кутта 4 порядка
Есть решение 18 задания: h=1; x=0:h:25; xi=0; yi=0; N=25/h; ...


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

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

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