MATLAB程序设计与应用第5章应用题（刘卫国第3版）答案：
1.
% 绘制方程 (1)
figure
subplot(2,2,1)
x = -10:0.1:10;
y1 = 100./(1 + x.^2);
plot(x, y1);
title('y = 100/(1 + x^2)');
% 绘制方程 (2)
subplot(2,2,2)
x = -10:0.1:10;
y2 = 1/(2*pi) * exp(-x.^2/2);
plot(x, y2);
title('y = 1/(2*pi) * exp(-x^2/2)');
% 绘制方程 (3)
subplot(2,2,3)
theta = 0:0.01:2*pi;
x3 = cos(theta);
y3 = sin(theta);
plot(x3, y3);
title('x^2 + y^2 = 1');
% 绘制方程 (4)
subplot(2,2,4)
x = 0:0.1:10;
y4 = 5.*x.^(3/2);
plot(x, y4);
title('y = 5*x^{3/2}');
2.区别：使用fplot绘制y = sin(1/x)的更加密集，更符合实际
x = linspace(-10, 10, 1000); % 生成-10到10之间的1000个点，避免包括0
x(x==0) = []; % 移除x中的0值，因为sin(1/0)是未定义的
y = sin(1./x); % 计算对应的y值
% 绘制曲线
plot(x, y)
grid on % 显示网格
title('Using plot function for y = sin(1/x)')
xlabel('x')
ylabel('y')%fplot函数专门用于绘制数学函数的图形，它接受一个函数句柄作为输入，并自动在指定的x范围内计算并绘制函数值。
%这对于绘制复杂函数或不需要预先计算所有y值的情况特别有用。
% 使用fplot绘制y = sin(1/x)
figure
f = @(x) sin(1./x); % 定义函数句柄
fplot(f, [-10 10]) % 在-10到10的范围内绘制函数
grid on % 显示网格
title('Using fplot function for y = sin(1/x)')
xlabel('x')
ylabel('y')3.
figure
subplot(2,2,1)
t=0:0.01:2*pi;
r=5.*cos(t)+4;
polar(t,r)
subplot(2,2,2)
t1=0.01:0.01:2*pi;
r1=12./(sqrt(t1));
polar(t1,r1)
subplot(2,2,3)
r2=5.*cos(t)./-7;
polar(t,r2)
subplot(2,2,4)
r3=(pi.*(t.^2))/3;
polar(t,r3)4.
t=0:0.01:pi;
x=sin(3*t).*cos(t);
y=sin(3*t).*sin(t);
plot(x,y);
hold on
x1=-1:0.01:1;
y1=2*x1-0.5;
plot(x1,y1);
%%交点
[xi,yi]=polyxpoly(x,y,x1,y1)
plot(xi,yi,'*','MarkerSize',10,'MarkerFaceColor','y')
grid on5.
figure
subplot(2,2,1)
t=0:0.01:4*pi;
x=cos(t);
y=sin(t);
z=t;
plot3(x,y,z)
grid on
subplot(2,2,2)
u=0:0.01:4*pi;
v=0:0.01:4*pi;
x1=(1+cos(u)).*cos(v);
y1=(1+cos(u)).*sin(v);
z1=sin(u);
plot3(x1,y1,z1)
grid on
subplot(2,2,3)
% 定义X和Y的范围
x = linspace(-10, 10, 100); % 你可以根据需要调整范围和点数
y = linspace(-10, 10, 100);
% 使用meshgrid创建网格
[X, Y] = meshgrid(x, y);
% 设置Z为常数5
Z = 5 * ones(size(X));
% 使用surf绘制平面
surf(X, Y, Z, 'EdgeColor', 'none'); % 'EdgeColor', 'none' 用于移除网格线
shading interp; % 使颜色更加平滑
subplot(2,2,4)
[X,Y,Z]=sphere(100)
R=10;
surf(R*X,R*Y,R*Z,'EdgeColor', 'none');
shading interp;
axis equal