矩形膜的格林函数

最新推荐文章于 2024-04-24 20:18:20 发布

西了小了贝

最新推荐文章于 2024-04-24 20:18:20 发布

阅读量331

点赞数

分类专栏：笔记

本文链接：https://blog.csdn.net/qq_39912914/article/details/102964027

版权

笔记专栏收录该内容

6 篇文章 1 订阅

订阅专栏

研究生理论声学课程学习，备忘

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);
%公式计算
for t=1:length(T)
for x=1:length(X)
   for y=1:length(Y)
    for M=1:2:m
       for N=1:2:n
         g0(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));
       end
    end
    g(x,y)=sum(sum(g0));
   end
%    G{t}=g;
end
figure(1)
imagesc(x,y,g);
end