matlab 第三章作業(yè)

1、第三章作業(yè)第一題：format longsyms xc1=int(x4*(1-x)4)/(1+x2),0,1); %符號工具箱int計算c1=eval(c1); %指令eval一般應(yīng)用于int后結(jié)果較復(fù)雜時，進(jìn)一步求數(shù)值解。f=(x)(x.4.*(1-x).4)./(1+x.2); c2=quad(f,0,1); %數(shù)值積分指令quad要求積分限為常數(shù) 第二題：format longf=(x)(2/sqrt(pi)*exp(-x.2); %定義誤差函數(shù)for k=0.1:0.1:1.0 c1=quad(f,0,k); c2=erf(k); c=abs(c1-c2); bar(k,c) hold

2、 onend 第三題：vp=1578;f=(v)(4/sqrt(pi)*(v.2)/vp3).*exp(-v.2)/vp2);c1=quad(f,0,vp)c2=quad(f,0,3.3*vp)c3=quad(f,3e4,3e8)結(jié)果：c1 =0.427593285448571c2 =0.999927444001040 在03.3Vp間的分子數(shù)接近為全部c3 =3.520658577050803e-151 速度很大的分子數(shù)很少第四題：clear;R=1; y=-4:0.11:4; z=-4:0.11:4 ;the=-0.5*pi : pi/400 : 0.5*pi;Y, Z,T=meshgri

3、d(y,z,the);r=sqrt(cos(T).2+Z.2+(Y-sin(T).2);dv = 1./r ;v=trapz(dv,3);ey,ez=gradient(-v,0.5);figure(1)subplot(2,2,1)contour(Y(:,:,1),Z(:,:,1),v,10)subplot(2,2,2)SY,SZ=meshgrid(-4:0.2:4,0.1,-0.1);streamline(Y(:,:,1),Z(:,:,1),ey,ez,SY,SZ);SY,SZ=meshgrid(-4:0.2:4,-0.005,0.005);streamline(Y(:,:,1),Z(:,:,

4、1),ey,ez,SY,SZ);R=1; the= 0: pi/20 :2* pi;x= -0.2:0.11:3; y= -0.2:0.11:3; z= -0.2:0.11:3 ;X,Y, Z,T=ndgrid(y,x,z,the);r=sqrt(X-cos(T).2+(Y-sin(T).2+Z.2);dv = 1./r ;v=trapz(dv,4);ex,ey,ez=gradient(-v,0.5);XX,YY,ZZ=meshgrid( -0.2:0.11:3);SX,SY,SZ=meshgrid(0:0.5:3, 0:0.5:3, 0.05);subplot(2,2,3)Zmin=min(

5、Z(:); Zmax=max(Z(:);streamslice(XX,YY,ZZ,ex,ey,ez,0.01)subplot(2,2,4)x=cos(the);y=sin(the);z=zeros(1,length(the);plot3(x,y,z,LineWidth,3,Color,r);hold onh1=streamline(XX,YY,ZZ,ex,ey,ez,SX,SY,SZ) ;h2=copyobj(h1,gca); rotate(h2,1,0,0,180,0 0 0);h3=copyobj(allchild(gca),gca);rotate(h3,0,1,0,180,0 0 0);

6、view(-28,8)第五題：%亥姆霍茲線圈組成clear allM0=4*pi*1e-7; %初始化, 給定環(huán)半徑、電流、圖形I0=5.0; Rh=1;C0=M0 /(4* pi)* I0;Nx=21; Ny=21; % 設(shè)定觀測點(diǎn)網(wǎng)格數(shù)x=linspace(-Rh,Rh,Nx); %設(shè)定觀測點(diǎn)范圍及數(shù)組y=linspace(-Rh,Rh,Ny);Nh=20; %電流環(huán)分段T0= linspace( 0, 2*pi, Nh+1); % 環(huán)的圓周角分段T1= T0( 1: Nh) ;yb= Rh* cos( T1);zb= Rh* sin(T1); % 設(shè)置環(huán)各段向量的起點(diǎn)坐標(biāo)yb, zbT2

7、= T0( 2: Nh+ 1) ;ye= Rh* cos( T2) ;ze= Rh* sin( T2) ; %設(shè)置環(huán)各段向量的終點(diǎn)坐標(biāo)ye, zedlx=0; dly=ye- yb; dlz=ze- zb; % 計算環(huán)各段向量dl的三個長度分量xc=0; yc=( yb+ ye) /2; zc= ( zb+ ze) /2; %計算環(huán)各段向量中點(diǎn)的三個坐標(biāo)分量for i=1:Ny %循環(huán)計算各網(wǎng)格點(diǎn)上的B( x, y)值 for j=1:Nx rx=x(j) - xc; ry= y(i) - yc; rz= 0- zc; %觀測點(diǎn)在z= 0平面上 r3= sqrt( rx.2+ ry.2+ rz.2).3; dlXr_x=dly.*rz- dlz.*ry; %計算叉乘積dlX r的x 和y的分量 dlXr_y=dlz.*rx- dlx.*rz; Bx(i,j)=sum(C0*dlXr_x./r3); %把環(huán)各段產(chǎn)生的磁場分量累加 By(i,j)=sum(C0*dlXr_y./r3); endendBax= Bx(:,11:21) + Bx(:,1:11) ;Bay= By(:,11:21) + By(:,1:11) ;subpl