Printable View
من باید برای درس محاسبات عددی یک برنامه با matlab برای حل عددی معادلات دیفرانسیل مرتبه اول به 4 روش: اولر، تیلور،رونگه کوتای مرتبه 2 و 4 بنویسم. به این صورت که برنامه تابع و نوع حل و شرط اولیه و طول تقسیمات و نقطه ی دلخواه A رو از کاربر دریافت کنه و جواب رو در نقطه ی A محاسبه کنه.
من با C کار کردم و متاسفانه آشنایی کافی با matlab ندارم.
لطفا بگید چطوری باید از کاربر تابع دو متغیره تحویل بگیرم و چطوری باید مشتقات جزیی رو برای اون تابع به دست بیارم.(لطفا یه جوری باشه که من بفهمم!!)
ممنون
این نمونه کد رو نگاه بنداز ، اول یک M فایل ایجاد کن و این برنامه رو توش ذخیره کن :نقل قول:
و بعد برای اجرا این کد رو بنویس :کد:% ورودي : ابتدا و انتهاي بازه و آرايه يک بعدي حامل معادله
% خروجي : ريشه معادله در بازه ورودي در صورتي که حاصلضرب مقدار تابع به ازاي
% مقدار ابتدا و انتهاي بازه منفي شود به عبارتي لااقل يک ريشه در بازه ورودي
% وجود داشته باشد
% عملکرد : پيدا کردن ريشه معادله به روش نابجايي
function NabejaeeOutput = Nabejaee(A , B , F)
clc;
% پاک کرد صفحه
Max_Count=20;
% حداکثر تعداد دفعاتي که حلقه تکرار شود
%در صورتي که مقدار تابع در ابتدا و انتهاي تابع تغيير علامت ندد در مورد
%وجود ريشه ها نميتوان اظهار نظر کرد
if polyval(F,A)*polyval(F,B)>0
fprintf('\nF(%f) * F(%f) > 0 \nF(%f)=%f , F(%f)=%f\n',A,B,A,polyval(F,A),B,polyval(F,B));
return;
end
n=1;
%مقدار اوليه به شمارشگر کنرل حداکثر تکرار
X= ( (A*polyval(F,B) - (B*polyval(F,A)) ) / ( polyval(F,B) - polyval(F,A) ));
hold on;
%شروع به رسم نمودار
X_Temp=[A:0.1:B];
Y_Temp=polyval(F,X_Temp);
plot(X_Temp,Y_Temp);
while ( (polyval(F,X)~=0) & (n<=Max_Count) )
% رسم خطوط قرمز رنگ
Y_Temp=[polyval(F,A) polyval(F,B)];
X_Temp=[A B];
plot(X_Temp,Y_Temp,'Red');
X= ( (A*polyval(F,B) - (B*polyval(F,A)) ) / ( polyval(F,B) - polyval(F,A) ));
fprintf('%2d) [%f,%f] , X=%f , F(%f)=%f , F(%f)=%f , F(%f)=%f\n',n,A,B,X,A,polyval(F,A),B,polyval(F,B),X,polyval(F,X))
if polyval(F,A)*polyval(F,X)<0
B=X;
else
A=X;
end
n=n+1;
end
hold off;
% رسم نمودار
NabejaeeOutput = X;
این برنامه بخشی از پروژه آنالیز عددی 1 ام بود و سه تا برنامه دیگه هم هست که روشهای دیگه حل معادلات هست . فکر کنم تو نوشتن سینتکس کمکت میکنه اگه اونا رو هم خواستی میلت رو بده واست میل کنمکد:clc;
f=[1 0 -2];
Nabejaee(1,2,f)
برنامه بالا در مورد دریافت بهت کمک میکنه و برای مشتق هم از تابع polyder استفاده کن
این هم نمونه ای که با مشتق کار میکنه :
کد:% روش نيوتن رافسون
% ورودي مقدار اوليه طول و آرايه ضرايب چند جمله اي
% خروجي : صفر تابع
function NROutput = NR( X0 , F )
clc;
FPerim=polyder(F);
n=1;
X= X0;
hold on;
%Stop Show Drawing
%Draw Function F
X_Temp= X-5:0.01:X+5;
Y_Temp= polyval(F,X_Temp);
plot(X_Temp,Y_Temp);
while ( (n <= 20) & (round(polyval(F,X)*10000) ~= 0) )
X0=X;
X= X - (polyval(F,X))/(polyval(FPerim,X));
%Calcute X(n+1)
fprintf('\n%2d) X%2d=%f F(%f)=%f',n,n,X,X,polyval(F,X));
X_Temp=[X X0];
Y_Temp=[0 polyval(F,X0)];
plot(X_Temp,Y_Temp,'Red')
n = n+1;
end
hold off;
%Show Drawing
NROutput=X;