zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

modified to work with octave-3.2.4

pull/1/head
Sebastian Held 2012-08-19 00:12:50 +02:00
parent ec1f3ec387
commit 3dc6cbebbd
1 changed files with 29 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
matlab/examples/microstrip/directional_coupler.m
View File

@ -13,7 +13,7 @@ function directional_coupler
% %
% Tested with % Tested with
% - Matlab 2010b % - Matlab 2010b
% - Octave 3.3.54 % - Octave 3.2.4
% - openEMS v0.0.17 % - openEMS v0.0.17
% %
% (C) 2010 Sebastian Held <sebastian.held@gmx.de> % (C) 2010 Sebastian Held <sebastian.held@gmx.de>
@ -157,9 +157,9 @@ FDTD = SetBoundaryCond( FDTD, BC );
%% Write openEMS compatible xml-file %% Write openEMS compatible xml-file
if runSimulation if runSimulation
[~,~,~] = rmdir(Sim_Path,'s'); [dummy,dummy,dummy] = rmdir(Sim_Path,'s');
end end
[~,~,~] = mkdir(Sim_Path); [dummy,dummy,dummy] = mkdir(Sim_Path);
WriteOpenEMS([Sim_Path '/' Sim_CSX],FDTD,CSX); WriteOpenEMS([Sim_Path '/' Sim_CSX],FDTD,CSX);
if showStructure if showStructure
@ -233,3 +233,29 @@ else
axis equal axis equal
axis off axis off
end end
function s = y2s(y, ZL)
% S = y2s(Y, ZL)
%
% Admittance to Scattering transformation
% for square matrices at multiple frequencies
%
% ZL defaults to 50 Ohm
if nargin < 2
ZL = 50;
end
if size(size(y),2) > 2
nF = size(y,3);
else
nF = 1;
end
I = diag(ones(1, size(y,2)))/ZL;
for i=1:nF
%s(:,:,i) = inv(I+y(:,:,i)) * (I-y(:,:,i));
s(:,:,i) = (I+y(:,:,i)) \ (I-y(:,:,i));
end