MATLAB Answers

Is it possible to plot this solution for a vibrating drum in a polar plot and possibly make a movie/gif for time 0 to 2 seconds?

2 ビュー (過去 30 日間)
Anik Faisal
Anik Faisal 2018 年 12 月 7 日
r=linspace(0,1,100);
theta=linspace(0,2*pi,100);
tdrum=linspace(0,2,100);
fun01=@(a) -besselj(0,2.2048*a);
A1=(1/(sqrt(2)*2.4048*(besselj(1,2.2048))^2))*integral(fun01,0,0.05);
fun02=@(a) -besselj(0,5.5201*a);
A2=(1/(sqrt(2)*5.5201*(besselj(1,5.5201))^2))*integral(fun02,0,0.05);
fun03=@(a) -besselj(0,8.6537*a);
A3=(1/(sqrt(2)*8.6537*(besselj(1,8.6537))^2))*integral(fun03,0,0.05);
fun04=@(a) -besselj(0,11.7915*a);
A4=(1/(sqrt(2)*11.7915*(besselj(1,11.7915))^2))*integral(fun04,0,0.05);
fun05=@(a) -besselj(0,14.9309*a);
A5=(1/(sqrt(2)*14.9309*(besselj(1,14.9309))^2))*integral(fun05,0,0.05);
for m=1:length(r)
for n=1:length(tdrum)
udrum1(m,n)=A1*sin(2.2048*sqrt(2)*t(n))*besselj(0,2.2048*r(m));
udrum2(m,n)=A2*sin(5.5201*sqrt(2)*t(n))*besselj(0,5.5201*r(m));
udrum3(m,n)=A3*sin(8.6537*sqrt(2)*t(n))*besselj(0,8.6537*r(m));
udrum4(m,n)=A4*sin(11.7915*sqrt(2)*t(n))*besselj(0,11.7915*r(m));
udrum5(m,n)=A5*sin(14.9309*sqrt(2)*t(n))*besselj(0,14.9309*r(m));
end
end
udrum=udrum1+udrum2+udrum3+udrum4+udrum5;

  0 件のコメント

サインインしてコメントする。

回答 (0 件)

タグ

製品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by