حل معادله غیر خطی با تابع fsolve در متلب

حل معادله غیر خطی با تابع fsolve در متلب

من می خواستم معادله غیر خطی زیر رو با تابع fsolve متلب حل کنم ولی بهم خطا میده .

معادله غیر خطی این هست :

خودم هم توی متلب این جوری نوشتم :
PHP Code:

function out=myfun(E)

h=(6.63*10^-34)/(2*pi);

m=9.1*10^-31 ;

a=10^-10 ;

v0=1 ;

me1=sqrt(2*m*E / h ) ;

me2=sqrt(2*m*(v0-E)) / h ;

out=[me1*tan(me1*a)-me2 ];
آخرین ویرایش توسط skyzare; 2011/11/03 در تاریخ 23:20 انجام شده است.

دستور حل دستگاه معادلات غیرخطی در متلب دستور fsolve می باشد. ابتدا نحوه نوشتن این دستور را توضیح می دهم:
fsolve ( function , initial_guess)

ه جای کلمه function شما باید اسم تابع خودتان را به همراه @ بنویسید.

به جاب initial_guess هم باید بردار حدس اولیه برای متغیرهایتان را وارد کنید. فرض کنید که بردار حدس اولیه برای متغیرهایمان به صورت زیر است:

[x0 y0]=[2 -3]

حال می توانیم دستور fsolve را به کار ببریم:
>> M=fsolve(@myfun , [2 -3])

Equation solved.

fsolve completed because the vector of function values is near zero

as measured by the default value of the function tolerance, and

the problem appears regular as measured by the gradient.

<stopping criteria details>

M =

3.1074 -4.1460

همانطور که میبینید متلب یک بردار با دو عضو را به شما خواهد داد که عضو اول همان x و عضو دوم همان y می باشد. مسلما اگر دستگاه 4 معادله و 4 مجهول داشته باشید M دارای 4 مولفه خواهد بود.

نکات مهم : 1- روشهای تکراری مثل نیوتن بسیار به حدس اولیه حساس می باشند. چنانچه حدس اولیه خیلی پرت زده شود ممکن است الگوریتم نتواند جوابی برای معادلات بدست آورد!

2- حتما در نوشتن function دقت کنید که خروجی تابع (در اینجا out) باید ستونی باشد.