بدست آوردن متغیرهای حالت یک سیستم در mfile متلب

[قاسمی2020]

با سلام خدمت دوستان گرامی
اگر سیستمی به صورت فضای حالت داده شده باشد،و تنها مجهول ما متغیرهای حالت آن باشد،با چه دستوراتی می توان آن را در mfile متلب بدست آورد؟
اگر بخوام از طریق لاپلاس آن را حل کنم از چه دستوراتی باید استفاده کنم؟
با تشکر

[BlueAbble]

پاسخ : بدست آوردن متغیرهای حالت یک سیستم در mfile متلب

سلام
معادلات حالت رو دارین، معنیش اینه که ماتریسهای A,B,C,D معلومند. ورودی u هم که قاعدتا معلومه. مقدار اولیه x که شرایط اولیه هست هم معلومه، حالا هدف اینه که x و y رو بدست بیارین.
.
x = Ax + Bu
y = Cx + Du
با یک حلقه for میشه اینکارو کرد. در این حلقه محاسبات زیر انجام میشه:
1- ابتدا از معادله اول xdot محاسبه میشه.
2- از معادله دوم y محاسبه میشه.
3- حالا قبل از شروع دور بعدی، x جدید از روی xdot و x قبلی محاسبه میشه:
xnew = xold + xdot*dt
که در اون dt پله های زمانی است که طی اون قرار است x محاسبه شود. هر چه dt کوچکتر باشد، رزولوشن زمانی بهتر است.

اگر هدفتون صرفا بدست آوردن y می بود میتونستید از دستور ss2tf هم استفاده کنید. این دستور ماتریسهای فضای حالت رو از شما میگیره و تابع تبدیل در فضای لاپلاس رو میده. با این تابع تبدیل میشه خروجیهارو محاسبه کرد.

[قاسمی2020]

پاسخ : بدست آوردن متغیرهای حالت یک سیستم در mfile متلب

با سلام و خسته نباشید خدمت شما
ضمن تشکر از پاسخ شما،از جنابعالی خواهشمندم در صورت امکان برنامه ای که اشاره کردید که می شه با حلقه for و با استفاده از Xnew=Xold+Xdot نوشت که با استفاده از برنامه مذکور متغیرهای حالت[X(t)] را می توان بدست آورد.آن را بنویسید؟!
با تشکر فراوان