محدودیت‌ها


 

20.1 مقدمه

20.1.1 گرادیان‌ها برای راه‌حل‌های بهینه با قیود

به ندرت «ناهار رایگان»ی وجود دارد. اگر بخواهیم کمیتی را بیشینه کنیم، اغلب مجبوریم با قیود کار کنیم. موانع ممکن است ما را از تغییر دلخواه پارامترها بازدارند. گرادیان همچنان می‌تواند به عنوان یک اصل راهنما استفاده شود. در حالی که نمی‌توانیم به صفر بودن f دست یابیم، می‌توانیم به دنبال نقاطی باشیم که گرادیان بر قید عمود است. این یک نقطه بهینه تحت محدودیت به ما می‌دهد. اگر در کوهستان در مسیری پیاده‌روی کنید، اغلب بدون اینکه در قله کوه باشید به یک حداکثر محلی می‌رسید. آنچه در چنین نقاطی x رخ می‌دهد این است که f ( x ) بر منحنی عمود است، به این معنی که f ( x ) با g ( x ) موازی است.

شکل 1. وضعیتی که در آن یک تابع f ( x , y ) در امتداد یک منحنی g ( x , y ) = c بهینه می‌شود، چارچوبی است که می‌توان با لاگرانژ به آن پرداخت. شرط حداکثر بودن به این معنی است که گرادیان f بر منحنی عمود است. این بدان معناست که گرادیان‌های f و g موازی هستند. f = λ g .

20.1.2 جادوی لاگرانژ: قیود متعدد، یک راه‌حل

روش لاگرانژ بسیار کلی‌تر است. می‌توانیم با تعداد دلخواه قید کار کنیم و همچنان از همان اصل استفاده کنیم. گرادیان f : n سپس بر سطح قید عمود است، به این معنی که ترکیبی خطی از گرادیان‌های همه m قید است: اینها n معادله f = j = 1 m λ j g j هستند زیرا بردارها دارای n مؤلفه هستند. همراه با m معادله g j = c j ، n + m معادله برای n + n متغیر x 1 , , x n ، λ 1 , , λ m داریم.

20.2 سخنرانی

20.2.1 یافتن حداکثر در فضاهای محدود

اگر بخواهیم یک تابع f : m را روی قید S = { x m g ( x ) = c } بیشینه کنیم، آنگاه هر دو گرادیان f و g اهمیت دارند. دو بردار v و w را موازی می‌نامیم اگر v = λ w یا w = λ v برای یک λ حقیقی. بردار صفر با همه چیز موازی است. در اینجا گونه‌ای از فرما آمده است:

قضیه 1. اگر x 0 یک حداکثر از f تحت قید g = c باشد، آنگاه f ( x 0 ) و g ( x 0 ) موازی هستند.

اثبات. با برهان خلف: فرض کنید f ( x 0 ) و g ( x 0 ) موازی نیستند و x 0 یک حداکثر محلی است. فرض کنید T صفحه مماس بر S = { g = c } در x 0 باشد. از آنجا که f ( x 0 ) بر T عمود نیست، می‌توانیم آن را روی T تصویر کنیم تا یک بردار غیرصفر v در T به دست آوریم که بر f عمود نیست. در واقع زاویه بین f و v تند است به طوری که cos ( α ) > 0 . یک منحنی r ( t ) در S با r ( 0 ) = x 0 و در نظر بگیرید. داریم با تقریب خطی، می‌دانیم که f ( r ( t ) ) > f ( r ( 0 ) ) برای t > 0 به اندازه کافی کوچک. این تناقضی است با این واقعیت که f در x 0 = r ( 0 ) روی S حداکثر بود. ◻

شکل 2. یک مسئله لاگرانژ

20.2.2 بررسی ضریب‌های لاگرانژ و شرایط لازم

این بلافاصله دلالت دارد: (بین g 0 و g = 0 تمایز قائل شوید)

قضیه 2. برای یک حداکثر از f روی S = { g = c } ، یا معادلات لاگرانژ f ( x 0 ) = λ g ( x 0 ) ، g = c برقرارند، یا g ( x 0 ) = 0 ، g = c .

برای توابع f ( x , y ) , g ( x , y ) از دو متغیر، این بدان معناست که باید یک دستگاه با سه معادله و سه مجهول را حل کنیم:

20.2.3 یافتن حداکثر واقعی

