How to solve for 4 variables from 5 equations that I have derived.

Question

Shaumik 2012-10-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/49606-how-to-solve-for-4-variables-from-5-equations-that-i-have-derived

Hi, I have a problem in MATLAB. I am doing research in Optics ... I have to find 4 values from 5 equations and having some trouble. Kindly help.

I have tried reg this as below: I would be grateful to have any suggestions. Thanks in advance

clear

T = 0.89; R = 0.11; lamda = 0.98; d = 1.65; pi = 3.14; phit = 0.6; phir = 0.31; %((-(a))*d) == (23.03)*k;

% sym a b n k; % a and b are "alpha" and "beta"

a = sym('a'); b = sym('b'); n = sym('n'); k = sym('k');

eq1 = T == '({(1-R)^2+4*R*(sin(phir))^2}*exp((-a)*d))/((1-R*exp((-a)*d))^2+4*R*exp((-a)*d)*(sin(b+phir))^2)';

eq2 = (phit) == '((2*(pi)*n*d)/lamda)-(atan((k(n^2+k^2-1))/((k^2+n^2)*(2+n)*n)))+(atan((R*exp((-a)*d))*(sin2(b+phir))/(1-R*exp((-a)*d)*cos2(b+phir))))';

eq3 = (phir) == 'atan(((-2)*k)/(1-((n)^2)-((k)^2)))';

eq4 = (b) == '((2*(pi)*n*d)/lamda)';

eq5 = (a) == '((4*(pi)*k)/lamda)';

%eq6 = (e) == '((-(a))*d) == (23.03)*k';

[a b n k] = solve(eq1,eq2,eq3,eq4,eq5,a,b,n,k);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Azzi Abdelmalek 2012-10-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/49606-how-to-solve-for-4-variables-from-5-equations-that-i-have-derived#answer_60628

编辑：Azzi Abdelmalek 2012-10-2

在 MATLAB Online 中打开

Corrected code

syms a b n k;
eq1 ='T==(((1-R)^2+4*R*(sin(phir))^2)*exp((-a)*d))/((1-R*exp((-a)*d))^2+4*R*exp((-a)*d)*(sin(b+phir))^2)';
eq2 = 'phit==((2*(pi)*n*d)/lamda)-(atan((k(n^2+k^2-1))/((k^2+n^2)*(2+n)*n)))+(atan((R*exp((-a)*d))*(sin2(b+phir))/(1-R*exp((-a)*d)*cos2(b+phir))))';
eq3 = 'phir==atan(((-2)*k)/(1-((n)^2)-((k)^2)))';
eq4 = 'b==((2*(pi)*n*d)/lamda)';
eq5 =  'a==((4*(pi)*k)/lamda)';
sol= solve(eq1,eq2,eq3,eq4,eq5,a,b,n,k)

but it seems that you system does'nt have solutions. Notice that you have more equations then unknown variables

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Shaumik 2012-10-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/49606-how-to-solve-for-4-variables-from-5-equations-that-i-have-derived#answer_60649

Thanks a lot for your help Azzi.. I will work on this and get back.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 3

Shaumik 2012-10-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/49606-how-to-solve-for-4-variables-from-5-equations-that-i-have-derived#answer_60841

hello

I just put in the value of 'a' in the equations (program given below the error message), but getting this error. Any suggestions would be very helpful. Thanks

________________________________________

??? Error using ==> solve>getEqns at 182 ' 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)) ' is not a valid expression or equation.

Error in ==> solve at 67 [eqns,vars] = getEqns(varargin{:});

Error in ==> RI at 19 sol = solve(eq1,eq2,eq3,eq4,b,n,k);

________________________________________________________

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))-(atan((k(n^2+k^2-1))/((k^2+n^2)*(2+n)*n)))+(atan((R*exp(((4*(pi)*k)/(lamda))*d))*(sin2(b+phir))/(1-R*exp(((4*(pi)*k)/(lamda))*d)*cos2(b+phir))))';

eq4 = 'b==(((2)*(pi)*(n)*(d))/(lamda))';

%eq5 = 'a==((4*(pi)*k)/(lamda))';

sol = solve(eq1,eq2,eq3,eq4,b,n,k);