طريقة نيوتن لحل المعادلات الغيرخطية

برقاوي نجد

عضو جديد
إنضم
8 سبتمبر 2007
المشاركات
16
مجموع الإعجابات
0
النقاط
0
البرنامج التالي يحل النظام(معادلتين) الغير الخطي بطريقة نيوتن، ولكن وجدت بعض الصعوبة في تعميمه الى ثلاث معادلات أو أكثر ، نبي مساعدتكم لحل المشكلة ............

البرنامج

function sol=newton3(fn5,jac,x0,tol)
old=x0+1;while max(abs(x0-old))>tol;old=x0;
f=feval(fn5,old);f1=f(1);f2=f(2);
J=feval(jac,old);f1x=J(1,1);f1y=J(1,2);f2x=J(2,1);f2y=J(2,2);
D=f1x*f2y-f1y*f2x;
h=(f2*f1y-f1*f2y)/D;
k=(f1*f2x-f2*f1x)/D;
x0=old+[h,k]
end;sol=x0;
 

نور سالم

عضو جديد
إنضم
1 ديسمبر 2007
المشاركات
28
مجموع الإعجابات
0
النقاط
0
يابرقاوي نجد برنامجك هذا سلمك الله
الله بالخير اولا كتابتك للجاكوبين بهذه الطريقة يجعل برنامجك غير جيد لانك ستضطر لكتابة جميع عناصر الجاكوبين بمعنى اذا كان لديك ثلاث معادلات يجب ان تكتب التسعة عناصر والمفروض ان تجعل برنامجك يحسب الجاكوبين لا ان يكون مدخل
ايضا هناك بعض الجمل تحتاج تعديل
استطيع اعدلك برنامجك ولكن دعني اسالك سؤال لماذا لا تستخدم امر fsolve
 

برقاوي نجد

عضو جديد
إنضم
8 سبتمبر 2007
المشاركات
16
مجموع الإعجابات
0
النقاط
0
أخوي نور سالم
البرنامج جربته وناجح في حل نظام غير خطي بمعادلتين ولا فيه مشكلة
لكن المشكلة مثل ماذكرت في حل ثلاث معادلات طبعاً أكيد بتصعب لأن بيكون لدي تسع عناصر وفيه صعوبة بتحديد D , h ,k
يمكن فيه طرق ثانية أسهل من هذا بكثير لكن ليس عندي خلفية عنها
اذا كانت لديك فإرجو نشرها في المنتدى لنستفيد منكم يعطيك 1000 عافيه أخوي
 

برقاوي نجد

عضو جديد
إنضم
8 سبتمبر 2007
المشاركات
16
مجموع الإعجابات
0
النقاط
0
أخوي نور سالم
البرنامج جربتة وناجح في حل نظام غيرخطي في معادلتين ولايوجد فيه مشكلة
ولكن المشكلة في حل ثلاث معادلات بيصبح لدي مثل ماذكرت تسع عناصر وبالتالي تصعب العملية في إيجاد D ,h , k
يمكن فيه برنامج أسهل وأفضل من هذا البرنامج ولكن ليس لدي معرفة به
أرجو إذا كان لديك أعلام به
يعطيك 1000عافيه وسلمك من كل مكروه
 
التعديل الأخير:

نور سالم

عضو جديد
إنضم
1 ديسمبر 2007
المشاركات
28
مجموع الإعجابات
0
النقاط
0
اخي الكريم
ارسلي معادلاتك على *****ي وابشر بالحل
وااكون لك مقدر اذا تنشرها هنا
لاني والله ليس لدي وقت

XXXXXXXXXXXXXXXX
 

برقاوي نجد

عضو جديد
إنضم
8 سبتمبر 2007
المشاركات
16
مجموع الإعجابات
0
النقاط
0
اولاً : أشكر نور المنتدى وملك الماتلاب الاخ/ نور سالم
على جهده ومساعدتة لي في هذا الموضوع
ثانياً : استئذن ألاخ نور لطرح الموضوع حل نظام غير خطي بشكل عام هنا

