Solving Transcendental equation

6 views (last 30 days)
tuan
tuan on 22 Sep 2011
Hi all,
I try to solve these system equations in matlab. I try to use fzero and fsolve but it did not work (R2010a). The system equation is
Re(λ)= μ=-αe^(-μ) cos(ω)
Im(λ)=ω=-αe^(-μ) sin(ω)
I appreciate for your help.
TN
  5 Comments
tuan
tuan on 31 Oct 2011
It just 2 variable mu and w, alpha is just a parameter, we can chose in any interval.
Walter Roberson
Walter Roberson on 31 Oct 2011
lambda and alpha are fixed values for any one problem? And you are looking for mu and w values that satisfy
lambda = alpha * exp(-mu + I*w)
Then
lambda/alpha = exp(-mu + I*w)
ln(lambda/alpha) = -mu + I*w
If mu and w are constrained to be real-valued then this would appear to have a single solution (unless alpha is 0).

Sign in to comment.

Answers (2)

Walter Roberson
Walter Roberson on 4 Oct 2011
You will need to indicate which variable(s) you are trying to solve for.
lambda = alpha * exp(-mu + I*w)
If you are trying to solve for alpha = 0, then that happens if either lambda = 0 or mu or w are infinite.
If you are trying to solve for lambda = 0, then that happens if either alpha = 0 or else mu and w are both zero.
You can also solve for mu or w being 0 without difficulty.

tuan
tuan on 31 Oct 2011
I don't think matlab can solve an equation with 2 variables like Walter said. I have the code that matlab can solve with 1 iteration but can not solve for 19 iteration. Any suggestion let me know thanks.
[a,u,v] = solve('-a*exp(-u)*cos(v)=u', 'v=a*exp(-u)*sin(v)', 'a=1')
it work well but with for loop
for a=0:0.1:1.8
[a,u,v] = solve('-a*exp(-u)*cos(v)=u', 'v=a*exp(-u)*sin(v)','a')
plot(a,u)
end
and it came out with this
Warning: 2 equations in 1 variables.
Warning: Explicit solution could not be found.
> In solve at 81
In Untitled2 at 3
a =
[ empty sym ]
u =
[]
v =
[]

Categories

Find more on Loops and Conditional Statements in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!