平面波极化

泸州月

已于 2022-05-27 16:21:06 修改

阅读量3.3k

点赞数 2

分类专栏：计算电磁文章标签：平面波极化线极化圆极化椭圆极化

于 2020-12-24 19:18:32 首次发布

本文为博主原创文章，未经博主允许不得转载。

本文链接：https://blog.csdn.net/qq_17119267/article/details/111646918

版权

计算电磁专栏收录该内容

2 篇文章 2 订阅

订阅专栏

平面波的极化

对于平面波，设传播方向为矢量 $\overrightarrow{z}$ ,电场E垂直与波的传播方向，极化状态按照E的矢量端点轨迹定义，根据场的矢端轨迹，分为线极化、圆极化、椭圆极化三类。根据情况我们对三种极化状态进行简单分析。

考虑到以下几种特殊情况的极化状态，如图3-2-1所示：

在这里插入图片描述

假设 $\overrightarrow{E}(z,t)=\overrightarrow{e}_x \overrightarrow{E}_{x0}(\omega t- \varphi_1)+\overrightarrow{e}_y \overrightarrow{E}_{y0}(\omega t- \varphi_2)$ ，其极化类型取决于 $E_{x0},E_{y0}$ 以及 $\varphi_1,\varphi_2$ 。

情形一：当 $\varphi_2-\varphi_1=2m\pi$ 或者 $\varphi_2-\varphi_1=(2m+1)\pi$ 为线极化，如 3-2-1（a）和（b）所示。

情形二：若 $E_{x}=E_{y}$ 为圆极化， $\varphi_2-\varphi_1=\pm\pi/2$ ，其中，"+“为右旋圆极化，”-"为左旋圆极化，如图3-2-1（c）和（d）所示。

情形二：若 $E_{x}\ne E_{y}$ 为圆极化， $\varphi_2-\varphi_1=\pm\pi/2$ ，其中，"+“为右旋椭圆极化，”-"为左旋椭圆极化，如图3-2-1(e)和(f)所示。

Matlab极化波实现:

线性极化波

%% 沿x方向phase =0传输时的线性极化波
x=(0:0.3:30);                   	%传输距离
beta=0.8;
Eym=1;
Ezm=2;
l=zeros(size(x));                  %设置l为与x尺寸相同的零向量
t=0;                            	%时间变量
for i=1:1500                     	%帧数
    omega=2*pi;
    Ey=Eym*cos(omega*t-beta*x);   %电场横向y分量
    Ez=Ezm*cos(omega*t-beta*x);   %电场横向z分量
    quiver3(x,l,l,l,Ey,Ez); 
    %以(x,0,0)为起点画出传输方向上每一点的电场矢量图
    axis([0,30,-4,4,-4,4]); view(20,40); 	%观察范围
    pause(0.01);                   	%矢量图显示0.01s；
    mov(i)=getframe(gcf);         %捕捉当前图像作为一帧；
    t=t+0.01;                     	%时间变量变化微小量
end;
hold off; 
VideoWriter('mov.avi','Motion JPEG AVI');    	
%和前面的循环配合，实现动画输出为avi文件

在这里插入图片描述

左旋圆极化波

%% 沿x方向 phase =-1/2传输时的左旋圆极化波
x=(0:0.3:30);                   	%传输距离
l=zeros(size(x)); 
t=0;                          	%时间变量
for i=1:1500                    	%帧数
omega=2*pi;
Ey=cos(omega*t-0.8*x);        	%电场横向y分量
Ez=cos(omega*t-0.8*x+pi/2);    	%电场横向z分量
quiver3(x,l,l,l,Ey,Ez);          	
%以(x,0,0)为起点画出传输方向上每一点的电场矢量图
axis([0,30,-4,4,-4,4]); view(20,40); 	%观察范围
pause(0.01);                   	%矢量图显示0.01s；
mov(i)=getframe(gcf);          	%捕捉当前图像作为一帧；
t=t+0.01;                     	%时间变量变化微小量
end;
hold off;
VideoWriter('mov.avi','Motion JPEG AVI'); 	
%和前面的循环配合，实现动画输出为avi文件

