File:Drum vibration mode13.gif
Aus Wikimedia Commons, dem freien Medienarchiv
Zur Navigation springen
Zur Suche springen
Drum_vibration_mode13.gif (250 × 130 Pixel, Dateigröße: 137 KB, MIME-Typ: image/gif, Endlosschleife, 19 Bilder, 1,9 s)
Dateiinformationen
Strukturierte Daten
Bildtexte
BeschreibungDrum vibration mode13.gif | Illustration of vibrations of a drum. |
Datum | (UTC) |
Quelle | self-made with MATLAB |
Urheber | Oleg Alexandrov |
Andere Versionen |
Abgeleitete Werke dieser Datei: |
Public domainPublic domainfalsefalse |
Ich, der Urheberrechtsinhaber dieses Werkes, veröffentliche es als gemeinfrei. Dies gilt weltweit. In manchen Staaten könnte dies rechtlich nicht möglich sein. Sofern dies der Fall ist: Ich gewähre jedem das bedingungslose Recht, dieses Werk für jedweden Zweck zu nutzen, es sei denn, Bedingungen sind gesetzlich erforderlich. |
Source code (MATLAB)
[Bearbeiten]function main()
k = 1; % k-th asimuthal number and bessel function
p = 3; % p-th bessel root
q=find_pth_bessel_root(k, p);
N=20; % used for plotting
% Get a grid
R1=linspace(0.0, 1.0, N);
Theta1=linspace(0.0, 2*pi, N);
[R, Theta]=meshgrid(R1, Theta1);
X=R.*cos(Theta);
Y=R.*sin(Theta);
T=linspace(0.0, 2*pi/q, N); T=T(1:(N-1));
for iter=1:length(T);
t = T(iter);
Z=sin(q*t)*besselj(k, q*R).*cos(k*Theta);
figure(1); clf;
surf(X, Y, Z);
caxis([-1, 1]);
shading faceted;
colormap autumn;
% viewing angle
view(108, 42);
axis([-1, 1, -1, 1, -1, 1]);
axis off;
file=sprintf('Frame_mode%d%d_%d.png', k, p, 1000+iter);
disp(sprintf('Saving to %s', file));
print('-dpng', '-zbuffer', '-r100', file);
pause(0.1);
end
% converted to gif with the command
% convert -antialias -loop 10000 -delay 10 -scale 50% Frame_mode13* Drum_vibration_mode13.gif
function r = find_pth_bessel_root(k, p)
% a dummy way of finding the root, just get a small interval where the root is
X=0.5:0.5:(10*p+1); Y = besselj(k, X);
[a, b] = find_nthroot(X, Y, p);
X=a:0.01:b; Y = besselj(k, X);
[a, b] = find_nthroot(X, Y, 1);
X=a:0.0001:b; Y = besselj(k, X);
[a, b] = find_nthroot(X, Y, 1);
r=(a+b)/2;
function [a, b] = find_nthroot(X, Y, n)
l=0;
m=length(X);
for i=1:(m-1)
if ( Y(i) >= 0 & Y(i+1) <= 0 ) | ( Y(i) <= 0 & Y(i+1) >= 0 )
l=l+1;
end
if l==n
a=X(i); b=X(i+1);
%disp(sprintf('Error in finding the root %0.9g', b-a));
return;
end
end
disp('Root not found!');
Dateiversionen
Klicke auf einen Zeitpunkt, um diese Version zu laden.
Version vom | Vorschaubild | Maße | Benutzer | Kommentar | |
---|---|---|---|---|---|
aktuell | 21:00, 8. Jun. 2008 | 250 × 130 (137 KB) | Oleg Alexandrov (Diskussion | Beiträge) | {{Information |Description={{en|1=x}} |Source=Own work by uploader |Author=Oleg Alexandrov |Date=x |Permission=x |other_versions=x }} x {{ImageUpload|full}} |
Du kannst diese Datei nicht überschreiben.
Dateiverwendung
Keine Seiten verwenden diese Datei.
Globale Dateiverwendung
Die nachfolgenden anderen Wikis verwenden diese Datei:
- Verwendung auf ca.wikipedia.org
- Verwendung auf de.wikipedia.org
- Verwendung auf en.wikipedia.org
- Verwendung auf en.wikiversity.org
- Verwendung auf es.wikipedia.org
- Verwendung auf fr.wikipedia.org
- Verwendung auf ko.wikipedia.org
- Verwendung auf sq.wikipedia.org
- Verwendung auf sr.wikipedia.org
- Verwendung auf te.wikipedia.org