function sweepDpTickle(opt, sigAC, sigDC, prbI, prbQ) %% test=1; %% if test prbI = pr.REFL_1I; prbQ = pr.REFL_1Q; end %% driveNames = {'ETMX','ETMY','ITMX','ITMY','PRM','SRM'}; cDrvNames = {'DARM','CARM','MICH','PRC','SRC'}; mDrv = [1,-1,0,0,0,0; %DARM 1,1,0,0,0,0; %CARM 1,-1,1,-1,0,0;%MICH 0,0,0,0,1,0; %PRC 0,0,0,0,0,1]; %SRC %% dp sweep Ndrv=length(cDrvNames); iPrb=[prbI, prbQ]; iDrv=getDriveNumbers(opt,driveNames); th=linspace(0,2*pi,20); sigI = zeros(numel(th), Ndrv); sigQ = zeros(numel(th), Ndrv); ofsI = zeros(numel(th),1); ofsQ = zeros(numel(th),1); [fDC, sigDC, sigAC]=tickle(opt,[],0.0001); sigIc = sigAC(prbI,iDrv)*mDrv'; phaseOrigin = conj(sigIc/abs(sigIc)); for ii=1:numel(th) opt=setDemodPhase(opt, prbI, th(ii)*180/pi); [fDC, sigDC, sigAC]=tickle(opt,[],0.0001); sigIc = sigAC(prbI,iDrv)*mDrv'; sigQc = sigAC(prbQ,iDrv)*mDrv'; sigI(ii,:) = real(sigIc .* phaseOrigin); sigQ(ii,:) = real(sigQc .* phaseOrigin); ofsI(ii) = sigDC(prbI); ofsQ(ii) = sigDC(prbQ); end %% Plot [colors, styles]=pickLineStyle(Ndrv,[1;10],[1]); figure(1) hold on for ii=1:Ndrv h=plot(th, sigI(:,ii)); %legend(h,cDrvNames(ii)); set(h,'Color', colors(ii,:),'LineStyle',styles{ii}, ... 'DisplayName', cDrvNames{ii}); end legend('show') figure(2) hold on h=plot(th, ofsI); legend(h,'Offset');