記得有篇報道中,有個記者說到大功率激光要從傳統光學來,初步想法有了,等元器件回來,立即安裝調試,資金不夠了,先合成10個10w的激光源100w,實驗如能成功,後期籌集資金合成500-1000個,合成10*1000=10kw 激光,初步具有武器功能2公里內能擊穿0.5cm鋼板。
激光模塊已做好,頻率調整模塊已做好,Arduino2560對主路激光光電探測後,對每個激光源進行微分計算後輸出TTL電平動態調整頻率、相位,使每個激光源之間發生光的相干。
不用摻鐿光纖可以嗎?,我的激光源不在900nm
今晚就把主幹光纖上光電探測出的信號用號傅里葉變化分出 兩路 信號的Arduiono 類C 代碼寫好。
1、方法1直接探測分路上激光源的相位,用同一激光源不同的頻率發射兩次,取值後計算,實現方法簡單,分別得到每一路激光源的相位,再調整每路相位,使兩路激光相干。
2、方法2CCD探測出的信號用號傅里葉變化分出 兩路 信號的相位,再調整每路相位,使兩路激光相干,實現方法稍難。
3、方法3CD探測出的信號,生成圖片,直接對圖片 光相干 分析,調整每路相位,使光相干最強,實現方法簡單。
-------------------你眼中看似落葉紛飛變化無常的世界,實際只是躺在上帝懷中一份早已譜好的樂章。
人類世界現在還不存在複雜的量子計算機和生物計算機,哈哈哈
資料稍後上傳.....
平臺已搭建好,可以分別做上述3種實驗,看那種方法的相干效果好,下週我請了年休,封閉開發一週....,稍後上傳代碼和資料
所有的技術問題都要解決了,我馬上要見到光明瞭,還記得《珊瑚島上的激光》這部科幻片嗎?
今晚在Matlab或OpenCV上 做模擬測試
代碼如下:
參考1、
hold off
clear all
f=100;%信號頻率Hz
fs=10000;%採樣頻率Hz
N=400;%採樣點數
t=(0:N-1)/fs;%採樣時間s
for j = 175 :225
x=sin(2*pi*f*t);%信號採樣值
D = j;
x1 = x(1:D);
x2 = x(D+1:N);
for i = 1 : D
x3(i) = x1(D+1-i);
end
plot(x2);
hold on
plot(x3);
hold off
pause(0.5);
end
參考2、
MATLAB
方法一:
N=200;dt=0.001;n=1:200;
x=3*sin(2*pi*10*n*dt)+3*sin(2*pi*30*n*dt)+sin(2*pi*40*n*dt)+sin(2*pi*50*n*dt)+sin(2*pi*60*n*dt)+6*sin(2*pi*80*n*dt);%建立時間序列
X=zeros(1,200); %給X一個預先的內存空間,提高運行速度
figure(1),plot(n,x); %畫出時間系列圖像
%傅里葉變換
for k=1:200
for n=1:200
X(k)=X(k)+x(n)*exp(-i*2*pi*n*k/N);
end
end
f=abs(X); %對傅里葉變換後的圖像取正數部分
figure(2),plot(f); %畫出傅里葉變換後的圖像
%濾波
H=ones(200);H(8:14)=0;H(186:193)=0; %建立40、50、60Hz的時間序列所在空間域數值爲0的一維矩陣
for k=1:200
Y(k)=X(k)*H(k);
end
k=1:200;
figure(3);
plot(k,abs(Y)); %畫出濾波之後的時間序列
%傅里葉逆變換
y=zeros(1,200)% 給y一個預先的內存空間,提高運行速度
for n=1:200
for k=1:200
y(n)=Y(k)*exp(i*2*pi*n*k/N)+y(n);
end
y(n)=y(n)/N;
end
figure(4) ;
n=1:200;
plot(n,y,'r-',n,x,'b-');% 畫出原時間序列和濾波後的時間序列,可見濾波之後的時間序列振幅明顯變小
方法二:
N=200;dt=0.001;n=1:200;
x=3*sin(2*pi*10*n*dt)+3*sin(2*pi*30*n*dt)+sin(2*pi*40*n*dt)+sin(2*pi*50*n*dt)+sin(2*pi*60*n*dt)+6*sin(2*pi*80*n*dt);
X=zeros(1,200);
figure(1),plot(n,x);
X=fft(x,N);
f=abs(X);
figure(2),plot(f);
H=ones(200);H(8:14)=0;H(186:194)=0;
%濾波
for k=1:200
Y(k)=X(k)*H(k);
end
k=1:200;
figure(3);
plot(k,abs(Y));
%傅里葉逆變換
y=zeros(1,200);
y=ifft(Y,N) ;
figure(4) ;
n=1:200;
plot(n,y,'r-',n,x,'b-');