MATLAB Answers

0

SciLab to Matlab transfer impedance code issue

Asked by Phillip Rauscher on 29 Jan 2019
Latest activity Commented on by Sean de Wolski
on 29 Jan 2019
I was provided he attached Scilab code for a continuing ed lesson and have gone back and forth with the instructor as to why my Matlab ode wont work to find the transfer impedance like his provided scilab code does:
k_1=2000;
k_2=5000;
k_3=1000;
m_1=2.5;
m_2=5.0;
f=linspace(1,25,100);
omega= 2*pi*f;
while i<20
R_m= input('R_m for both elastic elements =');
Z_1=R_m+i*(omega*m_1-(k_1+k_3)./omega);
Z_2=R_m+i*(omega*m_2-(k_2+k_3)./omega);
M=i*k_3./omega;
Z_12=(Z_1.*Z_2-M.^2)./M;
Level=20*log10(abs(Z_12));
plot(f,Level)
i=i+1;
An=input('Another?, Yes=1, No=0');
if An==0
abort
end
end

  0 Comments

Sign in to comment.

Products


Release

R2018a

1 Answer

Answer by Sean de Wolski
on 29 Jan 2019

abort
should be
break

  2 Comments

Thank you for the quick reply- but it continues with the same error. If it is selected to try a new R_m value the chart continues to distort without resetting like the scilab code allows for.

Sign in to comment.