IgorPro實(shí)用教程代碼

1、Igor Pro 實(shí)用教程代碼第 1 章Variable v1,v2String s1=”hello ”,s2Make/N=100 data1Print 2+sqrt(3)+53+60.3Print sin(1),ln(3),exp(2)Make/O/N=100 data1Setscale/I x,0,2*pi,data1Data1=sin(x)Fft data1Make/O/N=100 data1Data1=x+gnoise(1)Curvefit line data1/DDisplay data1Edit data1Display data1ModifyGraph mode(data1)=3

2、,marker(data1)=19Sin(1) /errorPrint sin(1)Make/O/N=100 dataSetscale/I x,0,1,dataData=xEdit data /只顯示 data內(nèi)容Edit data.xy /顯示 data內(nèi)容和 x 坐標(biāo)信息Make/O/N=200 sinx;Setscale/I x,0,2*pi,sinx;sinx=sin(x);Display sinx;Make/O/N=200/C cmpwave;Make/O/T/N=200 twave;Setdatafolder root;/設(shè)定當(dāng)前目錄Newdatafolder tmp;/創(chuàng)建一個(gè)名

3、為tmp 的目錄Make/O/N=200 root:tmp:cmpwave;setscale/P x,0,2*pi/199,”s”,sinxsetscale/I x,0,2*pi,”s”,sinxDuplicate/O oldwave,newwave;dupplicate/O/R=(x1,x2) oldwave,newwave;dupplicate/O/R=p1,p2 oldwave,newwave;Make/O data=1,2,3Make/O dataSetscale/P x,0,0.001,dataVariable tau=0.01Data=exp(-x/0.01)Make/O data

4、1Data1=datapVariable v1; /申明一個(gè)普通變量Variable/C v1; /申明一個(gè)復(fù)數(shù)變量Variable/G v1; /申明一個(gè)全局變量String str; /申明一個(gè)字符串變量Function func()String curr=getdatafolder(1)Setdatafloder mydestfd /mydestfd/some thing to be doneSetdatafolder curr/恢復(fù)當(dāng)前目錄/ 獲取當(dāng)前目錄應(yīng)該存在，存放了要處理的數(shù)據(jù)EndIntegrate1D(UserFunctionName, min_x, max_x , opti

5、ons , count )Integrate /DIM = d /METH=m /P/TtypeFlags yWaveA /X = xWaveA /D = destWaveA , yWaveB /X = xWaveB /D = destWaveB , .第 2 章make/N=100/O gaussfun,lorfunsetscale/I x,-1,1,gaussfun,lorfungaussfun=exp(-x*x/0.01)lorfun=1/(x*x)+0.04)make/N=10 XY_Y=4,6,8,4,6,7,9,1,4,3make/N=10 XY_X=3,5,1,4,7,5,8,6

6、,9,2Display gaussfun/ 繪制 gaussfun 曲線DisplayXY_Y vs XY_X as “mygraph”/ 繪制 XY 曲線并且設(shè)定 graph 標(biāo)題為 mygraph Display/L=newLeft/B=NewBottomgaussfun as "mygraph"/ 以新建自由坐標(biāo)軸 newLeft 和 NewBottom 繪制 gaussfun 曲線。Preferences 0;Display/L=newLeft/B=NewBottom gaussfun as "mygraph"Graph1Style();Pref

7、erences 1;/ 以新建坐標(biāo)軸 newLeft 和 NewBottom 為坐標(biāo)軸繪制 gaussfun 曲線，并采用預(yù)定義樣式 Grap h1Style 。AppendToGraph/L=newLeft/B=newBottomlorfun/ 向當(dāng)前 graph 添加 lorfun，坐標(biāo)軸為 newLeft 和 newBottom。AppendToGraph/W=Graph0 lorfun/向名為“ Graph0 ”的窗口添加 lorfun曲線2.2.2 Make/N=50 sinx,cosx;/ 創(chuàng)建兩個(gè) waveSetscale/I x,0,2*pi,sinx,cosx;/設(shè)置兩個(gè) w

