MATLAB Answers

SciLab to Matlab transfer impedance code issue

4 views (last 30 days)
Phillip Rauscher
Phillip Rauscher on 29 Jan 2019
Commented: 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.

Answers (1)

Sean de Wolski
Sean de Wolski on 29 Jan 2019
abort
should be
break

  2 Comments

Phillip Rauscher
Phillip Rauscher on 29 Jan 2019
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.

Sign in to answer this question.

Products


Release

R2018a