حل معادله با تابع كتاب خانه اي ()fzero
fzero سعي ميكند كه مقدار ريشه واقعي را حول و حوش يك حدس اوليه يا نقطه شروع starting point كه دستي وارد ميشود، پيدا كند ()fzero در واقع محل تغيير علامت تابع را پيدا ميكند، لذا بهتر است براي توابع پيوسته به كار رود. اين تابع كه به اين شكل نوشته ميشود (fzero(f,x0 دو آرگومان اصلي دارد. تابع f كه ميتواند به صورت رشته، گيره، يا خط فرمان وارد شود. حدس اوليه x0 ، مقداري است حتي الامكان نزديك به ريشه و دستي وارد ميشود.
پارامترهاي ديگر:
P كه به اين صورت وارد ميشود (fzero(f,X0,P . در صورت مراجعه به گيره تابع، P مقدار ورودي تابع f خواهد بود.
options كه به اين صورت وارد ميشود (fzero(f,X0,options براي درك بهتر options به مثال زير يا help optimset مراجعه كنيد.
مثال ها: ارسال به صورت رشته
>> y = ‘cos(x)+sin(x)+log(x)’;
>> X0 = fzero(y,1) X0 = 0.2885
ريشه را امتحان ميكنيم:
>> cos(X0) + sin(X0) + log(X0)
ans = 0
ارسال به صورت گيره
function y = fz(q)
y = exp(q)*sin(q)+log(q);
>> x0 = 1;
>> rt = fzero(@fz,x0)
rt = 0.4771
اگر مايل باشيم تعداد دفعات تكرار براي رسيدن به نتيجه نمايش داده شود:
>> options = optimset(‘Display’,’iter’);
x0 = 1;
>> rt = fzero(@fz,x0,options)
نتيجه را خودتان امتحان و مشاهده كنيد.
حل دستگاه معادلات غير خطي با()fsolve (جعبه ابزار بهينه سازي Optimization Toolbox)
فرم كلي اين تابع به شكل زير است:
x = fsolve(fun,x0,options,P1,P2, … ) fun
تابع مخصوصي است كه فرم معادلات غيرخطي خود را در آن ميگذاريم، و ويژه كاربرد در آرگومان ()fsolve يا نظائرش مانند ()ode45 است. بقيه پارامترها مشابه fzero هستند. اگر سمت چپ تساوي [x,fval] قرار دهيم، جواب معادلات دستگاه نيز در ازاي ريشه هاي به دست آمده چاپ خواهند شد.
مثال: دستگاه زير را حل ميكنيم.
function F = fsl(x)
F(1) = x(2)-x(1) ^ -x(2);
F(2) = 2*x(1)- x(2) ^ -x(1); % The following form of writing F is also accepted
% F = [x(2)-x(1) ^ -x(2); % 2*x(1)- x(2) ^ -x(1)];
>> format bank
>> x0 = [2 2];
>> [a,fv] = fsolve(@fsl,x0);
a = 2.98 0.55
fv = -0.00 -0.00
نوشتن fv به عنوان پارامتر مشابه عمل زير است:
>> fsl(a)
ans = -0.00 -0.00
نویسنده: مصطفی همت آبادی
مقاله متلب,مطلب,متلب,مقاله برق,مقاله قدرت,مقاله مطلب,مقاله سیمولینک,دانلود متلب,دانلود مقاله متلب,مقالهmatlab ,آموزش متلب,مطلب,متلب,آموزش برق,آموزش قدرت,آموزش مطلب,آموزش سیمولینک,دانلود متلب,دانلود آموزش متلب,آموزشmatlab ,پروژه متلب,مطلب,متلب,پروژه برق,پروژه قدرت,پروژه مطلب,پروژه سیمولینک,دانلود متلب,دانلود پروژه متلب,پروژهmatlab ,