【物理应用】基于matlab二维边界单元法计算腐蚀电位【含Matlab源码 523期】

一、简介

基于matlab二维边界单元法计算腐蚀电位

二、源代码

clear;

fid=fopen('input.dat','r');
indat=fscanf(fid,'%g%g%d%g',[4,inf]);
indat=indat';
xb=indat(:,1);
yb=indat(:,2);
bt=indat(:,3);
bv=indat(:,4);
n=length(xb)-1;

for i=1:n
    xm(i)=0.5d0*(xb(i)+xb(i+1));
    ym(i)=0.5d0*(yb(i)+yb(i+1));
    lm(i)=sqrt((xb(i+1)-xb(i))^2d0+(yb(i+1)-yb(i))^2d0);
    nx(i)=(yb(i+1)-yb(i))/lm(i);
    ny(i)=(xb(i)-xb(i+1))/lm(i);
end

for m=1:n
    b(m)=0d0;
    for k=1:n
        if(k==m)
            G=0.0;
            F=lm(k)/(2.0*pi)*(log(lm(k)/2.0)-1.0);
            del=1.0;
        else
            [F,G]=findfg(xm(m),ym(m),xb(k),yb(k),nx(k),ny(k),lm(k));
            del=0.0;
        end
        if(bt(k)==0)
            A(m,k)=-F;
            b(m)=b(m)+bv(k)*(-G+0.5d0*del);
        else
            A(m,k)=G-0.5d0*del;
            b(m)=b(m)+bv(k)*F;
        end
    end
end
z=A\b';

for m=1:n
   u(m)=(1-bt(m))*bv(m)+bt(m)*z(m);
   q(m)=(1-bt(m))*z(m)+bt(m)*bv(m);
end

for j=1:9
   y(j)=0.1*j;
   for i=1:9
      x(i)=0.1*i;
      s(j,i)=0d0;
      for k=1:n
         [F,G]=findfg(x(i),y(j),xb(k),yb(k),nx(k),ny(k),lm(k));
         s(j,i)=s(j,i)+u(k)*G-q(k)*F;
      end
   end
end

三、运行结果

四、备注

完整代码或者代写添加QQ 1564658423
往期回顾>>>>>>
【物理应用】基于matlab杨氏双缝干涉实验【含Matlab源码 001期】
【物理应用】基于matlab井筒多相流【含Matlab源码 055期】
【物理应用】基于matlab 二维对流扩散温度场【含Matlab源码 056期】
【物理应用】基于matlab内弹道【含Matlab源码 057期】
【物理应用】基于matlab两个光纤光栅的新光子微波发生器的仿真分析【含Matlab源码 058期】
【物理应用】基于matlab计算油气井井底压力【含Matlab源码 059期】
【物理应用】基于matlab光流场【含Matlab源码 060期】
【物理应用】基于matlab 带电粒子在放射状电场和均匀磁场下的混沌运动模拟【含Matlab源码 061期】
【物理应用】基于matlab异构交通流【含Matlab源码 062期】
【物理应用】基于matlab模拟斜抛运动含GUI界面【含Matlab源码 063期】
【物理应用】基于matlab计算晶体结构的x射线衍射图谱【含Matlab源码 064期】
【物理应用】基于matlab衍射光栅教学版【含Matlab源码 065期】
【物理应用】基于matlab电荷电场线和电势【含Matlab源码 066期】
【物理应用】基于matlab速度场找到漩涡【含Matlab源码 067期】
【物理应用】基于matlab车道检测仿真【含Matlab源码 068期】
【物理应用】基于matalb GUI界面的电偶极子【含Matlab源码 069期】
【物理应用】基于matlab车载自组织网络中路边性能及防碰撞算法研究【含Matlab源码 263期】
【物理应用】基于matlab Q学习的无线体域网路由方法【含Matlab源码 264期】
【物理应用】基于matlab GUI家庭室内温湿度控制【含Matlab源码 326期】
【电路仿真】基于simulink三相相控变流器的设计【含Matlab源码 327期】
【物理应用】基于matlab的导线网坐标计算【含Matlab源码 328期】
【物理应用】基于matlab的功率谱估计【含Matlab源码 329期】
【物理应用】基于动态矩阵控制方法的热水器出口温度控制【含Matlab源码 330期】
【电路仿真】基于matlab simulink三相半波可控整流电路【含Matlab源码 331期】
【电路仿真】基于matlab simulink三相桥式全控整流电路【含Matlab源码 332期】
【电路仿真】基于simulink模拟锁相环设计【含Matlab源码 339期】
【电路仿真】基于simulink全数字锁相环设计【含Matlab源码 340期】
【物理应用】基于matlab直角坐标的电力系统潮流计算【含Matlab源码 341期】
【光学】基于matlab单缝衍射【含Matlab源码 342期】
【物理应用】基于matlab RBF神经网络算法控制卫星轨道和姿态【含Matlab源码 377期】
【物理应用】基于matlab红绿灯车辆数量识别【含Matlab源码 378期】
【电路仿真】基于simulink 四旋翼的PID控制【含Matlab源码 454期】
【光学】 基于matlab GUI光的双缝干涉【含Matlab源码 455期】
【电路仿真】基于matlab光伏电池simulink仿真模型【含Matlab源码 486期】
【物理应用】 基于matlab PIC模型的太阳风粒子模拟【含Matlab源码 493期】
【光学】基于matlab光栅衍射仿真【含Matlab源码 502期】
【物理应用】基于matlab波数谱计算【含Matlab源码 508期】
【物理应用】基于matlab 仿真圆孔的菲涅尔衍射【含Matlab源码 522期】