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?

Question

Anik Faisal 2018-12-7

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/434392-is-it-possible-to-plot-this-solution-for-a-vibrating-drum-in-a-polar-plot-and-possibly-make-a-movie

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;