ghost imging.zip_finallydnl_ghost_ghost imaging_ghost-imaging_op

%object:****** %function:****** clc;clear all;close all; %array initialization CS=zeros(1,1000); sum_RAD=zeros(1,1000); sum_SAMD=zeros(1,1000); sum_RAS_SAMD=zeros(1,1000); delta_RAD=zeros(1,1000); delta_SAMD=zeros(1,1000); sum_2delta_RAD=zeros(1,1000); sum_2delta_SAMD=zeros(1,1000); mean_sum_RAS_SAMD=zeros(1,1000); mean_sum_2delta_RAD=zeros(1,1000); mean_sum_2delta_SAMD=zeros(1,1000); ghost_image=zeros(1,1000); SNR=zeros(1,1000); %carrier signal for p=1:1000 if p<=100 CS(p)=1; elseif (100<p)&&(p<=200) CS(p)=0; elseif (200<p)&&(p<=300) CS(p)=1; elseif (300<p)&&(p<=400) CS(p)=0; elseif (400<p)&&(p<=500) CS(p)=1; elseif (500<p)&&(p<=600) CS(p)=0; elseif (600<p)&&(p<=700) CS(p)=1; elseif (700<p)&&(p<=800) CS(p)=0; elseif (800<p)&&(p<=900) CS(p)=1; elseif (900<p)&&(p<=1000) CS(p)=0; end end;%carrier signal for n=1:5000 %random light intensity of a laser in 10ns LI=normrnd(0.5,0.01,[1 1000]);%light intensity per 10ps %chect whether the minimum is minus number minLI=min(LI);%(meet the need) maxLI=max(LI);% j=1:1000; hold on; subplot(221); plot(j,LI,'g-'); axis([300 900 0 1]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('primary signal'); %split beam RA=LI; SA=LI; %modulation SAM=SA.*CS; %detecting %detecting the modulated signal %unable to resolve the modulated signal MMS=ones(1,500);%band:5ns SAMD=imfilter(SAM,MMS,'replicate'); %SAMD=SAMD+randn(1,1000); %detecting the reference signal %be quality to resolve the primary signal MMQ=ones(1,5);%band:50ps RAD=imfilter(RA,MMQ,'replicate'); %RAD=RAD+randn(1,1000); %coresponding for m=1:1000 sum_SAMD(m)=sum_SAMD(m)+SAMD(m); sum_RAD(m)=sum_RAD(m)+RAD(m); mean_RAD=sum_RAD./n; max_mean_RAD=max(mean_RAD); mean_SAMD=sum_SAMD./n; max_mean_SAMD=max(mean_SAMD); delta_RAD(m)=RAD(m)-mean_RAD(m); delta_SAMD(m)=SAMD(m)-mean_SAMD(m); sum_RAS_SAMD(m)=sum_RAS_SAMD(m)+(delta_RAD(m)*delta_SAMD(m)); mean_sum_RAS_SAMD(m)=sum_RAS_SAMD(m)./n; sum_2delta_RAD(m)=sum_2delta_RAD(m)+(delta_RAD(m)*delta_RAD(m)); mean_sum_2delta_RAD(m)=sum_2delta_RAD(m)./n; sum_2delta_SAMD(m)=sum_2delta_SAMD(m)+(delta_SAMD(m)*delta_SAMD(m)); mean_sum_2delta_SAMD(m)=sum_2delta_SAMD(m)./n; ghost_image(m)=abs(mean_sum_RAS_SAMD(m))/sqrt(mean_sum_2delta_RAD(m)*mean_sum_2delta_SAMD(m)); SNR(m)=ghost_image(m)*sqrt(n)/(1-(ghost_image(m)*ghost_image(m))); end end %above image show j=1:1000; subplot(221); plot(j,LI,'r-'); axis([300 900 0 1]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('primary signal'); subplot(222); plot(j,CS,'b-'); axis([300 900 0 2]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('carrier signal'); subplot(223); plot(j,SAM,'b-'); axis([300 900 0 2]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('modulated signal of signal arm'); subplot(224); plot(j,SAMD,'g-'); axis([300 900 -10 300]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('detecting modulated signal'); figure; subplot(221); plot(j,RAD,'g-'); axis([300 900 -10 50]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('detecting reference signal'); subplot(222); plot(j,ghost_image,'r-'); axis([300 900 -1 1.5]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('ghost image'); subplot(223); plot(j,CS,'b-'); axis([300 900 0 1.5]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('carrier signal'); subplot(224); plot(j,SNR,'b-'); axis([300 900 -1 1.5]); xlabel('Time index£¨ps£©');ylabel('Itensity'); title('SNR');