Dٹ DC} D D5 D9 D Df D D D- D: E& E E=J Ek E
% E E- Ev E EE E E2! E E G E"u E% E(8 E*M E-j E0 E3l} E6V E9MM E E7 Eٲ E; E E{ E1 Ea E Ei E EF E/ F F' F8 Ff F
> F F( F|Q F} F@~ F} F, F F B F"ܲ F%& F(2 F*A F-E F0 F3f% F6Pb F9F F F F9 F/ F< F3 F Fb FM) FCy FF FUD Fq% Fř F F: Fd F F/| F٪l F3 F Fs F)R F Fq F FY F` F$ G G| G3 Gb G
W G G# Gw, GB G;. G G'( G G J J J J +j J"ū J%j J( J* J- J0n J3L J66 J9, J<. J?= JBY JE JH JK JOJs JR JV JY} J]x J`' JdV Jh Jk Jo Js Jw J{= J J J% J#0 JQ( J0 Jo J Je, J J( J3 J J J( J Jg J J J Jk JI J3b J)G J+| J:5 JU J~ Jȳ J^ JF JҤ J2 Jً J J. JR J~ Jͬ J JV J|U J J Kߔ K K KN K
K K Kb K[ K% . ¬ Angle ? ]o ]o B
Bq$BEBIBPBB7BaB1.B B BBiB6*BB&BB^B&BB;ByB>BB0BBFB7BBB?SBBBo'B'BBRBJBBrBc1BB&BqVB MB
BtB BźBlgBBBXBBB8BBqB *B0B#aB?8lB@KBA]:BBlBCyBDBEZBF.BG9BHYBIgBJ@BKBLvBMf0BNQBO9BP wBP)BQءBRBSiBTPBUBUޓBVXBWY$BXBXrBYjBZBZZB[LqB[B\qB\ B]bB]B^wB^2B_VB_^B` B`x B`BaBaaBa"BaBbBb>BbeBbBbBbTBbIBbǜBbKBbOBbBbLBb=BbayBb8Bb BaBa'BaSBa B`B`a2B`B_B_2GB^B^FB]VB]AB\EB\ B[BZBZ?BY:BXwBX$WBWcBV*BU5BTBT#BSDxBR_ BQsBPBOBNBMVBLBBKvBJcBIK'BH-qBG
BEBD~BCBBLBA[B?MB>B=CB;B::B9OB7B6bB58B3CB2lB1{B/KB."B, B+7*B)B(AB&ÌB%CB#B"<6B B .?B B YBBBsBABVBKB61BBiB B
B dB8BF2BB*;BB AvAAAfA7AGAAAAAAA[A A=_A` A͇AʴIAA (AYAsAA/AAA99ASAXAwRAKAiKA[AsAA A0ADAx A$AA1AP|AAWAAjAgA| AxiAtHAp3&Al)CAh*Ad8A`PA\tnAXlAT}AQ"AMriAI
AF2QABA? CA;A8/HA4A1j_A.A*~A'A$TA!&A AA}AAÅA GA"AA7A f
@8@@<@ө@y?@,@M@ջ3@Е@}V@q@p@|9@.@W@o@,@\K@@@^.@@8@@6@I@v@m@e%@\x@Se@KGt@B3@:I#@1ۿ@)y@!!@c@v@U@ #[???߬?u?KLP? u=>P
>6>|=8% 5(s6 "_>^#~0nz꿯@ Ef _l!i~)t2*D:qCK T8\{enOw 1vQQz9Gw ^ %RPی*Ӎ("~y ~ tv d&\H]f|xP ! % (_+.2_$5Ȧ9<kik7m.Un/phr su&vx>yw{Hv|~CugP_m\ ̷
y$hatm^u64g$ rBϐZabjsupeܞR9R_\0X^hC Q9nF)
F §§C§¨ ¨Pc¨¨¨ܵ©
©6©b©©C© ª Sª3ª[ªªªª@« «>e«b«>«««^¬ V¬,¬L¬l%¬¬¬x¬ :VDq{C3® ®$>®\ TransformAxis .s@>\ rescaleMin rescaleMax tickDensity @ tickSep @ LastTabUsed ModifyingTransformAxis TargetName DZ.s@ #Target: \K(1,16019,65535)\f01Graph0S_name ame DZ.s@ Panel0
WindowCoordinates s@>\ e b @ @ W_windowCoordinates ? ? ? ? 4 NewTransformAxisPanel4603029245730 ! "
PSFontInfo inates s@>\ C & * TTtoPS ? ? ? ? 5 HelveticaArialHelvetica-NarrowArial NarrowPalatinoBook AntiquaBookmanBookman Old StyleAvant GardeCentury GothicNew Century SchlbkCentury SchoolbookCourierCourier NewZapf ChanceryMonotype CorsivaZapf DingbatsMonotype SortsSymbolTimesTimes New RomanHelveticaHelveticaHelvetica-NarrowHelvetica-NarrowPalatinoPalatinoBookmanBookmanAvantGardeAvantGardeNewCenturySchlbkNewCenturySchlbkCourierCourierZapfChanceryZapfChanceryZapfDingbatsZapfDingbatsSymbolTimesTimes * 2 > E V a o ) 1 9 @ G Q [ k { v & 7 TTPSFNames ? ? ? ? 6 HelveticaHelvetica-NarrowPalatinoBookmanAvantGardeNewCenturySchlbkCourierZapfChanceryZapfDingbatsSymbolTimesHelveticaHelvetica-NarrowPalatino-RomanBookman-LightAvantGarde-BookNewCenturySchlbk-RomanCourierZapfChancery-MediumItalicZapfDingbatsSymbolTimes-RomanHelvetica-BoldHelvetica-Narrow-BoldPalatino-BoldBookman-DemiAvantGarde-DemiNewCenturySchlbk-BoldCourier-BoldZapfChancery-MediumItalicZapfDingbatsSymbolTimes-BoldHelvetica-ObliqueHelvetica-Narrow-ObliquePalatino-ItalicBookman-LightItalicAvantGarde-BookObliqueNewCenturySchlbk-ItalicCourier-ObliqueZapfChancery-MediumItalicZapfDingbatsSymbolTimes-ItalicHelvetica-BoldObliqueHelvetica-Narrow-BoldObliquePalatino-BoldItalicBookman-DemiItalicAvantGarde-DemiObliqueNewCenturySchlbk-BoldItalicCourier-BoldObliqueZapfChancery-MediumItalicZapfDingbatsSymbolTimes-BoldItalic ! ( 2 B I U a g l u % 2 > M b n * C O U a v & 6 z R & PSSymbolNames ? ? ? ? 7 SymbolZapfDingbats
) d WMT0IGR0 q Instructions Notebook0 Notebook0:Instructions H D 9 H H (Fg( H H ( d h 6 C 6 C $ Normal
<
B H H H H $ $ Z r Z f f q q q q q q q q Amplitude waves are Usually Functionx. Only names starting with func appear in popup menus. Function is the default wave. jw is the complex frequency wave.
p Waves can be plotted by entering expression of the form of function=cabs(200000/(1+jw/150)) in the command line
m The function Resize sets the scale of amplitude and phase windows based on the width of the amplitude window
G MakeNewPlot sets up the axes, and sets function to the magnitude of jw
PolesAndZeroes make a function from real poles and zeros. These are specifed as strings, with a null string indicating no singularity.
PairOfComplexPoless sets waves to the amplitude and phase of a complex pair. To combine with real poles, multiply waves Using Multiply macro. Z v 8 Normal+ #
Normal+ "
Geneva // Platform=Macintosh, IGORVersion=6.010, architecture=PowerPC
Silent 101 // use | as bitwise or -- not comment.
DefaultFont "Helvetica"
MoveWindow/P 241,409,959,726
OpenNotebook/N=Notebook0/W=(146,328,691,580)/J=87098 "Instructions"
MoveWindow/C 18,590,986,721
Graph2()
Amplitude()
Window Amplitude() : Graph
PauseUpdate; Silent 1 // building window...
Display /W=(53,77,989,480) Amplitude vs Omega as "Amplitude"
ModifyGraph grid=1
ModifyGraph log=1
ModifyGraph logLabel(bottom)=4
Label left "Gain"
Label bottom "Angular frequency \\F'Symbol'w \\F'Helvetica'(rads/s)"
SetAxis left 0.547109665450513,248.762359619141
ControlBar 25
SetVariable MinFreq,pos={16,6},size={128,15},proc=SetScales,title="Min omega = 10^"
SetVariable MinFreq,value= MinFreqP
SetVariable MaxFreq,pos={158,5},size={128,15},proc=SetScales,title="Max omega = 10^"
SetVariable MaxFreq,value= MaxFreqP
SetVariable p1,pos={400,5},size={50,15},proc=Calculate,title="P1",value= p1
SetVariable p2,pos={455,5},size={50,15},proc=Calculate,title="P2",value= p2
SetVariable p3,pos={512,5},size={50,15},proc=Calculate,title="P3",value= p3
SetVariable p4,pos={568,5},size={50,15},proc=Calculate,title="P4",value= p4
SetVariable z1,pos={624,5},size={50,15},proc=Calculate,title="Z1",value= z1
SetVariable z2,pos={680,5},size={50,15},proc=Calculate,title="Z2",value= z2
SetVariable z3,pos={736,5},size={50,15},proc=Calculate,title="Z3",value= z3
SetVariable z4,pos={792,5},size={50,15},proc=Calculate,title="Z4",value= z4
SetVariable Gain,pos={330,4},size={50,15},proc=Calculate,title="Gain"
SetVariable Gain,format="%g",value= Gain
EndMacro
Window Graph2() : Graph
PauseUpdate; Silent 1 // building window...
Display /W=(54,517,996,747) Angle vs Omega as "Phase"
ModifyGraph grid=1
ModifyGraph log(bottom)=1
ModifyGraph minor(bottom)=1
ModifyGraph logLabel(bottom)=4
ModifyGraph manTick(left)={0,30,0,0},manMinor(left)={2,0}
Label left "Phase (degrees)"
Label bottom "Angular frequency \\F'Symbol'w \\F'Helvetica'(rads/s)"
SetAxis/A/N=1 left
EndMacro
:#pragma rtGlobals=1 // Use modern global access method.
Function SetScales(ctrlName,varNum,varStr,varName) : SetVariableControl
String ctrlName
Variable varNum
String varStr
String varName
NVAR MinFreqP, MaxFreqP //Power of minimum and maximum frequency being plotted
Variable Points=1000
make /o /n=(points) Omega // Omega is the magnitude of the frequency
Omega=10^MinFreqP*10^(p*(MaxFreqP-MinFreqP)/points) //Distribute points logarithmically over required range
Duplicate/o Omega,Amplitude, Angle
make /o /c /n=(points) jw
jw=cmplx(0,omega) //jw is the complex wave with omega as its imag part
Calculate("",0,"","") //Calculate the current function and set vertical axis
End
Function Calculate(ctrlName,varNum,varStr,varName) : SetVariableControl
String ctrlName
Variable varNum
String varStr
String varName
SVAR gain, p1,p2,p3,p4,z1,z2,z3,z4
string SFunc=gain //Make string expression for execution. This allows leaving out terms for non existent singularities
if (strlen(z1)>0)
SFunc+="*(jw+"+z1+")"
Endif
if (strlen(z2)>0)
SFunc+="*(jw+"+z2+")"
Endif
if (strlen(z3)>0)
SFunc+="*(jw+"+z3+")"
Endif
if (strlen(p1)>0)
SFunc+="/(jw+"+p1+")"
Endif
if (strlen(p2)>0)
SFunc+="/(jw+"+p2+")"
Endif
if (strlen(p3)>0)
SFunc+="/(jw+"+p3+")"
Endif
Print SFunc
String Command="Amplitude=cabs("+sFunc+"); Duplicate /o Amplitude, Angle; Angle=360/2/pi*atan2(imag("+sFunc+"),real("+sfunc+"))"
execute Command
SetVertAx() //Set the vertical axis
End
Function SetVertAx() //Sets vertical scale to match decades to horizontal
NVAR MinFreqP, MaxFreqP
WAVE Amplitude
GetWindow Amplitude psize // reads plot area in points to get length of x axis
SetAxis left waveMax(Amplitude)/(10^((MaxFreqP-MinFreqP)/(V_right-V_Left)*(V_bottom-V_top))),waveMax(Amplitude)
End