Уравнение 3-го порядка явным методом эйлера

@kolkor · Регистрация: 18.10.2016

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

Добрый день, уважаемые форумчане, помогите написать программу для решения уравнения 3-го порядка явным методом эйлера.
6y'''(t) + 35y''(t) + 1,7y'(t) + 0,5y(t) =5
y''(0) = y'(0) = y(0) =0
t = [0;140]
Со вторым порядком, как мне казалось, разобралась, но стоило взглянуть на третий, как тут же впала в ступор.

Добавлено через 22 часа 17 минут

Matlab M

clear, clc
h = 0.5;
T = 0.01:h:140;
f = @(x,y,dy,ddy) 1/12*(10-y-70*ddy-3.4*dy);
 
Y = zeros(size(T)); 
dY = zeros(size(T));
ddY = zeros(size(T));
 
Y(1) = 0; 
dY(1) = 0;
ddY(1) = 0;
 
for j=1:length(T)-1
    ddY(j+1)= ddY(j) + h*f(T(j),Y(j),dY(j),ddY(j)); % считаем производную
    dY(j+1)=dY(j)+ h*ddY(j);
    Y(j+1)=Y(j)+ h*dY(j); % считаем ф-цию
end 
plot(T,Y,T,dY)
legend('y(x)','dy/dx')

помогите, пожалуйста, найти ошибку х___х
он строит совсем не то, что получается при решении с помощью ode45, вот даже не близко.

@Krasme · 25.11.2016, 22:52

а где решение с ode45?

Зосима · 25.11.2016, 23:39

kolkor, а если шаг уменьшить?

Не по теме:

*выражение f верное, но ИМХО "(5-35*ddy-1.7*dy-0.5*y)/6" было бы понятней

Зосима · 25.11.2016, 23:54

Так и есть! Уже при h=0.3 имеем сходимость численного решения

Новые блоги и статьи
Как написать микросервис на Go/Golang InfoMaster 14.01.2025 Определение микросервиса, преимущества использования Go/ Golang Микросервис – это архитектурный подход к разработке программного обеспечения, при котором приложение состоит из небольших, независимо. . .	Как написать микросервис с нуля на C# InfoMaster 14.01.2025 В современном мире разработки программного обеспечения микросервисная архитектура стала стандартом де-факто для создания масштабируемых и гибких приложений. Этот архитектурный подход предполагает. . .	Как создать интернет-магазин на PHP и JavaScript InfoMaster 14.01.2025 В современном мире электронная коммерция стала неотъемлемой частью бизнеса. Создание собственного интернет-магазина открывает широкие возможности для предпринимателей, позволяя достичь большей. . .	Как написать Тетрис на Ассемблере InfoMaster 14.01.2025 Тетрис – одна из самых узнаваемых и популярных компьютерных игр, созданная в 1984 году советским программистом Алексеем Пажитновым. За прошедшие десятилетия она завоевала симпатии миллионы людей по. . .	Как создать игру "Танчики" на Unity3d и C# InfoMaster 14.01.2025 Разработка игр – это увлекательный процесс, сочетающий в себе творчество и технические навыки. В этой статье мы рассмотрим создание классической игры "Танчики" с использованием Unity3D и языка. . .	Организую платный онлайн микро-курс по доработке Android-клиента Telegram _Ivana 14.01.2025 Официальная версия и распространенные форки не полностью устраивают? Сделай свою кастомную версию клиента! 4 занятия по 2 часа (2 недели пн, ср 19:00-21:00 по Москве). Первое вводное занятие. . .
Как создать приложение для фитнеса для iOS/iPhone на Kotlin InfoMaster 14.01.2025 Создание собственного фитнес-приложения — это не только захватывающий, но и полезный процесс, ведь оно может стать вашим верным помощником на пути к здоровому и активному образу жизни. В современных. . .	Как создать приложение магазина для iOS/iPhone на Swift InfoMaster 14.01.2025 Введение в разработку iOS-приложений Разработка приложений для iPhone и других устройств на базе iOS открывает огромные возможности для создания инновационных мобильных решений. В данной статье мы. . .	Это работает. Скорость асинхронной логики велика. Вопрос видимо останется в стабильности. Плата - огонь! Hrethgir 13.01.2025 По прошлому проекту в Logisim Evolution https:/ / www. cyberforum. ru/ blogs/ 223907/ blog8781. html прилагаю файл архива проекта в Gowin Eda. Восьмибитный счётчик из сумматора+ генератор сигнала. . .	UserScript для подсветки кнопок языков программирования в зависимости от текущего раздела volvo 13.01.2025 В результате работы этого скрипта подсвечиваются нужные кнопки не только в форме быстрого ответа, но и при редактировании сообщения: / / ==UserScript== / / @name CF_DefaultLangSelect / / . . .	Введение в модели и алгоритмы машинного обучения InfoMaster 12.01.2025 Машинное обучение представляет собой одну из наиболее динамично развивающихся областей искусственного интеллекта, которая фокусируется на разработке алгоритмов и методов, позволяющих компьютерам. . .	Как на Python создать нейросеть для решения задач InfoMaster 12.01.2025 В контексте стремительного развития современных технологий особое внимание уделяется таким инструментам, как нейросети. Эти структуры, вдохновленные биологическими нейронными сетями, используются для. . .

@Krasme 6962 / 4991 / 2108 Регистрация: 02.02.2014 Сообщений: 13,249
	25.11.2016, 22:52	2
	а где решение с ode45? 0

Зосима 5243 / 3571 / 379 Регистрация: 02.04.2012 Сообщений: 6,473 Записей в блоге: 17
	25.11.2016, 23:39	3
	kolkor, а если шаг уменьшить? Не по теме: выражение f верное, но ИМХО "(5-35ddy-1.7dy-0.5y)/6" было бы понятней 0

Зосима 5243 / 3571 / 379 Регистрация: 02.04.2012 Сообщений: 6,473 Записей в блоге: 17
	25.11.2016, 23:54	4
	Так и есть! Уже при h=0.3 имеем сходимость численного решения 0