Printable View
در ابتدا اگه مجموعه تهی باشه صفر بر می گردونهنقل قول:
نوشته شده توسط coronaa
سپس در S می گردیم اگر تعداد اعضای s کمتر از از اون Smimi بود همین طور شمارش می کنیم تا برسیم به جایی که تعداد اعضای S با Smimi برابر باشه با استفاده از ismember عضوی از S که با Smini برابره پیدا کرده با شماره های اندکس قبلی که جمع کرده ایم جمع می کنیم در نتیجه اندکس مورد نظر به دست می آید بعدش هم v متناظر با اون اندکس بدست میاریم
منظور از Smini چیه؟نقل قول:
نوشته شده توسط rahnema1
و rows در این خط چیکار میکنه؟
[val,dx]=ismember(Smini, S{i},'rows');
smini همون s-i هست و rows باعث میشه مقایسه به جای اینکه عنصر به عنصر باشه، سطر به سطر باشهنقل قول:
نوشته شده توسط coronaa
مهندس ببخشید که من هی سوال می پرسم.منظور از [val,dx] چیه؟
val که کاری باهاش نداریم اما dx شماره سطری هست که smini درسلول مورد نظر s پیدا میشهنقل قول:
نوشته شده توسط coronaa
ret رو بدست میاریم که چی بشه؟
چرا یه عدد اعشاریه؟
در ضمن، Smini رو که از قبل معرفی نکردی چطوری مساوی صفر قرار دادی؟ توی خط :
if (numel(Smini)==0)
چرا ما تابع تعریف می کنیم؟نقل قول:
نوشته شده توسط coronaa
برای اینکه یک کار تکراری که چند بار می خواهید انجام بدیدم یک بار کد اون را بنویسیم و تنها با صدا زدن نام اون و ورد چند تا پارامتر اون کد ( یعنی تابع) را اجرا کنیم
وقتی وسط برنامه این تابع داره صدا زده میشه
indexV(setdiff( S{i}(j,:) ,a(z)),S,v)
سه تا آرگیومنت داره بهت وارد میشه که با کاما جاد شدند. اولین مورد یعنی
setdiff( S{i}(j,:) ,a(z))
همون Smini هست
دومین و سومین مورد هم اینجا نامهاشون با اسامی وافعی متغیرها یکی قرار داده شده در صورتی که می تونستیم یک اسم دیگه بذاریم مثلا تابع را به این صورت زیر تعریف می کردیم باز هم فرقی نمی کرد
function ret= indexV(Smini,SSS,vvv)
ret=0;
if (numel(Smini)==0)
return;
end
for i=1:numel(SSS)
if(size(SSS{i},2)<numel(Smini))
ret=ret+size(SSS{i},1);
else
[val,dx]=ismember(Smini, SSS{i},'rows');
ret=ret+dx;
break;
end
end
ret=vvv(ret);
end
واینکه چرا اعشاری هست ما در پایان نوشتیم
ret=v(ret);
تا قبل از این دستور معنای ret اندکسی از v بود که متناظر بود با s-i
وقتی که این دستور را زدیم یعنی اون V با این اندکس را استخراج کن و دوباره اون را به جای ret بذار. دراین حالت ret دیگه اندکس نیست بلکه v مورد نظر هست
خدا عمرت بده. همون اول اینطوری توضیح میدادی :لبخندساده:
توضیحت عالی بود. ممنونم.
مهندس میشه بگی هر یک از حلقه ها چیکار میکنن؟
for i=1:numel(S)
for j=1:size(S{i},1)
for z=1:numel(n)
if (ismember(n(z),S{i}(j,:)))
phi(z)=phi(z)+factorial(numel(S{i}(j,:))-1).* factorial(numel(n) -numel(S{i} (j,:))).* (v(k)-indexV(setdiff( S{i}(j,:) ,n(z)),S,v))./ factorial(numel(n));
end
end
k=k+1;
end
end
ُS یک سلوله با 3 تا خونه که داخل هر خونه هم یک ماتریس قرار داره (اگه چاپ کنی معلوم میشه)
دو تا for اولی واسه پیمایش S هست
for سومی هم واسه پیمایش phi هست چون phi 3 تا خونه داره
این معادلات جواب ندارن یا دستورش اشتباهه؟
h=solve('x1+x2=39815.1','x1+x3=35644.2','x2+x3=124 49.0','x1+x2+x3=53413.2')
h=[h.x1 h.x2 h.x3]
میگم می خواهی nucleolus را حساب کنی؟نقل قول:
نوشته شده توسط coronaa
نه. بخش 3.2 (هسته)
معادلات 15 تا 18
توی این سایت یه سری فایل زیپ مخصوص تئوری بازی ها در متلب هست که shapely و core و nucleolus و چیزهای مختلف را حساب می کنه دیگه فکر کنم دردسر نداشته باشه
http://webs.uvigo.es/mmiras/TUGlab
ممنونم. این خیلی عالیهنقل قول:
نوشته شده توسط rahnema1
phi رو با مقداری که خودمون بدست آوریم تست کردم. درست جواب داد. ولی سر هم کردنش سخته.
بعضی فایلهاش اسپانیاییه. ولی راهنمای اصلیش انگلیسیه. فعلن دارم کار میکنم.
مرسی
مهندس این برای قسمت nucleolus خطا میده بصورت زیر :
??? Undefined function or method 'colmmd' for input arguments of
type 'double'.
مفهومش چیه؟ یعنی فانکشن colmmd رو نمیشناسه؟ توی اسکریپتاش پیداش نکردم!
این یک دستور در نسخه های قدیمی متلب هست ( توی یک سایت توضیح داده بود) این را بذاری درست میشهنقل قول:
نوشته شده توسط coronaa
colmmd=@(xxx) colamd(xxx);
کجای برنامه باید بذارم؟
آقا حل شد. ممنونم.نقل قول:
نوشته شده توسط coronaa
چطوری یه دستوری که معرفی نشده رو باید معرفی کرد مثل همین؟
قبلنم چند بار انجام داده بودی!
یک روش ساده همین روشی هست که گذاشتم اینجا 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 تا) داریم و تعداد معلوماتمون کمه. به کجا می خواهیم برسیم؟