基于matlab的圆偏振光/两相交流电模拟-导出gif

最新推荐文章于 2023-06-30 22:15:18 发布

置顶范范TT西西

最新推荐文章于 2023-06-30 22:15:18 发布

阅读量1.8k

点赞数 8

分类专栏： Wave 文章标签：偏振 matlab gif

本文链接：https://blog.csdn.net/qq_24694761/article/details/103578224

版权

Wave 专栏收录该内容

0 篇文章 1 订阅

订阅专栏

光学和电学具有诸多相似的地方，本文针对两个平面偏振光的合成进行讨论。

从波动光学中我们可以得到：

相同频率、相同振幅、相位也完全相同的两电磁波，合成的依然是平面偏振光；
相同频率、相同振幅、相位差为pi/2的两电磁波，合成的是右旋圆偏振光；
相同频率、相同振幅、相位差为-pi/2的两电磁波，合成的是左旋圆偏振光；
相同频率、相同振幅、相位差为0~pi/2的两电磁波，合成的是椭圆偏振光;

上述四条对于掌握圆二色谱的基本原理有一定帮助。

同时我们知道：

空间相差90度，相位相差90度（利用电容）的两相交流电可以合成一个旋转电磁场；

本matlab程序包含两个子函数：Polarization.m和OutGif.m，前者作图并输出gif文件，后者用于设定gif的参数。

可以在这里下载

%% 偏振光演示
function [] = Polarization(A,f,Tend,Tint,Angle)
if nargin == 0
    A = [1,2];    % 幅值，可改参数哈
    f = [1,1];    % 频率，可改参数哈
    Tint = 0.01;  % 每帧的时长，可改参数哈
    Tend = 1;     % 运行一秒，可改参数哈
    Angle = pi/2;   % 相位，可改参数哈，注意单位
end
%% 开始作图
figure('color',[1,1,1]);
w = 2*pi*f;
iter = 0;
x = [];
y = [];
for t = 0:Tint:Tend
    iter = iter+1;
    % 线偏振光
    X = A(1)*sin(w(1)*t);
    plot([0,X],[0,0],'r','LineWIdth',1.5);
    hold on;
    plot([X,X],[0,0],'ro','LineWIdth',1.5,'MarkerFaceColor','r');
    % 另外一个线偏振光
    Y = A(2)*sin(w(2)*t+Angle);
    plot([0,0],[0,Y],'b','LineWIdth',1.5);
    hold on;
    plot([0,0],[Y,Y],'bo','LineWIdth',1.5,'MarkerFaceColor','b');
    % 合成
    plot([0,X],[0,Y],'k','LineWIdth',1.5);
    hold on;
    plot([0,X],[0,Y],'ko','LineWIdth',1.5,'MarkerFaceColor','k');
    % 轨迹
    x = [x X];
    y = [y Y];
    plot(x,y,'k--','LineWIdth',0.5);
    % 其他作图设置
    hold off
    axis equal
    xlim([-A(1),A(1)]);
    ylim([-A(2),A(2)]);
    axis off
    M = getframe();
    OutGif(M,'POLAR',iter,Tint)
end

end

%% 输出gif文件(简易版)
function Status = OutGif(GetFrameOutput,FileName,NumFrame,LoopTime,Path)
if nargin == 3
    LoopTime = 0.3;
    Path = 'C:\Users\ZLY\Desktop';
elseif nargin == 4
    Path = 'C:\Users\ZLY\Desktop';
end
% Status = PathCheck(Path);
imind = frame2im(GetFrameOutput);
[imind,cm] = rgb2ind(imind,256);
T = datestr(now,1);
if NumFrame == 1
    imwrite(imind,cm,strcat(Path,'\',[T,FileName],'.gif'),'gif','LoopCount',inf,'DelayTime',LoopTime);
else
    imwrite(imind,cm,strcat(Path,'\',[T,FileName],'.gif'),'gif','WriteMode','append','DelayTime',LoopTime);
end
end

运行结果：

不同参数对比
A = [1,2]; % 幅值 f = [1,1]; % 频率 Tint = 0.01; % 每帧的时长 Tend = 2; % 运行时间 Angle = pi/2; % 相位
A = [1,1]; % 幅值 f = [1,1]; % 频率 Tint = 0.01; % 每帧的时长 Tend = 2; % 运行时间 Angle = pi/2; % 相位
A = [1,1]; % 幅值 f = [1,1]; % 频率 Tint = 0.01; % 每帧的时长 Tend = 2; % 运行时间 Angle = -pi/2; % 相位
A = [1,1]; % 幅值 f = [1,1]; % 频率 Tint = 0.01; % 每帧的时长 Tend = 2; % 运行时间 Angle = pi/3; % 相位
A = [1,1]; % 幅值 f = [1,2]; % 频率 Tint = 0.01; % 每帧的时长 Tend = 2; % 运行时间 Angle = 0; % 相位
A = [1,1]; % 幅值 f = [1,2]; % 频率 Tint = 0.01; % 每帧的时长 Tend = 2; % 运行时间 Angle = pi/3; % 相位
A = [1,1]; % 幅值 f = [1,2]; % 频率 Tint = 0.01; % 每帧的时长 Tend = 2; % 运行时间 Angle = pi/2; % 相位