matlab实验-旋转体

求曲线 y=x^(1/4)*sin⁡x^2(0≤x≤π^(1/2))与 x 轴所围成图形分别绕x 轴、y 轴旋转所成的旋转体体积，并画出两个旋转体的图形。

为了计算曲线 y=x^1/4*sin(x^2) 在区间 [0,π​] 与 x 轴所围成的图形绕 x 轴和 y 轴旋转所形成的旋转体的体积，一般使用 MATLAB 中的数值积分方法。旋转体体积可以通过定积分来计算。

当图形绕 x 轴旋转时，体积的计算公式为：

而绕 y 轴旋转时，体积的计算公式为：

其中 a=0 和 b=​。

MATLAB 代码：

% 定义函数 y = x^(1/4) * sin(x^2)
y = @(x) x.^(1/4) .* sin(x.^2);

% 定义积分区间
a = 0;
b = sqrt(pi);

% 计算绕 x 轴旋转的旋转体体积
Vx = pi * integral(@(x) y(x).^2, a, b);

% 计算绕 y 轴旋转的旋转体体积
Vy = 2 * pi * integral(@(x) x.*y(x), a, b);

% 输出结果
fprintf('绕 x 轴旋转的旋转体体积：%.4f\n', Vx);
fprintf('绕 y 轴旋转的旋转体体积：%.4f\n', Vy);

% 绘制旋转体
fplot(@(x) y(x), [a, b]);
xlabel('x');
ylabel('y');
title('曲线 y = x^{1/4} * sin(x^2)');

% 创建旋转体的图形（绕 x 轴和 y 轴）
figure;
f = @(x,z) x.^(1/4) .* sin(x.^2);
fsurf(@(x,z) x, @(x,z) f(x,z).*cos(z), @(x,z) f(x,z).*sin(z), [a, b, 0, 2*pi]);
title('绕 x 轴旋转');

figure;
fsurf(@(x,z) f(x,z).*cos(z), @(x,z) f(x,z).*sin(z), @(x,z) x, [a, b, 0, 2*pi]);
title('绕 y 轴旋转');

输出图形和结果：