برای یافتن یک حداکثر، معادلات لاگرانژ را حل کنید و یک لیست از نقاط بحرانی g روی قید اضافه کنید. سپس نقطه‌ای را انتخاب کنید که f در بین همه نقاط حداکثر است. ما خود را با آزمون مشتق دوم مشغول نمی‌کنیم. اما در اینجا یک بیان ممکن است: d 2 d t 2 D t v D t v f ( x 0 ) | t = 0 < 0 برای همه v عمود بر g ( x 0 ) ، آنگاه x 0 یک حداکثر محلی است.

البته، حالت حداکثر و حداقل مشابه هستند. اگر f یک حداکثر روی g = c داشته باشد، آنگاه f یک حداقل در g = c دارد. می‌توانیم یک حداکثر از f تحت یک قید هموار S = { g = c } داشته باشیم بدون اینکه معادلات لاگرانژ برآورده شوند. یک مثال f ( x , y ) = x و g ( x , y ) = x 3 y 2 است که در شکل (20.3) نشان داده شده است.

شکل 3. مثالی از یک تابع که در آن معادلات لاگرانژ حداقل را نمی‌دهند، در اینجا ( 0 , 0 ) . این حالتی است که g = 0 .

20.2.4 صعود لاگرانژ: بیشینه‌سازی با قیود متعدد

روش لاگرانژ می‌تواند توابع f را تحت چندین قید بیشینه کند. بیایید این را در مورد یک تابع f ( x , y , z ) از سه متغیر و دو قید g ( x , y , z ) = c و h ( x , y , z ) = d نشان دهیم. مشابه اصل فرما این است که در یک حداکثر از f ، گرادیان f در صفحه ای است که توسط g و h گسترده شده است. این منجر به معادلات لاگرانژ برای 5 مجهول x , y , z , λ , μ می‌شود.

به عنوان مثال، اگر آنگاه نقاطی روی بیضی g = 1 ، h = 4 با حداقل یا حداکثر فاصله تا 0 پیدا می‌کنیم.

شکل 4. وضعیتی را می‌بینیم که در آن سعی می‌کنیم یک تابع f را تحت دو قید بیشینه کنیم. در این مورد، تقاطع g = c ، h = d یک بیضی است.

20.3 مثال‌ها

مثال 1. مسئله: f ( x , y ) = x 2 + 2 y 2 را تحت قید g ( x , y ) = x + y 2 = 1 کمینه کنید.
راه‌حل: معادلات لاگرانژ 2 x = λ ، 4 y = λ 2 y هستند. اگر y = 0 آنگاه x = 1 . اگر y 0 می‌توانیم معادله دوم را بر y تقسیم کنیم و 2 x = λ ، 4 = λ 2 به دست آوریم که دوباره x = 1 را نشان می‌دهد. نقطه x = 1 ، y = 0 تنها راه‌حل است.

مثال 2. مسئله: کدام قوطی نوشابه استوانه‌ای با ارتفاع h و شعاع r دارای حداقل سطح A برای حجم ثابت V است؟
راه‌حل: داریم V ( r , h ) = h π r 2 = 1 و A ( r , h ) = 2 π r h + 2 π r 2 . با x = h π ، y = r ، باید f ( x , y ) = 2 x y + 2 π y 2 را تحت قید g ( x , y ) = x y 2 = 1 بهینه کنید. این کار را در کلاس انجام خواهیم داد.

مثال 3. مسئله: اگر 0 p k 1 احتمال این باشد که یک تاس عدد k را نشان دهد، آنگاه داریم g ( p ) = p 1 + p 2 + + p 6 = 1 . این بردار p یک توزیع احتمال نامیده می‌شود. آنتروپی شانون p به صورت تعریف می‌شود.

توزیع p را پیدا کنید که آنتروپی S را بیشینه می‌کند.
راه‌حل: f = ( 1 log ( p 1 ) , , 1 log ( p n ) ) , g = ( 1 , , 1 ) . معادلات لاگرانژ عبارتند از 1 log ( p i ) = λ , p 1 + + p 6 = 1 , که از آن p i = e ( λ + 1 ) به دست می‌آید. معادله آخر 1 = i exp ( ( λ + 1 ) ) = 6 exp ( ( λ + 1 ) ) λ = log ( 1 / 6 ) 1 را ثابت می‌کند به طوری که p 1 = p 2 = = p 6 = 1 / 6 . این تاس عادلانه است که دارای حداکثر آنتروپی است. حداکثر آنتروپی به معنای کمترین محتوای اطلاعاتی است.

