Zaloguj się
Blog na Matlab.pl
Forum polskich użytkowników
 
UŻYTKOWNICY GRUPY PROFIL Zaloguj się, by sprawdzić wiadomości FAQ
 



Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

SciLab Rozwiązanie równania różniczkowego 2 stopnia
Forum MATLAB Strona Główna-> Klony Matlaba (Scilab, Octave, ...)
Post Wysłany: 27 Styczenia 2016, Sro 11:40 am Temat postu: SciLab Rozwiązanie równania różniczkowego 2 stopnia Odpowiedz z cytatem
 
AUTOR:
dzikidarek



Dołączył: 27 Sty 2016
Posty: 1


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam. Mam do rozwiązania równanie:
y"+3y'+2y=3e^(-2t)
Muszę je rozwiązać w SciLabie w XCosie oraz w konsoli.
Przekształciłem to równanie do postaci:
y=x(1)
y'=x'(1)=x(2)
y"=x'(2)
czyli
x'(2)+3*x(2)+2*x(1)=3*e^(-2*t)

x'(2)=-3*x(2)-2*x(1)+3*e^(-2*t)

W konsoli spróbowałem to rozwiązać tak:
-->function [xprim]=pochodna(t,x)
-->global c;
-->global d;
-->xprim(1)=x(2);
-->xprim(2)=-c*x(2)-d*x(1)+3*exp(-2*t);
-->endfunction
-->global c
-->c=3;
-->global d
-->d=2;
-->x0=[-5;5]; t=0:.01:50;
-->[x]=ode(x0,0,t,pochodna);
-->plot(t,x(1,:),'r')

Czy tak powinno to wyglądać? Co ewentualnie robię źle?
Niestety w Xcos'ie nie wiem jak i gdzie wprowadzić e^(-2t)...
W Xcos mam takie coś jak na załączonym rysyunku.
Proszę o pomoc i z góry dzięki.
Rysunek xcos:
https://goo.gl/photos/zstbPJmkF31TB6g5A


 

Post Wysłany: 12 Lutego 2016, Pią 6:07 pm Temat postu: Prośba o pomoc - zadanie Odpowiedz z cytatem
 
AUTOR:
Mariuszek



Dołączył: 12 Lut 2016
Posty: 4


Ogląda profil użytkownika Wyślij prywatną wiadomość
Serdecznie witam,

mam spory problem z rozwiązaniem zadania następującej treści:

Rozwiąż co najmniej dwoma metodami (Eulera i Rungego-Kutty 2 i 4 rzędu) równanie różniczkowe:
(d^2 x(t))/(dt^2 )-8 dx(t)/dt+16x(t)=0
x(0)=1,(dx(0))/dt=1
w przedziale t∈[0,7].
Wyznacz i porównaj błędy obliczeń.


Zrobiłem 2 skrypty: do Eulera i RK II rzędu - wklejam:

function x1=f1(t, x1,x2)
x1 = x2;
endfunction
function x2=f2(t, x1,x2)
x2 = 8*x2-16*x1; !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
endfunction
function [t, x1, x2] = Euler2(t0, tk, h, x0, x00)
N = (tk - t0)/h;
t(0) = t0;
x1(1) = x0;
x2(1) = x00;
for n = 1:N
t(n+1) = t(n) + h;
x2(n+1) = x2(n) + h * f2(t(n), x1(n), x2(n));
x1(n+1) = x1(n) + h * f1(t(n), x1(n), x2(n));

end
endfunction
t0=0; !!!!!!!!!!!!!!
tk=7;
h=0.01;
x0=0;
x00=0;
[t, x1, x2] = Euler2(t0, tk, h, x0, x00);
plot2d(t,x1,-1)
plot2d(t,x2,3)

bl=f-x1;
plot2d(f,bl,-1)


RK II rzędu:

function x1=f1(t, x1,x2)
x1 = x2;
endfunction
function x2=f2(t, x1,x2)
x2 = 8*x2-16*x1;
endfunction
function [t, x1, x2] = RK2(t0, tk, h, x0, x00)
N = (tk - t0)/h;
t(0) = t0;
x1(1) = x0;
x2(1) = x00;
for n = 1:N
t(n+1) = t(n) + h;
k12 = f2(t(n),y1(n), y2(n));
k22 = f2(t(n) + h/2,y1(n), y2(n)+ h/2 * k1);
k32 = f2(t(n) + h/2,y1(n), y2(n)+ h/2 * k2);
k42 = f2(t(n) + h ,y1(N), y2(n)+ h * k3);
x2(n+1) = x2(n) + h/6 * (k12 + 2*k22 + 2*k32 + k42);
k1 = f1(t(n), y1(n),y2(n));
k2 = f1(t(n) + h/2, y1(n)+ h/2 * k1,y2(n));
k3 = f1(t(n) + h/2, y1(n)+ h/2 * k2,y2(n));
k4 = f1(t(n) + h , y1(n)+ h * k3,y2(n));
x1(n+1) = x1(n) + h/6 * (k1 + 2*k2 + 2*k3 + k4);
end
endfunction
t0=0;
tk=7;
h=0.01;
x0=0;
x00=0;
[t, x1, x2] = RK2(t0, tk, h, x0, x00);
plot2d(t,x2,-1)
plot2d(t,x1,3)

bl=f-x1;
plot2d(f,bl,-1)

Jeszcze raz, bardzo proszę o pomoc!


 

Forum MATLAB Strona Główna-> Klony Matlaba (Scilab, Octave, ...)
Wyświetl posty z ostatnich:   

Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Wszystkie czasy w strefie CET (Europa)

Skocz do:  

Statystyki forum:



Od dnia 08.06.2006 forum odwiedzano 37634451
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
Faktoring  (22 Sierpnia 2017, Wto 1:25 pm)
Nauka profesjonalnego montażu filmów  (22 Sierpnia 2017, Wto 1:24 pm)
Simulink/ xPC target - błąd polecenia slrtexplr  (21 Sierpnia 2017, Pon 7:41 pm)
Niedoparki - 2017 - Cały Film Online/Pobierz  (19 Sierpnia 2017, Sob 6:24 pm)
Fatima. Ostatnia tajemnica- 2017 - Cały Film Online/Pobierz  (19 Sierpnia 2017, Sob 6:24 pm)
Mamuśki mają wychodne- 2017 - Cały Film Online/Pobierz [CDA]  (19 Sierpnia 2017, Sob 6:23 pm)
Gang Wiewióra 2 - 2017 - Cały Film Online/Pobierz [CDA]  (19 Sierpnia 2017, Sob 6:22 pm)
Dom wygranych - 2017 - Cały Film Online/Pobierz [CDA]  (19 Sierpnia 2017, Sob 6:22 pm)
Agentka specjalnej troski - 2017 - Film Online/Pobierz  (19 Sierpnia 2017, Sob 6:21 pm)
Królewicz Olch - 2017 - Cały Film Online/Pobierz [CDA]  (19 Sierpnia 2017, Sob 6:20 pm)
Twoje prawa:
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Nie możesz ściągać plików na tym forum