Scribd
Upload
47 views

'R' 'r1/21' 'R - P' 'Pressure Reflection Coefficient Vs r1/r1'

The document contains MATLAB code that calculates and plots acoustic reflection and transmission coefficients as functions of the ratio of acoustic impedances. It defines functions for pressure reflection coefficient, pressure transmission coefficient, intensity reflection coefficient, and intensity transmission coefficient. It then plots these coefficients versus r1/r2 and calculates their limits as r1/r2 approaches 0, 1, and infinity. It also contains code that calculates the angle of incidence for minimum reflection of a sound wave, plots the normalized imaginary impedance versus k1L, and solves additional problems regarding acoustic reflection coefficients.

Uploaded by

Adam Foltz
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
47 views

'R' 'r1/21' 'R - P' 'Pressure Reflection Coefficient Vs r1/r1'

The document contains MATLAB code that calculates and plots acoustic reflection and transmission coefficients as functions of the ratio of acoustic impedances. It defines functions for pressure reflection coefficient, pressure transmission coefficient, intensity reflection coefficient, and intensity transmission coefficient. It then plots these coefficients versus r1/r2 and calculates their limits as r1/r2 approaches 0, 1, and infinity. It also contains code that calculates the angle of incidence for minimum reflection of a sound wave, plots the normalized imaginary impedance versus k1L, and solves additional problems regarding acoustic reflection coefficients.

Uploaded by

Adam Foltz
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

6.2.

6c)
%% Plots
x = .1:.1:10;
xi = 1./x;
% Pressure Reflection Coefficient
R_p = (xi-1)./(xi+1);
subplot(2,2,1)
plot(x,R_p,'r')
xlabel('r1/21')
ylabel('R_p')
title('Pressure Reflection Coefficient vs r1/r1')
% Pressure Transmission Coefficient
T_p = (2*xi)./(xi+1);
subplot(2,2,2)
plot(x,T_p,'g')
xlabel('r1/21')
ylabel('T_p')
title('Pressure Transmission Coefficient vs r1/r1')
% Intensity Reflection Coefficient
R_i = ((xi-1)./(xi+1)).^2;
subplot(2,2,3)
plot(x,R_i)
xlabel('r1/21')
ylabel('R_i')
title('Intensity Reflection Coefficient vs r1/r1')
% Intensity Transmission Coefficient
T_i = (4*xi)./(xi+1).^2;
subplot(2,2,4)
plot(x,T_i,'k')
xlabel('r1/21')
ylabel('T_i')
title('Intensity Transmission Coefficient vs r1/r1')
%% Limits
syms x
xi = 1/x;
% Pressure Reflection Coefficient
R_p = (xi-1)/(xi+1);
R_p_0 = limit(R_p,x,0)
R_p_1 = limit(R_p,x,1)
R_p_inf = limit(R_p,x,inf)
% Pressure Transmission Coefficient
T_p = (2*xi)/(xi+1);
T_p_0 = limit(T_p,x,0)
T_p_1 = limit(T_p,x,1)
T_p_inf = limit(T_p,x,inf)
% Intensity Reflection Coefficient
R_i = ((xi-1)/(xi+1))^2;
R_i_0 = limit(R_i,x,0)
R_i_1 = limit(R_i,x,1)
R_i_inf = limit(R_i,x,inf)
% Intensity Transmission Coefficient
T_i = (4*xi)/(xi+1)^2;
T_i_0 = limit(T_i,x,0)
T_i_1 = limit(T_i,x,1)
T_i_inf = limit(T_i,x,inf)

Pressure Transmission Coefficient vs r1/r1


2

0.5

1.5

Pressure Reflection Coefficient vs r1/r1


1

-0.5

-1

0.5

10

Intensity Reflection Coefficient vs r1/r1

0.6

0.9

0.5

0.8

0.4

0.7

0.3

0.6

0.2

0.5

0.1

0.4
0

6
r1/21

10

6
r1/21

R_p_0 = 1

R_i_0 = 1
T_i_0 = 0

** When r1/r2 goes to 0, the intensity wave is completely reflected, while the
pressure wave reflects as the same magnitude of the incident wave, and the
transmitted pressure wave doubles in magnitude.

R_p_1 = 0
T_p_1 = 1
R_i_1 = 0

** When r1/r2 goes to 1, the intensity wave is completely transmitted, as is the


pressure wave.

T_i_1 = 1

R_p_inf = -1
T_p_inf = 0
R_i_inf = 1
T_i_inf = 0

10

Intensity Transmission Coefficient vs r1/r1

0.7

T_p_0 = 2

r1/21

r1/21

** When r1/r2 goes to infinity, the intensity wave is completely reflected, as is


the pressure wave, but the pressure wave is out of phase.

10

6.6.1)
%% Part (a)
% angle of incidence
syms x
% normal specific acoustic resistance
rn = 900;
% normal specific acoustic reactance
xn = -1200;
% acoustic impedance of air
r1 = 415;
% power reflection coefficient
Rii = ((rn-r1/cos(x))^2+xn^2)/((rn+r1/cos(x))^2+xn^2);
% differentiate with respect to the angle of incidence
par1 = diff(Rii);
% set to 0 and solve for the angle of incidence
par2 = solve(par1 == 0);
% convert from symbolic to double
par3 = double(par2);
% convert from radians to degrees and find the positive angle
ans_a = max(180/pi*par3)
%% Part (b)
% convert symbolic expression to a MATLAB function
Rii = matlabFunction(Rii);
% find the value for the power reflection coefficient at 80 degrees
ans_b = Rii(80*pi/180)
%% Part (b)
% find the value for the power reflection coefficient at 0 degrees
ans_c = Rii(0)
ans_a = 73.9386 deg.
ans_b = 0.2984
ans_c = 0.5286

Add. Prob ii)


% k1L
x = 0:.01:4*pi;
% normalized imaginary impedance
y = -cot(x);
% plot
plot(x,y)
axis([0,4*pi,-10,10])
title('Normalized Imaginary Impedance vs k1L')
xlabel('k1L')
ylabel('Normalized Imaginary Impedance')

Normalized Imaginary Impedance vs k1L


10

Normalized Imaginary Impedance

8
6
4
2
0
-2
-4
-6
-8
-10

6
k1L

10

12

You might also like