0 = 0 /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ /* Comparison of our formulation - 21 December 2019: */ /* https://groups.io/g/nanovna-users/message/8459 */ /* */ /* with the one by Erik, PD0EK - 21 December 2019: */ /* https://groups.io/g/nanovna-users/message/8482 */ /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ /* Version : 20191221 - First - for 64-bit Maxima : 5.43.0 */ /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ /* (C) Copyright 2019 Nikolitsa I. Giannopoulou, OE3ZGN/SV7DMC, gin, and */ /* Petros E. Zimourtopoulos, OE3ZZP/SV7BAX, pez, gin.pez.arg@gmail.com, */ /* Antennas Research Group ARG, Informal Association of Individuals IAOI,*/ /* No Finance Involved NFI, ARG IAOI NFI, www.op4.eu, Austria/Hellas, EU */ /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ /* Permission is hereby granted, free of charge, to any person obtaining */ /* a copy of this software and associated documentation files (the */ /* "Software"), to deal in the Software without restriction, including */ /* distribute, sublicense, and/or sell copies of the Software, and to */ /* permit persons to whom the Software is furnished to do so, subject to */ /* the following conditions: */ /* */ /* The above copyright notice and this permission notice shall be */ /* included in all copies or substantial portions of the Software. */ /* */ /* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, */ /* EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF */ /* MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.*/ /* IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY */ /* CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, */ /* TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE */ /* SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. */ /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ /* /F/L/O/S/S/ Free Libre Open Source Software under "The MIT License" */ /* approved by OSI Open Source Initiative www.opensource.org/licenses/MIT*/ /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */; 1=1 /* ------------------------------------------------------------------- Using Erik's symbols and relations : b <- e00 c <- e11 a <- Delta-e Gamma <- G GammaM <- g and replacing: a -> aa b -> bb c -> cc Gn -> numerator of G Gd -> denominator of G we put: GnGPA -> numerator of our G formulation GdGPA -> denominator of our G formulation for Loads: S <- SHORT L <- Load O <- Open G <- DUT and for Measurements: s <- SHORT l <- Load o <- Open g <- DUT ------------------------------------------------------------------- */; linsolve ([b+s*S*c-a*S=s ,b+l*L*c-a*L=l,b+o*O*c-a*O=o], [a,b,c]); aa:(-1)*((l-o)*s*S+s*(o*O-l*L)-l*o*O+l*o*L)/((o*O+s*(L-O)-l*L)*S+L*(l*O-o*O)); bb:((s*(o*L-l*O)+l*o*O-l*o*L)*S+s*L*(l*O-o*O))/((o*O+s*(L-O)-l*L)*S+L*(l*O-o*O)); cc:(-1)*((l-o)*S+s*(O-L)-l*O+o*L)/((o*O+s*(L-O)-l*L)*S+L*(l*O-o*O)); linsolve([bb+g*G*cc-aa*G=g], [G]); Gn:expand((-1)*((((l-g)*s+(g-l)*o)*O+((g-o)*s+l*o-g*l)*L)*S+((o-l)*s-g*o+g*l)*L*O)); Gd:expand((((o-l)*s-g*o+g*l)*S+((g-o)*s+l*o-g*l)*O+((l-g)*s+(g-l)*o)*L)); GdGPA:expand((S-L)*(g*o+s*l)+(L-O)*(g*s+l*o)+(O-S)*(g*l+o*s)); GnGPA:expand(O*(S-L)*(g*o+s*l)+S*(L-O)*(g*s+l*o)+L*(O-S)*(g*l+o*s)); is(equal(Gn,-GnGPA)); is(equal(Gd,-GdGPA));