طراحی orbital propagator در سیمولینک

خبببین صورتم مسئله چیه ما یه propagator مدار می خوایم که توش اثرات j2 و j4 لحاظ شده باشه اونم همون سکولار effectهاش .تو گام اولم در واقع حل تحلیلی باشه حل تحلیلی باشه حالا حل عددی هم می تونه باشه چون عددی که یه دونه ode45 خروجش م.حالا فعلا همون حل تحلیلی روش تمرکز کنید ورودی ها رو چی می گیره ورودی که می گه یکی epoch time تو دستگاه utc که همون تو timescale که همون بحث سال و ماه و در واقع روز و ساعت دقیق و ثانیه و این رو تبدیل کنه به julian date برد دیگه ای که باید بگیره شیش تا در واقع المان مداری کپلریه. حالا چیزی که به نظرم جالب برای اینکه کار در واقع منظمی هم بشه یه ورودی یه بخش بشه ورودی ها رو بگیره خب می دونی که تو simulink یه بخشش در واقع کدنویسی یه بخشیش simulink یعنی شما هر مدل simulink که ببینی یه چند تا در واقع اسکریپ هم داره . یه input ها رو بگیره یه بخش تبدیل می کنه ورودی ها رو مثلا julian date و این چیزارو می تونه تبدیل کنه البته کلی دارم می گم بستگی به سلیقه و در واقع نظر کسی که سیمولینکو اون داره انجام می ده داره و بعد هم می شه بلوک propagator . و خب تو بلوک propagator هم که این ورودی ها رو می گیره در واقع شیش تا المان مداری رو می گیره و julian date رو می گیره و step time اینا مهمه و time simulation . استپت تایم هم که قابل تغییر باشه و معمولا همین چیزی پیچیده نیست و زمان simulation رو هم بگیر حالا اینم که باز محدودیت ایجاد نمی کنه برای کار . propagation انجام می شه 6 المان مداری در واقع توش ذخیره می شه بعد یه بلوک دیگری باشه که این شیش تا المان مداری رو تبدیل کنه به در واقع مولفه های مختصات اون پوزیشن و سرعت vx vy vz سرعت و پوزیشن x y z تو دستگاه eci. پس این شد در واقع بلوک دوم باز می گم بلوک دو می تونه توی بلوک هم باشه ولی خب حالا جدا هم باشه بهتره چون که براِ partهای مداری رو برای جای دیگه نیاز باشه درواقع المان های کپلری که propagate می شه و خروجی نتایجم در واقع شیش تا valadation داریم تو شیش تا پارامتر مدار کپلرین که تو شیش تا گراف جداگانه می تونه به صورت sub plot باشه و مولفه های سرعت و مولفه های پوزیشن تو دستگاه eci که اونم توی در واقع می شه شیش تا گراف که توی یه sub plot باشه در واقع ما دو تا پلات داریم که هر پلات شامل شش تا گرافه و یه رد زمینی و دو انیمیشن یکی رد زمینی یکی هم حرکت مدار حول کره زمین این هم در واقع تو کره زمین هم در واقع تو چیز توی این عادی اگه بنویسی کره زمین ثابته کره زمین که در حال چرخشه سنگم شده باشه به صورت واقعی دیگه یعنی نوع مدارک که داره حرکت می کنه تو اون زمان مثلا تو بیست چهار ساعت کره زمین هم یه دوره رو خودش می چرخه یعنی فقط یه عکس کره زمین به جای چیز نباشه به جای یه دونه کره مثلا بی رنگ عادی متلب . اینم می شه در واقع خروج ها ..همین حالا اینکه چجوری باشه شما مثلا بلوک هایی که برام مهمه وجود داشته باشه یک بلوک های اصلی جلوتر می گم ما تستمون خب یه ورودی بهش می دی مثلا تو aerospace بعد حتما متلبتون aerospace toolbox و aerospace blockset اگه اشتباه نکنم داشته باشه چون که مثلا خیلی از اینا یصورت ماژولار هست مثلا شما می خواید تبدیلش کنید به julian date اصلا یه دونه در واقع ایکون هست از اون می تونه استفاده کنه حتی می تونه هم کد بنویس چون کدشم باز هست (استاد آماده) فقط این بلاکس ها و این توول باکس رو نصب داشته باشید ای نیست همین دیگه...............ببین شما چون که خب فرمول ها که در واقع مشخصه یعنی خروج از مرکزیت و شیب و semi axis که در واقع ثابتن در اون سه تا پارامتر زاویه ای در واقع arguman prigee و آنومالی و رایتن شن تغییر می کنن ناشی از این اثرات j2 j4 سکولار خب یعنی معادلات که دم دستتون و چیز آسونه فقط بشه تمرکز روی سیمولینک باشه و این بلوک ها مثلا به صورت ارتفاع بگیره ورودی رو مثلا تو گفتم یه بلوکی باشه که حالا می تونید اسمش رو بزنید حالا اینکه می گن بلوک می تونه در واقع چیز باشه می تونه یه چرخ script باشه ولی حالا یه ورودی رو بگیره که برای کاربر خیلی آسون باشه مثل زمان در واقع شروط هایی simulation و پارامترهای مدنی مثل ارتفاع و فلان و فلان و فلان و بعد اینا رو تبدیل کنه به Semijor axis و نمی دونم حالا درجه ها رو می خواد بکنه رادیان رو یا تایم اوقات بکنه تو جولی اندیت و این چیزا یه بلوک که در واقع یه اسکریپت بلوک برای تبدیل یه propagator یه تبدیل کپلر به Eci و در واقع رسم شکل ها