研究生理論聲學(xué)課程學(xué)習(xí)，備忘

clf clc close all clear all %% 初始化 a=0.5; b=0.5; c=3;%薄膜中的聲速 m=3; n=3; X=linspace(0,a,100); Y=linspace(0,b,100); T=linspace(0,1,10); %公式計(jì)算 for t=1:length(T) for x=1:length(X)for y=1:length(Y)for M=1:2:mfor N=1:2:ng0(M,N)=((-1).^((M+N-2)./2))*(4./(a.*b.*(c*pi*sqrt((M./a).^2+(N./b).^2))))*sin((M.*pi.*X(x))./a).* sin((N.*pi.*Y(y))./b)*sin((c*pi*sqrt((M./a).^2+(N./b).^2)).*T(t));endendg(x,y)=sum(sum(g0));end % G{t}=g; end figure(1) imagesc(x,y,g); end

總結(jié)

以上是生活随笔為你收集整理的矩形膜的格林函数的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。