MOM

                Never    
%---------FUNZIONE-------------
%Copiare la funzione in un file a parte chiamato integrando.m
function f_int=integrando(delta_z,k,a,z_mp,zm,x)
        f_int=(exp(-1i*k*sqrt(a.^2+(z_mp-x*delta_z-zm)^2))/(sqrt(a.^2+(z_mp-x*delta_z-zm)^2)))
    end
%----------------------
clear
close all
clc

%Definizioni parametri elettromagnetici
freq=300E6;
c=3E8;
lambda=c/freq;
mi=4*pi*1E-7;
k=2*pi/lambda;
l=lambda/4;
a=lambda/100;
epsilon=8.854E-12;
omega=2*pi*freq;
Vg=1;

%Discretizzazione incognite
delta_z=lambda/30;
N=2*l/delta_z;

zm=linspace(-l,l-delta_z,N)
z_mp=zm+delta_z/2;
w1=0.5;
w2=w1;
x1=(-1/sqrt(3)+1)*0.5;
x2=(1/sqrt(3)+1)*0.5;
      
%Riempo la matrice
A=zeros(N);
for zi=1:N
    for zj=1:N
        if zj==zi
      
          A(zi,zj)= (mi/(4*pi))*log((sqrt(1+((4*a^2)/delta_z^2))+1)/(sqrt(1+(4*a^2)/delta_z^2)-1))-((1i*k*delta_z)/(4*pi));
        else
          A(zi,zj) =mi/(4*pi)*((w1*delta_z*integrando(delta_z,k,a,z_mp(zj),zm(zi),x1))+ w2*delta_z*(integrando(delta_z,k,a,z_mp(zj),zm(zi),x2)));
         
        end
    end
end

%Calcolo il termine noto
u = zeros(1, N);    
u(1) = 1;
u(N) = 1;

c=(cos(k*z_mp))'
s=(-1i*((2*pi)*freq)*epsilon*Vg)/(2*k)*sin(k*abs(zm))';
C=-(u*inv(A)*s)/(u*inv(A)*c)
I=inv(A)*C*c+inv(A)*s
plot(abs(I));

Raw Text