楼主试试下述语句,其中nTmax就是最大的T对应的n值(近似值):
n=linspace(600,4000,1001);
T=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;
[Tmax ind]=max(T); %求T的最大值(近似值)
nTmax = n(ind); %求T的最大值对应的n值(近似值)
plot(n,T,'linewidth',2) %绘曲线T=T(n)
hold on
plot(nTmax,Tmax,'p','markersize',8,'color','r')%绘最大值对应的点
首先这样编程
n=600:0.6pi:4000
T=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4
从而可以确定n值