《matlab答案第五章.doc》由会员分享，可在线阅读，更多相关《matlab答案第五章.doc（8页珍藏版）》请在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;el
4、seif 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)8), ( ) ( ; )线 ),) ) ( + ) =) _ _ -)( : )_ ( ;, , )( ) 0*) . )0,0 .0循用 , 0 )(*=( .)( 数 )+)(. ()=(- ( (.0)0 0 二 ( - ( 0 0= )0 0 一数函答线的数列绘 和 用 ) (
5、 ( ) ) ,) ( = ,( :=0 :形六和圆方绘分形中窗图 球的径半% , ; = 球单%下 ,线格网0分的，球面 出就“度纬,经 _ = _ & % . : 话 数面面来出 ) % = ) , ;) *0: ( = ,)()( . ) ) *0/ 0 ( ) *0 : 函 （球为径 形图下绘 ,( . .0 .0)= ).- /0 0 ,)+ .0:. - 5. 在同一坐标轴中绘制下列两条曲线并标注两曲线的交叉点。（1）（2）答: 相邻两点x1、x2 曲线（1）相应的y值为y11、y12, 曲线（2）y相应的y值y21、y22, 若(y11-y21)与(y12-y22)符号相反或等于
6、0，则x1或x2就是一个交叉点。方法一:t=linspace(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);