مثال 4. فرض کنید احتمال اینکه یک سیستم فیزیکی یا شیمیایی در حالت k باشد p k است و انرژی حالت k E k است. طبیعت انرژی آزاد F ( p 1 , , p n ) = i [ p i log ( p i ) E i p i ] را کمینه می‌کند اگر انرژی‌های E i ثابت باشند. توزیع احتمال p i که i p i = 1 را برآورده می‌کند و انرژی آزاد را کمینه می‌کند، یک توزیع گیبس نامیده می‌شود. این توزیع را به طور کلی اگر E i داده شده باشند، پیدا کنید.
راه‌حل: f = ( 1 log ( p 1 ) E 1 , , 1 log ( p n ) E n ) , g = ( 1 , , 1 ) . معادله لاگرانژ log ( p i ) = 1 λ E i ، یا p i = exp ( E i ) C است، که در آن C = exp ( 1 λ ) . قید p 1 + + p n = 1 C ( i exp ( E i ) ) = 1 را می‌دهد به طوری که C = 1 / ( i e E i ) . راه‌حل گیبس p k = exp ( E k ) / i exp ( E i ) است.1

مثال 5. اگر f یک تابع درجه دوم روی m باشد و g خطی باشد، یعنی f ( x ) = B x x / 2 با B M ( m , m ) و اگر قید g ( x ) = A x = c خطی باشد A M ( 1 , m ) ، آنگاه f ( x ) = B x و g ( x ) = A T . بیایید b = A T M ( m , 1 ) m بنامیم. معادلات لاگرانژ سپس B x = λ b ، A x = c هستند. به طور کلی می‌بینیم که برای f درجه دوم و g خطی، به یک دستگاه معادلات خطی می‌رسیم.

مثال 6. مرتبط با نکته قبلی، مشاهده زیر است. اغلب می‌توان مسئله لاگرانژ را به یک مسئله بدون قید تقلیل داد. این دیدگاهی است که اغلب در اقتصاد اتخاذ می‌شود. بیایید به آن در بعد 2 نگاه کنیم، جایی که f ( x , y ) را تحت قید g ( x , y ) = 0 حدی می‌یابیم. F ( x , y , λ ) = f ( x , y ) λ g ( x , y ) را تعریف کنید. معادلات لاگرانژ برای f , g اکنون معادل F ( x , y , λ ) = 0 در سه بعد هستند.

تمرین‌ها

تمرین ۱. سبد استوانهای را بیابید که از بالا باز است و بیشترین حجم را برای مساحت ثابت π دارد. اگر x شعاع و y ارتفاع باشد، باید f ( x , y ) = π x 2 y را تحت قید g ( x , y ) = 2 π x y + π x 2 = π بیشینه کنیم. از روش ضریب‌های لاگرانژ استفاده کنید.

تمرین ۲. با فرض یک ماتریس متقارن n × n به نام B ، تابع f ( x ) = x B x را در نظر می‌گیریم و به دنبال نقاط اکسترمم f تحت قید g ( x ) = x x = 1 می‌گردیم. این به معادله B x = λ x منجر می‌شود. یک جواب x را بردار ویژه می‌نامند. ثابت لاگرانژ λ یک مقدار ویژه است. جواب‌های B x = λ x ، | x | = 1 را بیابید اگر B یک ماتریس 2 × 2 باشد، که در آن f ( x , y ) = a x 2 + ( b + c ) x y + d y 2 , g ( x , y ) = x 2 + y 2 . سپس مسئله را با a = 4 ، b = 1 ، c = 1 ، d = 4 حل کنید.

تمرین ۳. کدام هرم به ارتفاع h بر روی یک مربع [ a , a ] × [ a , a ] با مساحت سطح 4 a h 2 + a 2 + 4 a 2 = 4 دارای حداکثر حجم V ( h , a ) = 4 h a 2 / 3 است؟ با استفاده از متغیرهای جدید ( x , y ) و ضرب V در یک ثابت، به مسئله معادل بیشینه‌سازی f ( x , y ) = y x 2 تحت قید g ( x , y ) = x y 2 + x 2 + x 2 = 1 می‌رسیم. از متغیرهای اخیر استفاده کنید.

