N=1;
for t=0:0.01:50
h(N)=3*sin((5*pi/32)*t);
N=N+1;
end
N=1;
for t=0:0.01:50
c(N)=(0.4/(1+log(1/(30000*h(N)))))^2;
u=(375/8*pi)*cos((5*pi/32)*t);
y(N)=1025*c(N)*u^2;
N=N+1;
end
aa=max(c);
bb=min(c);
cc=max(y);
dd=min(y);
%%[aa,bb]是c的定义域 [cc,dd]是y的定义域
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024