استفسار عن رسم بياني مهم جددددددددا


(ثلجة) #1

السلام عليكم وكل عام وانتم بألف خير

لو سمحتوا الرجاء الرد بأسرع وقت ممكن اذا كان عندي رسمه فيها 3منحنيات و

لكن المشكله اذا كانت المنحنيات مع بعض تطلع الرسمه غير صحيحه واذا كان كل

منحنى لوحده تكون الرسمه صحيحه المطلوب ان تكون المنحنيات الثلاث مع بعض

وبالصوره الصحيحه

ملف m-fail المستخدم

[LEFT]

function dM=alm(t,M)
a=1.8^-5;b=1.9^-4;c=2.27^6;
dM=[-aM(1);aM(1)-bM(2);bM(2)-cM(3);-cM(3)][/left]

ولكم جزيل الشكر


(ايجيبت) #2

عن طريق plotyy يمكنك رسم منحنيين لهم نفس الْْx ولكن مختلفيين في مدي ال y وانا حاضاوارلك علي امر اخر لو كانو 3 بس ياريت ال mfile


(ايجيبت) #3

[LEFT]this function plotyyy allow you to plot three equation have major differance in range of axis y

put the file after extract in work folder and used the form >>plotyyy(x1,y1,x2,y2,x3,y3)

ex:
x=0:360
x1=x
x2=x
y=90.*sind(x
y1=20.*cosd(x
y2=sind(x

plotyyy(x,y,x1,y1,x2,y2[/left]

i hope i understand your problem


(ثلجة) #4

شكرا جزيلا على اهتمامك بالموضوع و اتمنى ان توضح لي أكثر لأني مبتدئه و بالنسبه لـm-fail مكتوب في الطلب


(ايجيبت) #5

ممكن توضحي اكتر الثلاث معادلات المطلوب رسمهم


(ثلجة) #6

المعادله الأولى dM1\dt=-a*M1

المعادله الثانيه dM2\dt=aM1-bM2

المعادله الثالثه dM3\dt=bM2-cM3

a=1.59^-18 ,b=3.86^-9 ,c=3.14^-5

وشكرا لك


(ايجيبت) #7

المفروض في initial condition

m1,m2,m3 at initial ?


(ثلجة) #8

انا باشرحلك اللي عملته كتبت في ملف mfail التالي


function dM=alm(t,M)
a=1.8^-5;b=1.9^-4;c=2.27^6;
dM=[-aM(1);aM(1)-bM(2);bM(2)-cM(3);-cM(3)]

[RIGHT]
وحفظته بعدين استدعيته باستخدام الأمر ([t,M]=ode45(@new,[0 1^9] ,[1000 0 0]

حيث ان M01=1000,M02=0,M03=0,t=0:1^9

وطلعت لي قيم كثيره بعدها كتبت أمر الرسم التالي[/right]

plot(t,M(:,1)’-’,t,M(:,2),’-.’,t,M(:,3),’.’)

وطلعت لي الرسمه

اتمنى أن تكون وضحت المسئله لديك


(ايجيبت) #9

الطريقه دي سليمه جدا ومتجربه كتيير قبل كده

او ممكن عن طريق السميولينك لوحبيتي بس لو ممكن تحطي صوره للرسمه الناتجه عشان لوطلعت زي اللي معايا نجربها سيميوليك
بس الطريقه اللس انتي شغاله بيها اكيد صح جدا لو فيها غلط قد يكون في كتابه الكود


(ايجيبت) #10

[CENTER]

[/center]


(ايجيبت) #11

[CENTER]


بعد استخدام السيميولينك تاكد استحاله رسمهم مع بعض في ظل وجود تلك الثوابت

في المرفقات موجود ملف السيميولينك ممكن تغيري الثوابي منه[/center]


(ثلجة) #12

طبقت نفس اللي قلتلي عليه لكن تطلعلي رسالة error

[FONT=Times New Roman][COLOR=#000000][FONT=Times New Roman]
[LEFT]
Error in ==> rigid at 3
y(1)= y(2)*y(3);
[/left]

وبالنسبه للرسمه حاولت احطها بس لما احطها في الرد ما تطلع او تكون كبيره مره

[/font][/color][/font]


(ايجيبت) #13

ده الملف بتاع السيميولينك


(ايجيبت) #14

[CENTER]

يارب تكون صح لوصح حبعتلك الكود

[/center]

note: interval 0:12

(ثلجة) #15

ايوه كذا صحيحه و شكرا جزيلا ايجبت وما قصرت


(ايجيبت) #16