Info

此问题已关闭。 请重新打开它进行编辑或回答。

Symbolic calculation and to solve out a variable

1 次查看(过去 30 天)
I have the equation :
I_1=2.15 ; I_2=8.625;
n2=0.039; n1=0.1;
L_1=16.76*10^-10
L_2=15.56*10^-10;
I tried with this code but it is wrong:
clc
clear all
syms I1 I2 n1 n2 L1 L2 d
I_1=2.15 ; I_2=8.625;
n2=0.039; n1=0.1;
L_1=16.76*10^-10
L_2=15.56*10^-10;
I1/I2=((n1*L1)/(n2*L2))*(1-exp(-d/L1))/(exp(-d/L2));
I want to solve d

回答(2 个)

Steven Lord
Steven Lord 2020-3-11
Write your last line as an equation and use solve.

Ameer Hamza
Ameer Hamza 2020-3-11
编辑:Ameer Hamza 2020-3-11
First, your syntax is incorrect, as pointed out in Steven's answer. However, if you correct the syntax, the solve will still not be able to provide the solution because an analytical solution does not exist for such equations.
Your other choice is to use a numerical solver, such as fsolve. However, because of the scale of L_1 and L_2, and their presence in the exponent, directly using fsolve will get you nowhere. Therefore, you first need to scale the exponent and then apply fsolve. See the code
I_1=2.15 ; I_2=8.625;
n2=0.039; n1=0.1;
L_1=16.76*10^-10;
L_2=15.56*10^-10;
f = @(d) I_2/I_1*((n1*L_1)/(n2*L_2))*(1-exp(-d/L_1*10^-10))/(exp(-d/L_2*10^-10)) - 1;
d = fsolve(f, 1)*10^-10;
Here, I first scaled the exponents by a factor of 10^-10 and corrected the final solution by scaling again.

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by