From: Jesse M. <jes...@gm...> - 2014-05-20 21:52:15
|
Hello, I am trying to perform a Parameter Sweep of Transient AC RLC Circuit using Qucs. What I want to do is plot the max current, and max voltage across each component for each value of the swept parameter, the inductor value. I don't believe I can attach anything to these messages so I copied and pasted the text of my schematic to the end of this message. My problem is that I can not seem to plot something that is dependent on time and the swept parameter. I tried to add an equation to the schematic that was Vmax=max(Vl.Vt), where Vl is the voltage from a voltage probe across my inductor. When I plot this there is only one value. I thought there would be one value for each value of the swept parameter. Also I can see that something like my Vl.Vt parameter is dependent on time and the swept variable. When I go to plot this in the Cartesian plot the x axis is time and I can not figure out how to plot with the other variable as the x axis variable. I tried typing "Vl.Vt(INDvalue)" in the Graph Input but it was invalid. I am not sure if I am suffering from a fundamental misunderstanding or the software just doesn't support what I am trying to do. Thank you for any help, Jesse <Qucs Schematic 0.0.18> <Properties> <View=0,-60,800,844,1,0,60> <Grid=10,10,1> <DataSet=EGR330_HW_1.dat> <DataDisplay=EGR330_HW_1.dpl> <OpenDisplay=1> <Script=EGR330_HW_1.m> <RunScript=0> <showFrame=0> <FrameText0=Title> <FrameText1=Drawn By:> <FrameText2=Date:> <FrameText3=Revision:> </Properties> <Symbol> </Symbol> <Components> <C C1 1 600 250 17 -26 0 1 "1 pF" 1 "" 0 "neutral" 0> <GND * 1 340 390 0 0 0 0> <L L1 1 460 200 -26 10 0 0 "INDvalue" 1 "" 0> <R R1 1 310 200 -26 15 0 0 "100 Ohm" 1 "26.85" 0 "0.0" 0 "0.0" 0 "26.85" 0 "US" 0> <.TR TR1 1 140 440 0 77 0 0 "lin" 1 "0" 1 "100 ms" 1 "101" 0 "Trapezoidal" 0 "2" 0 "1 ns" 0 "1e-16" 0 "150" 0 "0.001" 0 "1 pA" 0 "1 uV" 0 "26.85" 0 "1e-3" 0 "1e-6" 0 "1" 0 "CroutLU" 0 "no" 0 "yes" 0 "0" 0> <VProbe Vr 1 310 130 28 -31 0 0> <VProbe Vl 1 460 130 28 -31 0 0> <VProbe Vc 1 730 250 -16 28 0 3> <Vac V1 1 80 250 18 -26 0 1 "1 V" 1 "60 Hz" 1 "0" 0 "0" 0> <IProbe I 1 170 200 -26 16 0 0> <.SW SWL 1 300 440 0 77 0 0 "TR1" 1 "lin" 1 "INDvalue" 1 "1 nH" 1 "10 uH" 1 "100" 1> <Eqn MAX 1 510 450 -23 16 0 0 "MAX_Vl=max(Vl.Vt)" 1 "yes" 0> </Components> <Wires> <340 200 360 200 "" 0 0 0 ""> <490 200 510 200 "" 0 0 0 ""> <600 200 600 220 "" 0 0 0 ""> <600 280 600 340 "" 0 0 0 ""> <340 340 600 340 "" 0 0 0 ""> <340 340 340 390 "" 0 0 0 ""> <250 200 280 200 "" 0 0 0 ""> <250 170 250 200 "" 0 0 0 ""> <250 170 300 170 "" 0 0 0 ""> <300 150 300 170 "" 0 0 0 ""> <320 150 320 170 "" 0 0 0 ""> <320 170 360 170 "" 0 0 0 ""> <360 200 420 200 "" 0 0 0 ""> <360 170 360 200 "" 0 0 0 ""> <420 200 430 200 "" 0 0 0 ""> <420 170 420 200 "" 0 0 0 ""> <420 170 450 170 "" 0 0 0 ""> <450 150 450 170 "" 0 0 0 ""> <470 150 470 170 "" 0 0 0 ""> <470 170 510 170 "" 0 0 0 ""> <510 200 600 200 "" 0 0 0 ""> <510 170 510 200 "" 0 0 0 ""> <710 200 710 240 "" 0 0 0 ""> <600 200 710 200 "" 0 0 0 ""> <710 260 710 340 "" 0 0 0 ""> <600 340 710 340 "" 0 0 0 ""> <80 340 340 340 "" 0 0 0 ""> <80 280 80 340 "" 0 0 0 ""> <80 200 80 220 "" 0 0 0 ""> <80 200 140 200 "" 0 0 0 ""> <200 200 250 200 "" 0 0 0 ""> </Wires> <Diagrams> </Diagrams> <Paintings> </Paintings> |
From: Guilherme B. T. <gui...@gm...> - 2014-05-22 07:50:24
|
Hi, As far as I can tell the equations only look on the last set of the swept simulations results. Fixing this behavior is a nice feature request. You can always use Octave or Python to read the result file and compute things from there. Regards, Guilherme On 5/20/14, 11:52 PM, Jesse Millwood wrote: > Hello, > > I am trying to perform a Parameter Sweep of Transient AC RLC Circuit > using Qucs. What I want to do is plot the max current, and max voltage > across each component for each value of the swept parameter, the > inductor value. > > I don't believe I can attach anything to these messages so I copied and > pasted the text of my schematic to the end of this message. > > My problem is that I can not seem to plot something that is dependent > on time and the swept parameter. I tried to add an equation to the > schematic that was Vmax=max(Vl.Vt), where Vl is the voltage from a > voltage probe across my inductor. When I plot this there is only one > value. I thought there would be one value for each value of the swept > parameter. > > Also I can see that something like my Vl.Vt parameter is dependent on > time and the swept variable. When I go to plot this in the Cartesian > plot the x axis is time and I can not figure out how to plot with the > other variable as the x axis variable. I tried typing "Vl.Vt(INDvalue)" > in the Graph Input but it was invalid. > > I am not sure if I am suffering from a fundamental misunderstanding or > the software just doesn't support what I am trying to do. > > Thank you for any help, > Jesse > > > <Qucs Schematic 0.0.18> > <Properties> > <View=0,-60,800,844,1,0,60> > <Grid=10,10,1> > <DataSet=EGR330_HW_1.dat> > <DataDisplay=EGR330_HW_1.dpl> > <OpenDisplay=1> > <Script=EGR330_HW_1.m> > <RunScript=0> > <showFrame=0> > <FrameText0=Title> > <FrameText1=Drawn By:> > <FrameText2=Date:> > <FrameText3=Revision:> > </Properties> > <Symbol> > </Symbol> > <Components> > <C C1 1 600 250 17 -26 0 1 "1 pF" 1 "" 0 "neutral" 0> > <GND * 1 340 390 0 0 0 0> > <L L1 1 460 200 -26 10 0 0 "INDvalue" 1 "" 0> > <R R1 1 310 200 -26 15 0 0 "100 Ohm" 1 "26.85" 0 "0.0" 0 "0.0" 0 > "26.85" 0 "US" 0> > <.TR TR1 1 140 440 0 77 0 0 "lin" 1 "0" 1 "100 ms" 1 "101" 0 > "Trapezoidal" 0 "2" 0 "1 ns" 0 "1e-16" 0 "150" 0 "0.001" 0 "1 pA" 0 "1 > uV" 0 "26.85" 0 "1e-3" 0 "1e-6" 0 "1" 0 "CroutLU" 0 "no" 0 "yes" 0 "0" > 0> > <VProbe Vr 1 310 130 28 -31 0 0> > <VProbe Vl 1 460 130 28 -31 0 0> > <VProbe Vc 1 730 250 -16 28 0 3> > <Vac V1 1 80 250 18 -26 0 1 "1 V" 1 "60 Hz" 1 "0" 0 "0" 0> > <IProbe I 1 170 200 -26 16 0 0> > <.SW SWL 1 300 440 0 77 0 0 "TR1" 1 "lin" 1 "INDvalue" 1 "1 nH" 1 "10 > uH" 1 "100" 1> > <Eqn MAX 1 510 450 -23 16 0 0 "MAX_Vl=max(Vl.Vt)" 1 "yes" 0> > </Components> > <Wires> > <340 200 360 200 "" 0 0 0 ""> > <490 200 510 200 "" 0 0 0 ""> > <600 200 600 220 "" 0 0 0 ""> > <600 280 600 340 "" 0 0 0 ""> > <340 340 600 340 "" 0 0 0 ""> > <340 340 340 390 "" 0 0 0 ""> > <250 200 280 200 "" 0 0 0 ""> > <250 170 250 200 "" 0 0 0 ""> > <250 170 300 170 "" 0 0 0 ""> > <300 150 300 170 "" 0 0 0 ""> > <320 150 320 170 "" 0 0 0 ""> > <320 170 360 170 "" 0 0 0 ""> > <360 200 420 200 "" 0 0 0 ""> > <360 170 360 200 "" 0 0 0 ""> > <420 200 430 200 "" 0 0 0 ""> > <420 170 420 200 "" 0 0 0 ""> > <420 170 450 170 "" 0 0 0 ""> > <450 150 450 170 "" 0 0 0 ""> > <470 150 470 170 "" 0 0 0 ""> > <470 170 510 170 "" 0 0 0 ""> > <510 200 600 200 "" 0 0 0 ""> > <510 170 510 200 "" 0 0 0 ""> > <710 200 710 240 "" 0 0 0 ""> > <600 200 710 200 "" 0 0 0 ""> > <710 260 710 340 "" 0 0 0 ""> > <600 340 710 340 "" 0 0 0 ""> > <80 340 340 340 "" 0 0 0 ""> > <80 280 80 340 "" 0 0 0 ""> > <80 200 80 220 "" 0 0 0 ""> > <80 200 140 200 "" 0 0 0 ""> > <200 200 250 200 "" 0 0 0 ""> > </Wires> > <Diagrams> > </Diagrams> > <Paintings> > </Paintings> > > > ------------------------------------------------------------------------------ > "Accelerate Dev Cycles with Automated Cross-Browser Testing - For FREE > Instantly run your Selenium tests across 300+ browser/OS combos. > Get unparalleled scalability from the best Selenium testing platform available > Simple to use. Nothing to install. Get started now for free." > http://p.sf.net/sfu/SauceLabs > _______________________________________________ > Qucs-help mailing list > Quc...@li... > https://lists.sourceforge.net/lists/listinfo/qucs-help |
From: Guilherme B. T. <gui...@gm...> - 2014-05-22 16:19:47
|
Hi Jesse, Frankly I haven't used Octave. Below you can find a snippet that uses the `parse_result.py`: https://github.com/Qucs/qucs/blob/master/qucs/qucs/python/parse_result.py ~~~~ import parse_result as pa data = pa.parse_file('divider.dat') # <= result file print data.keys() sweep=data['Vdiv.Vt'] # <= node or variable for node in sweep: print max(node) ~~~~ Hope it helps! Regards, Guilherme On 5/22/14, 10:03 AM, Jesse Millwood wrote: > Guilherme, > > Thank you, using Python or Octave was my next step. I just wasn't sure > if I was missing something obvious or not. I am looking at my .dat > file from my simulation and I only see variables dependent on time. Is > there a way to extract the data from the simulation of each value of > the swept parameter? > > Thank you, > Jesse > > On Thu, May 22, 2014 at 3:50 AM, Guilherme Brondani Torri > <gui...@gm...> wrote: >> Hi, As far as I can tell the equations only look on the last set of >> the swept simulations results. Fixing this behavior is a nice feature >> request. You can always use Octave or Python to read the result file >> and compute things from there. Regards, Guilherme On 5/20/14, 11:52 >> PM, Jesse Millwood wrote: >> >> Hello, I am trying to perform a Parameter Sweep of Transient AC >> RLC Circuit using Qucs. What I want to do is plot the max >> current, and max voltage across each component for each value of >> the swept parameter, the inductor value. I don't believe I can >> attach anything to these messages so I copied and pasted the text >> of my schematic to the end of this message. My problem is that I >> can not seem to plot something that is dependent on time and the >> swept parameter. I tried to add an equation to the schematic that >> was Vmax=max(Vl.Vt), where Vl is the voltage from a voltage probe >> across my inductor. When I plot this there is only one value. I >> thought there would be one value for each value of the swept >> parameter. Also I can see that something like my Vl.Vt parameter >> is dependent on time and the swept variable. When I go to plot >> this in the Cartesian plot the x axis is time and I can not >> figure out how to plot with the other variable as the x axis >> variable. I tried typing "Vl.Vt(INDvalue)" in the Graph Input but >> it was invalid. I am not sure if I am suffering from a >> fundamental misunderstanding or the software just doesn't support >> what I am trying to do. Thank you for any help, Jesse <Qucs >> Schematic 0.0.18> <Properties> <View=0,-60,800,844,1,0,60> >> <Grid=10,10,1> <DataSet=EGR330_HW_1.dat> >> <DataDisplay=EGR330_HW_1.dpl> <OpenDisplay=1> >> <Script=EGR330_HW_1.m> <RunScript=0> <showFrame=0> >> <FrameText0=Title> <FrameText1=Drawn By:> <FrameText2=Date:> >> <FrameText3=Revision:> </Properties> <Symbol> </Symbol> >> <Components> <C C1 1 600 250 17 -26 0 1 "1 pF" 1 "" 0 "neutral" >> 0> <GND * 1 340 390 0 0 0 0> <L L1 1 460 200 -26 10 0 0 >> "INDvalue" 1 "" 0> <R R1 1 310 200 -26 15 0 0 "100 Ohm" 1 "26.85" >> 0 "0.0" 0 "0.0" 0 "26.85" 0 "US" 0> <.TR TR1 1 140 440 0 77 0 0 >> "lin" 1 "0" 1 "100 ms" 1 "101" 0 "Trapezoidal" 0 "2" 0 "1 ns" 0 >> "1e-16" 0 "150" 0 "0.001" 0 "1 pA" 0 "1 uV" 0 "26.85" 0 "1e-3" 0 >> "1e-6" 0 "1" 0 "CroutLU" 0 "no" 0 "yes" 0 "0" 0> <VProbe Vr 1 310 >> 130 28 -31 0 0> <VProbe Vl 1 460 130 28 -31 0 0> <VProbe Vc 1 730 >> 250 -16 28 0 3> <Vac V1 1 80 250 18 -26 0 1 "1 V" 1 "60 Hz" 1 "0" >> 0 "0" 0> <IProbe I 1 170 200 -26 16 0 0> <.SW SWL 1 300 440 0 77 >> 0 0 "TR1" 1 "lin" 1 "INDvalue" 1 "1 nH" 1 "10 uH" 1 "100" 1> <Eqn >> MAX 1 510 450 -23 16 0 0 "MAX_Vl=max(Vl.Vt)" 1 "yes" 0> >> </Components> <Wires> <340 200 360 200 "" 0 0 0 ""> <490 200 510 >> 200 "" 0 0 0 ""> <600 200 600 220 "" 0 0 0 ""> <600 280 600 340 >> "" 0 0 0 ""> <340 340 600 340 "" 0 0 0 ""> <340 340 340 390 "" 0 >> 0 0 ""> <250 200 280 200 "" 0 0 0 ""> <250 170 250 200 "" 0 0 0 >> ""> <250 170 300 170 "" 0 0 0 ""> <300 150 300 170 "" 0 0 0 ""> >> <320 150 320 170 "" 0 0 0 ""> <320 170 360 170 "" 0 0 0 ""> <360 >> 200 420 200 "" 0 0 0 ""> <360 170 360 200 "" 0 0 0 ""> <420 200 >> 430 200 "" 0 0 0 ""> <420 170 420 200 "" 0 0 0 ""> <420 170 450 >> 170 "" 0 0 0 ""> <450 150 450 170 "" 0 0 0 ""> <470 150 470 170 >> "" 0 0 0 ""> <470 170 510 170 "" 0 0 0 ""> <510 200 600 200 "" 0 >> 0 0 ""> <510 170 510 200 "" 0 0 0 ""> <710 200 710 240 "" 0 0 0 >> ""> <600 200 710 200 "" 0 0 0 ""> <710 260 710 340 "" 0 0 0 ""> >> <600 340 710 340 "" 0 0 0 ""> <80 340 340 340 "" 0 0 0 ""> <80 >> 280 80 340 "" 0 0 0 ""> <80 200 80 220 "" 0 0 0 ""> <80 200 140 >> 200 "" 0 0 0 ""> <200 200 250 200 "" 0 0 0 ""> </Wires> >> <Diagrams> </Diagrams> <Paintings> </Paintings> >> ------------------------------------------------------------------------------ >> "Accelerate Dev Cycles with Automated Cross-Browser Testing - For >> FREE Instantly run your Selenium tests across 300+ browser/OS >> combos. Get unparalleled scalability from the best Selenium >> testing platform available Simple to use. Nothing to install. Get >> started now for free." http://p.sf.net/sfu/SauceLabs >> _______________________________________________ Qucs-help mailing >> list Quc...@li... >> https://lists.sourceforge.net/lists/listinfo/qucs-help >> |