8、ave 的 x坐標(biāo)為 0 到 2 Sinx=sin(x);/設(shè)置 wave sinx為它的橫坐標(biāo)的正弦Cosx=cos(x);/設(shè)置 wavecosx為它的橫坐標(biāo)的余弦Make/N=(50,2) arrowinfo;/ 創(chuàng)建 Arrowinfo信息 waveSetscale/I x,0,2*pi,arrowinfo;/設(shè)置 Arrowinfo wave的 x 坐標(biāo)為 0到 2 ，代表一圈Arrowinfo0=10;/設(shè)置所有的箭頭長(zhǎng)度都是10 ；Arrowinfo1=x;/設(shè)置隨著角度的增加箭頭偏轉(zhuǎn)的方向Display sinxvscosx;/顯示 sinx，并且以 cosx為橫坐標(biāo)Modif

9、yGrapharrow Marker(sinx)=arrowinfo,1,5,0.5,1/ 設(shè)置 Marker 為 Arrow 。這一步可以也可以通過(guò)曲線/ 外觀設(shè)置對(duì)話框設(shè)置箭頭 marker 完成，注意 Arrow/Info要選擇 arrowinfoMake/N=100 datax,datay,datazdatax=enoise(2); datay=enoise(2); dataz=exp(-(datax2+datay2)Display datay vs datax; ModifyGraph mode=3,marker=8ModifyGraph zmrkSize(datay)=dataz,

10、*,*,1,10make/O datasetscale/P x 0,1,datadata=exp(-x/60)display dataModifyGraph manTick(bottom)=0,30,0,0,manMinor(bottom)=3,2 make/O InverseTemperature= 30,20,14.2857,10,5,3.0303,2.22222,1.25 make/O Mobility= 0.211521,0.451599,0.612956,0.691259,0.886406,0.893136,0.921083,1Display Mobility vs InverseT

11、emperature/ 繪制曲線，使用默認(rèn)坐標(biāo)軸ModifyGraph mode=3,marker=8make/N=5/T TickLable=“20”, ”30”, ”50”, ”100”, ”400”/ 刻度線標(biāo)簽make/N=5 TickPositionTickPosition=450/str2num(TickLable)ModifyGraph userticks(bottom)=TickPosition,TickLablemake/O gaussdatasetscale/I x,-1,1,gaussdatagaussdata=gauss(x,0,0.2)display gaussdat

12、aModifyGraph zColor(gaussdata)=gaussdata,*,*,Rainbow256,0ModifyGraph mode=4,marker=19ModifyGraph mskip=5TextBox/C/N=text0/D=2,2,0/A=MC "JCZ15F'Times New Roman'Gauss Funct ionreS(-xS2MZ14S/2F'symbol'sF'Times New Roman'S2MZ14S)"Twice K0 is K0*2"Twice K0 is %g, an

13、d today is %s", 2*K0, date() “x value=%g ”,tagval(2) /輸出 x 的大小，每一個(gè)參數(shù)表示一種預(yù)定義類型 “mean value=%g”,mean(tagwaveref(),-inf,+inf /輸出引用 wave 的平均值， tagwaveref返回 wave 的引用。make/O/N=(128,128) gaussdata=gauss(x,64,10,y,64,10)function xyzcontour()make/O/N=(100,3) xyzwavevariable ifor(i=0;i<100;i+=1)xyzwav

14、ei0=gnoise(1)xyzwavei1=gnoise(1)xyzwavei2=gnoise(1)endfordisplay;appendxyzcontour xyzwaveendMake/O/N=(200,30) fdfun;Setscale/I x,-1,0.2,fdfun;fdfun=1/(exp(x*1.6*10000/1.38/30)+1);duplicate/O fdfun,lorfun1,lorfun2;lorfun1=0.005/(x+0.1)2+0.01);lorfun2=0.008/(x+0.3)2+0.01);fdfun=fdfun+lorfun1+lorfun2ne

