• Support PF! Buy your school textbooks, materials and every day products Here!

Matlab: Finding a quicker method than simple iteration

  • Thread starter rewmck
  • Start date
  • #1
3
0

Homework Statement



I am trying to speed up a function in matlab that solves the equation by simple iteration:

λ= [itex]\frac{gT^{2}\frac{\text{tanh}(2\pi h)}{λ}}{2\pi}[/itex]

Where λ is wavelength, h is wave height, T is period and g=9.81.

Also here where it's a bit easier to read.

Homework Equations





The Attempt at a Solution



Currently I'm using simple iteration so I'm guessing lambda then putting the equation in a while loop where the new value of lambda is given by putting the old value in the above equation. The while terminates when the values of lambda converge.

This is one of the slowest parts of my code and thousands of iterations are done when I run my code.

Does anyone know of a method which converges faster or would be computationally faster? I tried the Newton Raphson method and rearranged the equation by making one side zero to get f(x) and differentiated for f(x). This converges quicker but often gives the wrong answer.

Thanks
 
Last edited:

Answers and Replies

  • #2
gneill
Mentor
20,795
2,773
Matlab probably has built-in solvers that may behave better. Investigate "roots of scalar functions", for example.
 
  • #3
3
0
Thanks for the push in the right direction. Matlab has a built in function fzero which gets the roots. Iteration is quicker for 99% of my calculations but the final 1% is much faster with fzero so I have found the point where it becomes beneficial to use fzero and put an if in my function. So now the function is about 50 times faster. Thanks!
 
  • #4
rewmck - what are you trying to solve for? Also, is that formulae stated correctly? I read your equation as follows:

[tex]
\lambda = \frac{g \cdot T^2}{2\, \pi} \cdot \frac{\tanh(2\, \pi \, h)}{\lambda}
[/tex]
 
  • #5
3
0
Hi Theo, solving for lambda and you're right the equation is wrong, it should be:
[tex]
\lambda = \frac{g \cdot T^2}{2\, \pi} \cdot \tanh\left(\frac{2\, \pi \, h}{\lambda}\right)
[/tex]
Sorry about that.

So I solve for lambda over a range of T values. I discovered that the fzero function which works much quicker than simple iteration for large T values is actually giving completely the wrong answer but at smaller T values it gives the same as simple iteration. I'm not sure what is happening as I haven't had the chance to look closely at the fzero code but I guess it is round off errors or diverging. I think it is probably diverging as it is out by a factor of 10^16 for the value I checked.


rewmck - what are you trying to solve for? Also, is that formulae stated correctly? I read your equation as follows:
...
 
  • #6
gneill
Mentor
20,795
2,773
Have you plotted the function f(λ) = 0 to investigate why it might be 'tricky' in some regions of its domain?
 
  • #7
20,230
4,264
Have you plotted the function f(λ) = 0 to investigate why it might be 'tricky' in some regions of its domain?
Derive analytic expressions for what lambda approaches at very large T and at very small T.
 

Related Threads on Matlab: Finding a quicker method than simple iteration

Replies
17
Views
11K
  • Last Post
Replies
2
Views
1K
Replies
9
Views
473
  • Last Post
Replies
20
Views
16K
Replies
0
Views
3K
  • Last Post
Replies
6
Views
29K
  • Last Post
Replies
4
Views
1K
Replies
1
Views
679
Top