搜索
©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.()JournalofChangchunNormalUniversity(NaturalScience)20102Feb.2010MATLAB,(,723003)[],MATLAB,[];;MATLAB;;[]O4-39[]A[]1008-178X(2010)01-0020-04[]2009-08-18[](1975-),,,,,,,MATLAB,,111,,Idl_PdB_=04Idl_e_rr2,[1]r_r_PIdl_,r_Idl_r_=r_+r_,r_=xi_+yj_+zk_,r_=R(cosi_+sinj_)(R),dl_=Rdcos(+2)i_+sin(+2)j_=Rd(-sini_+cosj_),1z,,P(x,0,z):dB_=04Idl_e_rr2=0IRd4r3zcosi_+zsinj_+(R-xcos)k_.dBx=0IRd4r3zcos,Bx=dBx=200IR4r3zcosd.(1)dBy=0IRd4r3zsin,By=dBy=200IR4r3zsind.(2)dBz=0IRd4r3(R-xcos),Bz=dBz=200IR4r3(R-xcos)d.(3)r=x2+z2+R2-2Rxcos.2MATLAB02©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.(1)(2)(3),MATLAB[2],symssitaxzR%sitaxzRR=1;%R1mf=R3z3cos(sita)/((R.^2+x.^2+z.^2-23R3x3cos(sita)).^1.5);g=R3z3sin(sita)/((R.^2+x.^2+z.^2-23R3x3cos(sita)).^1.5);h=R3(R-x3cos(sita))/((R.^2+x.^2+z.^2-23R3x3cos(sita)).^1.5);Bx=int(f,sita,0,23pi);By=int(g,sita,0,23pi);Bz=int(h,sita,0,23pi);%,0I4,:Bx=-23(EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))3x^2-EllipticE(23(1/(1+x^2+z^2+23x)3x)^(1/2))3x^2-23EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))3x+EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))3z^2-EllipticE(23(1/(1+x^2+z^2+23x)3x)^(1/2))+EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))-EllipticE(23(1/(1+x^2+z^2+23x)3x)^(1/2))3z^2)3((1+x^2+z^2-23x)/(1+x^2+z^2+23x))^(1/2)3z/(1+x^2+z^2-23x)^(3/2)/x;By=0;Bz=23(EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))3x^2-EllipticE(23(1/(1+x^2+z^2+23x)3x)^(1/2))3x^2-23EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))3x-EllipticE(23(1/(1+x^2+z^2+23x)3x)^(1/2))3z^2+EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))+EllipticK(23(1/(1+x^2+z^2+23x)3x)^(1/2))3z^2+EllipticE(23(1/(1+x^2+z^2+23x)3x)^(1/2)))3((1+x^2+z^2-23x)/(1+x^2+z^2+23x))^(1/2)/(1+x^2+z^2-23x)^(3/2);212MATLABEllipticE(x)EllipticK(x),Matlab,mfun(EllipticE,x)mfun(EllipticK,x)[3]surfl(x,z,Bx)surfl(x,z,Bz),2323:[x,z]=meshgrid(0:0.03:2,-1:0.03:1);12©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.=2.3(mfun(EllipticK,(2.3(1./(1+x.^2+z.^2+2.3x).3x).^(1/2))).3x.^2-mfun(EllipticE,(23(1./(1+x.^2+z.^2+2.3x).3x).^(1/2))).3x.^2-23mfun(EllipticK,(23(1./(1+x.^2+z.^2+2.3x).3x).^(1/2))).3x-mfun(EllipticE,(23(1./(1+x.^2+z.^2+2.3x).3x).^(1/2))).3z.^2+mfun(El2lipticK,(23(1./(1+x.^2+z.^2+2.3x).3x).^(1/2)))+mfun(EllipticK,(23(1./(1+x.^2+z.^2+2.3x).3x).^(1/2))).3z.^2+mfun(EllipticE,(23(1./(1+x.^2+z.^2+23x).3x).^(1/2)))).3((1+x.^2+z.^2-2.3x)./(1+x.^2+z.^2+2.3x)).^(1/2)./(1+x.^2+z.^2-2.3x).^(3/2);[4]surfl(x,z,Bz);xlabel(x);ylabel(z);zlabel(Bz);3311,(x,0,z)yBy=0,(xoz),(,x)31223,,z=0,x=1,,,,z=0,x=11,1z=0,x=1zxBx(103)Bz(103)-0.00200.99800.99901.00001.00101.0020-0.5005-0.8004-1.0000-0.7996-0.49950.50750.40740.0073-0.3926-0.4926-0.00100.99800.99901.00001.00101.0020-0.4004-1.0005-2.0000-0.9995-0.39960.80801.00810.0080-0.9919-0.792000.99800.99901.00001.00101.002000NaN001.00832.0090NaN-1.9910-0.99170.00100.99800.99901.00001.00101.00200.40041.00052.00000.99950.39960.80801.00810.0080-0.9919-0.79200.00200.99800.99901.00001.00101.00200.50050.80041.00000.79960.49950.50750.40740.0073-0.3926-0.49261,(z=0,x=110000),Br=031322©1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreserved.=0,(),22Bz(103)xBzxBzxBzxBz00.00630.9910.22901.001-1.99101.100-0.01590.2000.00650.9920.25701.002-0.99171.200-0.00670.3000.00670.9930.29281.003-0.65881.300-0.00380.4000.00720.9940.34061.004-0.49241.400-0.00250.5000.00780.9950.40741.005-0.39261.500-0.00180.6000.00880.9960.50761.006-0.27871.600-0.00140.7000.01060.9970.67461.007-0.24311.700-0.00100.8000.01420.9981.00841.008-0.21551.800-0.00080.9000.02460.9992.00901.009-0.19341.900-0.00070.9800.10601.000NAN1.010-0.17532.000-0.00052(x11000)Bz0,(x11000)Bz0,,3,,,x=01999,B210090;x=01998,B110084;x=01980,B0110600199901991,,019900,,,,,x=11001,B-119910;x=11002,B-019917;x=11010,B-0117531100111009,x11009,,x(x=11000),(x=01998x=11002),,,,,4,,MAT2LAB,,,()[5];MATLAB[][1].()[M].:,2001.[2]..Matlab7.0[M].:,2004.[3].MATLAB6[M].:,2001.[4].Matlab[M].:,2004:292-295.[5],.[J].,2004(3).AnalysisontheMagneticFieldDistributionofRingElectricCurrentwithMATLABWANGYu-mei,SUNQing-long(PhysicsDepartment,ShaanxiUniversityofTechnology,Hanzhong723003,China)Abstract:ThispaperderivestheintegralrepresentationofthemagneticfielddistributionofcircularcurrentaccordingtoBiot-SavartLaw,givesthecomputedresultsusingthesymbolicintegrationofMATLAB,anddrawsthethree-dimensionalcurveofthemagneticfielddistribution.Somerepresentativepointsinthecomputedresultsareselectedtodiscussthedis2tributedruleofthemagneticfield.Keywords:ringelectriccurrent;magneticfield;MATLAB;symbolicintegration;three-dimensionalcartography32