تبدیل z ، تبدیل لاپلاس و عکس لاپلاس و تفکیک کسرها در Matlab
تبدیلات لاپلاس ، Z و معکوس تبدیل لاپلاس در حل بسیاری از مسائل مهندسی کاربرد دارد. در ادامه این مطلب ،کد برنامه تبدیل z ، تبدیل لاپلاس و عکس لاپلاس ، تفکیک کسرها در Matlab ارائه شده است.
function x=ilaplace my(B,A)
% To find the inverse Laplace transform of B(s)/A(s) using residue()
if ˜isnumeric(B)
[B,A]=numden(simple(B)); B=sym2poly(B); A=sym2poly(A);
end
[r,p,k]= residue(B,A); EPS = 1e-4;
N= length(r); x=[]; n=1;
while n<=N
if n>1, x = [x ’ + ’]; end
if nEPS & abs(sum(imag(p([n n+1]))))
sigma=real(p(n)); w=imag(p(n)); Kc=2*real(r(n)); Ks=-2*imag(r(n));
sigma_=num2str(sigma); w_=num2str(w); Kc_=num2str(Kc); Ks_=num2str(Ks);
if abs(sigma)>EPS
x = [x ’exp(’ sigma_ ’*t).*’];
if abs(Kc)>EPS&abs(Ks)>EPS
x = [x ’(’ Kc_ ’*cos(’ w_ ’*t) + ’ Ks_ ’*sin(’ w_ ’*t))’];
elseif abs(Kc)>EPS, x = [x Kc_ ’*cos(’ w_ ’*t)’];
else x = [x Ks_ ’*sin(’ w_ ’*t)’];
end
end
n = n+2;
elseif n<=N & abs(imag(r(n)))
if abs(p(n))>EPS, x = [x num2str(r(n)) ’*exp(’ num2str(p(n)) ’*t)’];
else x = [x num2str(r(n))];
end
n = n+1;
end
end
if ˜isempty(k), x = [x ’ + ’ num2str(k(end)) ’*dirac(t)’]; end
% To find the inverse Laplace transform of B(s)/A(s) using residue()
if ˜isnumeric(B)
[B,A]=numden(simple(B)); B=sym2poly(B); A=sym2poly(A);
end
[r,p,k]= residue(B,A); EPS = 1e-4;
N= length(r); x=[]; n=1;
while n<=N
if n>1, x = [x ’ + ’]; end
if nEPS & abs(sum(imag(p([n n+1]))))
sigma=real(p(n)); w=imag(p(n)); Kc=2*real(r(n)); Ks=-2*imag(r(n));
sigma_=num2str(sigma); w_=num2str(w); Kc_=num2str(Kc); Ks_=num2str(Ks);
if abs(sigma)>EPS
x = [x ’exp(’ sigma_ ’*t).*’];
if abs(Kc)>EPS&abs(Ks)>EPS
x = [x ’(’ Kc_ ’*cos(’ w_ ’*t) + ’ Ks_ ’*sin(’ w_ ’*t))’];
elseif abs(Kc)>EPS, x = [x Kc_ ’*cos(’ w_ ’*t)’];
else x = [x Ks_ ’*sin(’ w_ ’*t)’];
end
end
n = n+2;
elseif n<=N & abs(imag(r(n)))
if abs(p(n))>EPS, x = [x num2str(r(n)) ’*exp(’ num2str(p(n)) ’*t)’];
else x = [x num2str(r(n))];
end
n = n+1;
end
end
if ˜isempty(k), x = [x ’ + ’ num2str(k(end)) ’*dirac(t)’]; end
از ارسال دیدگاه های نا مرتبط با متن خبر، تکرار نظر دیگران، توهین به سایر کاربران و ارسال متن های طولانی خودداری نمایید.
لطفا نظرات بدون بی احترامی، افترا و توهین به مسئولان، اقلیت ها، قومیت ها و ... باشد و به طور کلی مغایرتی با اصول اخلاقی و قوانین کشور نداشته باشد.
در غیر این صورت، «برق نیوز» مطلب مورد نظر را رد یا بنا به تشخیص خود با ممیزی منتشر خواهد کرد.