|
本帖最后由 fanxiaoduo 于 2019-11-25 10:42 编辑
记得有篇报道中,有个记者说到大功率激光要从传统光学来,初步想法有了,等元器件回来,立即安装调试,资金不够了,先合成10个10w的激光源100w,实验如能成功,后期筹集资金合成500-1000个,合成10*1000=10kw 激光,初步具有武器功能2公里内能击穿0.5cm钢板。


激光模块已做好,频率调整模块已做好,Arduino2560对主路激光光电探测后,对每个激光源进行微分计算后输出TTL电平动态调整频率、相位,使每个激光源之间发生光的相干。
不用掺镱光纤可以吗?,我的激光源不在900nm
今晚就把主干光纤上光电探测出的信号用号傅里叶变化分出 两路 信号的Arduiono 类C 代码写好。
1、方法1 直接探测分路上激光源的相位,用同一激光源发射,反射两次,得到相位差后取值后通过微分计算出相位,实现方法简单,分别得到每一路激光源的相位,再调整每路相位,使两路激光相干。
2、方法2 CCD探测出的信号用号傅里叶变化分出 两路 信号的相位,再调整每路相位,使两路激光相干,实现方法稍难。
3、方法3 CCD探测出的信号,生成图片,直接对图片 光相干 分析,调整每路相位,使光相干最强,实现方法简单。
-------------------你眼中看似落叶纷飞变化无常的世界,实际只是躺在上帝怀中一份早已谱好的乐章。
人类世界现在还不存在复杂的量子计算机和生物计算机,哈哈哈
资料稍后上传.....
平台已搭建好,可以分别做上述3种实验,看那种方法的相干效果好,下周我请了年休,封闭开发一周....,稍后上传代码和资料

所有的技术问题都要解决了,我马上要见到光明了,还记得《珊瑚岛上的激光》这部科幻片吗?


自行研发的激光相位调整高频震荡装置见下图(说明,还有两种装置1、电磁直接震荡光纤,2、压电陶瓷震荡光纤,稍后实验后上图。)

今晚在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 );
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-');
|
|