15、wwaterfall fdfun;/繪制 waterfall圖Newwaterfall/Host=graphname wavenameDuplicate/O fdfun fdcolor;Fdcolor=y;ModifyGraph zColor(fdfun)=fdcolor,*,*,Rainbow256,0;Modifywaterfall angle=45Make/N=(100,100) gaussfun;Setscale/I x,-0.5,0.5,gaussfun;Gaussfun=exp(-x*x/0.04-y*y/0.04);make/O/N=(100,100) datasetscale/

16、I x,-pi,pi,datasetscale/I x,-pi,pi,datadata=sin(x)*cos(y)ModifyObject=surface0 property= frontColor,0.250004,0.make/O/N=(20,20,3) dspheresetscale/I x,0,pi,dspheresetscale/I y,0,2*pi,dspheredsphere0=1*sin(x)*cos(y)dsphere1=1*sin(x)*sin(y)dsphere2=1*cos(x)NewGizmo/N=Gizmo0/T="Gizmo0"ModifyGi

17、zmo startRecMacroAppendToGizmo Surface=root:dsphere,name=surface0ModifyGizmo ModifyObject=surface0 property= surfaceColorType,1ModifyGizmo ModifyObject=surface0 property= fillMode,3ModifyGizmo ModifyObject=surface0 property= srcMode,4ModifyGizmo996109,0.250004,1ModifyGizmo ModifyObject=surface0 prop

18、erty= backColor,0.996109,0.664073,0,1ModifyGizmo setDisplayList=0, object=surface0ModifyGizmo SETQUATERNION=0.488416,0.302182,0.455614,0.680110Make/O/N=(20,3) data=gnoise(5)data2=2*datap0-3*datap1+datap02+gnoise(0.05)make/N=(100,100)/O datasetscale/I x,-2,2,datasetscale/I y,-2,2,datadata=gauss(x,-1,

19、0.5,y,-1,0.5)+gauss(x,-1,0.5,y,1,0.5)data=data+gauss(x,1,0.5,y,-1,0.5)+gauss(x,1,0.5,y,1,0.5)第 3 章Make/O/N=100 trialwave;Setscale/I x,-1,3,trialwave;Trialwave=exp(-0.1*x)+gnoise(0.02);Display trialwave;ModifyGraph mode=3;function fitfun(coef,x):Fitfuncwave coefvariable x<Exressions>endmake/O c

20、oef=a1,a2,a3,funcfit fitfun,coef,ydata/Dmake/O LorFsetscale/I x,-10,10,LorFsetrandomseed 0.5LorF=1/(x*x+1)+gnoise(0.05)display LorFModifyGraph mode=3,marker=8Make/O/N=(200,200) w2dModifyGraph lsize(fit_LorF)=2,rgb(fit_LorF)=(0,0,65280)Make wf(x)=a*x*x+b*x+cMake/O/N=200 trialdataSetscale/P x,-2,0.02,

21、trialdataSetrandomseed 0.5Trialdata=x*x-2*x-1+gnoise(0.2)Display trialdataModifyGraph mode=3,marker=8,rgb=(0,0,65280)Function f(w,x):FitFuncWave wVariable xEndFunction f(w,x,y):FitFuncWave wVariable xVariable y/ function bodyEndCurveFit flags fitType, kwCWave=coefWaveName , waveName flag parameters

22、CurveFit fitType wavenameMake/O/N=20 trialdataSetscale/I x,0,10,trialdataSetrandomseed 0.5Trialdata=exp(-4*(x-5)2)Display trialdataModifygraph mode=3,marker=8Variable V_FitOptions=8/control fittingCurvefit gauss trialdataEdit M_iteratesK0=1;K1=2;K2=3;K3=4CurveFit/G gauss trialdataMake/O/N=200 triald

23、ataSetscale/I x 0,5,trialdataSetrandomseed 0.5Trialdata=exp(-(x-2)/2)+gnoise(0.1)Display trialdataModifygraph mode=3,marker=8Curvefit/TBOX=768 exp_Xoffset trialdata /DCurvefit/K=2 /TBOX=768 exp_Xoffset trialdata /DSetAxis left 0,4SetAxis bottom -1,8Curvefit/X/TBOX=896 exp_Xoffset trialdata /Dmake/O/

