matlab生成数字散斑

生成数字散斑Matlab代码一采用一阶形函数，要求在水平和竖直方向的位移均以非均匀线性变化的方式进行位移u v采用函数方程形式；代码二采用二阶形函数，同样要求在在水平和竖直方向的位移均以非均匀线性变化的方式进u v行位移采用函数方程形式，对于其余个参数只考虑位移为数值的形式，不将其设置为函数方10程，可在以下程序中改编该程序生成的散斑图位移、均为数值，不是函数方程u v%deformationgradient measurementusing digitalimage/speckle〈散%input forimage sizeNX’NY%=========================================================================斑图大小像素＞prompt={Enter hprozontalxresolution for simulatedimagefpixels]Enter verticalyresolution forsimulated image{pixels}:;dlg_title=Input forimage size;num_lines=l;def={‘512512};创建并打开输入对话框answer=inputdlgprompt,dlg_title,num_lines,def;%输入的水平分辨率的值的NX=str2doublecell2mat answer1,1;%=川「；输入的竖直分辨率的值5”26101|131©“©1120131311562,1%〈散斑数量%input fornumble of specklesS prompt={enter totalnumber ofspeckles;dlg_title=input numberofspeckles;num_lines=l;def={*1200*;创建并打开输入对话框answer=inputdig prompt,dlg_t itle,num_lines,def;%%输入的散斑数量S=str2doublecell2mat answer1,1;%＜散斑大小%input forspeckle sizeaprompt=enter speckle size};dlg_title=input forspecklesize,;num_lines=l;def={＜4,};answer=inputdlg prompt,,dlg_t itle^num_lines,def;输入的散斑大小2doublecell2mat answer1,1;%＜散斑峰值强度%input forpeak intensityof eachspeckle10prompt={enter peakintensityof eachspeckle};dlg_title=input forpeak intensityof eachspeckle;num_lines=l;def={,l,};answer=inputdlgpromptdlg_title^num_lines^def;I0=str2doublcccll2matanswer1^1;＜散斑图位移该位移均为数值，将其改变为方程的形式，使%input fordisplacementUX,UY其产生非均匀线性变化prompt=enter verticalydisplacementfor simulateimagefpixels]:.enter horizontalxdisplacementforsimulated image[pixels]:1};dlg_title=input forsub-pixel diaplacement,;num_lines=l;def={,0,/0,};answer=inputdlgprompt dlg_title num_lines,def;J JUX=str2doublecell2matanswer1^1;UY=str2doublecell2matanswer2^1;▽丫八位移梯度分量》%input fordeformation gradientux,uy,vx,prompt={enter deformationgradientux forsimulated image[pixels]:.enter deformationgradientuy forsimulatedimage[pixels]:

1.,enter deformationgradientvx forsimulated imagefpixels]:

1.enter deformationgradientvx forsimulated imagefpixels]:1};dlg_title=input forsub-pixel displacement;num_lines=l;def={,

0.002,/0,/0,/0,};answer=inputdlg prompt,,dlg_t itle^num_lines,def;ux=str2doublecell2matanswer1^1;uy=str2doublecell2matanswer2^1;vx=str2doublecell2matanswer3^1;vy=str2doublecell2matanswer41;主程序%main program创建一个的全矩阵1=zerosNX,NY;%512*5120%rand*state;产生阶均匀分布的随机数矩阵X=NX*randS,l;%S*10-1Y=NY*randS l;J初始散斑图的生成%generation ofundeformed speckleimage fori=l:NXfor j=l:NYIi j=

10.*pi./

4.*a.A

2.*sumerfi-Xl:end./a-erfi+1-4Xl:end./a.*erfj-Yl:end./a-erfj+l-Yl:end./a;endend二工是将读入的图像工的数据转换为类型的数据A doubleD;%double uint8double是一个计算机函数，功能是实现图像矩阵的归一化操作所谓”归一G=mat2grayA;%mat2gray化”就是使矩阵的每个元素的值都在和之间01生成数据为名称为的格式图像imwriteG^undeformed.tif%G1tifM=[l-ux-uy;-vx1-vy];3=detM;变形后的散斑图生成%generation ofspeckles imageafter deformationbfori=l:NXfor j=l:NYIi,,j=

10.*pi./

4.*a.A2,*sumerf i-i*ux-j*uy-UX-Xl:end./a-erfi+1-i+l*ux-j+l*uy-UX-Xl:end./a.*erfj-j*vy-i*vx-UY-Yl:end./a-erfj+l-j+l*vy-i+l*vx-UY-Yl:end./a;end endB=doubleI;D=mat2grayB;保存变形后的散斑图imwriteCD^deformed.tif%。