《判别分析实际应用优秀课件.ppt》由会员分享，可在线阅读，更多相关《判别分析实际应用优秀课件.ppt（50页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、判别分析实际应用第1页，本讲稿共50页 共设计六个变量：、，取f1=f2=0，引入四个变量，得生油岩演化阶段判别函数式如下：判别函数中变量的引入顺序及各演化阶段的正判率分别见表4.4和表4.5。第2页，本讲稿共50页表4.3 我国生油岩的演化参数表演化阶段 序 号 地区热演化参数未成熟1 松辽盆地（青2+3）337 125.005 10002 松辽盆地（青1）328 123 10003 松辽盆地（姚2）334 125 17504 岐口凹陷 347 33.75 16805 泌阳凹陷 342 27.999 14606 湖北（二叠系）338 242 32007 潜江凹陷 343 34.33 1600
2、8 高邮凹陷 348 15 19809 惠民凹陷 351 8.04 135010 沾化凹陷 344 10.76 160011 东明凹陷 352 8.5 225012 松辽盆地（姚2）228.5 127 1480成熟13 松辽盆地（青2+3）341 127 155014 松辽盆地（青2+1）367 127 219915 松辽盆地（青）362 124.2 179916 松辽盆地（姚2）364 120.4 185117 松辽盆地（姚2）370 120.4 197018 岐口凹陷 359 35 200119 泌阳凹陷 353 31 170020 泌阳凹陷 367 31 209921辽河凹陷 364 4
3、8.5 200122 东台凹陷（阜宁组）364 34 220123 湖北（二叠系）348 257.6 345024 高邮凹陷（集宁组）356 17.5 220025 高邮凹陷（集宁组）360 18 270026 沾化凹陷 367 10.7 230027 沾化凹陷 364 12.8 220028 沾化凹陷 368 13 250029 东明凹陷 362 11.6 3500第3页，本讲稿共50页演化阶段 序 号 地区 热演化参数高成熟30 松辽盆地（青2+3）380 129 310031 松辽盆地（青1）377 127.88 235032 松辽盆地（青1）402 127.88 249933 松辽盆地
4、（姚2）383 124.64 229934 松辽盆地（姚2）377 121.80 215035 松辽盆地（姚2）389 125.008 240136 岐口凹陷 419 36 370137 湖北（二叠系）392 217.25 405038 泌阳凹陷 376 34.01 230039 泌阳凹陷 392 34.05 279940 东台凹陷（阜宁组）410 54.25 320141 东台凹陷（阜宁组）398 37.005 290042 东台凹陷（阜宁组）403 35.5 370143 湖北（二叠系）399 271.25 400144 高邮凹陷（集宁组）380 23.5 300045 高邮凹陷（集宁组）
5、382 24 340046 惠民凹陷 383 13 280047 惠民凹陷 386 13.5 310048高邮凹陷（集宁组）399 25 410049 沾化凹陷 381 16 270050 沾化凹陷 385 16.5 300051 沾化凹陷 387 16.5 310052 东明凹陷 410 15.05 355053 东明凹陷 412 18 355054 松辽盆地（青2+3）400 130 3400第4页，本讲稿共50页演化阶段 序 号 地区 热演化参数过成熟55 松辽盆地（姚2）434 123.205 355556 江汉盆地（潜江组）440 37.005 430057 东台凹陷(阜宁组)444
6、 60.005 410058 东台凹陷(阜宁组)433 57.005 510059 湖北（二叠系）469 285 570160 湖北（二叠系）572 285.5 719961 高邮凹陷（集宁组）403 26 400162 东明凹陷 435 20 370063 东明凹陷 422 18 380064 泌阳凹陷 433 36.05 414065 湖北（二叠系）454 271.29 520066 湖北（二叠系）440 271.29 4900第5页，本讲稿共50页判别函数中变量的引入顺序及各演化阶段的正判率分别见表4.4和表4.5。第6页，本讲稿共50页表4.4变量引入顺序及检验量引入顺序变量号变量 检
7、验量1 90.10782 29.07253 0.433940.2979第7页，本讲稿共50页表4.5生油岩演化阶段正判率演化阶段正判率%演化阶段正判率%未成熟83高成熟96成熟94过成熟92第8页，本讲稿共50页n 各演化阶段的正判率都超过了80%，故可把上述判别函数用于判别生油岩演化阶段。如珠江口盆地第三系生油岩为中新世至晚渐新世，地层绝对年龄约为1630百万年，埋藏深度为2200米，地层温度104。按上述判别函数式计算，绝对年龄以25百万年计，判别函数值为第9页，本讲稿共50页 最大值为，判定属于高成熟阶段。n 用Matlab程序运行结果如下：n 调用程序zupanbie.m运行结果如下：
8、n 解 liti3.m中的数据为：%T+273 t H 1/H ln(T+273)1/(T+273)337 125.005 1000 1.0000000e003 6.4134590e+000 1.6393443e003 328 123 1000 1.0000000e003 6.3985949e+000 1.6638935e003 334 125 1750 5.7142857e004 6.4085288e+000 1.6474465e003 347 33.75 1680 5.9523810e004 6.4297195e+000 1.6129032e003 342 27.999 1460 6.84
9、93151e004 6.4216223e+000 1.6260163e003 338 242 3200 3.1250000e004 6.4150970e+000 1.6366612e003 343 34.33 1600 6.2500000e004 6.4232470e+000 1.6233766e003 348 15 1980 5.0505051e004 6.4313311e+000 1.6103060e003 351 8.04 1350 7.4074074e004 6.4361504e+000 1.6025641e003第10页，本讲稿共50页 344 10.76 1600 6.250000
10、0e004 6.4248690e+000 1.6207455e003 352 8.5 2250 4.4444444e004 6.4377516e+000 1.6000000e003 228.5 127 1480 6.7567568e004 6.2176036e+000 1.9940179e003 341 127 1550 6.4516129e004 6.4199949e+000 1.6286645e003 367 127 2199 4.5475216e004 6.4614682e+000 1.5625000e003 362 124.2 1799 5.5586437e004 6.4536250e
11、+000 1.5748031e003 364 120.4 1851 5.4024851e004 6.4567697e+000 1.5698587e003 370 120.4 1970 5.0761421e004 6.4661447e+000 1.5552100e003 359 35 2001 4.9975012e004 6.4488894e+000 1.5822785e003 353 31 1700 5.8823529e004 6.4393504e+000 1.5974441e003 367 31 2099 4.7641734e004 6.4614682e+000 1.5625000e003
12、364 48.5 2001 4.9975012e004 6.4567697e+000 1.5698587e003 364 34 2201 4.5433894e004 6.4567697e+000 1.5698587e003 348 257.6 3450 2.8985507e004 6.4313311e+000 1.6103060e003 356 17.5 2200 4.5454545e004 6.4441313e+000 1.5898251e003 360 18 2700 3.7037037e 004 6.4504704e+000 1.5797788e003第1 1页，本讲稿共50页 367
13、10.7 2300 4.3478261e004 6.4614682e+000 1.5625000e003 364 12.8 2200 4.5454545e004 6.4567697e+000 1.5698587e003 368 13 2500 4.0000000e004 6.4630295e+000 1.5600624e003 362 11.6 3500 2.8571429e004 6.4536250e+000 1.5748031e003 380 129 3100 3.2258065e004 6.4815771e+000 1.5313936e003 377 127.88 2350 4.2553
14、191e004 6.4769724e+000 1.5384615e003402 127.88 2499 4.0016006e004 6.5147127e+000 1.4814815e003 383 124.64 2299 4.3497173e004 6.4861608e+000 1.5243902e003 377 121.8 2150 4.6511628e004 6.4769724e+000 1.5384615e003 389 125.008 2401 4.1649313e004 6.4952656e+000 1.5105740e003 419 36 3701 2.7019724e004 6.
15、5395860e+000 1.4450867e003 392 217.25 4050 2.4691358e004 6.4997870e+000 1.5037594e003 376 34.01 2300 4.3478261e004 6.4754327e+000 1.5408320e003 392 34.05 2799 3.5727045e004 6.4997870e+000 1.5037594e003 410 54.25 3201 3.1240237e004 6.5264949e+000 1.4641288e003 398 37.005 2900 3.4482759e004 6.5087691e
16、+000 1.4903130e003第12页，本讲稿共50页 403 35.5 3701 2.7019724e004 6.5161931e+000 1.4792899e003 399 271.25 4001 2.4993752e004 6.5102583e+000 1.4880952e003 380 23.5 3000 3.3333333e004 6.4815771e+000 1.5313936e003 382 24 3400 2.9411765e004 6.4846352e+000 1.5267176e003 383 13 2800 3.5714286e004 6.4861608e+000
17、1.5243902e003 386 13.5 3100 3.2258065e004 6.4907235e+000 1.5174507e003 399 25 4100 2.4390244e004 6.5102583e+000 1.4880952e003 381 16 2700 3.7037037e004 6.4831074e+000 1.5290520e003 385 16.5 3000 3.3333333e004 6.4892049e+000 1.5197568e003 387 16.5 3100 3.2258065e004 6.4922398e+000 1.5151515e003 410 1
18、5.05 3.550 2.8169014e004 6.5264949e+000 1.4641288e003 412 18 3550 2.8169014e004 6.5294188e+000 1.4598540e003 400 130 3400 2.9411765e004 6.5117453e+000 1.4858841e003 434 123.205 3555 2.8129395e004 6.5610307e+000 1.4144272e003 440 37.005 4300 2.3255814e004 6.5694814e+000 1.4025245e003 444 60.005 4100
19、2.4390244e004 6.5750758e+000 1.3947001e003 433 57.005 5100 1.9607843e004 6.5596152e+000 1.4164306e003第13页，本讲稿共50页 469 285 5701 1.7540782e004 6.6093492e+000 1.3477089e003 572 285.5 7199 1.3890818e004 6.7393366e+000 1.1834320e003 403 26 4001 2.4993752e004 6.5161931e+000 1.4792899e003 435 20 3700 2.702
20、7027e004 6.5624441e+000 1.4124294e003 422 18 3800 2.6315789e004 6.5439118e+000 1.4388489e003 433 36.05 4140 2.4154589e004 6.5596152e+000 1.4164306e003 454 271.29 5200 1.9230769e004 6.5889265e+000 1.3755158e003 440 271.29 4900 2.0408163e004 6.5694814e+000 1.4025245e003 第14页，本讲稿共50页调用程序zupanbie.m运行结果如
21、下：*逐步判别分析的程序实现*请输入需判别样品的组数样品种类：4请输入需判别样品的总数：66请输入判别样品的变量数：6请输入引入变量f的检验临界值：0请输入剔除变量f的检验临界值：0第15页，本讲稿共50页请输入第1组样品数：12请输入第2组样品数：17请输入第3组样品数：15请输入第4组样品数：12第16页，本讲稿共50页n 每个变量的总均值为：xe(1)=383.4773xe(2)=80.2577xe(3)=2900.2879xe(4)=0.00040557xe(5)=6.4847xe(6)=0.0015301第17页，本讲稿共50页n 各组内变量的均值为：gxe(1,1)=332.708
22、3gxe(1,2)=73.3653gxe(1,3)=1695.8333gxe(1,4)=0.00064833gxe(1,5)=6.4048gxe(1,6)=0.0016564gxe(2,1)=360.9412gxe(2,2)=67.0412 第18页，本讲稿共50页gxe(2,3)=2248.2941gxe(2,4)=0.00046541gxe(2,5)=6.4519gxe(2,6)=0.0015777gxe(3,1)=391.8gxe(3,2)=99.9349gxe(3,3)=2963.4667gxe(3,4)=0.00035231gxe(3,5)=6.4993gxe(3,6)=0.0015048第19页，本讲稿共50页gxe(4,1)=399.9167gxe(4,2)=37.3133gxe(4,3)=3379.5833gxe(4,4)=0.00030128gxe(4,5)=6.5112gxe(4,6)=0.0014873log(q(1)=1.7047log(q(2)=1.3564log(q(3)=1.4816log(q(4)=1.7047第20页，本讲稿共50页