24、N=20 trialdatasetscale/I x,0,2,trialdatasetrandomseed 0.5trialdata=exp(-x)+gnoise(0.1)display trialdatamodifygraph mode=3,marker=8Curvefit/TBOX=768 exp trialdata /D/F=0.95,6,ErrorBarMake/O/N=200 trialdataSetscale/I x,0,2,trialdataSetrandomseed 0.5Trialdata=2*exp(-x)+gnoise(0.3)Display trialdataModif

25、ygraph mode=3,marker=8Curvefit/W=0 line trialdata /D/RMake/O/N=200 data1Setscale/I x,0,2*pi,data1Setrandomseed 0.5Data1=2*sin(x)+gnoise(0.2)Display data1Modifygraph mode=3,marker=8K1+K2 > 3K3/5 < 2*K1K1*K1 < 2。-100010000-100010000-100010000-10001K0,K1,K2,K31,1,0,3,0,2,1,1Edit M_FitConstrain

26、t,W_FitConstraintwaveFunction FitEllipse(w,x,y):FitFuncWave wVariable xVariable y/w0=a/w1=b/w2=x0/w3=y0/equition:x2/a2+y/b2=1Return (x-w2)2/w02+(y-w3)2/w12-1endMake/O/N=20 theta,ellipseX,ellipseYTheta=2*pi/20*pellipseX=2*cos(theta)+2ellipseY=3*sin(theta)+1Display ellipseY vs ellipseXModifygraph mode

27、=3, marker=8Modifygraph width= perUnit,72,bottom,height=perUnit,72,leftModifygraph width= perUnit,72,bottom,height=perUnit,72,leftSetrandomseed 0.5ellipseX += gnoise(.3)ellipseY += gnoise(.3)Duplicate/O ellipseY, ellipseYFit, ellipseXFitMake/D/O ellipseCoefs=3,2,1,2 / a, b, x0, y0FuncFit/ODR=3 FitEl

28、lipse, ellipseCoefs /X=ellipseX, ellipseY/XD=ellipse XFit,ellipseYFitAppendToGraph ellipseYFit vs ellipseXFitMake/O/N=(100,100) elllipseContourSetScale/I x -3,4.5,elllipseContourSetScale/I y -3, 5, elllipseContourelllipseContour = FitEllipse(ellipseCoefs, x, y)AppendMatrixContour elllipseContourModi

29、fyContour elllipseContour labels=0,autoLevels=*,*,0,moreLevels=0,mo reLevels=0ModifyContour elllipseContour rgbLines=(0,0,0)Make/O/N=20 data1Setscale/I x, -2,2,data1Setrandomseed 0.5Data1=abs(x)+gnoise(0.1)Display data1Modifygraph mode=3,marker=8Function myfit(w,x):FitFuncWave wVariable xIf(x>=0)

30、Return w0+w1*xEndifIf(x<0)Return w2-w3*xEndifEndMake/O/N=4 coCo=0,1,0,1FuncFit/W=0/NTHR=0 myfit,co,data1 /DFunction myfitfunc(pw,yw,xw):FitFuncWave pw,yw,xwYw=<expression involving pw and xw>Endduplicate yw tmp<do interpolation or extrapolation work>yw=tmp(x)Make/O/N=200 data1Setscale

31、/I x,-2,2,data1Setrandomseed 0.5Data1=1/(x*x+0.04)+gnoise(0.8)Display data1Modifygraph mode=3,marker=8Proc proc1()Make/O/N=101 gsfunSetscale/I x,-0.5,0.5,gsfunVariable width=0.2/2/ln(2)Gsfun=exp(-x*x/width2)Variable v=sum(gsfun,-inf,inf)Gsfun=gsfun/vDuplicate/O data1 data2Convolve/A gsfun data1Appen

32、dtograph data2EndFunction myfunc(pw,yw,xw):FitFuncWave pw,yw,xwWave w=gsfunYw=pw0+pw1/(xw-pw2)2+pw32)Convolve/A w ywEndMake/O/N=4 coCo=0,10,0,0.2FuncFit/W=0/NTHR=0 myfunc co data1 /DDuplicate/O data2 data3Data3=co0+co1/(x-co2)2+co32)Display data2Modifygraph mode=3,marker=8Appendtograph data3Legend/C

