matlab 计算结果为nan,matlab 计算 结果总是为Nan
本人刚刚接触matlab,对这些运算不是很懂 计算ni的位置(exp(-E_g./(2.*k.*T))+eps)这个值之前一直是0 加了eps后就有结果显示了 后面部分exp((alfa.*T)./(2.*k.*(T+beita)))一直为无穷大然后结果就为Nan 不知道怎么解决 求大神帮忙
q = 1.602*10^(-19);
k = 1.38*10^(-23);
fai_msf = -0.08;
fai_msb = 0.17;
kexi_Si = 1.05*10^(-12);
kexi_ox = 3.45*10^(-13);
d_fox = 10^(-6);
d_box = 3.8*10^(-5);
d_Si = 3*10^(-5);
%ni = 1.45*10^10;
Dn = 34;
%ni = 1.45*10^10;
%fai_f = 0.29;
fai_T = 0.026;
N_a = 10^15;
N_d = 10^20;
tao = 10^(-7);
s_0 = 10^4;
s_dSi = 10^4;
y = 10^(-7);
Q_oxf = 3*10^10*q;
Q_oxb = 10^10*q;
Q_dep = -q*N_a*d_Si;
Q_sb = 0;
V_BG = 0;
V_A = -1;
V_GK = 0.4;
E_g = 1.21;
alfa = 4.73*10^(-4);
beita = 636;
T = 250:10:350;
ni=3.39.*10^15.*T.^(3/2).*(exp(-E_g./(2.*k.*T))+eps).*exp((alfa.*T)./(2.*k.*(T+beita)));
fai_f = (k.*T./q).*log(N_a./ni);
fai_sb = fai_f;
C_Si = kexi_Si./d_Si;
C_ox = kexi_ox./d_fox;
C_box = kexi_ox./d_box;
A = -fai_msf+Q_oxf/C_ox+Q_dep/(2*C_ox)-(C_Si/C_ox)*(C_box./(C_box+C_Si))*(fai_msb-Q_oxb./C_box-Q_dep./(2*C_box));
B = 1+C_Si./C_ox-(C_Si./C_ox)*(C_box./(C_box+C_Si))*(C_Si./C_box);
C = V_BG-fai_msb+Q_oxb./C_box;
D = C_Si/C_box;
E = Q_dep/(2*C_box);
F = 1+C_Si/C_box;
G = q*N_a*y^2/(2*kexi_Si);
H = (q*N_a*d_Si)/(2*kexi_Si);
fai_sf = (V_GK+A)./B;
fai_sb = (C+D*fai_sf+E)/F;
fai_y = G+((fai_sf-fai_sb)./d_Si-H).*y+fai_sf;
% fai_sf1 = (V_GK2+A)./B;
% fai_sb1 = (C+D*fai_sf1+E)/F;
% fai_y1 = G+((fai_sf1-fai_sb1)./d_Si-H).*y+fai_sf1;
% n_VGK2 = ni*exp((fai_y1-fai_f)/fai_T);
n_y = ni.*exp((fai_y-fai_f)./fai_T);
p_y = ni./exp((fai_y-fai_f)/fai_T);
V_D = fai_T.*log(p_y.*N_d./(ni.^2));
V_D1 = fai_T.*log(N_d.*n_y/(ni.^2));
V_D0 = fai_T.*log(N_a.*N_d./(ni.^2));
V_D10 = fai_T*log(N_d/N_a);
fai_N = (V_D0+V_D10)/2-V_A;
V_A1 = -(fai_N-(fai_y-fai_f)-V_D);
Ld = (2.*kexi_Si.*(V_D-V_A1).*(N_a+N_d)./(q.*N_a.*N_d)).^(0.5);
xp = (2.*kexi_Si.*(V_D-V_A1).*N_d./(q.*N_a.*(N_a+N_d))).^(0.5);
L_d1 = (2.*kexi_Si.*(V_D+(1+V_A1)).*(N_a+N_d)/(q.*N_a.*N_d)).^(0.5);
xp1 = (2.*kexi_Si.*(V_D+(1-V_A1)).*N_d/(q.*N_a.*(N_a+N_d))).^(0.5);
tao_eff = 1/(p_y/(tao*ni)+(s_0+s_dSi)/d_Si);
% tao_eff = tao*ni/p_y;
Ln = (Dn*tao_eff)^(0.5);
L = 2*10^(-3);
L1 = L-xp-xp1;
x = 0:L1/100:L1;
nx = (N_d*exp(q*V_A1/(k*T))-n_y)*sinh(x/Ln)/sinh(L1/Ln)+n_y;
Jnx = q*Dn*(N_d*exp(q*V_A1/(k*T))-n_y)*cosh(x/Ln)/(sinh(L1/Ln)*Ln);
figure(1),plot(x,nx);
figure(2),plot(x,Jnx);
matlab 计算结果为nan,matlab 计算 结果总是为Nan相关推荐
- matlab求和,取均值以及计算标准差的时候数据有nan(nansum,nanmean,nanstd)
matlab求和,取均值以及计算标准差的时候数据有nan使用下面三个特定的函数 (nansum,nanmean,nanstd)
- float gpu 加速_(总结篇)使用 MATLAB GPU 加速计算|MATLAB 并行计算与分布式服务器|MATLAB技术论坛...
本帖最后由 蓝云风翼 于 2013-12-18 17:28 编辑 注: 利用gpu加速有一下工具 1.JACKET 可从帖子中寻找 2.MATLAB a.并行计算工具箱 gpuArray,查看支持gp ...
- matlab计算 var,matlab 在险价值 VaR 的计算
matlab 在险价值 VaR 的计算 matlab 在险价值 VaR 的计算 VaR 模型 数据获取 历史模拟法 蒙特卡罗模拟法 参数模型法 代码和数据下载 VaR 模型 Value at Risk ...
- 科学计算:Python VS. MATLAB (1)----给我一个理由先
科学计算:Python VS. MATLAB (1)----给我一个理由先 MATLAB 是一种用于算法开发.数据可视化.数据分析以及数值计算的高级技术计算语言和交互式环境.使用 MA ...
- 科学计算:Python VS. MATLAB(3)----线性代数基础
科学计算:Python VS. MATLAB(3)----线性代数基础 按:在介绍工具之前先对理论基础进行必要的回顾是很必要的.没有理论的基础,讲再多的应用都是空中楼阁.本文主要设涉及线性代数和矩阵论 ...
- 科学计算:Python VS. MATLAB(4)----图形系统简介
科学计算:Python VS. MATLAB(4)----图形系统简介 一.一般概念 图形系统就是用来实现所谓的可视化的.在学习具体的命令之前,先了解一个可视化的图形具有哪些特征,以及这些特征具有什么 ...
- matlab中计算不等式的解,matlab解不等式
MATLAB求解非线性规划good (1)_理学_高等教育_教育专区.数学建模课... x4 ? x5 MATLAB 程序如下: f=[1; 2; 3; 1; 1]; intcon=[1 2 3 4 ...
- matlab 电力系统分析 毕设,基于MATLAB的P-Q分解法电力系统潮流计算毕业设计
基于MATLAB的P-Q分解法电力系统潮流计算毕业设计 摘 要 电力系统潮流计算是研究电力系统稳定运行情况的一种重要的计算,在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用它来定量地分析比 ...
- matlab将常值函数转换为变量,MATLAB与科学计算期末复习题题库15.11.12
MATLAB 与科学计算期末复习题题库(第一部分) 一.填空 1.MATLAB 的主界面是一个高度集成的工作环境,有四个不同职责分工的窗口,分别 为 . . .和 窗口. 2.MATLAB 的值,应在 ...
- matlab求三维图形表面积,基于MATLAB的不规则平面立体表面积计算
·12· 设计与研究 机械 2012 年第 2 期 总第 39 卷 --------------- 收稿日期:2011-06-10 作者简介:和蕴锋(1988-),女,河南漯河人,硕士研究生,主要研究 ...
最新文章
- 贾扬清、Alex Smola、Julia创始人等大咖齐聚,WAIC开发者日共话AI未来
- select子句顺序
- 前端笔记之NodeJS(三)Expressejs模板引擎请求识别
- 金融领域下的数据挖掘算法应用:智能推荐算法模型
- imperial college application status check portal
- php使用redis做缓存,php使用redis做缓存和使用redis保存session
- arm-none-eabi-gcc install
- 【LeetCode笔记】98. 验证二叉搜索树(Java、dfs、中序遍历、二叉树)
- web开发程序员有几种_每个程序员都应该知道的10种流行的Web开发工具
- linux中的管道和重定向
- 2019级C语言大作业 - 十步万度
- CS安装卸载测试总结
- 阿里大牛程序员的Java问题排查工具单
- VGA信号一些分辨率的参数(摘抄)
- 智哪儿头条 萤石网络冲刺科创板 深圳购买全屋智能有补贴
- php计算股票均线,php写智能选股,股票回测系统之--MACD指标计算方法及实例
- oracle分区表备份恢复
- 解决Skype一台电脑登陆多个账号的问题
- zxr10交换机配置手册vlan_中兴ZXR10配置说明.doc
- Ardunio开发实例-MAX30102脉搏血氧饱和度和心率监测传感器