প্রকৌশল শব্দবিজ্ঞান/গাড়ির মাফলার: অ্যানিমেশন: কোড

% প্রস্থ-বর্ধিত নালী (Duct with enlargement) - সর্বশেষ পরিবর্তন ১১/০৩/৯৯

% প্রফেসর এল. মঁজো-র সৌজন্যে পুনঃপ্রকাশিত

clear all;

clc;

close all;

Nsteps=30;

f=86*2; % ফ্রিকোয়েন্সি

rho=1.15;

c=344.;

rhoc=rho*c;

l=1; % নালীর দৈর্ঘ্য

a1=0.0762/2;

S1=pi*a1^2;

a2=0.1524/2;

S2=pi*a2^2;

m=S2/S1;

w=2*pi*f;

k=w/c;

t=linspace(2*pi/Nsteps/w,2*pi/w,Nsteps); % প্রতি পর্যায়ে Nsteps সময় ধাপ

alpha1=2.93e-5/a1.*sqrt(f);

alpha2=2.93e-5/a2.*sqrt(f);

khat=k-i*alpha1;

kl=khat*l;

U0=0.01;

R=(1-m)/(m+1);

A=rho*c*U0./(exp(i*kl)-R*exp(-i*kl));

B=R*A;

C=A+B;

x=linspace(0,2,100); % গণনার ক্ষেত্র

y=linspace(-0.2,0.2,32);

for m=1:50

for k=8:24

pt(k,m)=A*exp(i*khat*(l-x(m)))-B*exp(-i*khat*(l-x(m)));

end

end

for m=51:100

for k=1:32

pt(k,m)=C*exp(i*khat*(l-x(m)));

end

end

pmag=abs(pt);

pphase=angle(pt);

clim=max(max(pmag));

V=[-clim clim]; % রঙের সীমা নির্ধারণ

pt=pmag.*exp(i*pphase); % লাল: সর্বোচ্চ ধনাত্মক চাপ

% নীল: সর্বনিম্ন ঋণাত্মক চাপ

figure(1)

pt2=real(pt);

H=pcolor(x,y,pt2);

axis equal % বিকৃতি রোধে

axis([0 2 -.2 .2])

M=moviein(Nsteps); % মুভি অ্যানিমেশনের জন্য

for k=1:Nsteps

pt2=real(pt*exp(i*w*t(k)));

H=pcolor(x,y,pt2);

axis equal

xlabel('x');

ylabel('y');

title('চাপ')

shading interp;

caxis(V);

axis([0 2 -.2 .2])

M(:,k)=getframe;

% ওয়েবে প্রকাশের জন্য চিত্র সংরক্ষণ (বিটম্যাপ)

% filnme=strcat('spher',int2str(k));

% eval(['print -dbitmap ' filnme ' -f']);

end

movie(M,20)