[数学建模]三维图像

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/qq_37360631/article/details/86623787

function y=mfun(v0)
pi=3.14159;
a=18;
g=9.8;      
m=75;
u=0.02;
k=0.5;
v0=60;
w1=0.5;
w2=0.5;
u_front=0.8;
u_back=0.01;
%a=18;
L=2;
S=5;
r=2.8;
R=0.4;
qq=exp(-pi.*u);
h_min=1;            %最低是1m
h_max=8;
a_min=180;          %基础动作180    
a_max=1260;         %高难度1260
E0=1./2.*m.*v0.*v0;

%plot(v0,h,'b','LineWidth',2);
%h,L,r
r=1:1:10;
L=1:1:10;
[xx,yy]=meshgrid(r,L);                      %构成格点矩阵

qqq=m.*g.*xx.*cosd(a);


Eb=(E0.*qq)-qqq.*(u.*qq-2.*u.*u-1)./(4.*u.*u+1);
E_edge=(1-2.*u.*yy./r).*((Eb-u.*m.*g.*S).*qq-((2.*u.*u+1).*qq+u).*qqq./(4.*u.*u+1))-m.*g.*yy.*cosd(a);
h=(1/2*k)*log(1+2*k*E_edge/m/g)+1;

z=h;
surf(xx,yy,z);
shading interp;             %使曲线变得光滑

xlabel('r');
ylabel('L');
ylabel('Maximum vertical height(h)');

end