# matrix loop doesnt work

1 view (last 30 days)
samuel herrera on 13 Nov 2021
Commented: Star Strider on 13 Nov 2021
how can i make a loop with the matrix a, because its not working
G=10*10-9;%
l=1*10^-3;
c=10*10^-9;
u=4*pi*10^-7;
e=(l*c)/u;
g=(u*G)/l;
vs=1;
f=60;
zs=50;
long=10:10:100;
w=2*pi*f;
z0=((j*w*l)/(g+j*w*c))^(1/2);
prop=(j*w*l*(g+j*w*c))^(1/2);
zl=30+10j;
%problem%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:length(long)
a=[(z0+zs) (z0-zs);(z0-zl)*exp(-prop*long(i)) (z0+zl)*exp(prop*long(i))];
b=[z0*vs;0];
ain=inv(a);
v=ain*b;
vp(i) = v(1);
vn(i) = v(2);
vt(i)=vp(i)+vn(i);
end
plot(long,vt)
Warning: Imaginary parts of complex X and/or Y arguments ignored.
##### 3 CommentsShow 1 older commentHide 1 older comment
Jan on 13 Nov 2021
Edited: Jan on 13 Nov 2021
"It is not working" is not a useful description of the problem. Yu know already, what's going wrong. Then it helps to solve your problem, if you share this knowledge.
Your code runs fine. So maybe the result differ from your expectations?
Star Strider on 13 Nov 2021
See if the expm function helps to produce the desired result.
It looks as though the code is attempting to model a state-space system of some sort, although it is not clear to me what it does.
.
.