您现在的位置:首页 >> 技术文章 >> MATLAB技术 >> 内容

通过有限差分对平行板电容器进行差分

时间:2020-10-13 20:38:37 点击:

  核心提示:通过有限差分对平行板电容器进行差分...
m=100;n=100;n00=1e8;N0=1e14;aa=2-pi*sqrt(2)*sqrt(1/m^2+1/n^2); 
dt=1e-10;p=760;e=1.6e-19;e0=8.85e-14;U=2600;zl=0.5;rl=0.5;zh=zl/m;rh=rl/n; 
 me=2.9e5/p;mi=2.6e3/p;DL=1800;DT=2190;h=zh;dz=zh;dr=rh; 
 for j=2:1:m 
    for i=2:1:n 
        z(i,j)=j*h-h; 
    end 
end 
for i=2:1:n 
    for j=2:1:m 
        r(i,j)=i*h-h; 
    end 
end 
 
 
%��ʼ�ռ�����ܶȷֲ��� 
for i=2:1:n 
        for j=2:1:m 
                
                 ne(i,j)=N0*exp(-(((i-50)/20)/0.25)^2-(((j-50)/20)/0.25)^2); 
            
               NNN(i,j)=ne(i,j); 
               ni(i,j)=0; 
            end 
end 
for j=2:1:m 
    for i=2:1:n 
        A(i,j)=(j-1)*U/m; 
        AA(i,j)=A(i,j); 
    end 
end 
%��ʼ�ռ������ܶȷֲ��� 
tt=3; 
while tt>0 
for i=3:n-1 
    for j=3:m-1 
        B(i,j)=A(i,j); 
        A(i,j)=1/4*((1+1/(2*(i-1)))*A(i+1,j)+A(i,j+1)+(1-1/(2*(i-1)))*A(i-1,j)+A(i,j-1)-h^2*(-e*(ni(i,j)-ne(i,j))/e0)); 
    end 
end 
 
for i=3:n-1 
    for j=3:m-1 
        A(i,j)=B(i,j)+aa*(A(i,j)-B(i,j)); 
        rr(i,j)=A(i,j)-B(i,j); 
    end 
     
end 
tt=tt-1; 
end 
%��ʼ�ռ�ĵ�ѹ�ֲ��� 
 for i=3:n-2 
    for j=3:m-2 
        Er(i,j)=(A(i+1,j)-A(i-1,j))/(2*rh); 
        Ez(i,j)=(A(i,j+1)-A(i,j-1))/(2*zh); 
    end 
        
 end 
 

联系:highspeedlogic

QQ :1224848052

微信:HuangL1121

邮箱:1224848052@qq.com

网站:http://www.mat7lab.com/

网站:http://www.hslogic.com/

微信扫一扫:

作者:有限差分 来源:有限差分
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做|python人工智能代做|FPGA项目合作(www.hslogic.com) © 2020 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment