function sweepDp(opt, sigAC, sigDC, prbI, prbQ) %% test=0; %% 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 %% I Q Signals iPrb=[prbI, prbQ]; iDrv=getDriveNumbers(opt,driveNames); sigIc = sigAC(prbI,iDrv)*mDrv'; sigQc = sigAC(prbQ,iDrv)*mDrv'; phaseOrigin = conj(sigIc/abs(sigIc)); sigI0 = real(sigIc .* phaseOrigin); sigQ0 = real(sigQc .* phaseOrigin); ofsI0 = sigDC(prbI); ofsQ0 = sigDC(prbQ); %% Plot Ndrv=length(cDrvNames); [colors, styles]=pickLineStyle(Ndrv,[1;10],[1]); th=linspace(0,2*pi,100); sigI = cos(th)' * sigI0 - sin(th)' * sigQ0; sigQ = sin(th)' * sigI0 + cos(th)' * sigQ0; ofsI = cos(th)' * ofsI0 - sin(th)' * ofsQ0; ofsQ = sin(th)' * ofsI0 + cos(th)' * ofsQ0; figure(1) hold on for ii=1:Ndrv h=plot(180*th/pi, sigI(:,ii)); %legend(h,cDrvNames(ii)); set(h,'Color', colors(ii,:),'LineStyle',styles{ii}, ... 'DisplayName', cDrvNames{ii}); end legend('show'); grid on; figure(2) hold on h=plot(180*th/pi, ofsI); legend(h,'Offset'); grid on;