সমীকরণের রৈখিক সিস্টেমগুলি সমাধানের জন্য ক্রিলোভ সাবস্পেস পদ্ধতির রূপান্তর করার পেছনের নীতিটি কী?


24

আমি যেমন এটি বুঝতে পারি, সমীকরণের রৈখিক সিস্টেমগুলি সমাধান করার জন্য দুটি প্রধান বিভাগের পুনরাবৃত্তি পদ্ধতি রয়েছে:

  1. নিশ্চল পদ্ধতি (জ্যাকোবি, গাউস-সিডেল, এসওআর, মাল্টিগ্রিড)
  2. ক্রিলোভ সাবস্পেস পদ্ধতি (কনজুগেট গ্রেডিয়েন্ট, জিএমআরইএস ইত্যাদি)

আমি বুঝতে পারি যে বেশিরভাগ স্থির পদ্ধতিগুলি ত্রুটিটির ফুরিয়ার মোডগুলি পুনরাবৃত্তভাবে শিথিল করে (স্মুথ করা) কাজ করে। আমি যেমন এটি বুঝতে পেরেছি, কনজুগেট গ্রেডিয়েন্ট পদ্ধতি (ক্রিলোভ সাবস্পেস পদ্ধতি) ম অবশিষ্টাংশের ক্ষেত্রে প্রয়োগিত ম্যাট্রিক্সের শক্তিগুলি থেকে অনুসন্ধানের দিকনির্দেশের সর্বোত্তম সেটের মাধ্যমে "পদক্ষেপ" দ্বারা কাজ করে । এই নীতিটি কি সমস্ত ক্রিলোভ সাবসপেস পদ্ধতিতে সাধারণ? যদি তা না হয় তবে আমরা কীভাবে সাধারণভাবে ক্রিলোভ সাবস্পেস পদ্ধতির রূপান্তর করার মূলনীতিটি চিহ্নিত করব?এন


2
নিশ্চল পদ্ধতিগুলির আপনার বিশ্লেষণ সাধারণ মডেল সমস্যাগুলির দ্বারা পক্ষপাতদুষ্ট, কারণ এগুলি ফুরিয়ার মোডের ক্ষেত্রে বিশ্লেষণ করা যেতে পারে। এটি বিকল্প দিকটি অন্তর্নিহিত (এডিআই) এবং অন্যান্য অনেক পদ্ধতি উপেক্ষা করে। বেশিরভাগ "স্টেশনারি মেথডস" এর মূল বিষয়টি হ'ল বহু সরল "আনুমানিক আংশিক" সমাধানকারীকে এক পুনরাবৃত্তীয় দ্রাবকের সাথে একত্রিত করা । ক্রিলোভ পদ্ধতির বিন্দুটি প্রদত্ত স্টেশনিয় রৈখিক পুনরাবৃত্তির রূপান্তরকে ত্বরান্বিত করা (বা এমনকি প্রয়োগ করা)।
টমাস ক্লিম্পেল

4
আমি মনে করি যে একটি কাগজ আপনার প্রশ্নের উত্তর দেওয়ার জন্য লেখা হয়েছিল তা হলেন ইপসেন এবং মায়ার, ক্রেলোভ পদ্ধতির পিছনে ধারণা, আমের। ম্যাথ। মাসিক 105 (1998) পিপি 889-899। এটি একটি আশ্চর্যজনকভাবে ভাল লিখিত এবং স্পষ্ট করার কাগজ, এখানে উপলব্ধ ।
অ্যান্ড্রু টি বার্কার

@ অ্যান্ড্রুটি.বার্কার: দুর্দান্ত! ধন্যবাদ অ্যান্ড্রু! :)
পল

উত্তর:


21

সাধারণভাবে, সমস্ত ক্রিলোভ পদ্ধতিগুলি ম্যাট্রিক্সের বর্ণালীতে মূল্যায়ন করার সময় মূলত একটি বহুবর্ষের সন্ধান করে। বিশেষ করে, ম একটি Krylov পদ্ধতি (শূন্য প্রাথমিক অনুমান সহ) এর অবশিষ্ট আকারে লেখা যেতে পারেএন

