¿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!