《matlab答案第五章(12页).doc》由会员分享,可在线阅读,更多相关《matlab答案第五章(12页).doc(12页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-第五章1.绘制下列曲线。(1) (2)(3) (4)答: (1) x=-10:0.1:10;y=100./(1+x.2);plot(x,y)(2) x=-10:0.1:10;y=1/(2*pi)*exp(-x.2/2);plot(x,y)(3) ezplot(x2+y2=1)(4) t=-10:0.1:10; x=t.2;y=5*t.3; plot(x,y)2. 绘制下列三维图形。(1)(2)(3)(4)半径为10的球面(提示:用sphere函数答: (1)t=0:pi/100:2*pi;x=cos(t);y=sin(t);z=t;plot3(x,y,z)(2)u=0:pi/100:2*pi
2、;v=0:pi/100:2*pi;x=(1+cos(u).*cos(v);y=(1+cos(u).*sin(v);z=sin(u);plot3(x,y,z)(3)x,y=meshgrid(-100:100);z=5*ones(size(x);mesh(x,y,z)(4) x,y,z = sphere; %等价于sphere(20) 20是画出来的球面的经纬分面数.20的话就是20个, 20个纬度,或者说“就是划出2020面的单位球体,把球的表面分成400个小网格”surf(x,y,z); %绘制单位球面x = 10*x;y = 10*y;z = 10*z;surf(x,y,z); %绘制半径为
3、10的球面axis equal3.在同一图形窗口中采用子图形式分别绘制正方形、圆、三角形和六边形。答: n=3 4 6 210for i=1:4 subplot(2,2,i) theta=linspace(pi/n(i),2*pi+pi/n(i),n(i)+1); plot(cos(theta),sin(theta); axis equal;end4. 分别用plot和fplot绘制下列分段函数的曲线。答:plot函数:方法一:x=linspace(-10,10,200);y=;for x0=xif x00y=y,x0.2+(1+x0).(1/4)+5;elseif x0=0y=y,0;els
4、eif x00y=y,x0.3+sqrt(1-x0)-5;endendplot(x,y)方法二:x=linspace(-10,10,200);y=(x0).*(x.2+(1+x).(1/4)+5);plot(x,y)用fplot函数:fplot(x0).*(x.2+(1+x).(1/4)+5),-10,10)第 12 页-5. 在同一坐标轴中绘制下列两条曲线并标注两曲线的交叉点。(1)(2)答: 相邻两点x1、x2 曲线(1)相应的y值为y11、y12, 曲线(2)y相应的y值y21、y22, 若(y11-y21)与(y12-y22)符号相反或等于0,则x1或x2就是一个交叉点。方法一:t=l
5、inspace(0,pi,300);x=sin(3*t).*cos(t);y1=2*x-0.5;y2=sin(3*t).*sin(t);len=length(t);p=y1-y2;x3=x(find(p(1:len-1).*p(2:len)=0);y3=2*x3-0.5;plot(x,y1,x,y2,x3,y3,rd)方法二:用循环%t=0:0.01:pi;t=linspace(0,pi,100);x=sin(3*t).*cos(t);y1=2*x-0.5;y2=sin(3*t).*sin(t);plot(x,y1,x,y2);grid;delta_y=y1-y2;% figure;stem(delta_y);j=1;for i=1:length(t)-1 if delta_y(i)*delta_y(i+1)=0 position(j)=i+1; j=j+1; endendhold on;plot(x(position),y1(position),ro);title(两条曲线相交);xlabel(x);ylabel(y);h = legend(y1,y2);