Rএন=পিএন(একজন)

যেখানে হল ডিগ্রি এর কিছু বহুপদী । এনপিএনএন

যদি , আমাদের= ভি Λ ভি - 1একজনএকজন=ভীΛভী-1

rnVPn(Λ)V1b=κ(V)Pn(Λ)b.

যে পরিস্থিতি এ স্বাভাবিক (যেমন, প্রতিসম বা একক) আমরা জানি যে জিআরইআরএস আর্নল্ডি পুনরাবৃত্তির মাধ্যমে এই জাতীয় বহুবিন্দু তৈরি করে, অন্যদিকে সিজি একটি পৃথক অভ্যন্তরীণ পণ্য ব্যবহার করে বহুপদী গঠন করে ( বিশদগুলির জন্য এই উত্তরটি দেখুন) )। একইভাবে, বিসিজি সংক্ষিপ্ত আকারের ল্যাঙ্কজোস প্রক্রিয়াটির মাধ্যমে তার বহুপদী গঠন করে, তবে চেবিশেভ পুনরাবৃত্তি বর্ণালী সম্পর্কে পূর্বের তথ্য ব্যবহার করে (সাধারণত প্রতিসাম্যিক নির্দিষ্ট ম্যাট্রিকগুলির জন্য বৃহত্তম এবং ক্ষুদ্রতম আইজেনিয়ুলের অনুমান)।κ ( ভি ) = 1।Aκ(V)=1.

একটি দুর্দান্ত উদাহরণ হিসাবে (ট্রেফেন + বাউ দ্বারা অনুপ্রাণিত), এমন একটি ম্যাট্রিক্স বিবেচনা করুন যার বর্ণালী এটি:

ম্যাট্রিক্সের স্পেকট্রাম

ম্যাটল্যাবে, আমি এটি দিয়ে এটি নির্মাণ করেছি:

A = rand(200,200);
[Q R] = qr(A);
A = (1/2)*Q + eye(200,200);

আমরা যদি জিএমআরইএস বিবেচনা করি, যেগুলি বহুবচন তৈরি করে যা প্রকৃতপক্ষে ডিগ্রি এর সমস্ত মনিক বহুবর্ষের চেয়ে অবশিষ্টাংশকে ন্যূনতম করে তোলে, আমরা প্রার্থী বহুবর্ষটি দেখে সহজেই অবশিষ্ট ইতিহাসের পূর্বাভাস দিতে পারিn

Pn(z)=(1z)n

আমাদের ক্ষেত্রে যা দেয়

|Pn(z)|=12n

জন্য বর্ণালী মধ্যে ।zA

এখন, যদি আমরা GMRES এলোমেলো আরএইচএসে চালিত করি এবং অবশিষ্ট ইতিহাসকে এই বহুবর্ষের সাথে তুলনা করি তবে তাদের অবশ্যই বেশিরভাগ মিল হওয়া উচিত (প্রার্থীর বহুপদী মানগুলি GMRES অবশিষ্টের চেয়ে ছোট কারণ because ):b2>1

অবশিষ্ট ইতিহাস


"ম্যাট্রিক্সের বর্ণালীতে ছোট" বলতে কী বোঝাতে চেয়েছেন আপনি কি তা পরিষ্কার করতে পারেন?
পল

2
একটি জটিল বহুপদী হিসেবে নিয়েছে, বহুপদী জটিল সমতল যার বর্ণালী অন্তর্ভুক্ত একটি অঞ্চলের ছোট মডুলাস হয়েছে । ইগুভ্যালুগুলির একটি বিক্ষিপ্ত প্লটের উপর সুপারম্পোজ করা একটি কনট্যুর প্লট কল্পনা করুন। ছোট কত ছোট? এটি সমস্যার উপর নির্ভর করে, সাধারণ কিনা এবং ডান-হাতযদিও মূল ধারণাটি হ'ল বহুবর্ষের ক্রম বর্ণালীতে ক্রমান্বয়ে আরও ছোট এবং ছোট হওয়ার চেষ্টা করে যাতে আমার অবশিষ্ট অনুমান । একজন একজন ( পি এন ) 0PnAAb.(Pn)0
রিড.এচচসন

