你好,可以发给我份2011数学建模a题matlab源程序吗?18954269189@163.com,非常感谢
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/27 17:01:27
![你好,可以发给我份2011数学建模a题matlab源程序吗?18954269189@163.com,非常感谢](/uploads/image/z/6206275-19-5.jpg?t=%E4%BD%A0%E5%A5%BD%2C%E5%8F%AF%E4%BB%A5%E5%8F%91%E7%BB%99%E6%88%91%E4%BB%BD2011%E6%95%B0%E5%AD%A6%E5%BB%BA%E6%A8%A1a%E9%A2%98matlab%E6%BA%90%E7%A8%8B%E5%BA%8F%E5%90%97%3F18954269189%40163.com%2C%E9%9D%9E%E5%B8%B8%E6%84%9F%E8%B0%A2)
你好,可以发给我份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