tianxian chapter【DOC精选】.docxVIP

源程序Rwg1%RWG1 Geometry calculations - all Chapters% Uses the structure mesh file, e.g. platefine.mat, % as an input.%% Creates the RWG edge element for every inner edge of % the structure. The total number of elements is EdgesTotal.% Outputs the following arrays:%% Edge first node number Edge_(1,1:EdgesTotal)% Edge second node number Edge_(2,1:EdgesTotal)% Plus triangle number TrianglePlus(1:EdgesTotal)% Minus triangle number TriangleMinus(1:EdgesTotal)% Edge length EdgeLength(1:EdgesTotal)% Edge element indicator EdgeIndicator(1:EdgesTotal)%% Also outputs areas and midpoints of separate triangles:% Triangle area Area(1:TrianglesTotal)% Triangle center Center(1:TrianglesTotal) % % This script may handle surfaces with T-junctions % including monopoles over various metal surfaces and % certain metal meshes%% The modification for Chapters 6 and 7 passes the array % Feed (see Chapter 6) %% Copyright 2002 AEMM. Revision 2002/03/14 Chapter 6clearalltic;load(mesh/strip2);Feed=[];[s1 s2]=size(p);if(s1==2)p(3,:)=0; %to convert 2D to 3Dend%Eliminate unnecessary trianglesRemove=find(t(4,:)1); t(:,Remove)=[]; TrianglesTotal=length(t);%Find areas of separate trianglesfor m=1:TrianglesTotal N=t(1:3,m); Vec1=p(:,N(1))-p(:,N(2)); Vec2=p(:,N(3))-p(:,N(2));Area(m) =norm(cross(Vec1,Vec2))/2;Center(:,m)=1/3*sum(p(:,N),2);end%Find all edge elements Edge_ with at least two %adjacent trianglesEdge_=[];n=0;for m=1:TrianglesTotal N=t(1:3,m);for k=m+1:TrianglesTotal M=t(1:3,k); a=1-all([N-M(1) N-M(2) N-M(3)]);if(sum(a)==2) %triangles m and k have common edge n=n+1; Edge_=[Edge_ M(find(a))]; TrianglePlus(n)=m;TriangleMinus(n)=k; end; endendEdgesTotal=length(Edge_);%This block is only meaningful for T junctions%It leaves only two edge elements at a junction Edge__=[Edge_(2,:); Edge_(1,