跳至內容

File:Drum vibration mode12.gif

頁面內容不支援其他語言。
這個檔案來自維基共享資源
維基百科,自由的百科全書

Drum_vibration_mode12.gif (249 × 170 像素,檔案大小:148 KB,MIME 類型:image/gif、​循環、​19 畫格、​1.9秒)


描述 Illustration of vibrations of a drum.
日期 (UTC)
來源 self-made with MATLAB
作者 Oleg Alexandrov
其他版本

此檔案衍生的作品:

  A raster version of this image is available. It should be used in place of this vector image when superior.

File:Harmonic partials on strings.svg → File:Drum vibration mode12.gif

In general, it is better to use a good SVG version.
 
本GIF 點陣圖使用MATLAB創作。
Public domain 我,此作品的版權所有人,釋出此作品至公共領域。此授權條款在全世界均適用。
這可能在某些國家不合法,如果是的話:
我授予任何人有權利使用此作品於任何用途,除受法律約束外,不受任何限制。

Source code (MATLAB)

function VibratingDrum()

k = 1; % k-th asimuthal number and bessel function
p = 2; % 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

% To save as a GIF comment out the next the 3 lines
%   file=sprintf('Frame%d.png', 1000+iter);
%   fprintf('Saving to %s\n', file)
%   print('-dpng',  '-opengl',  '-r100', file);

  pause(0.01)
end

end

   % converted to gif with the command (run in command shell)
   % convert -antialias -loop 10000 -delay 10  -scale 50% Frame10* Drum_vibration_mode12.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;
end
   
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!')

end

說明

添加單行說明來描述出檔案所代表的內容

在此檔案描寫的項目

描繪內容

檔案歷史

點選日期/時間以檢視該時間的檔案版本。

日期/時間縮⁠圖尺寸使用者備⁠註
目前2023年3月30日 (四) 03:52於 2023年3月30日 (四) 03:52 版本的縮圖249 × 170(148 KB)Dndnrmn1Reverted to version as of 07:10, 12 January 2008 (UTC)
2008年1月16日 (三) 04:45於 2008年1月16日 (三) 04:45 版本的縮圖248 × 130(138 KB)Oleg AlexandrovStrip caption, looks better like this
2008年1月16日 (三) 04:45於 2008年1月16日 (三) 04:45 版本的縮圖248 × 130(138 KB)Oleg AlexandrovStrip caption, looks better like this
2008年1月12日 (六) 07:10於 2008年1月12日 (六) 07:10 版本的縮圖249 × 170(148 KB)Oleg Alexandrov{{Information |Description=Illustration of vibrations of a drum. |Source=self-made with MATLAB |Date=06:21, 12 January 2008 (UTC) |Author= Oleg Alexandrov |Permission= |other_versions= }} {{PD-self

下列3個頁面有用到此檔案:

全域檔案使用狀況

以下其他 wiki 使用了這個檔案:

檢視此檔案的更多全域使用狀況