Könnte mir jemand helfen, erklären, wie ich eine Differentialgleichung der Form d'=A*d numerisch löse?
Wobei A eine 2x2 Matrix ist und d offensichtlich ein Zeilenvektor.
Bisher habe ich das nur mit Vektoren aus einem Vektorfeld und dann plotten lassen.
Und die numerische Lösung würde ich gerne als Vergleich mit in meinen bisherigen plot miteinbeziehen.
Ich habe schon versucht das mit ode45 zu lösen weiß aber nicht ob das dann so funktioniert wie ich das gerne hätte.
Code:
[X, Y] = meshgrid(-2:0.25:3);
x = reshape(X, 1, []);
y = reshape(Y, 1, []);
M=[6,9;-1,0];
[a,b] = ode45(@(a,b) M*[x;y],[0,5],0);
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Ich hatte mir das auch schon angeguckt aber das ist ja kein Numerisches Verfahren.
Ich habe den Code von dir von dem Link mal benutzt um das Verfahren mit ode45 hineinzuplotten jedoch glaube ich nicht, dass das so richtig aussieht.
Vllt könntest du mal rübergucken?
Code:
[X, Y] = meshgrid(-5:0.5:5);
x = reshape(X, 1, []);
y = reshape(Y, 1, []);
M=[-1,0;0,-2];
Dxy = M * [x;y];
[a,b] = ode45(@(t,d) M*d,[0,5],[5;5]) quiver(x,y,Dxy(1,:), Dxy(2,:)) hold on
plot(a,b(:,1))
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
RSS