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 - rownania rozniczkowe Idź do strony Poprzedni  1 2
Forum MATLAB Strona Główna-> Klony Matlaba (Scilab, Octave, ...)
Post Wysłany: 5 Czerwca 2009, Pią 12:02 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
r0g0
Może pisać książki


Dołączył: 06 Lis 2007
Posty: 138
Skąd: Białystok


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email
A x to jest jakieś wymuszenie? wejscie?


 

Post Wysłany: 5 Czerwca 2009, Pią 6:57 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
WeW



Dołączył: 05 Cze 2009
Posty: 1


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




Ostatnio zmieniony przez WeW dnia 13 Czerwca 2009, Sob 8:18 pm, w całości zmieniany 1 raz
 

Post Wysłany: 13 Czerwca 2009, Sob 2:05 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
tyGrys



Dołączył: 12 Cze 2009
Posty: 2


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email
witam,
mam dane równanie:
dX=1/2a^2Xdt-asqrt(1-X^2)dW
czy mógłby ktoś pomóc mi przetłumaczyć to równanie tak aby można było zastosować istrukcje 'ode'?
z góry dziękuje,
pozdrawiam


 

Post Wysłany: 17 Styczenia 2010, Nie 12:17 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
misz



Dołączył: 17 Sty 2010
Posty: 2


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam,
mam do rozwiazanie takie samo zadanie jak R2D2 - ale nie wiem co dalej zrobic z ta funkcja ode, tak aby uzyskac te wykresy..

juz kombinuje od kilku godzin, ale bez efektu..;/


 

Post Wysłany: 4 Lutego 2010, Czw 3:54 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
LuckyF



Dołączył: 04 Lut 2010
Posty: 1


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam.

Mam do przeanalizowania zachowanie oscylatora Duffinga. Zadanie prezentuje się tak:


Mam problem z zapisem mojej funkcji do procedury ODE. Sam stworzyłem coś takiego:
Kod:
function xdot=f(t,x)
 
xdot = [-1*x^2 -k; 1 0]*x + A*cos(omega*t)
endfunction

Jednak przy wykonaniu dostaję błąd:
Kod:
x=ode([0;0],0,t,f)
 !--
error 5
Inconsistent column
/row dimensions.

at line       2 of function f called by
x=ode([0;0],0,t,f)


Wiecie może jak odpowiednio zapisać tę funkcję, bo sam się troche pogubiłem przez ten x^3. I w którym wierszu dostanę x(t), a w którym v(t)?

Z góry dzięki!


 

Post Wysłany: 6 Marca 2010, Sob 6:51 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
cinmar



Dołączył: 06 Mar 2010
Posty: 2


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam wszystkich,
Mam równania opisujące zjawisko uderzenia hydraulicznego. Jest to układ dwóch rownań różniczkowych cząstkowych:
dU/dt + U*(dU/dx) + g*(dH/dx) + a*|U|*U = 0
dH/dt + U*(dH/dx) + b*(dU/dx) = 0

gdzie x - połączenie, t-czas, U - prędkość, H - ciśnienie, a,b,g - to są stałe.

Pozdrawiam


 

Post Wysłany: 7 Marca 2010, Nie 9:34 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
r0g0
Może pisać książki


Dołączył: 06 Lis 2007
Posty: 138
Skąd: Białystok


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email
No ładne równanie, tylko teraz jeszcze żebyś przedstawił czego dotyczy problem...


 

Post Wysłany: 7 Marca 2010, Nie 11:31 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
cinmar



Dołączył: 06 Mar 2010
Posty: 2


Ogląda profil użytkownika Wyślij prywatną wiadomość
Fakt Smile Widać tak się przejąłem że zapomniałem dopisać co za problem mam Smile

Otoż szukam pomocy w przełożeniu tego na scilaba.


 

Post Wysłany: 12 Marca 2010, Pią 10:07 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
r0g0
Może pisać książki


Dołączył: 06 Lis 2007
Posty: 138
Skąd: Białystok


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email
Zastanowił mnie ten problem... bo zajmuję się również równaniami opisanymi pochodnymi cząstkowymi...

Znalazłem zaskakujące stwierdzenie na oficjalnej stronie SCILABA:

http://wiki.scilab.org/Missing_Mathematical_features_in_Scilab

Cytat:
There is currently no Scilab component to solve PDEs


Wygląda na to, że Scilab nie dysponuje solverem do równań różniczkowych cząstkowych...
I nie wiem co z tym dalej zrobić...


 

Post Wysłany: 9 Czerwca 2010, Sro 8:45 pm Temat postu: układ Odpowiedz z cytatem
 
AUTOR:
vass



Dołączył: 09 Cze 2010
Posty: 1


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam wszystkich chodząc online znalazłem ten wątek i przypomniałem sobie o pewnym układzie równań (patrz załączony obrazek), który swego czasu próbowałem rozwiązać w SCILAB czy MATLAB, uprzejmie proszę o pomoc w implementacji problemu
rho(r,t); u(r,t); p(r,t); e(r,t)


 

Post Wysłany: 22 Czerwca 2010, Wto 10:49 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
razor87



Dołączył: 22 Cze 2010
Posty: 1


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam.
Mam do rozwiązanie w scilabie następujące równanie:
dy/dt+0.2y+0.1*y^3-1.1cos(0.75t)= 0
Nie wiem jak zapisac funkcje ydot. Problem stwarza ta trzecia potęga przy y.
Próbowałem juz ale nie wychodzi:
Kod:

function ydot=f(t,y)
 
ydot=(-0.1*y^3)-0.2*y+1.1*cos(0.75t)
endfunction


Przy wstawieniu y^2 działa.
Może ktoś pomóc ?


 

Post Wysłany: 12 Lutego 2016, Pią 6:05 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! Smile


 

Forum MATLAB Strona Główna-> Klony Matlaba (Scilab, Octave, ...)
Wyświetl posty z ostatnich:   
Idź do strony Poprzedni  1 2

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 36054949
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
Model sieci SN  (25 Maj 2017, Czw 4:56 pm)
Samochód z gwarancją  (24 Maj 2017, Sro 6:29 pm)
Serwis samochodowy Honda  (24 Maj 2017, Sro 6:01 pm)
Sterowanie wektorowe silnika PMSM  (24 Maj 2017, Sro 12:53 pm)
Zapis obrazu z scope  (23 Maj 2017, Wto 7:05 pm)
Brak funkcji tr3d  (23 Maj 2017, Wto 4:11 pm)
FFT Problem  (23 Maj 2017, Wto 1:12 pm)
Nienawidzę swojej pracy  (23 Maj 2017, Wto 9:52 am)
Modelowanie układu.  (21 Maj 2017, Nie 8:40 pm)
Małe poprawki  (21 Maj 2017, Nie 8:33 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