diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2013-09-17 16:01:08 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2013-09-17 16:01:08 -0400 |
| commit | 0c533a1250ec1a1c6966db7854beb2b8fe5ea78e (patch) | |
| tree | 867c5d8b5b07834de2ccab3572d56505a3f8b3ff | |
| parent | bb04fb5d646defa7bfc061ec11cebc042d320299 (diff) | |
| download | mode_match-0c533a1250ec1a1c6966db7854beb2b8fe5ea78e.tar.gz mode_match-0c533a1250ec1a1c6966db7854beb2b8fe5ea78e.zip | |
cavity stability finction
| -rw-r--r-- | isCavityStable.m | 19 |
1 files changed, 19 insertions, 0 deletions
diff --git a/isCavityStable.m b/isCavityStable.m new file mode 100644 index 0000000..86dbd99 --- /dev/null +++ b/isCavityStable.m @@ -0,0 +1,19 @@ +function ret=isCavityStable (abcd_cavity) + % finds if a cavity is stable baced on its abcd matrix + % see eq.4 for stability condition + % doi:10.1088/1464-4258/7/1/002 + % http://dx.doi.org/10.1088/1464-4258/7/1/002 + + A=abcd_cavity(1,1); + B=abcd_cavity(1,2); + C=abcd_cavity(2,1); + D=abcd_cavity(2,2); + + u= (A+D+2)/4 + if ( (0 < u) && (u < 1) ) + ret = true; + else + ret = false; + end + +end |