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일)
이전 댓글 표시
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개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Boundary Conditions에 대해 자세히 알아보기
제품
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!