Printable View
یک روش ساده همین روشی هست که گذاشتم اینجا xxx یعنی پارامتر ورودی تابع. ممکنه تابع چند تا ورودی داشته باشه که ما هم می تونیم چند تا بدیم مثل xxx,yyy
یک روش دیگه اینه که ما یک تابع بنویسیم با نام جدید و داخل اون از تابع استفاده کنیم
function ret=colmmd(xxx)
ret=colamd(xxx);
end
function [a,b]=tabejadia(x,y)
[a,b]=tabeghadim(x,y);
end
آقای من هنوز توی محاسبه تابع هزینه مشکلم بر طرف نشد. گفته بودی که شاید ضرایب a مشکل داشته باشه. بررسی کردم دیدم زیاد فرقی نمیکنه . چون اعداد خیلی کوچیکی هستن.
خروجی pi ها رو که دیدم، دیدم p1,p2,p3 با P1,P2,P3 (بزرگ) برابرند!
راه حلی ندارین؟ کد مشکل نداره؟
p های کوچیک و بزرگ که باید باهم برابرا باشند. توی مقاله جدول 3 هم نگاه کنید برابرند
مشکل ما این بود که اولا P را نداده بود و ثانیا شاید ضرایب a مشکل داشته باشه از اون فایل ورد که واسه محاسبات بود چیزی دستتون نیومد؟
نه متاسفانه. توی اون هم با محاسبات لاندا و اینا کار داره!
ضرایب a مشکلی نداره.
میگید که p بزرگ و کوچیک باید با هم برابر باشند، فقط برای سه تای اول دیگه!؟ چون بعد از ائتلاف دیگه pi ها متفاوت میشه. ضمن اینکه ما فقط سه تا مقدار برای P بزرگ داریم ولی 7 مقدار برای pi ها.
آره واسه همون سه تای اولی
این تابع هزینه ای که از برنامه بدست میآد خیلی بزرگه که سود پایین میآد. تابع هزینه رو چطوری میشه بهینه کرد؟ برنامه ای وجود داره واسه optimization؟
مشکل اینه که ما ورودی درست نداریم حالا اگه تابع را هم بخواهیم بهینه کنیم فکر کنم صورت مساله پاک می شهنقل قول:
نوشته شده توسط coronaa
کدوم ورودیهامون درست نیس؟ همه درسته دیگه!
یا P ها درست نیست یا abc چون اینها را داخل فرمول ساده aP2+bP+c میذاریم جواب درست در نمیادنقل قول:
نوشته شده توسط coronaa
برای P که همون مقدار جدول رو قرار میدیم که! (همون تقسیم Revenue بر ro). در ضمن برای تابع هزینه که از pi (کوچیک) استفاده میشه!نقل قول:
نوشته شده توسط rahnema1
اگه واقعن این piها رو اشتباه بدست میاریم پس چرا weel درست بدست میاد؟
نمیشه در جدول (3) ، از سطر سوم به بعد pi رو مجهول قرار داد تا ببینیم مقدارش چنده؟
اون سه تای اولی درست به دست بیاد بقیه اش صدقه سرمون . abc غلطهنقل قول:
نوشته شده توسط coronaa
الآن به من بگو pi سه تای اولی ممکنه اشتباه بدست اومده باشه؟ اگه درست باشه که یعنی پارامترهای b,c درست هستن.
میمونه a که من یه بار توی تابع هزینه بجای p پارامتر a رو مجهول قرار دادم بعد دوباره با a جدید اجرا کردم دیدم باز نتایج اشتباهه!
بعد مقاله ای که توی این ژورنال با پابلیشر Elsevier چاپ شده میشه اشتباه داشته باشه؟
Elsevier که وحی منزل نیست اون مقاله قبلی که در مورد otsu بود مگه مشکل نداشت؟نقل قول:
نوشته شده توسط coronaa
این راببین اصلا b و c ساده میشن
http://www.sharefile.ir/uploads/1406323933.png
این اثباتی که نشون دادی خیلی خوب بود. مرسی. وقتی خروجی pii رو توی برنامه هم می بینیم، مشاهده میشه که دقیقن همون Revenue تقسیم بر ro هست. (فقط برای سه تای اول اینطوریه)
کدوم مقاله رو میگی؟
مگه همون مقاله که در مورد تشخیص حرکت لب بود شما نفرستادی ؟نقل قول:
نوشته شده توسط coronaa
نه من نفرستادم!نقل قول:
نوشته شده توسط rahnema1
پس شما الآن میگین ممکنه a,b,c اشتباه باشه؟
P که از روی Revenue/ro بدست اومد دیگه اشتباه نیست دیگه؟
P که از روی revenu بدست اومده که مقادیر موجود در جدول گرد شده به خاطر همین نتایجی که ما در متلب به دست میاریم مخصوصا ردیف آخر مقداری با ردیف آخر جدول متفاوتهنقل قول:
نوشته شده توسط coronaa
سلام مهندس میشه ضرایب a,b,c رو توی تابع a*p^2+b*p+c=100 با برنامه خاصی بدست آورد؟ اگه p رو داشته باشیم!
من در جدول (3)، برای تابع هزینه a رو مجهول قرار میدم ، a رو بدست میارم ولی با a=[.0009 .0022 .0024] متفاوته!
اینا بدست میاد :
a=[.005534 .00118279 .000008]
عجیب نیست؟
من فرض گرفتم b,c درست باشن. پس تصمیم گرفتم a ها را به دست بیارم که این شد نتیجه ولی نمیشه
http://www.sharefile.ir/uploads/1406517119.zip
ولی از نتیجه دوباره به abc برسیم و دوباره اونها را بذاریم تا نتیجه بدست بیاد چه فایده داره؟
sqp چیه ؟نقل قول:
نوشته شده توسط rahnema1
??? Undefined function or method 'sqp' for input arguments of type
'function_handle'.
سوال : در جدول (3) در هر سطر، چرا مقدار p که از Revenue بدست میاد با مقدار p که از تابع هزینه بدست میاد یکی نیست؟ مگه نگفتیم توی سه تای اولی P,p با هم برابرند؟
توی برنامه خط آخر را با این عوض کننقل قول:
نوشته شده توسط coronaa
[x,fval]=fminsearch (@pppz,[.0009 .0022 .0024]')
مقدار P به دست میاد از Rev نه p
خب این منظورش الآن چیه ؟نقل قول:
نوشته شده توسط rahnema1
x =
0.5974
1.2274
1.0865
fval =
4.9029e+003
ـــــــــــــــــــــــــ ـــــــــــــــــــــــــ ــــــــــ
خب توی سه سطر اول که مقدار P بزرگ و p کوچیک که باید شبیه هم باشه! فرقی نمیکنه!
خب این x ها همون سه مقدار a هست که اشتباهه
با نفهمیدم این کار چه فایده ای داره
من که دقیقن توی برنامه متوجه نشدم که چیکار کردی ولی من هدفم این بود که با مقدار p بدست آمده از revenue که مسلماً باید با p تابع هزینه برابر باشه، بتونیم ضریب a رو بدست بیاریم ببینیم همونی هست که توی مقاله داده یا نه!نقل قول:
نوشته شده توسط rahnema1
نظرت در مورد حرفی که زدم چیه؟ :
توی سه سطر اول که مقدار P بزرگ و p کوچیک که باید شبیه هم باشه! فرقی نمیکنه!
در سه تای اولی p با P برابرهنقل قول:
نوشته شده توسط coronaa
پس با بدست آوردن p از revenue و قرار دادن اون در تابع هزینه و مجهول قرار دادن a میشه مقدار صحیح a رو بدست آورد؟نقل قول:
نوشته شده توسط rahnema1
این برنامه هم همینکارو میکنه؟
بله مثلا می خواهیم این کار را بکنیم ولی اصولا چه فایده ای داره از جواب به سوال برسیم تا دو باره به جواب برسیم؟ یعنی چی آخه؟!
خب بالا گفتم دیگه. میخواستم ببینم همون مقدار a بدست میاد که توی مقاله هست یا خیر که دیدیم اشتباهه. ینی نمیتونیم نتیجه بگیریم که مقدار a داخل مقاله اشتباهه؟ من اینو میخواستم بفهمم!نقل قول:
نوشته شده توسط rahnema1
فکر می کنم وارد یک دور تسلسل باطل شدیمنقل قول:
نوشته شده توسط coronaa
ببین الان 9 تا مجهول ( a,b,c هر کدوم 3 تا) داریم و تعداد معلوماتمون کمه. به کجا می خواهیم برسیم؟
من اصلن کاری به سطر دوم و سوم ندارم. فعلن میخوام همون I1 رو بررسی کنم تا ببینم مشکل من کجاست که تابع هزینه درست بدست نمیاد!
گفتیم توی این سطر باید P(بزرگ) و p(کوچیک) با هم برابر باشند که نیستن. این از مشکل اول!
چون P که از revenue بدست میاد با p که از تابع هزینه بدست میاد باید یکی باشه که نیست.
حالا اگه فرض کنیم ضریب a اشتباه هست، اگه pرو از روی revenue بدست بیاریم و بذاریم توی تابع هزینه و a رو بدست بیاریم باید مقدار صحیح a بدست بیاد.
با این تسلسل دیگه خودمم دارم دیوونه میشم!
میشه بگید چه طور این کار را کردید؟ (کد متلب)نقل قول:
نوشته شده توسط coronaa
اون اشتباهه. این مقدار درسته :نقل قول:
نوشته شده توسط rahnema1
a=[.000607 .000172 .000553]
کد :
a1=solve('1.8156*1.8156*10e8*a+1.8156*10e5+420=200 581*10e3')
a2=solve('1.8156*1.8156*10e8*a+1.8156*10e5+316=571 48.9*10e3')
a3=solve('1.8156*1.8156*10e8*a+1.8156*10e5+156=184 40*10e3')
آخ اشتباه شدم!نقل قول:
نوشته شده توسط coronaa
به نظر می رسه شما b و c را درست فرض کردی و با استفاده از اونها a را به دست آوردی در صورتی که b و c هم خودشون مشکل دارن
من اینو بدست آوردم :
a=[.6079 2.5663 8.4345]
این خط از برنامه چیکار میکنه؟
[x,fval]=fminsearch (@pppz,[.0009 .0022 .0024]')
و دستور fminsearch چیه؟
من هم فرض کردم b و c درسته اما اشتباههنقل قول:
نوشته شده توسط coronaa
این دستور برای پیدا کردن مینیمم یک تابع به کار میره
تابع مورد نظری که باید مینیمم اون بدست بیاد در واقع مقدار تابع هزینه ای که از محاسبات بدست اومده با مقدار واقعی را تفاضل به دست میاریو به توان 2 می رسونیم و با هم جمع می کنیم
ولی همون طور که گفتم 9 تا مجهول داریم این طور نمیشه
منحنی سطح زیر نمودار رو چطوری میشه بدست آورد؟
حالا هر نموداری باشه. فقط میخوام یاد بگیرم.
اگه نمودار به صورت عکس باشه، عکس را می بری داخل یه نرم افزار مثل paint یا فوتوشاپ. قسمتهای زیر منحنی را مثلا قرمز می کنی و بقیه را یک رنگ دیگه
بعد توی متلب یا هر جای دیگه هیستوگرام عکس را به دست میاری. تعداد پیکسلهای قرمز برابر خواهد بود با مساحت زیر منحنی . البته فراموش نکنید که ابعاد پیکسل را حتما در یک ضریب ضرب کنید تا به اندازه واحد محورهای منحنی بشه