একটি পৃথক ইউনিফর্ম বিতরণ প্রতিস্থাপন ছাড়াই টানা নমুনার মধ্যে সর্বাধিক ফাঁক


16

এই সমস্যাটি রোবোটিক কভারেজে আমার ল্যাবের গবেষণার সাথে সম্পর্কিত:

সেটটি থেকে এলোমেলোভাবে সংখ্যাগুলি আঁকুন replacement প্রতিস্থাপন ছাড়াই, এবং সংখ্যাটিকে ক্রমবর্ধমান ক্রমে সাজান। ।n n{ 1 , 2 , , মি } {1,2,,m}1 n মি1nm

numbers সংখ্যার এই সাজানো তালিকা থেকে একটানা সংখ্যা এবং সীমানার মধ্যে পার্থক্য তৈরি করুন: । এটি ফাঁক দেয় ।{ একটি ( 1 ) , একটি ( 2 ) , ... , একটি ( এন ) }{a(1),a(2),,a(n)}g = { a ( 1 ) , a ( 2 ) - a ( 1 ) , , a ( n ) - a ( n - 1 ) , m + 1 - a ( n ) } g={a(1),a(2)a(1),,a(n)a(n1),m+1a(n)}n + 1n+1

সর্বোচ্চ ব্যবধানের বিতরণ কী?

পি ( সর্বোচ্চ ( ) = কে ) = পি ( কে ; এম , এন ) = ?P(max(g)=k)=P(k;m,n)=?

এটি অর্ডার পরিসংখ্যান ব্যবহার করে ফ্রেম করা যেতে পারে : পি ( জি ( এন + 1 ) = কে ) = পি ( কে ; মি , এন ) = ?P(g(n+1)=k)=P(k;m,n)=?

ফাঁকগুলির বিতরণের জন্য লিঙ্কটি দেখুন , তবে এই প্রশ্নটি সর্বোচ্চ ফাঁকগুলির বিতরণ জিজ্ঞাসা করে ।

আমি গড় মান, \ mathbb {E} [g _ {(n + 1)}] দিয়ে সন্তুষ্ট হয়েছি E [ g ( n + 1 ) ]E[g(n+1)]

যদি n = মিn=m সমস্ত ফাঁক হয় তবে আকার 1 n যদি এন + 1 = মি 2n + 1 = মিn+1=m আকারের একটি ফাঁক হয় , এবং n + 1 সম্ভাব্য অবস্থানগুলি। সর্বাধিক গ্যাপের আকার এম-এন + 1 , এবং এই ফাঁকটি মোট এন + 1 সম্ভাব্য পজিশনের জন্য, যে কোনও এন সংখ্যার আগে বা পরে স্থাপন করা যেতে পারে । সর্বাধিকতম সর্বাধিক গ্যাপের আকার \ lceil \ frac {mn} {n + 1} ce rceil । কোনো সমন্বয় সম্ভাবনা নির্ধারণ টি = {M N চয়ন \} ^ {- 1}22n + 1n+1মি - এন + 1mn+1এনnn + 1n+1এম - এনএন + + 1mnn+1টি = ( মিn ) -1T=(mn)1

