Datei:Anharmonic oscillators solutions.png
Anharmonic_oscillators_solutions.png (623 × 567 Pixel, Dateigröße: 13 KB, MIME-Typ: image/png)
Diese Datei und die Informationen unter dem roten Trennstrich werden aus dem zentralen Medienarchiv Wikimedia Commons eingebunden.
Beschreibung
BeschreibungAnharmonic oscillators solutions.png | English: solution to several anharmonic oscillator problems |
Datum | |
Quelle | Eigenes Werk |
Urheber | Jkrieger |
Andere Versionen |
![]() |
PNG‑Erstellung InfoField | |
Quelltext InfoField | MATLAB codefunction [ dy ] = harmonic( t,y,m,k,l,c )
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-(c*y(2)+k*y(1))/m;
end
function [ dy ] = anharmonicA( t,y,m,k,l,c )
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-(c*y(2)+k*y(1)+l*y(1)^2)/m;
end
function [ dy ] = anharmonicB( t,y,m,k,l,c )
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-(c*y(2)+k*y(1)+l*y(1)^3)/m;
end
m=1;
k=1;
l=0.2;
c=0;
x0=1;
x01=2.4;
x02=3;
p0=0;
dy=@(t,y) anharmonicA(t,y,m,k,l,c);
dyB=@(t,y) anharmonicB(t,y,m,k,l,c);
dyh=@(t,y) harmonic(t,y,m,k,l,c);
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
PSRange=[-10 10];
subplot(3,2,1);
[T,Y] = ode45(dy,[0 40],[x0, p0],options);
[Th,Yh] = ode45(dyh,[0 40],[x0, p0],options);
[TB,YB] = ode45(dyB,[0 40],[x0, p0],options);
plot(Th,Yh(:,1), 'b-', T,Y(:,1),'r-', TB,YB(:,1),'g-');
ylim([-1.5 1.5]*x0);
xlabel('time t');
ylabel('position x(t)');
title(['solution for x_0=' num2str(x0) ', p_0=0']);
subplot(3,2,2);
for x0=-5:0.5:10
[Th,Yh] = ode45(dyh,[0 100],[x0, p0],options);
plot(Yh(:,1),Yh(:,2), 'b-');
if x0==0.5
hold on
end
end
hold off
xlim(PSRange);
ylim(PSRange);
title(['phase space plot, harmonic oscillator']);
subplot(3,2,3);
[T,Y] = ode45(dy,[0 40],[x01, p0],options);
[Th,Yh] = ode45(dyh,[0 40],[x01, p0],options);
[TB,YB] = ode45(dyB,[0 40],[x01, p0],options);
plot(Th,Yh(:,1), 'b-', T,Y(:,1),'r-', TB,YB(:,1),'g-');
xlabel('time t');
ylabel('position x(t)');
title(['solution for x_0=' num2str(x01) ', p_0=0']);
ylim([-2 2]*x01);
subplot(3,2,4);
for x0=-5:0.5:10
[T,Y] = ode45(dy,[0 100],[x0, p0],options);
plot(Y(:,1),Y(:,2), 'r-');
if x0==0.5
hold on
end
end
hold off
xlim(PSRange);
ylim(PSRange);
title(['phase space plot, asymmetric force law']);
subplot(3,2,5);
[T,Y] = ode45(dy,[0 40],[x02, p0],options);
[Th,Yh] = ode45(dyh,[0 40],[x02, p0],options);
[TB,YB] = ode45(dyB,[0 40],[x02, p0],options);
plot(Th,Yh(:,1), 'b-', T,Y(:,1),'r-', TB,YB(:,1),'g-');
xlabel('time t');
ylabel('position x(t)');
title(['solution for x_0=' num2str(x02) ', p_0=0']);
ylim([-3 3]*abs(x02));
subplot(3,2,6);
for x0=-5:0.5:10
[TB,YB] = ode45(dyB,[0 100],[x0, p0],options);
plot(YB(:,1),YB(:,2), 'g-');
if x0==0.5
hold on
end
end
hold off
xlim(PSRange);
ylim(PSRange);
xlabel('position x');
ylabel('impulse p');
title(['phase space plot, symmetric force law']);
|
Lizenz
Ich, der Urheber dieses Werkes, veröffentliche es unter der folgenden Lizenz:
![w:de:Creative Commons](https://upload.wikimedia.org/wikipedia/commons/thumb/7/79/CC_some_rights_reserved.svg/90px-CC_some_rights_reserved.svg.png)
![Namensnennung](https://upload.wikimedia.org/wikipedia/commons/thumb/1/11/Cc-by_new_white.svg/24px-Cc-by_new_white.svg.png)
![Weitergabe unter gleichen Bedingungen](https://upload.wikimedia.org/wikipedia/commons/thumb/d/df/Cc-sa_white.svg/24px-Cc-sa_white.svg.png)
Diese Datei ist unter der Creative-Commons-Lizenz „Namensnennung – Weitergabe unter gleichen Bedingungen 3.0 nicht portiert“ lizenziert.
- Dieses Werk darf von dir
- verbreitet werden – vervielfältigt, verbreitet und öffentlich zugänglich gemacht werden
- neu zusammengestellt werden – abgewandelt und bearbeitet werden
- Zu den folgenden Bedingungen:
- Namensnennung – Du musst angemessene Urheber- und Rechteangaben machen, einen Link zur Lizenz beifügen und angeben, ob Änderungen vorgenommen wurden. Diese Angaben dürfen in jeder angemessenen Art und Weise gemacht werden, allerdings nicht so, dass der Eindruck entsteht, der Lizenzgeber unterstütze gerade dich oder deine Nutzung besonders.
- Weitergabe unter gleichen Bedingungen – Wenn du das Material wiedermischst, transformierst oder darauf aufbaust, musst du deine Beiträge unter der gleichen oder einer kompatiblen Lizenz wie das Original verbreiten.
Kurzbeschreibungen
Ergänze eine einzeilige Erklärung, was diese Datei darstellt.
In dieser Datei abgebildete Objekte
Motiv
Urheber
Einige Werte ohne einen Wikidata-Eintrag
Autor (Text): Jkrieger
Wikimedia-Benutzername: Jkrieger
Urheberrechtsstatus
urheberrechtlich geschützt
Lizenz
Creative Commons Namensnennung – Weitergabe unter gleichen Bedingungen 3.0 Generisch
Datum der Gründung, Erstellung, Entstehung, Erbauung
31. Mai 2013
Quelle der Datei
durch Hochlader/in erstelltes Original
MIME-Typ
image/png
Prüfsumme
2643ad9e33f9eb4d890d27af665c4b9725f51dcd
Bestimmungsmethode: SHA-1
Dateigröße
13.194 Byte
Höhe
567 Pixel
Breite
623 Pixel
Dateiversionen
Klicke auf einen Zeitpunkt, um diese Version zu laden.
Version vom | Vorschaubild | Maße | Benutzer | Kommentar | |
---|---|---|---|---|---|
aktuell | 15:58, 31. Mai 2013 | ![]() | 623 × 567 (13 KB) | AnonMoos | adjust margins |
14:43, 31. Mai 2013 | ![]() | 710 × 614 (23 KB) | Jkrieger | {{Information |Description ={{en|1=solution to several anharmonic oscillator problems}} |Source ={{own}} |Author =Jkrieger |Date =2013-05-31 |Permission = |other_versions = }} [[Category:Theoretical p... |
Dateiverwendung
Die folgende Seite verwendet diese Datei: