Skip to content

New element has been added, the ideal transformer. The stamp can be f… #5

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
manelvelasco wants to merge 1 commit into
base: master
Choose a base branch
from
Open

New element has been added, the ideal transformer. The stamp can be f… #5

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
5 changes: 5 additions & 0 deletions exampleT.cir
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
V1 1 0 symbolic
R1 1 2 1000
T1 2 3 4 5 20
R2 4 5 1000
R3 3 0 1000
50 changes: 46 additions & 4 deletions scam.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,7 @@
fprintf('Netlist:');
type(fname)
fid = fopen(fname);
fileIn=textscan(fid,'%s %s %s %s %s %s'); % Read file (up to 6 items per line
fileIn=textscan(fid,'%s %s %s %s %s %s','CommentStyle','#'); % Read file (up to 6 items per line), lines starting with "#" are ingnored
% Split each line into 6 columns, the meaning of the last 3 columns will
% vary from element to element. The first item is always the name of the
% circuit element and the second and third items are always node numbers.
Expand All @@ -35,7 +35,9 @@
for k1=1:nLines % Check all lines to find voltage sources
switch Name{k1}(1)
case {'V', 'O', 'E', 'H'} % These are the circuit elements with
m = m + 1; % We have voltage source, increment m.
m = m + 1; % voltage source, increment m.
case {'T'} % Ideal Transformer equivalent to 2 sources
m = m + 2;
end
end

Expand Down Expand Up @@ -154,7 +156,39 @@
end

j{vsCnt}=['I_' Name{k1}]; % Add current through source to unknowns

% Ideal transformer
case 'T' % TRAFO

n3 = str2double(arg3{k1}); % Control voltage, pos side
n4 = str2double(arg4{k1}); % Control voltage, neg side

% Change B and C matrices as appropriate for output nodes.
% (if node is not ground)

if n1~=0
B{n1,vsCnt+1} = [B{n1,vsCnt+1} '+ 1'];
C{vsCnt+2,n1} = [C{vsCnt+2,n1} ' - ' Name{k1}];
end
if n2~=0
B{n2,vsCnt+1} = [B{n2,vsCnt+1} '- 1'];
C{vsCnt+2,n2}= [C{vsCnt+2,n2} ' + ' Name{k1}];
end
if n3~=0
B{n3,vsCnt+2} = [B{n3,vsCnt+2} ' + 1'];
C{vsCnt+2, n3} = [C{vsCnt+2, n3} ' + 1'];
end
if n4~=0
B{n4,vsCnt+2} = [B{n4,vsCnt+2} ' - 1'];
C{vsCnt+2, n4} = [C{vsCnt+2, n4} ' - 1'];
end
D{vsCnt+1,vsCnt+1}=[D{vsCnt+1,vsCnt+1} '+ 1'];
D{vsCnt+1,vsCnt+2}=[D{vsCnt+1,vsCnt+2} ' - ' Name{k1}];

j{vsCnt+1}=['I_' Name{k1} '_P']; % Add current through source to unknowns
j{vsCnt+2}=['I_' Name{k1} '_S']; % Add current through source to unknowns
vsCnt = vsCnt + 2; % Keep track of number of voltage sources

% Voltage controlled current source
case 'G' % VCCS GXXXXXXX N+ N- NC+ NC- VALUE
nc1 = str2double(arg3{k1}); % Control voltage, pos side
Expand Down Expand Up @@ -273,7 +307,7 @@
disp(z);

% Find all variables in matrices (symvar) and make them symbolic (syms)
syms([symvar(A), symvar(x), symvar(z)]);
syms([symvar(A), symvar(x), symvar(z)],'real');

% Displey the matrix equation
fprintf('\nThe matrix equation: \n');
Expand All @@ -289,6 +323,14 @@

fprintf('\nThe solution: \n');
disp(x==eval(x))
if rank(A)~=size(A(1,:))
NullA=null(A);
[m,n]=size(NullA);
fprintf('\n There is an infinite set of solutions spaned by the vector\n')
syms K [n 1] real
disp(x==null(A)*K)

end

%% Lastly, assign any numeric values to symbolic variables.
% Go through the elements a line at a time and see if the values are valid
Expand All @@ -304,7 +346,7 @@
case {'H', 'F'}
[num, status] = str2num(arg4{k1}); %#ok<ST2NM>
% Elements defined by five variables, arg5 is the value.
case {'E', 'G'}
case {'E', 'G', 'T'}
[num, status] = str2num(arg5{k1}); %#ok<ST2NM>
end
if status % status will be true if the argument was a valid number.
Expand Down