MATLAB Files - Electromagnetics - Lecture Notes | ENGR 367, Study notes of Guiding Electromagnetic Systems

Download MATLAB Files - Electromagnetics - Lecture Notes | ENGR 367 and more Study notes Guiding Electromagnetic Systems in PDF only on Docsity! Electromagnetics (ENGR367) MATLAB files function y=magvector(R) %Calculates the magnitude of a Cartesian vector R y=sqrt(R(1)^2 + R(2)^2 + R(3)^2); function y=unitvector(R) %Calculates the unit vector of a Cartesian vector R y=R/magvector(R); %M-File: MLP0310 %Find the magnetic field intensity at any observation point %resulting from a ring of radius a and current I, %in the aphi direction centered in the x-y plane. df=1; %increment in degrees a=1; %ring radius in m I=1; %current in A Ro=input('vector location of observation point, in the form [x y z]: '); for j=1:df:360; %increments phi (f) from 1 to 360 in df steps Fr=j*pi/180; %converts phi (Fr) into radians Rs=[a*cos(Fr) a*sin(Fr) 0]; %defines src pt vector as ftn of phi (Fr) as=unitvector(Rs); %finds src pt unit vector dL=a*df*(pi/180)*cross([0 0 1],as); %finds dL vector (mag. & dir.) Rso=Ro-Rs; %finds vector (Rso) from src to obs pt aso=unitvector(Rso); %unitizes Rso dH=I*cross(dL,aso)/(4*pi*(magvector(Rso))^2); %finds dH by Biot-Savart dHx(j)=dH(1); %isolates x-component dHy(j)=dH(2); %isolates y-component dHz(j)=dH(3); %isolates z-component end %terminates phi loop H=[sum(dHx) sum(dHy) sum(dHz)] %sums 360 contrib's to each comp. of H %at the specified observation point