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

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

البرنامج

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;

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

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

اخي الكريم
ارسلي معادلاتك على *****ي وابشر بالحل
وااكون لك مقدر اذا تنشرها هنا
لاني والله ليس لدي وقت

XXXXXXXXXXXXXXXX

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

المطلوب اولاً فتح ملف جديد 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]

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

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

mshkooooooooooor

مشكوووووووووووووووووووووووووووووووووووور