今天中午，看到你还有这个提问，随便继续昨天晚上太晚未完成的工作，三维曲线图效果如下（图片只给你上传了H/D=5的情况）：
计算代码如下（附件便于你自己验证，复制好像不自动分段）：
%第二问：画出区间300<=Ra<=10000上H/D=5，10，15，20的Sopt图形
clear
clc
HD=[5 10 15 20];%表示H/D
options=optimset('display','off');
dRa=100;%Ra每次增加的大小
n1=length(300:dRa:10000);
n2=length(1:0.5:10);
Reslut_Sopt=zeros(n1,n2);
for i=1:4
k=0;
for D=1:0.5:10
flag=0;
k=k+1;
Hi=D*HD;
H=Hi(i);%分别计算H/D的五种情况
for Ra=300:dRa:10000
flag=flag+1;
fx=@(Sopt)Sopt./D*(2+Sopt/D)./(1+Sopt/D)^(2/3)/D-2.75*(H/D).^(1/3).*Ra.^(-1/4);
x0=1;%初始点
[sopt,fval]=fsolve(fx,x0,options);
Reslut_Sopt(flag,k)=sopt;
end
end
figure
Dx=1:0.5:10;
Ray=300:dRa:10000;
[X,Y] =meshgrid(Dx,Ray);
hold on
surf(X,Y,Reslut_Sopt)%三维曲面
xlabel('D');
ylabel('Ra');
zlabel('Sopt');
title(['H/D=',num2str(HD(i))])
grid on
end
随意问一下，这两天给你们做的这些题，是不是你们几个都是一个班的？那个同学，还这么为你们操心，到处找人帮你们做。
希望对你有帮助
