matlabbiancheng.rar_ansys 电磁_电磁场_电磁场 matlab_磁场

clear;clc; fid=fopen('ELEM.txt','r'); ELEM=fscanf(fid,'%g %g %g %g %g %g %g ',[7 inf]); fclose(fid); fid=fopen('n1.txt','r'); BOUNDARY=fscanf(fid,'%g %g %g ',[3 inf]); fclose(fid); a=0.005;b=0.004; z1=(a^2+b^2)/(3*a*b)/(pi*4e-4);z2=(a^2-2*b^2)/(6*a*b)/(pi*4e-4); z3=-(a^2+b^2)/(6*a*b)/(pi*4e-4);z4=(b^2-2*a^2)/(6*a*b)/(pi*4e-4); i=zeros(1000,1);j=zeros(1000,1);m=zeros(1000,1);q=zeros(1000,1);elem=ELEM'; i=elem(:,2);j=elem(:,3);m=elem(:,4);q=elem(:,5);K=zeros(1071,1071);R=zeros(1071,1); for n=1:1000 K(i(n),i(n))=K(i(n),i(n))+z1; K(i(n),j(n))=K(i(n),j(n))+z2; K(i(n),m(n))=K(i(n),m(n))+z3; K(i(n),q(n))=K(i(n),q(n))+z4; K(j(n),j(n))=K(j(n),j(n))+z1; K(j(n),m(n))=K(j(n),m(n))+z4; K(j(n),q(n))=K(j(n),q(n))+z3; K(m(n),m(n))=K(m(n),m(n))+z1; K(m(n),q(n))=K(m(n),q(n))+z2; K(q(n),q(n))=K(q(n),q(n))+z1; K(j(n),i(n))=K(i(n),j(n)); K(m(n),i(n))=K(i(n),m(n)); K(q(n),i(n))=K(i(n),q(n)); K(m(n),j(n))=K(j(n),m(n)); K(q(n),j(n))=K(j(n),q(n)); K(q(n),m(n))=K(m(n),q(n)); R(i(n))=250*a*b+R(i(n)); R(j(n))=250*a*b+R(j(n)); R(m(n))=250*a*b+R(m(n)); R(q(n))=250*a*b+R(q(n)); end; o=zeros(1071,1);boundary=BOUNDARY'; b1=zeros(140,1);b1=boundary(:,1); for m=1:140 K(b1(m),:)=o'; K(:,b1(m))=o; K(b1(m),b1(m))=1; R(b1(m))=0; end; Jie=zeros(1071,1); T=sparse(K); Jie=K\R; a=zeros(1071,2); for n=1:1071 a(n,1)=n; a(n,2)=Jie(n); end; d=a'; fprintf( 'DNSOL, %g, AZ,, %g\n',d);