2011数学建模A题答案

1 413 359 2 403 343 3 383.5 351 4 381 377.5 5 339 376 6 335 383 7 317 362 8 334.5 353.5 9 333 342 10 282 325 11 247 301 12 219 316 13 225 270 14 280 292 15 290 335 16 337 328 17 415 335 18 432 371 19 418 374 20 444 394 21 251 277 22 234 271 23 225 265 24 212 290 25 227 300 26 256 301 27 250.5 306 28 243 328 29 246 337 30 314 367 31 315 351 32 326 355 33 327 350 34 328 342.5 35 336 339 36 336 334 37 331 335 38 371 330 39 371 333 40 388.5 330.5 41 411 327.5 42 419 344 43 411 343 44 394 346 45 342 342 46 342 348 47 325 372 48 315 374 49 342 372 50 345 382 51 348.5 380.5 52 351 377 53 348 369 54 370 363 55 371 353 56 354 374 57 363 382.5 58 357 387 59 351 382 60 369 388 61 335 395 62 381 381 63 391 375 64 392 366 65 395 361 66 398 362 67 401 359 68 405 360 69 410 355 70 408 350 71 415 351 72 418 347 73 422 354 74 418.5 356 75 405.5 364.5 76 405 368 77 409 370 78 417 364 79 420 370 80 424 372 81 438 368 82 438.5 373 83 434 376 84 438 385 85 440 392 86 447 392 87 448 381 88 444.5 383 89 441 385 90 440.5 381.5 91 445 380 92 444 360 clear,clc count=319;%count±íʾȡÑùµãµÄ¸öÊý %fidheng=fopen('hengzuobiao.txt','r'); %heng=zeros(1,count); %[heng,count]=fread(fidheng,inf); x=load('xzuobiao.txt'); x=x'; y=load('yzuobiao.txt'); y=y'; z=load('zzuobiao.txt'); z=z'; as=load('as.txt'); as=as'; cd=load('cd.txt'); cd=cd'; cr=load('cr.txt'); cr=cr'; cu=load('cu.txt'); cu=cu'; hg=load('hg.txt'); hg=hg'; ni=load('ni.txt'); ni=ni'; pb=load('pb.txt'); pb=pb'; zn=load('zn.txt'); zn=zn'; gongneng=load('gongnengqu.txt'); gongneng=gongneng'; %plot3(x,y,as); hold on As=ones(size(as'))*as; [X,Y]=meshgrid(x,y); x1=linspace(min(x),max(x)); y1=linspace(min(y),max(y)); [X1,Y1]=meshgrid(x1,y1); %As1=griddata(X,Y,As,X1,Y1) %As1=interp2(x,y,As,X,Y); bianhao=1:1:319; %ԭʼÊý¾Ý´¦Àí maxas=max(as); minas=min(as); as=(as-minas)/(maxas-minas); maxcd=max(cd); mincd=min(cd); cd=(cd-mincd)/(maxcd-mincd); maxcr=max(cr); mincr=min(cr); cr=(cr-mincr)/(maxcr-mincr); maxcu=max(cu); mincu=min(cu); cu=(cu-mincu)/(maxcu-mincu); maxhg=max(hg); minhg=min(hg); hg=(hg-minhg)/(maxhg-minhg); maxni=max(ni); minni=min(ni); ni=(ni-minni)/(maxni-minni); maxpb=max(pb); minpb=min(pb); pb=(pb-minpb)/(maxpb-minpb); maxzn=max(zn) minzn=min(zn) zn=(zn-minzn)/(maxzn-minzn); %plot(bianhao,as,bianhao,cd,bianhao,cr,bianhao,cu,bianhao,hg,bianhao,ni,bianhao,pb,bianhao,zn) %mesh(X1,Y1,As1); subplot(4,2,1); plot(bianhao,as); subplot(4,2,2); plot(bianhao,cd); subplot(4,2,3); plot(bianhao,cr); subplot(4,2,4); plot(bianhao,cu); subplot(4,2,5); plot(bianhao,hg); subplot(4,2,6); plot(bianhao,ni); subplot(4,2,7); plot(bianhao,pb); subplot(4,2,8); plot(bianhao,zn); data=[x;y;as]'; data1=[x;y;gongneng]'; x=data(:,1); y=data(:,2); as=data(:,3); gongz=data1(:,3); %scatter(x,y,z);%É¢µãͼ figure [xx,yy,zz]=griddata(x,y,as,linspace(min(x),max(x))',linspace(min(y),max(y)),'v4');%²åÖµ [xx,yy,zz]=griddata(x,y,gongz,linspace(min(x),max(x))',linspace(min(y),max(y)),'v4');%²åÖµ daxiao=size(zz); zz=zz-3.6; for i=1:100 for j=1:100 if zz(i,j)<0 zz(i,j)=0; end end end 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);%²åÖµ gongnengi=griddata(x,y,gongneng,xi,yi); marker={'^','*','s','p','o'}; color={'r','b','k','c','y'}; mat={'As','Cd','Cr','Cu','Hg','Ni','Pb','Zn'}; str={'µÈ¸ßÏß','Éú»îÇø','¹¤ÒµÇø','ɽÁÖÇø','½»Í¨Çø','Â̵ØÇø'}; hold on figure %contourf(xi,yi,zi,0:10:500,'b-'); %figure,surf(xx,yy,zz)