33、/N=text0/A=MC/ 普通擬合函數(shù)結(jié)構(gòu)體變量Structure myBasicFitStructWave coefwVariable xEndStructure/all at once擬合函數(shù)結(jié)構(gòu)變量Structure myAllAtOnceFitStructWave coefwWave ywWave xwEndStructurefunction myfitfunc(s):Fitfuncstruct mystruct &s /mystruct為預(yù)先定義的結(jié)構(gòu)體類型<do something>endMake/D/O/N=100 expData,expDataXexpD

34、ataX = enoise(0.5)+100.5expData = 1.5+2*exp(-(expDataX-100)/0.2) + gnoise(.05)Display expData vs expDataXModifyGraph mode=3,marker=8/ The structure definition Structure expFitStructWave coefw / required coefficient wave Variable x / required X value input Variable x0 / constantEndStructure/ The fitt

35、ing functionFunction fitExpUsingStruct(s) : FitFuncStruct expFitStruct &sreturn s.coefw0 + s.coefw1*exp(-(s.x-s.x0)/s.coefw2)End/ The driver function that calls FuncFit: Function expStructFitDriver(pw, yw, xw, xOff)Wave pw /coefficient wave Wave yw /dataWave xw /x value Variable xOff Variable do

36、ODR STRUCT expFitStruct fs fs.x0 = xOffFuncFit fitExpUsingStruct, pw, yw /X=xw /D /STRC=fs print pwWave W_sigma print W_sigmaEndMake/D/O expStructCoefs = 1.5, 2, .2expStructFitDriver(expStructCoefs, expData, expDataX, 100)Structure expFitStructWave coefw /擬合參數(shù) wave ，必須Variable x/x值，必須STRUCT WMFitInf

37、oStruct fi / WMFitInfoStruct結(jié)構(gòu)成員，可選Variable x0 /常量，可選EndStructure/ 普通擬合函數(shù)：Structure My2DFitStructWave coefwVariable x2EndStructure/All at once擬合函數(shù)：Structure My2DAllAtOnceFitStructWave coefwWave ywWave xw2EndStructureFunction FitManyGaussian(w, x) : FitFuncWAVE wVariable xVariable returnValue = w0Var

38、iable iVariable numPeaks = floor(numpnts(w)-1)/3)Variable cfifor (i = 0; i < numPeaks; i += 1)cfi = 3*i+1returnValue += wcfi*exp(-(x-wcfi+1)/wcfi+2)2)endforreturn returnValueEndSetscale/I x,0,10,line1;Line1=x+enoise(0.1);Display line1;ModifyGraph marker=16;Display newtonring vs ringnum;ModifyGrap

39、h mode=3,marker=16;Display resistance vs temperature第 4 章Make/O/N=100 wSetscale/I x,0,2*pi,ww=sin(x)Print w(0.1),sin(0.1)Duplicate w1 wdestSetscale/I x,x1,x2,wdestWdest=w1(x)+w1(2*x0-x)Make/O data=3,2,5,3,1,5,8,9Make/O xw=1,2.1,3,3.9,5,5.8,7,8Print interp(2,xw,data)Make/O/N=10 xData, yData/ Make sou

40、rce dataxData = p; yData = 3 + 4*xData + gnoise(2)/ Put sample data in themDisplay yData vs xData/ Make a graphModify mode=2, lsize=3/ Display source data as dotsInterpolate2/N=200/T=2/E=2/Y=ydata_CS xData, ydataAppendToGraph yData_CS; Modify rgb(yData_CS)=(0,0,65535)Redimension/N=500 xData, yDataxD