المطلوب اولاً فتح ملف جديد m.file ثم تكتب معادلاتك على سبيل المثال

(function F=EQU(x

F(1) = 3*x(1)-cos(x(2)*x(3))-0.5;
F(2) = x(1)^2-81*(x(2)+0.1)^2 + sin(x(3))+1.06;
F(3) = exp(-x(1)*x(2))+ 20*x(3)+ ((10*pi)-3)/3;

ثم حفظ الملف بالاسم الموجود لديك

ثم بعد ذلك تستخدم الأمر fsolve في نافذة الاوامر command بالطريقة التالية

% call for equations
%[x, feval]=fsolve('EQU',[1 1 1])
حيث اننا لابد ان نستخدم قيمة إبتدائة في حل النظام وضعنا هنا بالقيمة [1 1 1]
 

رائد غسان

عضو جديد
إنضم
9 فبراير 2006
المشاركات
111
مجموع الإعجابات
1
النقاط
0
اليك الحل

الملف ادناه بصيغة m-file يحل ن من المعادلات غير الخطية بطريقة نيوتن واذا اردت طرق اخرى فهي متوفرة لدي ان شاء الله .......
 

المرفقات

  • ALG101.zip
    1.9 KB · المشاهدات: 2,864

برقاوي نجد

عضو جديد
إنضم
8 سبتمبر 2007
المشاركات
16
مجموع الإعجابات
0
النقاط
0
السلام عليكم ، أسف على التاخير
أخ رائد شكرأ على الملف الموجود ... (في طور التجريب الان فيه صعوبة في التطبيق !)
بالنسبة للطرق الأخرى نعم يوجد طرق أخرى للغير الخطية مثل
(nowten Method ; Fixd point; Quasi-newton ..)
وكل طريقة لها عملها وسلبياتها وتقاربها
لو تضعها في المنتدى ليستفيد منها الجميع نكون من الشاكرين
أخيرأً شكراً لمشاركتك الفعالة معنا
 

tariqsamer

عضو جديد
إنضم
6 مارس 2005
المشاركات
274
مجموع الإعجابات
1
النقاط
0
مشكوووووووووووووووووووووووووووووووووووور
 

khaled 88

عضو جديد
إنضم
10 أغسطس 2010
المشاركات
8
مجموع الإعجابات
0
النقاط
0
السلام عليكم ورحمه الله وبركاته ..

انا عضو جديد ومآصدقت لقيت هـ الموضوع , ويالبرقاوي تكفى طلبتك تراني ولد عمك وضبطني !!
انا عندي بكره اخر موعد لتسليم هوم ورك طالبه الدكتور علي وهو حل اي مسأله على newton او newton modified ..
اي معادله تخطر ببالك بس اللهم انها تنحل عن طريق الماتلاب !! والمشكله هـ الماتلاب ماحصلته دورته لين داخ راسي بس مالقيتو !!

ظروري يالغآلي حتى لو مسأله محلوله وجاهزه عندك اتمنى ترسلي اياها , لأن اصلاً مو مطلوب اني ادرس الماتلاب بس مجرد فلسفه من الدكتور يبغا يعطي درجات ومو عارف كيف !!

اتمنى اذا كان عند احد جواب يضعه هنا

تكفى يالبرقاوووي او اي احد لاتخلوني وربي محتاجها ظروري ظروري ظروري ,.
 
التعديل الأخير:

khaled 88

عضو جديد
إنضم
10 أغسطس 2010
المشاركات
8
مجموع الإعجابات
0
النقاط
0
احيييييييييه آخر رد كان قبل سنتين ؟!!

برْب ؟!!!
 

برقاوي نجد

عضو جديد
إنضم
8 سبتمبر 2007
المشاركات
16
مجموع الإعجابات
0
النقاط
0
أأخوي خاااااااااالد موجود يااااااالغاااااالي
توووني دااااخل المنتدى وشفت الموووضوع صدفة

لي رجعة أن شاء الله أأنتظرررك
 
أعلى