数学建模第四次优秀论文研读报告
- 格式:doc
- 大小:4.14 MB
- 文档页数:21
特等奖论文1、用到知识clearload shuju0.txtload shuju1.txtn1=load('1.txt') %生活区横纵坐标数据n2=load('2.txt') %工业区横纵坐标数据n3=load('3.txt') %山区横纵坐标数据n4=load('4.txt') %交通区横纵坐标数据n5=load('5.txt') %公园绿地区横纵坐标数据plot(n1(:,1),n1(:,2),'*') %生活区位置图hold onplot(n2(:,1),n2(:,2),'.') %工业区位置图hold onplot(n3(:,1),n3(:,2),'^') %山区位置图hold onplot(n4(:,1),n4(:,2),'+') %交通区位置图hold onplot(n5(:,1),n5(:,2),'o') %公园绿地区位置图legend('生活区','工业区','山区','交通区','公园绿地区')hold oncx=shuju0(:,2); %x轴坐标cy=shuju0(:,3); %y轴坐标cz=shuju1(:,5); %所有重金属浓度数据(As Cd Cr Cu Hg Ni Pb Zn)x=0:73:28654; %插值横坐标y1=0:73:18449;y=y1'; %插值纵坐标z1=griddata(cx,cy,cz,x,y,'cubic') %插值后的重金属浓度值 z=log2(z1/(1.5*35)); %计算插值后的地累积指数 [x,y]=meshgrid(x,y); %自动生成网格[c,h]=contour(x,y,z) %插值后的地累积指数等高线图 c 是等高线矩阵,h 是等高线句柄 set(h,'ShowText','on') %把h 值写上去00.511.522.53x 1040.20.40.60.811.21.41.61.82x 104-6-4-4-4-4-4-4-4-4-4-4-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-22-2-2-2-2-2-2-2-2-2-2-2-2-2-2-200000000000000000000000000000000000022222222222222222222222222222222224444444444444466688生活区工业区山区交通区公园绿地区cx=shuju0(:,2); %x 轴坐标 cy=shuju0(:,3); %y 轴坐标cz=shuju1(:,5); %所有重金属浓度数据(As Cd Cr Cu Hg Ni Pb Zn ) x=0:73:28654; %插值横坐标 y1=0:73:18449; y=y1'; %插值纵坐标 %z=300.*x+350.*y-385000;z=griddata(cx,cy,cz,x,y,'cubic') %插值后的重金属浓度值 mesh(x,y,z) h=shuju0(:,4)xlabel('x'),ylabel('y'); title('浓度');123x 1040.511.52x 104-500005000100001500020000x浓度y2、区域污染负荷指数计算(1)函数 cumprod() ——累积求积函数:cumprod(n:m) 输出一个向量[n n*(n+1) n(n+1)(n+2) …n(n+1)(n+2)…(m-1)m]cumprod(A) 若A为矩阵:输出同维数的矩阵,按列累积求积cumprod(A,dim) A为矩阵,dim=1或2,dim=1,默认,与上面一样;dim=2,按行累积求积。
(2)clcA=load('2.txt');As=A(:,1)./3.6;Cd=A(:,2)./130;Cr=A(:,3)./31;Cu=A(:,4)./13.2;Hg=A(:,5)./35;Ni=A(:,6)./12.3;Pb=A(:,7)./31;Zn=A(:,8)./69;B=[];B=[As Cd Cr Cu Hg Ni Pb Zn];%for i=1:length(As)P=cumprod(B,2); %A为矩阵,dim=1或2,dim=1,默认,列累积;dim=2,按行累积求积。
%endPL=P(:,8); %行累积后结果PLI=PL.^(1./8); %开根号PLIz=cumprod(PLI); %列累积PLIzero=PLIz.^(1./length(As));PLIzero(length(As)) %取最后一个数ans =1.8336ans =2.3484ans =1.0602ans =1.9209ans =1.57803、地累积指数相关性检验clcA=load('1.txt');As=A(:,1)./3.6.*1.5;Cd=A(:,2)./130.*1.5;Cr=A(:,3)./31.*1.5;Cu=A(:,4)./13.2.*1.5;Hg=A(:,5)./35.*1.5;Ni=A(:,6)./12.3.*1.5;Pb=A(:,7)./31.*1.5;Zn=A(:,8)./69.*1.5;B=[];B=[As Cd Cr Cu Hg Ni Pb Zn]把Matlab左侧的Workspace中的B矩阵的数据打开,然后复制到txt文件中。
相关系数As Cd Cr Cu Hg Ni Pb ZnSpearman 的rho As 相关系数 1.000 .132*.508**.410**.382**.594**.215**.239** Sig.(双侧). .019 .000 .000 .000 .000 .000 .000 N 319 319 319 319 319 319 319 319 Cd 相关系数.132* 1.000 .384**.705**.537**.292**.811**.829** Sig.(双侧).019 . .000 .000 .000 .000 .000 .000 N 319 319 319 319 319 319 319 319 Cr 相关系数.508**.384** 1.000 .573**.398**.792**.406**.512** Sig.(双侧).000 .000 . .000 .000 .000 .000 .000 N 319 319 319 319 319 319 319 319 Cu 相关系数.410**.705**.573** 1.000 .639**.490**.764**.839** Sig.(双侧).000 .000 .000 . .000 .000 .000 .000 N 319 319 319 319 319 319 319 319Hg 相关系数.382**.537**.398**.639** 1.000 .359**.621**.647** Sig.(双.000 .000 .000 .000 . .000 .000 .000 侧)N 319 319 319 319 319 319 319 319 Ni 相关系数.594**.292**.792**.490**.359** 1.000 .287**.438**.000 .000 .000 .000 .000 . .000 .000 Sig.(双侧)N 319 319 319 319 319 319 319 319 Pb 相关系数.215**.811**.406**.764**.621**.287** 1.000 .861**.000 .000 .000 .000 .000 .000 . .000 Sig.(双侧)N 319 319 319 319 319 319 319 319 Zn 相关系数.239**.829**.512**.839**.647**.438**.861** 1.000.000 .000 .000 .000 .000 .000 .000 .Sig.(双侧)N 319 319 319 319 319 319 319 319 *. 在置信度(双测)为 0.05 时,相关性是显著的。
**. 在置信度(双测)为 0.01 时,相关性是显著的。
4、标准化数据相关性检验clcA=load('1.txt');As=A(:,1);Cd=A(:,2);Cr=A(:,3);Cu=A(:,4);Hg=A(:,5);Ni=A(:,6);Pb=A(:,7);Zn=A(:,8);AS=(As-mean(As))./var(As);CD=(Cd-mean(Cd))./var(Cd);CR=(Cr-mean(Cr))./var(Cr);CU=(Cu-mean(Cu))./var(Cu);HG=(Hg-mean(Hg))./var(Hg);NI=(Ni-mean(Ni))./var(Ni);PB=(Pb-mean(Pb))./var(Pb);ZN=(Zn-mean(Zn))./var(Zn);B=[];B=[AS CD CR CU HG NI PB ZN]相关系数As Cd Cr Cu Ni Pb ZnSpearman 的rho As 相关系数 1.000 .132*.511**.417**.594**.217**.247** Sig.(双侧). .020 .000 .000 .000 .000 .000 N 319 311 316 314 319 318 305 Cd 相关系数.132* 1.000 .378**.704**.289**.815**.834** Sig.(双侧).020 . .000 .000 .000 .000 .000 N 311 311 309 307 311 311 297 Cr 相关系数.511**.378** 1.000 .571**.797**.402**.508** Sig.(双侧).000 .000 . .000 .000 .000 .000 N 316 309 316 312 316 315 302 Cu 相关系数.417**.704**.571** 1.000 .498**.765**.837** Sig.(双侧).000 .000 .000 . .000 .000 .000 N 314 307 312 314 314 313 300 Ni 相关系数.594**.289**.797**.498** 1.000 .288**.440** Sig.(双侧).000 .000 .000 .000 . .000 .000 N 319 311 316 314 319 318 305 Pb 相关系数.217**.815**.402**.765**.288** 1.000 .860** Sig.(双侧).000 .000 .000 .000 .000 . .000 N 318 311 315 313 318 318 304Zn 相关系数.247**.834**.508**.837**.440**.860** 1.000 Sig.(双侧).000 .000 .000 .000 .000 .000 .N 305 297 302 300 305 304 305 *. 在置信度(双测)为 0.05 时,相关性是显著的。