光栅莫尔条纹 matlab程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
global dt %dt是鼠标每次移动绕旋转中心转动角度的增量
N = 25; %number of lines
h_f = figure;
h_a = axes('parent',h_f);
set(h_f,...
'windowbuttonupfcn',@stopdragfcn,...
'windowbuttonmotionfcn','') %先将figure的鼠标移动的callback设为空，鼠标（按下后）弹起的callback设为stopdragfcn函数
set(h_a,'buttondownfcn',@startdragfcn) %将axes的鼠标按下的callback设为startdragfcn函数
set(h_a,'dataaspectratio',[1 1 1],...
'xlim',[-2,2],'ylim',[-2,2])
%mouse_down = false;
xy_pre = [NaN,NaN]; %存储鼠标移动前上一次的坐标
rt_pt = [0 0]; %旋转中心
rt_mat = repmat(rt_pt',1,2*N); %计算变换后坐标时要用
t = 0; %t为旋转角度
x0 = repmat(linspace(-1,1,N),2,1); 在图上不动的N条线的x坐标
y0 = repmat([-1;1],1,N);
x1 = x0; %旋转运动的N条线的x坐标
y1 = y0;
l1 = zeros(1,N); %运动线条的句柄
for k = 1:N
line(x0(:,k),y0(:,k),'color',[0 0 0],'linestyle','-',...
'linewidth',2);
l1(k) = line(x1(:,k),y1(:,k),'color',[0 0 1],'linestyle','-',...
'linewidth',2);
end %产生2N条线