2.4 0.9 0.6; 2.4 1.5 0.9; 2.4 2.1 1.1; 2.4 2.7 1.1; 5 0.3 0.4; 5 0.9 0.4; 5 1.5 0.5; 5 2.1 0.6; 5 2.7 0.6; 7.2 0.3 0.4; 7.2 0.9 0.2; 7.2 1.5 0.2; 7.2 2.1 0.2; 7.2 2.7 0.2;];
x=A(:,1);y=A(:,2);z=A(:,3); scatter(x,y,5,z)%散点图 figure
[X,Y,Z]=griddata(x,y,z,linspace(2.4,7.2)',linspace(0.3,2.7),'v4');%插值 pcolor(X,Y,Z);shading interp%伪彩色图 figure,contourf(X,Y,Z) %等高线图 figure,surf(X,Y,Z)%三维曲面 ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)]
clear%根据附件一做出的通道三流场分布图 clc
A=[2.4 0.3 0.4; 2.4 0.9 0.6; 2.4 1.5 0.7; 2.4 2.1 0.8; 2.4 2.7 0.9; 5 0.3 0.4; 5 0.9 0.5; 5 1.5 0.6; 5 2.1 0.7; 5 2.7 0.6; 7.2 0.3 0.4; 7.2 0.9 0.6; 7.2 1.5 0.6; 7.2 2.1 0.6; 7.2 2.7 0.5;];
x=A(:,1);y=A(:,2);z=A(:,3); scatter(x,y,5,z)%散点图
20
figure
[X,Y,Z]=griddata(x,y,z,linspace(2.4,7.2)',linspace(0.3,2.7),'v4');%插值 pcolor(X,Y,Z);shading interp%伪彩色图 figure,contourf(X,Y,Z) %等高线图 figure,surf(X,Y,Z)%三维曲面 ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] 附录Ⅱ：
clf%做出每个机柜任务量为0.5时每个高度温度分布图 clear
A=load('data41.txt');%%读入矩阵 for n=4:1:7
x=A(:,1);y=A(:,2);z=A(:,n);
[X,Y,Z]=griddata(x,y,z,linspace(0,9.2)',linspace(2.1,7.1),'v4');%插值 pcolor(X,Y,Z);
shading interp%伪彩色图 figure,surf(X,Y,Z) figure
ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] end
clf%做出每个机柜任务量为0.2时每个高度温度分布图 clear
A=load('data49.txt');%%读入矩阵 for n=4:1:7
x=A(:,1);y=A(:,2);z=A(:,n);
[X,Y,Z]=griddata(x,y,z,linspace(0,9.2)',linspace(2.1,7.1),'v4');%插值 pcolor(X,Y,Z);
shading interp%伪彩色图 figure,surf(X,Y,Z) figure
ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] end
clf%做出每个机柜任务量为0.3时每个高度温度分布图 clear
A=load('data412.txt');%%读入矩阵 for n=4:1:7
x=A(:,1);y=A(:,2);z=A(:,n);
21
[X,Y,Z]=griddata(x,y,z,linspace(0,9.2)',linspace(2.1,7.1),'v4');%插值 pcolor(X,Y,Z);
shading interp%伪彩色图 figure,surf(X,Y,Z) figure
ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] end
附录Ⅲ：
clf%作出温度极值点所在通道温度分布图 clear
A=load('data51.txt');%%读入矩阵 for n=4:1:6
x=A(:,2);y=A(:,3);z=A(:,n);
[X,Y,Z]=griddata(x,y,z,linspace(2.1,7.1)',linspace(0.2,2.6),'v4');%插值 pcolor(X,Y,Z);
shading interp%伪彩色图 figure,surf(X,Y,Z) figure
ii=find(imregionalmax(Z)); Imregionalmax=Z(ii) Points=[X(ii) Y(ii)] end
附录Ⅳ：
clf%温度极值点1温度与相关机柜任务量的拟合 clear clc
k=[0.2,0.3,0.3,0.3,0.3,0.5,0.5,0.5,0.5,0.5,0.8]';
y=[28.1,30.2,30,30.1,30.1,32.1,32.4,32.2,32.1,31.9,35.9]'; plot(k,y,'.') hold on
fun=inline('b(1).*b(2).^k','b','k'); b0=[1,1];
[b,r,j]=nlinfit(k,y,fun,b0); b
R=sum(r.^2) k=0:0.1:1;
plot(k,b(1).*b(2).^k,'-')
clf%温度极值点2温度与相关机柜任务量的拟合 clear
22
clc
k=[0.5,0.5,0.8,0.5,0.5,0.5,0.5,0.2,0.2,0.2,0.3]';
y=[32.0,32.0,35.9,32.0,32.0,32.0,32.0,28.1,28.0,28.0,30.0]'; plot(k,y,'.') hold on
fun=inline('b(1).*b(2).^k','b','k'); b0=[1,1];
[b,r,j]=nlinfit(k,y,fun,b0); b
R=sum(r.^2) k=0:0.1:1;
plot(k,b(1).*b(2).^k,'-')
clf%温度极值点3温度与相关机柜任务量的拟合 clear clc
k=[0.5,0.5,0.5,0.5,0.8,0.5,0.5,0.5,0.2,0.2,0.3,0.3]';
y=[32.2,32.2,32.2,32.6,35.8,32.2,32.2,32.2,28.1,28.1,30.2,30.2]'; plot(k,y,'.') hold on
fun=inline('b(1).*b(2).^k','b','k'); b0=[1,1];
[b,r,j]=nlinfit(k,y,fun,b0); b
R=sum(r.^2) k=0:0.1:1;
plot(k,b(1).*b(2).^k,'-')
附录Ⅴ：
clear%任务量为0.8时线性规划 clc
x1=1;x2=1;x3=1;x4=1; x11=0;x12=0;x13=0;x14=0;
A=[26.7333*1.4483^x1,26.0368*1.5055^x2,26.4856*1.4760^x4]; maxint=max(A); for x1=0:0.01:1 for x2=0:0.01:1 for x3=0:0.01:1 x4=3.2-x1-x2-x3;
A1=[26.7333*1.4483^x1,26.0368*1.5055^x2,26.4856*1.4760^x4]; if x4<1&x4>0
if max(A1) 23 x12=x2; x13=x3; x14=x4; end end end end end maxint x11 x12 x13 x14 clear%任务量为0.5时线性规划 clc x1=1;x2=1;x3=1;x4=1; x11=0;x12=0;x13=0;x14=0; A=[26.7333*1.4483^x1,26.0368*1.5055^x2,26.4856*1.4760^x4]; maxint=max(A); for x1=0:0.01:1 for x2=0:0.01:1 for x3=0:0.01:1 x4=2-x1-x2-x3; A1=[26.7333*1.4483^x1,26.0368*1.5055^x2,26.4856*1.4760^x4]; if x4<1&x4>0 if max(A1) 附录Ⅵ： 24 clf clear clc A=load('data81.txt');%%读入矩阵 for n=3:1:14 x=A(:,1);y=A(:,2);z=A(:,n); [X,Y,Z]=griddata(x,y,z,linspace(1.1,8)',linspace(2.1,7.1),'v4');%插值 surf(X,Y,Z)%三维曲面 figure,shading interp%伪彩色图 end 25