Clear Filters
Clear Filters

matrix loop doesnt work

1 view (last 30 days)
samuel herrera
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 Comments
Jan
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
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.
.
.

Sign in to comment.

Answers (0)

Categories

Find more on Loops and Conditional Statements in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!