تمرین ۴. با الهام از فیلم دیزنی "Tangled"، می‌خواهیم یک بالون هوای گرم با یک شبکه مکعب‌مستطیلی به ابعاد x ، y ، z بسازیم که همراه با استحکامات بالا و پایین از سیم‌هایی با طول کل g ( x , y , z ) = 6 x + 6 y + 4 z = 32 استفاده می‌کند. بالونی با حداکثر حجم f ( x , y , z ) = x y z را بیابید.

تمرین ۵. یک گلوله جامد ساخته شده از یک نیم‌کره و یک استوانه دارای حجم V = 2 π r 3 / 3 + π r 2 h و مساحت سطح A = 2 π r 2 + 2 π r h + π r 2 است. دکتر منهتَن گلوله‌ای با حجم ثابت و حداقل مساحت طراحی می‌کند. با g = 3 V / π = 1 و f = A / π ، او بنابراین f ( h , r ) = 3 r 2 + 2 r h را تحت قید g ( h , r ) = 2 r 3 + 3 r 2 h = 1 کمینه می‌کند. از روش لاگرانژ برای یافتن یک کمینه موضعی f تحت قید g = 1 استفاده کنید.

پیوست: تصویرسازی داده‌ها: کاب-داگلاس

۲۰.۳.۱ کاب-داگلاس: فرمولی برای رشد اقتصادی

ریاضیدان و اقتصاددان چارلز دبلیو. کاب در کالج آمهرست و اقتصاددان و سیاستمدار پل اچ. داگلاس که در آمهرست نیز تدریس می‌کرد، در سال ۱۹۲۸ به طور تجربی فرمول F ( K , L ) = L α K β را یافتند که با تولید کل F یک سیستم اقتصادی به عنوان تابعی از سرمایه‌گذاری K و نیروی کار L مطابقت دارد. دو نویسنده از متغیرهای لگاریتمی استفاده کردند و خطی بودن را برای یافتن α , β فرض کردند. در زیر داده‌ها به گونه‌ای نرمال‌سازی شده‌اند که مقدار سال ۱۸۹۹ برابر 100 است.

سال K L P
۱۸۹۹۱۰۰۱۰۰۱۰۰
۱۹۰۰۱۰۷۱۰۵۱۰۱
۱۹۰۱۱۱۴۱۱۰۱۱۲
۱۹۰۲۱۲۲۱۱۸۱۲۲
۱۹۰۳۱۳۱۱۲۳۱۲۴
۱۹۰۴۱۳۸۱۱۶۱۲۲
۱۹۰۵۱۴۹۱۲۵۱۴۳
۱۹۰۶۱۶۳۱۳۳۱۵۲
۱۹۰۷۱۷۶۱۳۸۱۵۱
۱۹۰۸۱۸۵۱۲۱۱۲۶
۱۹۰۹۱۹۸۱۴۰۱۵۵
۱۹۱۰۲۰۸۱۴۴۱۵۹
۱۹۱۱۲۱۶۱۴۵۱۵۳
۱۹۱۲۲۲۶۱۵۲۱۷۷
۱۹۱۳۲۳۶۱۵۴۱۸۴
۱۹۱۴۲۴۴۱۴۹۱۶۹
۱۹۱۵۲۶۶۱۵۴۱۸۹
۱۹۱۶۲۹۸۱۸۲۲۲۵
۱۹۱۷۳۳۵۱۹۶۲۲۷
۱۹۱۸۳۶۶۲۰۰۲۲۳
۱۹۱۹۳۸۷۱۹۳۲۱۸
۱۹۲۰۴۰۷۱۹۳۲۳۱
۱۹۲۱۴۱۷۱۴۷۱۷۹
۱۹۲۲۴۳۱۱۶۱۲۴۰
شکل ۵. نمودار F ( L , K ) = L 3 / 4 K 1 / 4 به خوبی با آن مجموعه داده مطابقت دارد. می‌توانید در داده‌ها یک نقطه پرت را مشاهده کنید.

۲۰.۳.۲ تجسم محدودیت‌های تولید

فرض کنید که نیروی کار و سرمایه‌گذاری با قید اضافی G ( L , K ) = L 3 / 4 + K 1 / 4 = 50 محدود شده‌اند. (این تابع G با تابع F ( L , K ) مرتبط نیست زیرا ما در یک مسئله لاگرانژ هستیم.) تولید P تحت این قید در کجا بیشینه است؟ دو تابع F ( L , K ) و G ( L , K ) را رسم کنید.


  1. این مثال از روفوس بوون، یادداشت‌های درسی در ریاضیات، ۴۷۰، ۱۹۷۸ است.↩︎