MATLAB程序设计与应用刘卫国（第三版）课后实验答案——13

一

h&＃61;figure(&＃39;MenuBar&＃39;,&＃39;figure&＃39;,&＃39;color&＃39;,&＃39;r&＃39;,&＃39;WindowButtonDownFcn&＃39;,&＃39;disp("Left Button Pressed")&＃39;)


二

x&＃61;-2:0.01:2;
y&＃61;1/2.*log(x&＃43;sqrt(1&＃43;x.^2));
h&＃61;line(x,y);
set(h,&＃39;color&＃39;,&＃39;r&＃39;,&＃39;linestyle&＃39;,&＃39;:&＃39;,&＃39;linewidth&＃39;,2)
text(1,exp(2),&＃39;y&＃61;1/2*log(x&＃43;sqrt(1&＃43;x^2))&＃39;)t&＃61;-2:0.01:2;
x&＃61;t.^2;
y&＃61;5.*t.^3;
h1&＃61;line(t,x);
set(h1,&＃39;color&＃39;,&＃39;r&＃39;,&＃39;linestyle&＃39;,&＃39;:&＃39;,&＃39;linewidth&＃39;,2)
text(-1.5,0,&＃39;x&＃61;t.^2&＃39;)
hold on
h2&＃61;line(t,y);
set(h2,&＃39;color&＃39;,&＃39;b&＃39;,&＃39;linestyle&＃39;,&＃39;:&＃39;,&＃39;linewidth&＃39;,3)
text(1.5,20,&＃39;y&＃61;5.*t.^3&＃39;)x&＃61;-2:0.01:2;
y&＃61;(1&＃43;x.^2)./(1&＃43;x.^4);
h&＃61;line(x,y);
set(h,&＃39;color&＃39;,&＃39;r&＃39;,&＃39;linestyle&＃39;,&＃39;:&＃39;,&＃39;linewidth&＃39;,2)
text(1,1,&＃39;y&＃61;(1&＃43;x.^2)./(1&＃43;x.^4)&＃39;)x&＃61;-2:0.01:2;
y&＃61;x.^2.*exp(2*x);
h&＃61;line(x,y);
set(h,&＃39;color&＃39;,&＃39;r&＃39;,&＃39;linestyle&＃39;,&＃39;:&＃39;,&＃39;linewidth&＃39;,2)
text(1,exp(2),&＃39;y&＃61;x^2*exp(2*x)&＃39;)


三

t&＃61;-20*pi:0.1*pi:20*pi;
[t,~]&＃61;meshgrid(t);
x&＃61;t;
y&＃61;2.*t.*cos(t);
z&＃61;5.*t.*sin(t);
axes(&＃39;view&＃39;,[-37.5,30]);
surface(x,y,z)
axis equaltheta&＃61;0:0.01*pi:2*pi;
t&＃61;0:0.01*pi:2*pi;
[theta,t]&＃61;meshgrid(theta,t);
x&＃61;2.*cos(theta);
y&＃61;2.*sin(theta);
z&＃61;t;
axes(&＃39;view&＃39;,[-37.5,30]);
surface(x,y,z)t&＃61;0:0.00001:0.001;
[t,x]&＃61;meshgrid(t);
v&＃61;10*exp(-0.01*x).*sin(2000*pi*t-0.2*x&＃43;pi);
axes(&＃39;view&＃39;,[-37.5,30]);
h&＃61;surface(t,x,v);
title(&＃39;v&＃61;10*exp(-0.01*x).*sin(2000*pi*t-0.2*x&＃43;pi)&＃39;);
xlabel(&＃39;t&＃39;),ylabel(&＃39;x&＃39;),zlabel(&＃39;v&＃39;)y&＃61;0:0.01:20;
x&＃61;0:0.01:20;
[x,y]&＃61;meshgrid(x,y);
z&＃61;y.^3;
axes(&＃39;view&＃39;,[-37.5,30]);
h&＃61;surface(x,y,z);


四

x&＃61;0:0.01:2*pi;
y1&＃61;sin(x);
y2&＃61;cos(x);
y3&＃61;tan(x);
y4&＃61;cot(x);
subplot(2,2,1);
plot(x,y1);
subplot(2,2,2);
plot(x,y2);
subplot(2,2,3);
plot(x,y3);
subplot(2,2,4);
plot(x,y4);


五

t&＃61;linspace(0,2*pi,6);
k&＃61;t(1:length(t)-1);
patch(sin(k),cos(k),1:length(k),&＃39;facecolor&＃39;,&＃39;interp&＃39;);
colormap summer;
axis equal;


六

r&＃61;150;
a&＃61;150;
b&＃61;150;
para1&＃61;[a-r,b-r,2*r,2*r];
para2&＃61;[a-r&＃43;300,b-r,2*r,2*r]
rectangle(&＃39;Position&＃39;,para1,&＃39;Curvature&＃39;,[1 1],&＃39;EdgeColor&＃39;,&＃39;r&＃39;)
hold on
rectangle(&＃39;Position&＃39;,para2,&＃39;Curvature&＃39;,[1 1],&＃39;EdgeColor&＃39;,&＃39;r&＃39;)
axis equal


注&＃xff1a;考试范围仅到此&＃xff0c;所以simulaition仿真及之后的章节没有作答。由于上次挂了个网址被删啦&＃xff0c;若读者需要&＃xff0c;请自行查阅相关资料&＃xff08;话说csdn不给返回底稿的吗…&＃xff09;