سلام دوست من دستور زیر به اسانی پاسخ مسئله را به شما میدهد.
(x,y]=ode45(@(x,y)(1.5-.8*y)/(.5+.8*x),[-2,2],[0]);plot(x,y]
و شکل پاسخ به صورت زیر می شود همانگونه که از شکل پاسخ پیداست این معادلات در فاصله ی بین 1- و 0.5- مجانب دارد که به همین دلیل صورت درستی از پاسخ را نمایش نمی دهد.
مشاهده پیوست 188715
به همین دلیل شما میتوانید بازه را به دو قسمت تقیسم کنید یکی از 2- تا مجانب و دیگری از مجانب تا 2+ و از دستور بالا دوبار استفاده کنید.
اما برای اینکه متوجه شوید این دستور درست است من با فرض اینکه[ x=[0,2 است و مسئله مجانب ندارد میتوان نمودار را به دقت دید.
(x,y]=ode45(@(x,y)(1.5-.8*y)/(.5+.8*x),[0,2],[0]);plot(x,y]
مشاهده پیوست 188716
برای دیدن آموزشهای بیشتر در زمینه ی متلب لطفا به سایت تخصصی ما مراجعه نمایید.
مکانیک سافت سایت تخصصی مهندسی مکانیک و کنترل