% function T = cal_Trans_all(omi, rushe)
clear;
clc;
c=3e8;
h=6.626e-34;
e=1.6e-19;
ev=1.6e-19;
mev=1.6e-22; %%%%%
epsilon_0=8.85e-12;
mu_0=4*pi*1e-7;
z0=mu_0*c;
alfa=1/137;
v_F=1e6;
E_F=0.25*ev*1.0;
% h_bar_gama=2.5*mev;
h_bar=h/(2*pi);
sigma0=e^2/(4*h_bar);
i=complex(0,1);
% omega=2*pi*4.4064e+11;
omega=2*pi*3.3310e+11;
beta=(2*alfa*E_F/h_bar/omega)^2;
omega0=2*alfa*E_F/h_bar/sqrt(7)/(2*pi);
theta_0=0*pi/180;
fp=fopen('Iin_Iout.dat','w');
for j=1:1001
Iout(1)=0;
% cc=(e*v_F/omega/E_F)^2;
cc=(e*v_F/omega/(0.25*ev))^2;
delta_Iout=1.9164e+08;
Iout(j)= Iout(1)+(j-1)*0.000705/cc;
e1=1;
u1=1;
e2=1;
u2=1;
gama1=sigma0*4*E_F/(pi*h_bar*omega);
gama3=sigma0*9*e^2*v_F^2/(2*pi*h_bar*omega^3*E_F);
gama5=sigma0*25*e^4*v_F^4/(4*pi*h_bar*omega^5*E_F^3);
k1=(omega/c)*sqrt(e1)*sqrt(u1)*sqrt(1-sin(theta_0)^2/(e1*u1));
k2=(omega/c)*sqrt(e2)*sqrt(u2)*sqrt(1-sin(theta_0)^2/(e2*u2));
q1(j)=sqrt(e1)/sqrt(u1)*sqrt(1-sin(theta_0)^2/(e1*u1));
q2(j)=sqrt(e2)/sqrt(u2)*sqrt(1-sin(theta_0)^2/(e2*u2));
Eout(j)=sqrt(Iout(j));
Ein(j)=Eout(j)*((k1+k2)/(2*k1)+i*0.5*mu_0*omega/k1*(gama1-gama3*Iout(j)-gama5*Iout(j)^2));
Ere(j)=-1*Eout(j)*i*0.5*mu_0*omega/k1*(gama1-gama3*Iout(j)-gama5*Iout(j)^2);
Iin(j)=Iout(j)*(((k1+k2)/(2*k1))^2+(0.5*mu_0*omega/k1*(gama1-gama3*Iout(j)-gama5*Iout(j)^2))^2); % Iin=Ein^2;Iout=Eout^2;
% T(j)=q2/q1*Iout(j)/Iin(j);
t(j)=Eout(j)/Ein(j);
r(j)=Ere(j)/Ein(j);
% R(j)=abs(r(j))^2;
T(j)=q2/q1*abs(t(j))^2;
theta(j)=angle(t(j)); % arg(); real(complex),imag(complex),abs(complex),angle(complex),分别求复数的实部、虚部、模、幅角(弧度制);
% theta(j)=atan(imag(t(j))/real(t(j)));
Ire(j)=Iout(j)*(0.5*z0*(gama1-gama3*Iout(j)-gama5*Iout(j)^2))^2; %Ire=Ere^2;
R(j)=Ire(j)/Iin(j);
A(j)=1-R(j)-T(j);
fprintf(fp,'%10.5f %10.5f\n',Iin(j)*sqrt(e1)/z0/2/1e9, Iout(j)*sqrt(e2)/z0/2/1e9);
% fprintf(fp,'%10.5f %10.5f\n',Iin(j)*sqrt(e1)/z0/2/1e9, T(j));
% fprintf(fp,'%10.5f %10.5f\n',Iin(j)*sqrt(e1)/z0/2/1e9, theta(j)*180/pi);
end
fclose(fp);
figure(1)
plot(Iin(1:1000)*sqrt(e1)/z0/2/1e9, Iout(1:1000)*sqrt(e2)/z0/2/1e9)
% plot(Iin(1:1000)*cc, Iout(1:1000)*cc)
% plot(abs(Ein(1:1000)), Eout(1:1000))
% plot(sqrt(Iin(1:1000)), sqrt(Iout(1:1000)))
% % hold on
%axis tight
% figure(2)
% plot(Iin(1:1000)*sqrt(e1)/z0/2/1e9, T(1:1000)) % GW/M^2
% % plot(sqrt(Iin(1:1000)), T(1:1000))
% figure(3)
% plot(Iin(1:1000)*sqrt(e1)/z0/2/1e9, theta(1:1000)*180/pi)
没有合适的资源?快使用搜索试试~ 我知道了~
基于MATLAB实现的石墨烯在太赫兹频率的非线性光学特性的数值模拟+使用说明文档.zip
![preview](https://csdnimg.cn/release/downloadcmsfe/public/img/white-bg.ca8570fa.png)
共2个文件
md:1个
m:1个
![preview-icon](https://csdnimg.cn/release/downloadcmsfe/public/img/scale.ab9e0183.png)
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 56 浏览量
2024-05-22
22:44:37
上传
评论
收藏 14KB ZIP 举报
温馨提示
CSDN IT狂飙上传的代码均可运行,功能ok的情况下才上传的,直接替换数据即可使用,小白也能轻松上手 【资源说明】 基于MATLAB实现的石墨烯在太赫兹频率的非线性光学特性的数值模拟+使用说明文档.zip 1、代码压缩包内容 主函数:main.m; 调用函数:其他m文件;无需运行 运行结果效果图; 2、代码运行版本 Matlab 2020b;若运行有误,根据提示GPT修改;若不会,私信博主(问题描述要详细); 3、运行操作步骤 步骤一:将所有文件放到Matlab的当前文件夹中; 步骤二:双击打开main.m文件; 步骤三:点击运行,等程序运行完得到结果; 4、仿真咨询 如需其他服务,可后台私信博主; 4.1 期刊或参考文献复现 4.2 Matlab程序定制 4.3 科研合作 功率谱估计: 故障诊断分析: 雷达通信:雷达LFM、MIMO、成像、定位、干扰、检测、信号分析、脉冲压缩 滤波估计:SOC估计 目标定位:WSN定位、滤波跟踪、目标定位 生物电信号:肌电信号EMG、脑电信号EEG、心电信号ECG 通信系统:DOA估计、编码译码、变分模态分解、管道泄漏、滤波器、数字信号处理+传输+分析+去噪、数字信号调制、误码率、信号估计、DTMF、信号检测识别融合、LEACH协议、信号检测、水声通信 5、欢迎下载,沟通交流,互相学习,共同进步!
资源推荐
资源详情
资源评论
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
收起资源包目录
![package](https://csdnimg.cn/release/downloadcmsfe/public/img/package.f3fc750b.png)
![file-type](https://csdnimg.cn/release/download/static_files/pc/images/minetype/UNKNOWN.png)
![file-type](https://csdnimg.cn/release/download/static_files/pc/images/minetype/UNKNOWN.png)
共 2 条
- 1
资源评论
![avatar-default](https://csdnimg.cn/release/downloadcmsfe/public/img/lazyLogo2.1882d7f4.png)
![avatar](https://profile-avatar.csdnimg.cn/default.jpg!1)
IT狂飙
- 粉丝: 4778
- 资源: 2640
![benefits](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-1.c8e153b4.png)
下载权益
![privilege](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-2.ec46750a.png)
C知道特权
![article](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-3.fc5e5fb6.png)
VIP文章
![course-privilege](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-4.320a6894.png)
课程特权
![rights](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-icon.fe0226a8.png)
开通VIP
上传资源 快速赚钱
我的内容管理 展开
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助
![voice](https://csdnimg.cn/release/downloadcmsfe/public/img/voice.245cc511.png)
![center-task](https://csdnimg.cn/release/downloadcmsfe/public/img/center-task.c2eda91a.png)
安全验证
文档复制为VIP权益,开通VIP直接复制
![dialog-icon](https://csdnimg.cn/release/downloadcmsfe/public/img/green-success.6a4acb44.png)