第三章作业。
第一题：format long
syms x
c1=int((x^4*(1-x)^4)/(1+x^2),0,1); 符号工具箱int计算。
c1=eval(c1指令eval一般应用于int后结果较复杂时，进一步求数值解。
f=@(x)(x.^4.*(1-x).^4)./1+x.^2);
c2=quad(f,0,1数值积分指令quad要求积分限为常数
第二题：format long
f=@(x)((2/sqrt(pi))*exp(-x.^2定义误差函数。
for k=0.1:0.1:1.0
c1=quad(f,0,k);
c2=erf(k);
c=abs(c1-c2);
bar(k,c)
hold on
end 第三题：
vp=1578;
f=@(v)((4/sqrt(pi))*v.^2)/vp^3).*exp((-v.^2)/vp^2));
c1=quad(f,0,vp)
c2=quad(f,0,3.3*vp)
c3=quad(f,3e4,3e8)
结果：c1 =0.427593285448571
c2 =0.999927444001040 在0~3.3间的分子数接近为全部。
c3 =3.520658577050803e-151 速度很大的分子数很少。
第四题：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]=meshgrid(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(:,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(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 on
h1=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]);
view(-28,8)
第五题：亥姆霍兹线圈组成。
clear all
m0=4*pi*1e-7; %初始化， 给定环半径、电流、图形。
i0=5.0; rh=1;
c0=m0 /(4* pi)* i0;
nx=21; ny=21; %设定观测点网格数。
x=linspace(-rh,rh,nx); 设定观测点范围及数组。
y=linspace(-rh,rh,ny);
nh=20; %电流环分段。
t0= linspace( 0, 2*pi, nh+1); 环的圆周角分段。
t1= t0( 1: nh) ;
yb= rh* cos( t1);
zb= rh* sin(t1); 设置环各段向量的起点坐标yb, zb
t2= t0( 2: nh+ 1) ;
ye= rh* cos( t2) ;
ze= rh* sin( t2) ;设置环各段向量的终点坐标ye, ze
dlx=0; dly=ye- yb; dlz=ze- zb; %计算环各段向量dl的三个长度分量。
xc=0; yc=( yb+ ye) /2; zc= (zb+ ze) /2; %计算环各段向量中点的三个坐标分量。
for i=1:ny %循环计算各网格点上的b( x, y)值。
for j=1:nx
rx=x(j) -xc; ry= y(i) -yc; rz= 0- zc; %观测点在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); 把环各段产生的磁场分量累加。
by(i,j)=sum(c0*dlxr_y./r3);
endend
bax= bx(:,11:21) +bx(:,1:11) ;
bay= by(:,11:21) +by(:,1:11) ;
subplot(1,2,1);
mesh(x(11:21),y,bax); xlabel(x); ylabel(y) ;画出其b分布三维图。
subplot(1,2,2);
plot(y,bax), grid, xlabel(y); ylabel(bx);
第三章作业
v s 顺序执行下述两个动作 1.s值加1，即s s 1 2.如果s 0，则该进程继续运行 3.如果s 0，则唤醒等待信号量s阻塞队列中的头一个进程 把阻塞态改为就绪态 执行v操作的进程继续运行。procedure s var s semaphore begin s s 1 if s 0 then ...
第三章作业
1 顺序栈空 栈满条件2 链栈栈空 栈满条件。3 循环队列队空 队满条件，如何表示队列中数据元素的个数4 链队列队空 队满条件。5 以下运算实现在顺序栈上的进栈，请在 处用适当的语句予以填充。int push sqstacktp sq,datatype x if sp top sqstack max...
第三章作业
1.论述各类绿地的环境特点和树种的选择。一 高层建筑中的狭窄街巷绿地绿地内的环境特点 直射辐射量少，日照时间短 夏季气温偏低，冬季因受周围建筑物热辐射的影响，气温偏高 风速一般偏低，但有时会产生狭管效应，使风速增大。这些地方裸露土面极少，多为水泥铺装，严重阻碍了土壤与大气的水 气交换，且存在一定程度...