کد خبر: ۲۵۹۸۴
تاریخ انتشار : ۱۰:۴۵ - ۰۵ مهر ۱۳۹۶
تبدیلات لاپلاس ، Z و معکوس تبدیل لاپلاس در حل بسیاری از مسائل مهندسی کاربرد دارد. در ادامه این مطلب ،کد برنامه تبدیل 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
ارسال نظر قوانین ارسال نظر
لطفا از نوشتن با حروف لاتین (فینگلیش) خودداری نمایید.
از ارسال دیدگاه های نا مرتبط با متن خبر، تکرار نظر دیگران، توهین به سایر کاربران و ارسال متن های طولانی خودداری نمایید.
لطفا نظرات بدون بی احترامی، افترا و توهین به مسئولان، اقلیت ها، قومیت ها و ... باشد و به طور کلی مغایرتی با اصول اخلاقی و قوانین کشور نداشته باشد.
در غیر این صورت، «برق نیوز» مطلب مورد نظر را رد یا بنا به تشخیص خود با ممیزی منتشر خواهد کرد.
نتیجه عبارت زیر را وارد کنید
captcha =
وضعیت انتشار و پاسخ به ایمیل شما اطلاع رسانی میشود.
پربازدیدها
برق در شبکه های اجتماعی
اخبار عمومی برق نیوز
عکس و فیلم
پربحث ترین ها
آخرین اخبار