برنامه متلب محاسبه ژیزمان با استفاده از زاویه حامل

برنامه متلب (matlab) محاسبه ژیزمان با استفاده از زاویه حامل:

مثال:
X1=100 Y1=100 X2=150 Y2=200
V=Arc tan (150-100/200-100)=50/100=+∆x/+∆y
چون ∆x,∆y مثبت می باشد پس G=V خروجی:
X1=200 Y1=250 X2=300 Y2=150
V=Arc tan (300-200/150-250)=100/-100=+∆x/-∆y
چون ∆x مثبت و ∆y منفی می باشد پس G=180-V خروجی
X1=150 Y1=200 X2=100 Y2=100
V=Arc tan (100-150/100-200)=-50/-100=-∆x/-∆y
چون ∆x,∆y منفی می باشد پس G=180+V خروجی
X1=200 Y1=150 X2=150 Y2=200
V=Arc tan (150-200/200-150)=-50/50=-∆x/+∆y
چون ∆xمنفی و ∆y مثبت می باشد پس G=360-V خروجی
در این مسئله ابتدا برنامه باید تشخیص بدهد که وضعیت ∆x,∆y چه است (مثبت یا منفی ) سپس ژیزمان را حساب کند یعنی اگر من ( X1,Y1 ) و ( X2,Y2 ) را دارم برنامه تشخیص بدهد که در کدام ربع قرار دارد سپس ژیزمان را حساب کند.
کلیه ورودیها مختصات X,Y
خروجی ژیزمان G