@ রেড.এচচসন: খুব ভাল করে বলেছেন। আমি recommend লেখার প্রস্তাব দিতে পারি যেমন এবং উল্লেখ যে এটা স্বাভাবিক ম্যাট্রিক্স জন্য কেউ নেই? κ ( ভি )VV1κ(V)
জ্যাক পলসন

অনুকূল এসওআর দ্বারা পূর্বশর্তযুক্ত ল্যাপ্লেসিয়ার উদাহরণটি ম্যাট্রিক্সের সাথে খুব মিল রয়েছে। বিশদগুলি
জেড ব্রাউন

কড়া কথায় বলতে গেলে সিজিএনই বর্ণালী থেকে স্বতন্ত্র কারণ এটি কেবলমাত্র একক মানের উপর নির্ভর করে।
জেড ব্রাউন

17

নিয়ম অনুসারে

রেডের যোগসূত্র হিসাবে। অ্যাচচসনের উত্তরে আমি নিয়ম সম্পর্কিত কিছু বিষয় পরিষ্কার করতে চাই। এ পুনরাবৃত্তির, GMRES বহুপদী খুঁজে বের করে যে ছোট অবশিষ্ট এর -norm P n 2nthPn2

rn=Axnb=(Pn(A)1)bb=Pn(A)b.

ধরুন তাই হয়, এসপিডি হয় একটি আদর্শ সংঘটিত এবং তাই করে । তারপর- AAA1

rnA1=rnTA1rn=(Aeএন)টিএকজন-1একজনএন=এনটিএকজনএন=এনএকজন

যেখানে আমরা ত্রুটিটি ব্যবহার করেছি

এন=এক্সএন-এক্স*=এক্সএন-একজন-1=একজন-1Rএন

সুতরাং ত্রুটির -norm অবশিষ্টের আদর্শের সমান । কনজ্যুগেট গ্রেডিয়েন্টগুলি ত্রুটির -nnor হ্রাস করে যা এটি কম শক্তি মোডগুলি সমাধান করার ক্ষেত্রে তুলনামূলকভাবে আরও নির্ভুল করে তোলে। অবশিষ্ট এর -norm, যা ছোট GMRES, মত হল ত্রুটির -norm, এবং এইভাবে অর্থে দুর্বল যে কম শক্তি মোড কম ভাল-সমাধান করা হয়। নোট করুন যে অবশিষ্টাংশের ইউনির্ম মূলত অকেজো কারণ এটি স্বল্প-শক্তি মোডে এমনকি দুর্বল।একজনএকজন-1একজন2একজনটিএকজনএকজন

অভিমুখে সীমাবদ্ধতার তীক্ষ্ণতা

অবশেষে, জিআরইআরএস রূপান্তরটির বিভিন্ন ক্রাইলোভ পদ্ধতি এবং সূক্ষ্মতা সম্পর্কিত আকর্ষণীয় সাহিত্য রয়েছে, বিশেষত অ-সাধারণ অপারেটরদের জন্য।



11

সংক্ষেপে Iterative পদ্ধতি:

  1. একজনএক্স=সি

    এক্স=এক্স+ +সি-সিএকজনএক্স
    আমি-সিএকজন<1সিসি=ডি-1ডিএকজন
  2. ইউ,ভীসিএনএক্স~ইউ-একজনএক্স~ভীইউএকজনভীভী=ইউভী=একজনইউ

    ভীএক্স~ইউইউ

    ইউএকজনএক্স~

    এটি ইউটিফ সাদের বইটিতে পুনরুক্তি পদ্ধতিতে সুন্দরভাবে ব্যাখ্যা করা হয়েছে ।

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.