Simulação de Reflexões Múltiplas em MATLAB

f=input('Insira o valor da frequencia (Hz):');
Ex10p=input('Insira valor do Campo eletrico (V/m):'); n=input('Insira o numero de interfaces:');

l=(round(n-1));%bloco de declaração de vetores er=(round(n+1)); mr=(round(n+1)); v=(round(n+1)); b=(round(n+1)); eta=(round(n+1)); lambda=(round(n+1)); w =2*pi*f;

er=input('insira as permissividades relativas dos meios : '); %inserir 'er', 'mr' e 'l' em formato de vetor linha mr=input('insira as permeabilidades relativas dos meios : '); %para n interfaces, serão n+1 meios l=input('insira a espessura do meio em metros: '); %para n interfaces, serão n-1 espessuras

for k=1:1:n+1, %cálculo de velocidade de propagação, constante de defasamento, impedância intrínseca e comprimento de onda v(k)=3e8/sqrt(mr(k)*er(k)); b(k) = w*sqrt(mr(k)*er(k))/3e8; eta(k)=sqrt(mr(k)/er(k))*sqrt(4*pi*1e-7/(8.85e-12)); lambda(k)=2*pi/b(k); end for k=n:-1:2, %cálculo da impedância de entrada. A cada iteração, esta é armazenada no vetor das impedâncias intrínsecas eta(k)=eta(k)*(eta(k+1)*cos(b(k)*l(k-1))+1i*eta(k)*sin(b(k)*l(k-1)))/(eta(k)*cos(b(k)*l(k-1))+1i*eta(k+1)*sin(b(k)*l(k-1))); end gama=(eta(2)-eta(1))/(eta(2)+eta(1)); %cálculo dos coeficientes de reflexão e transmissão tau = 1+gama;

fprintf('\nImpedância do primeiro meio (Ohm) : ') %impressão da impedância intrínseca do 1º meio fprintf('%f',real(eta(1))) fprintf(' + i *') fprintf('%f',imag(eta(1))) fprintf('\nImpedância de entrada (Ohm) : ') %impressão da impedância de entrada das n interfaces fprintf('%f',real(eta(2))) fprintf(' + i *') fprintf('%f',imag(eta(2))) fprintf('\nCoeficiente de transmissao : ') %impressão de tau fprintf('%f',real(tau)) fprintf(' + i *') fprintf('%f',imag(tau)) fprintf('\nCoeficiente de reflexao : ') %impressão de gama fprintf('%f',real(gama)) fprintf(' + i *') fprintf('%f',imag(gama)) for t=1:20, x = -2*lambda(1):.025*lambda(1):0; %coordenadas para