41、ata = p/50; yData = 10*sin(xData) + gnoise(1.0)Modify lsize(yData)=1Interpolate2/N=200/T=2/E=2/Y=ydata_CS xdata, ydataAppendToGraph yData_CS; Modify rgb(yData_CS)=(0,0,65535)Make/N=(100,100)/O ydataSetscale/I x,-1,1,ydataSetscale/I y,-1,1,ydataYdata=x*x+y*yMake/O/N=(10,20,30) ddd=gnoise(10)Print int

42、erp3d(ddd,1,0,0)Make/O/N=(10,4) ddd=gnoise(10)Triangulate3D/OUT=1 dddPrint interp3d(ddd,1,0,0,M_3DVertexList)Interp3DPath 3dWave tripletPathWaveYdata_new=interp2d(ydata,M-1 (x,y)Make/O/N=100 xData = .01*x + gnoise(.01)Make/O/N=100 yData = 1.5 + 5*exp(-(xData-.5)/.1)2)Display yData vs xDataSort xdata

43、,xdata,ydataDuplicate/O yData, wDataSetScale/I x 0, 1, wDatawData = interp(x, xData, yData)AppendToGraph wDatainterp(x1, xwaveName, ywaveName )Interpolate2/T=2/N=200/E=2/Y=yData_CS xData, yDataMake/n=(tot,4) a3ddataMake/N=(n1,n2,n3) new3d/set new3d x,y and z scale/here needs some loops to get x0,y0,

44、z0/loop beginsTransformXYZ (x1,y1,z1,x0,y0,z0)A3ddatai0=x1A3ddatai1=y1A3ddatai2=z1A3ddatai3=v /v表示（ x0,y0,z0）處的測(cè)量值i=i+1/loop endsTriangulate3D a3ddataWave M_3DVertexListNew3d=interp3d(srcwave,x,y,z, M_3DVertexList)ImageInterpolate/S=-3,0.01,1,-2,0.01,2 Voronoi myTripletWaveImageInterpolate/STW/S=0,1

45、,0,1,1,1 Voronoi myTripletWaveImageInterpolate/PTW=W_TriangulationData/S=-3,0.02,1,-2,0.02,2 Voronoimy1DZWaveMy1DZWave= myTripletWave2pContourZ(graphNameStr, contourNameStr, instance, x, y)空字符串，can use ContourZ to convert XYZ triplet data into gridded data./ Make example XYZ triplet contour data Mak

46、e/O/N=50 wx,wy,wzwx= enoise(2)wy= enoise(2)wz= exp(-(wxp*wxp + wyp*wyp)/ x = -2 to 2/ y = -2 to 2/ XY gaussian, z= 0 to 1/ ContourZ requires a displayed contour data set Display; AppendXYZContour wz vs wx,wy;DelayUpdate ModifyContour wz autoLevels=*,*,0ModifyContour wz xymarkers=1/ show the X and Y

47、locations/ Set the null (out-of-XY domain) valueModifyContour wz nullValue=NaN/ default is min(wz) - 1/ Convert to grid: Make matrix that spans X and YMake/O/N=(30,30) matrixSetScale/I x, -2, 2, "", matrixSetScale/I y, -2, 2, "", matrixmatrix= ContourZ("","wz"

48、,0,x,y)/ or = ContourZ("","",0,x,y)AppendImage matrixBinarySearch(waveName,val)Make/O data = 1, 2, 3.3, 4.9/ Monotonic increasingPrint BinarySearch(data,3)BinarySearchInterp(waveName,val)Make/O data = 1, 2, 3.3, 4.9/ Monotonic increasingPrint BinarySearchInterp(data,3)Print data1

49、.76923duplicate w,w1w1=leftx(w)+deltax(w)*pvariable x0=interp(y0,w,w1)FindLevel /B=box /EDGE=e /P/Q/R=(startX,endX ) waveName, levelMake/O/N=100 data1Setscale/P x,0,0.01,data1Data1=2*sin(2*pi*3*x)+gnoise(0.1)Display data1Findlevel/B=5/R=(0.6,0.8) data1,1Make/O/N=2 ylevelSetscale/I x,0,1,ylevelYlevel=1Make/N=2/O xcross_x,xcross_y=-2,2Xcross_x=V_LevelXAppendtograph ylevelAppendtograph xcross