资源描述：
《电磁场生医实验2》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、《生物医学电磁场数值分析》实验报告一、实验目的熟悉二维有限元前处理技术，学握二维均匀介质平面场域和非均匀介质平面场域的剖分方法，为电磁场有限元数值计算打下基础。二、实验条件硕件：微型计算机。软件：MATLAB6.0软件。三、实验内容1、根据场域特点编制二维平面场域剖分程序。2、对所剖分场域进行图形显示。四、实验步骤(按照下面步骤中的小条贴上相应的程序)1、对二维均匀介质场域进行剖分，首先找到节点他标并编号，然后进行单元的剖分并编号。2、对二维非均匀介质场域进行剖分，方法同上。3、绘制二维场域节点分布图形。4、绘制二
2、维场域剖分单元图形。五、实验结果(将实验结果，如一些生成的图，贴在下面)举例：二维均匀介质场域的剖分结果如下所示。(a)二维均匀场域剖分结果图clearall;nx=l;ny=l;xl=[0z6];yi=[o,4]；nxl=6;ny1=4;num_nodex=nxl+l；num_nodey=nyl+l；total_node=num_nodex*num_nodey;total_element=nxl*nyl*2;stepx=(xl(2)-xl(1))/nxl;stepy=(yl(2)-yl(1))/nyl;p=l;f
3、ori=l:num_nodextempx=xl(1)+stepx*(i-1);forj=l:num_nodeynode(pz1)=tempx;node(p,2)=yl(1)+stepy*(j-1);P=P+1；endendfigure(1);plot(node(:f1)fnode(:z2)f1b*1);xlabel(1竹点分布图1);3.50.5-•0123456节点分布图p=l;fori=l:num_nodex-lforj=l:num_nodey-lK=num_nodey*(i-1)+j;N=K+1；M=N+nu
4、m_nodey;element(pz:)=[K,M,N];P=P+1；%匕类endforj=l:num_nodey-lK=num_nodey*(i-1)+j;M=K+num_nodey;N=M+1;element(p,:)=[K,M,N];P=P+1；endendfigure(2);forp=l:total_elementa=node(element(p,1)z1);b=node(element(p,2),1);c=node(element(p,3),1);d=node(element(p,1)z2);e=node
5、(element(pr2),2);f=node(element(p,3)z2);plot([a,b,c,a],[d,e,f,d]);holdon;endxlabel(1单位剖分图1)3.5弟位扁分啊(b)二维非均匀场域剖分结果图clearall;nx=3;ny=l;xl=[0,2,4z6];yl=[Of4];nxl=[4,2,2];nyl=4;num_nodex=sum(nxl)+1;num_nodey=sum(nyl)+1；total_node=num_nodex*num_nodey;total_element=
6、(num_nodex-l)*(num_nodey-l)*2;stepy=(yl(2)-yl(1))/nyl;p=l;fork=l:nxstepx(k)=(xl(k+1)-xl(k))/nxl(k);fori=l:nxl(k)tempx=xl(k)+stepx(k)*(i-1);forj=l:ny1+1node(pz1)=tempx;node(p,2)=yl(1)+stepy*(j-1);P=P+1；endendendforj=l:nyl+1node(p,1)=xl(nx+1)；node(p,2)=yl(1)+ste
7、py*(j-1);P=P+1；endfigure(1);plot(node(:f1)rnode(:f2)r匕*1);xlabel('非均匀介质节点图');4-04i44i♦卜012345非均匀介质节点图P=l；fori=l:num_nodexforj=l:num_nodey-lelement(pz1)=num_nodey*(i-1)+j;element(p,2)=element(p,1)+l+num_nodey;element(pz3)=element(p，l)+1;P=P+1；endforj=l:num_node
8、y-lelement(p,1)=num_nodey*(i-1)+j;element(pz2)=element(p,1)+num_nodey;element(p，3)=element(p,2)+1;P=P+1；endendfigure(2);forp=l:total_elementa=node(element(p,1),1);b=node(element(p,
