作业帮你好,可以发给我份2011数学建模a题matlab源程序吗?18954269189@163.com,非常感谢
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/06 04:27:30
你好,可以发给我份2011数学建模a题matlab源程序吗?18954269189@163.com,非常感谢
你好,可以发给我份2011数学建模a题matlab源程序吗?18954269189@163.com,非常感谢
你好,可以发给我份2011数学建模a题matlab源程序吗?18954269189@163.com,非常感谢
clear
close all
clc
C=xlsread('cumcm2011A附件_数据','附件2','B04:I322');
A=xlsread('cumcm2011A附件_数据','附件1','B04:E322');
X=A(:,1); Y=A(:,2);Z=A(:,3);F=A(:,4);
XI=linspace(min(X),max(X),100);
YI=linspace(min(Y),max(Y),100);
[XI,YI]=meshgrid(XI,YI);
ZI=griddata(X,Y,Z,XI,YI,'cubic');
FI=griddata(X,Y,F,XI,YI,'nearest');
%作该城市地势图
figure
meshz(XI,YI,ZI);
xlabel('X');
ylabel('Y');
zlabel('海拔');
%作该城市功能分区图
figure
h=surf(XI,YI,ZI);
set(h,'cdata',FI);
colormap cool
colorbar
axis([min(X),max(X),min(Y),max(Y)]);
for i=1:8
%在地势等高线上加各中金属的浓度,和所属分区
figure
grid on
hold on
color={'k','r','y','c','b'};
CI=griddata(X,Y,C(:,i),XI,YI,'cubic');
contour(XI,YI,ZI,[50]);
for j=1:5
27
loc=F==j;
plot(X(loc),Y(loc),'o','markerfacecolor',color{j},'MarkerEdgeColor',color{j
});
end
str={'地势等高线','生活区','工业区','山区','主干道路区','公园绿地区'};
legend(str,'location','best')
for k=1:length(X)
text(X(k)-200,Y(k)+200,num2str(C(k,i)),'fontsize',8);
end
xlabel('X');
ylabel('Y');
hold off
%插值做有效区间的浓度变化示意图
figure
grid on
mesh(XI,YI,CI);
xlabel('X');
ylabel('Y');
zlabel('浓度');
end