¿MATLAB en electromagnética? ¿Ejemplo? código fuente!

El primero: Fluctuación de la fuente de onda sinusoidal unidimensional

% Propagación de onda sinusoidal en el espacio libre

Borrar todo

Cerrar todo

tmax=200

t= (0:1.tmax-1);%timestep

frecuencia=1.0e+9; % Frecuencia del pulso incidente (1 GHz)

omiga=2.0*pi*frequency;% Calcular omiga

c= 2.998e+8;% velocidad de onda

fuente =sin(omiga*t);%Función de onda inicial: onda sinusoidal

Gráfico

plot(fuente)%Dibuja la onda inicial

título( 'Pulso de fuente: onda sinusoidal de 1 GHz')

pausa(1);

xmax = 200; % x número de muestras.

xcenter = fix(xmax/2); % centro de la cuadrícula

xcenter = fix(xmax/2); % centro de la cuadrícula

ex = ceros (1); ,xmax); % Inicializar campo Ex

hy = zeros(1,xmax % Inicializar campo Hy

c=1.0%, resolver el coeficiente 1D obtenido en la tabla

p>

c=1. En la tabla de soluciones

para t=1:tmax

para x=2:xmax

ex( 1, x) = ex(1,x) + (hy(1,x-1) - hy(1,x) % Ver la fórmula en la tabla de soluciones

hy = zeros(1, xmax); ); % Inicializar campo Hy

c=1.0%. 0% coeficientes derivados para la hoja de trabajo 1D

end

pulso = fuente; % pulso gaussiano

% Agregar pulso en el centro de la cuadrícula

ex(1,xcenter) = pulso(t);

para x=1: (xmax-1)

hy(1,x) = hy(1,x ) + (ex(1,x) - ex(1,x+1))

fin

[Emax] = max(fuente)

[Emin] = min(fuente)

imagesc(ex,[Emin Emax]);

title('propagación de ondas 1D: fuente sin')

plot(ex);

pause(0.01);

end

También tengo una fuente de onda sinusoidal tridimensional

1D y una fuente de onda gaussiana tridimensional Fluctuación

El código de publicación excedería el límite de caracteres. Si aún lo necesitas, ¡envíame un correo electrónico!