下段代码报错：数组的大小不兼容。如何修改：x1 = sqrt(R*R-r1.*r1).*cos(phai); y1 = sqrt(R*R-r1.*r1).*sin(phai);

这段代码中出现数组大小不兼容的问题，是因为 `r1` 和 `phai` 的长度不同。`r1` 的长度为 `(R/0.01)+1`，而 `phai` 的长度为 `(2*pi/0.01)+1`，两者长度不同，导致在执行 `cos(phai)` 和 `sin(phai)` 时出现了错误。

解决这个问题的方法是，将 `r1` 和 `phai` 的长度保持一致。由于 `phai` 的范围是从 `0` 到 `2*pi`，所以可以将 `r1` 的长度修改为 `length(phai)`，即：

r1 = linspace(0, R, length(phai));

这样，`r1` 的长度就和 `phai` 相同了，即 `(2*pi/0.01)+1`。修改后的代码如下：

clear all; 
clc; 

% 初始值
h = 8; 
R = 3; 
phai = 0:0.01:2*pi; 
r1 = linspace(0, R, length(phai)); % 修改此处 

% 球部分
x1 = sqrt(R*R-r1.*r1).*cos(phai); 
y1 = sqrt(R*R-r1.*r1).*sin(phai); 
z1 = r1+h; 

% 圆锥部分
r2 = 0:0.01:R*h/h; 
x2 = (R/h)*(r2+h).*cos(phai); 
y2 = (R/h)*(r2+h).*sin(phai); 
z2 = r2+h; 

% 绘制图像
figure 
plot3(x1,y1,z1,'LineWidth',2); 
hold on; 
plot3(x2,y2,z2,'LineWidth',2); 
xlabel('x'); 
ylabel('y'); 
zlabel('z'); 
title('降落伞');

% 添加背景色和光照效果
set(gcf,'Color',[0.5,0.5,0.5]); 
light('Position',[0,0,10],'Style','local'); 
lighting gouraud;

运行以上代码，即可得到正确的降落伞图像。