Trans_线极化波matlab_线极化转圆极化_

clear all; close all; clc; % load the experimental data of transmission matrix Back_te=load('1.txt'); %Background for nomorlization Back_tm=load('0.txt'); %Background for nomorlization Txx_exp=load('11.txt'); %载入x方向偏振输入，x方向偏振输出文件 Tyy_exp=load('22.txt'); %载入y方向偏振输入，y方向偏振输出文件 Tyx_exp=load('21.txt'); %载入x方向偏振输入，y方向偏振输出文件 Txy_exp=load('12.txt'); %载入y方向偏振输入，x方向偏振输出文件 Txx_dB=Txx_exp(:,2); %透射的幅度 Tyy_dB=Tyy_exp(:,2); %透射的幅度 Tyx_dB=Tyx_exp(:,2); %透射的幅度 Txy_dB=Txy_exp(:,2); %透射的幅度 Txx_phase=Txx_exp(:,3); %透射的相位 Tyy_phase=Tyy_exp(:,3); %透射的相位 Tyx_phase=Tyx_exp(:,3); %透射的相位 Txy_phase=Txy_exp(:,3)+ones(201,1).*50; %透射的相位 Txx_Linear = 10.^(Txx_dB./20); %由dB化为普通 Txy_Linear = 10.^(Txy_dB./20); %由dB化为普通 Tyx_Linear = 10.^(Tyx_dB./20); %由dB化为普通 Tyy_Linear = 10.^(Tyy_dB./20); %由dB化为普通 Txx = Txx_Linear.*exp(i*Txx_phase*pi./180); %转为复数形式 Txy = Txy_Linear.*exp(i*Txy_phase*pi./180); %转为复数形式 Tyx = Tyx_Linear.*exp(i*Txy_phase*pi./180); %转为复数形式 Tyy = Tyy_Linear.*exp(i*Tyy_phase*pi./180); %转为复数形式 T_RR = 0.5*(Txx+Tyy+1i*(Txy-Tyx)); %t_++ T_RL = 0.5*(Txx-Tyy-1i*(Txy+Tyx)); %t_+- T_LR = 0.5*(Txx-Tyy+1i*(Txy+Tyx)); %t_-+ T_LL = 0.5*(Txx+Tyy-1i*(Txy-Tyx)); %t_-- freq=Txx_exp(:,1); figure plot(freq,abs(T_RR)); hold on plot(freq,abs(T_LL)); hold on plot(freq,abs(T_RL)); hold on plot(freq,abs(T_LR)); % % figure % plot(freq,abs(Txx_Linear)); % hold on % plot(freq,abs(Tyy_Linear)); % hold on % plot(freq,abs(Txy_Linear)); % hold on % plot(freq,abs(Tyx_Linear)); % Back_te_complex = Back_te(:,2).*cos(Back_te(:,3)/180*pi)-j*Back_te(:,2).*sin(Back_te(:,3)/180*pi); % Back_tm_complex = Back_tm(:,2).*cos(Back_tm(:,3)/180*pi)-j*Back_tm(:,2).*sin(Back_tm(:,3)/180*pi); % Txx_exp_complex = Txx_exp(:,2).*cos(Txx_exp(:,3)/180*pi)-j*Txx_exp(:,2).*sin(Txx_exp(:,3)/180*pi); % Txy_exp_complex = Txy_exp(:,2).*cos(Txy_exp(:,3)/180*pi)-j*Txy_exp(:,2).*sin(Txy_exp(:,3)/180*pi); % Tyx_exp_complex = Tyx_exp(:,2).*cos(Tyx_exp(:,3)/180*pi)-j*Tyx_exp(:,2).*sin(Tyx_exp(:,3)/180*pi); % Tyy_exp_complex = Tyy_exp(:,2).*cos(Tyy_exp(:,3)/180*pi)-j*Tyy_exp(:,2).*sin(Tyy_exp(:,3)/180*pi); % % Txx=Txx_exp_complex./Back_te_complex; % Tyx=Tyx_exp_complex./Back_te_complex; % Tyy=Tyy_exp_complex./Back_te_complex; % Txy=Txy_exp_complex./Back_te_complex; % % % % Txx=Txx_exp_complex % % % Tyx=Tyx_exp_complex % % % Tyy=Tyy_exp_complex % % % Txy=Txy_exp_complex % % % % % T_RR = 0.5*(Txx+Tyy+1i*(Txy-Tyx)); %t_++ % T_RL = 0.5*(Txx-Tyy-1i*(Txy+Tyx)); %t_+- % T_LR = 0.5*(Txx-Tyy+1i*(Txy+Tyx)); %t_-+ % T_LL = 0.5*(Txx+Tyy-1i*(Txy-Tyx)); %t_-- % % % % figure % plot(freq,abs(T_RR)); % hold on % plot(freq,abs(T_LL)); % hold on % plot(freq,abs(T_RL)); % hold on % plot(freq,abs(T_LR)); %