画7/8圆球，matlab

zhuhuawe

close  all;   clc;      clear all;
% % % % % % % % % % % % % % %   下半圆
t=linspace(pi/2,pi,100);
p=linspace(0,2*pi,100);
[theta,phi]=meshgrid(t,p);
x=sin(theta).*sin(phi);
y=sin(theta).*cos(phi);
z=cos(theta);
[m,n]=size(z);
re=[255/256 255/256 0/256];%改数字 改变颜色、数字为0-1
colormap(re)
surf(x,y,z);
% shading interp
shading flat
axis equal;
% % % % % % % % % % % % % % % 上四分之三圆
hold on
t=linspace(0,pi/2,100);
p=linspace(pi/2,2*pi,100);
[theta,phi]=meshgrid(t,p);
x=sin(theta).*sin(phi);
y=sin(theta).*cos(phi);
z=cos(theta);
[m,n]=size(z);
re=[255/256 255/256 0/256];%改数字 改变颜色、数字为0-1
colormap(re)
surf(x,y,z);
% shading interp
shading flat
axis equal;
% % % % % % % % % % % % % % % % 内圆环
R = 1; r = 0.5;
theta=linspace(0,pi/2,90);ph=linspace(r,R,30);
[t,p]=meshgrid(theta,ph);
r=t*0;
[x,y,z]=pol2cart(t,p,r);
surf(x,z,y,'edgecolor','none');
re=[255/256 255/256 0/256];
colormap(re)
R = 1; r = 0.5;
theta=linspace(0,pi/2,90);ph=linspace(r,R,30);
[t,p]=meshgrid(theta,ph);
r=t*0;
[x,y,z]=pol2cart(t,p,r);
surf(x,y,z,'edgecolor','none');
re=[255/256 255/256 0/256];
colormap(re)
R = 1; r = 0.5;
theta=linspace(0,pi/2,90);ph=linspace(r,R,30);
[t,p]=meshgrid(theta,ph);
r=t*0;
[x,y,z]=pol2cart(t,p,r);
surf(z,y,x,'edgecolor','none');
re=[255/256 255/256 0/256];
colormap(re)
% % % % % % % % % % % % % % % %

hold on
t=linspace(0,1/2*pi,1000);
x=1/2*cos(t);
y=zeros(1,1000);
z=1/2*sin(t);
fill3([0 x],[0 y],[0 z],[0/256 199/256 255/256],'edgecolor','none')
re=[255/256 255/256 0/256];
colormap(re)
x=1/2*cos(t);
y=1/2*sin(t);
z=zeros(1,1000);
fill3([0 x],[0 y],[0 z],[0/256 199/256 255/256],'edgecolor','none')
x=zeros(1,1000);
y=1/2*sin(t);
z=1/2*cos(t);
fill3([0 x],[0 y],[0 z],[0/256 199/256 255/256],'edgecolor','none')
axis off
% % % % % % % % % % % % % % % % % % % % % %
light

zhuhuawe

huatiantian

very good thanks  useful