在这里插入图片描述

右旋圆极化波

%% 沿x方向phase =0传输时的线性极化波
x=(0:0.3:30);                   	%传输距离
beta=0.8;%% 沿x方向 phase =-1/2传输时的右旋圆极化波
x=(0:0.3:30);                   	%传输距离
l=zeros(size(x)); 
t=0;                          	%时间变量
for i=1:1500                    	%帧数
		omega=2*pi;
		Ey=cos(omega*t-0.8*x);        	%电场横向y分量
		Ez=cos(omega*t-0.8*x-pi/2);    	%电场横向z分量
		quiver3(x,l,l,l,Ey,Ez);          
		%以(x,0,0)为起点画出传输方向上每一点的电场矢量图
		axis([0,30,-4,4,-4,4]); view(20,40); 	%观察范围
		pause(0.01);                   	%矢量图显示0.01s；
		mov(i)=getframe(gcf);          	%捕捉当前图像作为一帧；
		t=t+0.01;                     	%时间变量变化微小量
end;
hold off;
VideoWriter('mov.avi','Motion JPEG AVI'); 	
%和前面的循环配合，实现动画输出为avi文件

Eym=1;
Ezm=2;
l=zeros(size(x));                  %设置l为与x尺寸相同的零向量
t=0;                            	%时间变量
for i=1:1500                     	%帧数
    omega=2*pi;
    Ey=Eym*cos(omega*t-beta*x);   %电场横向y分量
    Ez=Ezm*cos(omega*t-beta*x);   %电场横向z分量
    quiver3(x,l,l,l,Ey,Ez); 
    %以(x,0,0)为起点画出传输方向上每一点的电场矢量图
    axis([0,30,-4,4,-4,4]); view(20,40); 	%观察范围
    pause(0.01);                   	%矢量图显示0.01s；
    mov(i)=getframe(gcf);         %捕捉当前图像作为一帧；
    t=t+0.01;                     	%时间变量变化微小量
end;
hold off; 
VideoWriter('mov.avi','Motion JPEG AVI');    	
%和前面的循环配合，实现动画输出为avi文件

在这里插入图片描述

均匀平面电磁波的一般表达式

一般表达式：

$\overrightarrow{E}_y= \overrightarrow{E}_{ym}(\omega t- \beta x)$

$\overrightarrow{E}_z= \overrightarrow{E}_{zm}(\omega t- \beta x-\varphi)$

%% J45 可以在均匀平面电磁波的一般表达式
clc; 
clear;	                                 %清屏幕，清内存
t=0;                                	%时间变量t赋初值
x=0;                               	%沿传播方向取一个截面，x=0
Eym=1; Ezm=2; phi=-1*pi/4;      %两个电场分量的设置，修改可以得到不同极化
Scope=1.5*max([Eym Ezm]);	%设置显示范围
omega=2*pi;	%角频率
for i=1:200                         	%循环100次，获取100帧数据
		Ey=Eym*cos(omega*t-0.8*x);      	%电场横向分量Ey
		Ez=Ezm*cos(omega*t-0.8*x-phi);  	%电场横向分量Ez
		quiver(0,0,Ey,Ez,1);            	%以(0,0)为起点绘制电场矢量箭头图
		% plot(Ey,Ez,'*');	%另外一种展示方式
		axis equal;	%等比例显示
		% hold on;	%绘图保持模式打开
		axis([-Scope Scope -Scope Scope]); 	%设置显示范围
		pause(0.01);                    	%矢量图显示0.01s；
		mov(i)=getframe(gcf);               	%捕捉当前图像作为一帧
		t=t+0.01;                       	%时间变量变化微小量
end;
VideoWriter('mov.avi','Motion JPEG AVI');    	
%将所有帧的数据转换为avi文件并保存