Condiciones de contorno- Propagación de ondas planas

Parte 1. Propagación de ondas planas

1. Desarrolle una función en MATLAB que a partir de las características del medio (ε, μ, σ) y la frecuencia calcule la constante de atenuación, fase, velocidad y longitud de onda). Nota: la función debe indicar que tipo de medio es. (Espacio libre, dieléctrico, dieléctrico disipativo o buen conductor).

function [A, B, up, L, E]= OP(er,ur,g,f)

er=72; ur=1; g=4; f=(10^8)/(2*pi);

%Permitvidad y permeabilidad absoluta

e0=(1/(36*pi))*10^(-9)

u0=4*pi*10^(-7);

%Permitvidad y permeabilidad relativa

e=er*e0

u=ur*u0;

%Otros parámetros

w=2*pi*f

z=w*e

TP=g/(w.*e) %Tangente de pérdidas.

x1=(1+TP^2)^(1/2)

x2=w*(u*e/2)^(1/2)

%Costante de atenuación.

A=w*(((u*e)/2)*((1+TP^2)^(1/2)-1))^(1/2)

%Costante de beta.

B=w*(((u*e)/2)*((1+TP^2)^(1/2)+1))^(1/2)

%Velocidad de propagación

up=w/B

%Longitud de onda

L=(2*pi)/B

%Tipo de medio

if g==0

    if e==e0 & u==u0

         fprintf('El medio es el vacío.')

    elseif (e~=e0 && u==u0) || g < z

         fprintf('Dieléctricos sin pérdidas.')

    end

elseif g~=0

    if e==e0 & u==u0

    fprintf('Dieléctricos disipativos.')  

    elseif (g > 1000 && e==e0 && u~=u0) || g > z

    fprintf('Buenos conductores.')          

 end

end

2. Utilizando la función realizada en el inciso anterior para desarrolle otra función en MATLAB que dibuje la intensidad del campo eléctrico en función de la distancia, utilice t como constante, t=0. Nota: exprésela la distancia en longitudes de onda, y el campo en V/m  (dibuje hasta 3λ) debe solicitar al usuario E0 y las características del medio.

function [Grafica]= Intensidad_De_CampoE(E0,w, L, d)

...