Question

MATLAB程序求解

我需要这个问题的MATLLAB的程序





Answer 1

今天中午，看到你还有这个提问，随便继续昨天晚上太晚未完成的工作，三维曲线图效果如下（图片只给你上传了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

 随意问一下，这两天给你们做的这些题，是不是你们几个都是一个班的？那个同学，还这么为你们操心，到处找人帮你们做。

希望对你有帮助

Answer 2