12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849 |
- function [FE] = calc_emn_lam(no, alp, eps1, eps2)
- FE = zeros(2*no-1,2);
- te1 = eps1 - eps2;
- te2 = 1/eps1 - 1/eps2;
-
- FE(no,1) = eps1*alp + eps2*(1-alp);
- FE(no,2) = alp/eps1 + (1-alp)/eps2;
-
- ind = linspace(1,no-1,no-1);
- ifun = transpose(sin(ind*pi*alp)./(ind*pi));
- FE(no+1:2*no-1,1) = te1*ifun;
- FE(no+1:2*no-1,2) = te2*ifun;
- FE(no-1:-1:1,1) = FE(no+1:2*no-1,1);
- FE(no-1:-1:1,2) = FE(no+1:2*no-1,2);
- end
|