Solving Non-Linear Equations
2 ビュー (過去 30 日間)
古いコメントを表示
I am stuck up with the following program. I want to solve for 3 variables b, n, k. Please help.
for b,n,k I am getting <1x1 sym> how to get the values? Any help would be highly apprecaited.. Thanks
clear
T = 0.89;
R = 0.11;
lamda = 0.98;
d = 1.65;
pi = 3.14;
phit = 0.60;
phir = 0.31;
syms b n k;
eq1 = 'phir=atan(((-2)*k)/(1-((n)^2)-((k)^2)))';
eq2 = 'T=((1-R)^2)+((4)*(R)*(sin(phir)^2))*(exp((-((4*(pi)*k)/(lamda)))*(d)))/((1)-(R)*(exp(((4*(pi)*k)/(lamda))*(d)))^2)+((4)*(R)*(exp(((4*(pi)*k)/(lamda))*(d)))*(sin(b+phir)^2))';
eq3 = 'phit=((2)*(pi)*(n)*(d)/(lamda))+arctan(((k)*((n)^2+(k)^2-(1)))/(((k)^2+(n)^2)*(2+n)*(n)))+arctan((((R)*exp((4)*(pi)*(k))/(lamda))*(sin(b+phir)^2))/((1)-((R)*exp((4)*(pi)*(k)))/(lamda))*(cos(b+phir)^2))';
eq4 = 'b=(((2)*(pi)*(n)*(d))/(lamda))';
sol1 = solve(eq1);
sol2 = solve(eq2);
sol3 = solve(eq3);
sol4 = solve(eq4);
0 件のコメント
回答 (3 件)
Rizwana Junaid
2012 年 10 月 15 日
you can solve ur system by calling ode45 function. for that you will have to save all equations in one file and call the solver function from the other main file.
Walter Roberson
2012 年 10 月 15 日
You need to remove all of your apostrophes. If you are using R2012a or later, change your = that are in the equations to == but if you are using an earlier version, convert your A=B forms to (B)-(A) or (A)-(B) with no = or == in the definition of the expression. For example,
eq4 = (b)-((((2)*(pi)*(n)*(d))/(lamda)));
The result of your solve() is going to be a structure array with one field for each variable name. Access the field to see the value, such as
sol4.k
2 件のコメント
Sachin Ganjare
2012 年 10 月 16 日
Undefined function or method 'arctan' for input arguments of type 'sym'.
Matlab function name is 'atan' and not 'arctan'
Hope it helps!!!
3 件のコメント
Sachin Ganjare
2012 年 10 月 16 日
編集済み: Sachin Ganjare
2012 年 10 月 16 日
Probably you do not need quotes in your equation.
Check in your version of matlab documentation.
Hope it helps!!!
参考
カテゴリ
Help Center および File Exchange で Equation Solving についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!