আমি P (g _ {(n + 1)} = কে) = পি (কে; এম, এন) = \ আরম্ভ {কেস} 0 & কে <ce লিসিল \ ফ্র্যাক {এমএন} as হিসাবে সম্ভাব্য ভর কার্যকারিতা আংশিক সমাধান করেছি পি ( জি ( এন + 1 ) = কে ) = পি ( কে ; এম , এন ) = { 0 কে < এম - এনn + 11কে=মি-এনএন + + 1 1=1 (ঘটে যখন মি=)টি(এন+ +1)=2 (ঘটে যখন মি=+ +1)টি(এন+ +1)k=মি-(এন-1)এন ? মি-(এন-1)nkmn+1T(n+1)k=mn+10k>mn+1P(g(n+1)=k)=P(k;m,n)=011T(n+1)T(n+1)?T(n+1)0k<mnn+1k=mnn+1k=1 (occurs when m=n)k=2 (occurs when m=n+1)k=m(n1)nm(n1)nkmn+1k=mn+1k>mn+1(1)

বর্তমান কাজ (1): প্রথম ফাঁক, একটি _ {(1)} এর সমীকরণটি a(1)a(1)সোজা: P(a(1)=k)=P(k;m,n)=1(mn)mn+1k=1(mk1n1)

P(a(1)=k)=P(k;m,n)=1(mn)k=1mn+1(mk1n1)
প্রত্যাশিত মানটির একটি সাধারণ মান রয়েছে: E[P(a(1))]=1(mn)mn+1k=1(mk1n1)k=mn1+nE[P(a(1))]=1(mn)mn+1k=1(mk1n1)k=mn1+n । প্রতিসাম্য দ্বারা, আমি আশা করি যে সমস্ত nn ফাঁকগুলি এই বন্টন হবে। সম্ভবত এই বিতরণ থেকে nn বার অঙ্কন করে সমাধানটি সন্ধান করা যেতে পারে ।

বর্তমান কাজ (2): মন্টি কার্লো সিমুলেশনগুলি চালানো সহজ।

simMaxGap[m_, n_] := Max[Differences[Sort[Join[RandomSample[Range[m], n], {0, m+1}]]]];
m = 1000; n = 1; trials = 100000;
SmoothHistogram[Table[simMaxGap[m, n], {trials}], Filling -> Axis,
Frame -> {True, True, False, False},
FrameLabel -> {"k (Max gap)", "Probability"},
PlotLabel -> StringForm["m=``,n=``,smooth histogram of maximum map for `` trials", m, n, trials]][![enter image description here][1]][1]

1
এই শর্তগুলির সাথে আপনার অবশ্যই এন <= মি। আমার মনে হয় আপনি g = {a_ (1), a_ (2) -a_ (1), ..., a_ (n) -a_ (n-1) want চান} এলোমেলোভাবে নির্বাচিত হওয়ার অর্থ কি প্রথম অঙ্কে সম্ভাব্যতা 1 / মি সহ প্রতিটি সংখ্যা নির্বাচন করা? যেহেতু আপনি প্রতিস্থাপন করবেন না সম্ভাবনাটি দ্বিতীয়টিতে 1 / (মি -1) এবং এমথ অঙ্কনে 1-এ নেমে যদি n = m থাকে। যদি এন <এম এটি প্রথম ড্র হয় তবে সম্ভাব্যতা 1 / (এম- (এন -1)) এর সাথে থাকবে নবম ড্রতে stop
মাইকেল আর চেরনিক

2
আপনার মূল বিবরণটি কোনও অর্থহীন নয়, কারণ (আমি বিশ্বাস করি) আপনি দুটি সাবস্ক্রিপ্ট স্থানান্তরিত করেছেন। দয়া করে যাচাই করুন যে আমার সম্পাদনাটি আপনার অভিপ্রায় অনুসারে রয়েছে: বিশেষত, দয়া করে নিশ্চিত করুন যে আপনি সেখানে ফাঁক থাকতে চান, যার মধ্যে । প্রথম। ggnna(1)a(1)
হুড়হুড়ি

1
@ গুং আমি মনে করি এটি আত্ম-অধ্যয়নের চেয়ে গবেষণাটি
গ্লেন_বি -রিনস্টেট মনিকা 18

1
আমি মনে করি আপনার সর্বনিম্ন এবং সর্বাধিক গ্যাপের মাপ এবং হওয়া উচিত । একটানা পূর্ণসংখ্যা নির্বাচন করা হলে সর্বনিম্ন ব্যবধানের আকার হ'ল এবং আপনি যখন এবং প্রথম পূর্ণসংখ্যা (বা এবং ) নির্বাচন করেন তখন সর্বাধিক ফাঁক আকার হয়11mn+1mn+1mmn1n11,,n11,,n111mn+2,,mmn+2,,m
সম্ভাব্যতাব্লোগিক

1
ধন্যবাদ মাইকেল চেরনিক এবং সম্ভাব্যতা সম্পর্কিত, আপনার সংশোধন করা হয়েছে। সংশোধন করার জন্য আপনাকে @ শুভ ধন্যবাদ!
অ্যারোনবেকার

উত্তর:


9

যাক সুযোগ যে ন্যূনতম হতে , সমান ; যে, নমুনা নিয়ে গঠিত এবং এর -subset । সেখানে মতো সাবসেটগুলি সমানভাবে সম্ভাব্য বাইরে রয়েছেf(g;n,m)f(g;n,m)a(1)a(1)ggggn1n1{g+1,g+2,,m}{g+1,g+2,,m}(mgn1)(mgn1)(mn)(mn)

Pr(a(1)=g=f(g;n,m)=(mgn1)(mn).

Pr(a(1)=g=f(g;n,m)=(mgn1)(mn).

যোগ করার পদ্ধতি সব সম্ভাব্য মান জন্য তার চেয়ে অনেক বেশী বেঁচে থাকার ফাংশন উৎপাদf(k;n,m)f(k;n,m)kkgg

Pr(a(1)>g)=Q(g;n,m)=(mg)(mg1n1)n(mn).

Pr(a(1)>g)=Q(g;n,m)=(mg)(mg1n1)n(mn).

যাক হতে দৈব চলক বৃহত্তম ফাঁক কর্তৃক প্রদত্ত:Gn,mGn,m

Gn,m=max(a(1),a(2)a(1),,a(n)a(n1)).

Gn,m=max(a(1),a(2)a(1),,a(n)a(n1)).

(এটি এবং মধ্যে ফাঁক অন্তর্ভুক্ত করার আগে এটি সংশোধন করার আগে এই প্রশ্নের উত্তর দেয় )a(n)a(n)mm আমরা এর বেঁচে থাকার ক্রিয়া গণনা করব যা থেকে of এর সম্পূর্ণ বিতরণ সহজেই উদ্ভূত হয়। পদ্ধতিটি দিয়ে শুরু হওয়া একটি গতিশীল প্রোগ্রাম , যার জন্য এটি স্পষ্টP(g;n,m)=Pr(Gn,m>g),

P(g;n,m)=Pr(Gn,m>g),
Gn,mGn,mn=1n=1

P(g;1,m)=Pr(G1,m>1)=mgm, g=0,1,,m.

P(g;1,m)=Pr(G1,m>1)=mgm, g=0,1,,m.(1)

বৃহত্তর , নোট করুন যে ইভেন্ট হল ইভেন্টটির ইউনিয়নn>1n>1Gn,m>gGn,m>g

a1>g,

a1>g,

যার জন্য প্রথম ব্যবধানটি ছাড়িয়ে গেছে এবং পৃথক ইভেন্টgggg

a1=k and Gn1,mk>g, k=1,2,,g

a1=k and Gn1,mk>g, k=1,2,,g

যার জন্য প্রথম ব্যবধানটি সমান এবং চেয়ে বড় ব্যবধানটি নমুনায় পরে আসে। মোট সম্ভাবনার আইন এই ইভেন্টগুলির সম্ভাব্যতাগুলি যুক্ত করে, কোথাও যুক্ত করেkkgg

P(g;n,m)=Q(g;n,m)+gk=1f(k;n,m)P(g;n1,mk).

P(g;n,m)=Q(g;n,m)+k=1gf(k;n,m)P(g;n1,mk).(2)

ফিক্সিং এবং দ্বি-মুখী অ্যারে রাখুন এবং আমরা ব্যবহার করে গণনা করতে পারি এর প্রথম সারিটি পূরণ করতে এবং প্রতি সারি ক্রিয়াকলাপ ব্যবহার করে প্রতিটি ধারাবাহিক সারি পূরণ করতে । ফলস্বরূপ টেবিলটি অপারেশনে সমাপ্ত হতে পারে এবং মাধ্যমে জন্য সমস্ত টেবিলগুলি যেতে পারে ।ggi=1,2,,ni=1,2,,nj=1,2,,mj=1,2,,mP(g;n,m)P(g;n,m)(1)(1)(2)(2)O(gm)O(gm)O(gmn)O(gmn)g=1g=1g=mn+1g=mn+1O(m3n)

ব্যক্তিত্ব

এই গ্রাফগুলি জন্য থাকার ফাংশনটি । হিসাবে বেড়ে যায়, বাম দিকে গ্রাফ প্যাচসমূহ, বড় ফাঁক এর কমে সম্ভাবনা সংশ্লিষ্ট।gP(g;n,64)n=1,2,4,8,16,32,64n

বন্ধ সূত্রগুলি বিশেষত বড় জন্য অনেকগুলি বিশেষ ক্ষেত্রে পাওয়া যেতে পারে , তবে আমি সমস্ত প্রয়োগ করে এমন একটি বদ্ধ সূত্র পেতে সক্ষম হইনি । একটানা ইউনিফর্ম ভেরিয়েবলের জন্য অভিন্ন সমস্যাটির সাথে এই সমস্যাটি প্রতিস্থাপন করে ভাল আনুমানিকতা সহজেই উপলব্ধ।P(g;n,m)ng,n,m

শেষ , of এর প্রত্যাশা থেকে শুরু হয়ে বেঁচে থাকার ক্রিয়াটি প্রাপ্ত হবে :Gn,mg=0

E(Gn,m)=mn+1g=0P(g;n,m).

চিত্র 2: প্রত্যাশার কনট্যুর প্লট

প্রত্যাশার এই কনট্যুর প্লটটি অন্ধকার থেকে আলোর দিকে স্নাতক হওয়া এ সংক্ষিপ্তসার দেখায় ।2,4,6,,32


প্রস্তাবনা: লাইন 'আসুন দৈব চলক বৃহত্তম ফাঁক কর্তৃক প্রদত্ত হতে: ", শেষ ফাঁক যোগ করুন । আপনার প্রত্যাশার প্লটটি আমার মন্টি কার্লো সিমুলেশনের সাথে মেলে। Gn,mm+1an
অ্যারোনবেকার
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.