关于用模糊综合法评价水质

模¡ê糊y综á?合?评¤价?法¤¡§ private void btnevaluate3_click(object sender, routedeventargs e) { if (evaluatecombobox3.selecteditem == **) { **box.show("请?先¨¨选?择?测a站?名?称?！ê?"); return; } d**playwindow.title = "模¡ê糊y综á?合?评¨¤价?法¤¡§结¨¢果?展1示º?"; evaluate evaluate = new evaluate(); //记?录?每?种?物?质¨º的ì?所¨´有®d含?量¢? l**t<double> l**tmn = new l**t<double>(); l**t<double> l**tnh3n = new l**t<double>(); l**t<double> l**tcod = new l**t<double>(); l**t<double> l**to2 = new l**t<double>(); l**t<double> l**tcr = new l**t<double>(); l**t<double> l**tpb = new l**t<double>(); l**t<double> l**tcd = new l**t<double>(); l**t<double> l**tphenol = new l**t<double>(); l**t<double> l**tn = new l**t<double>(); l**t<double> l**tp = new l**t<double>(); l**tmn = weightcalculatedata("高?锰¨¬酸¨¢盐?"); l**tnh3n = weightcalculatedata("氨ã¡à氮ìa"); l**tcod = weightcalculatedata("生¦¨²化¡¥需¨¨氧?量¢?"); l**to2 = weightcalculatedata("溶¨¨¹解a氧?"); l**tcr = weightcalculatedata("六¢¨´价?铬?"); l**tpb = weightcalculatedata("铅|"); l**tcd = weightcalculatedata("镉¨®"); l**tphenol = weightcalculatedata("挥¨®发¤¡é酚¤¨®"); l**tn = weightcalculatedata("总ᨹ氮ìa"); l**tp = weightcalculatedata("总ᨹ磷¢¡á"); //记?录?每?种?物?质¨º含?量¢?的ì?所¨´有®d和¨ª double summn, sumnh3n, sumcod, sumo2, sumcr, sumpb, sumcd, sumphenol, sumn, sump; summn = evaluate.sumelement(l**tmn); sumnh3n = evaluate.sumelement(l**tnh3n); sumcod = evaluate.sumelement(l**tcod); sumo2 = evaluate.sumelement(l**to2); sumcr = evaluate.sumelement(l**tcr); sumpb = evaluate.sumelement(l**tpb); sumcd = evaluate.sumelement(l**tcd); sumphenol = evaluate.sumelement(l**tphenol); sumn = evaluate.sumelement(l**tn); sump = evaluate.sumelement(l**tp); //记?录?每?种?物?质¨º每?条¬?记?录?的ì?权¨¡§重? l**t<double> wimn = new l**t<double>(); l**t<double> winh3n = new l**t<double>(); l**t<double> wicod = new l**t<double>(); l**t<double> wio2 = new l**t<double>(); l**t<double> wicr = new l**t<double>(); l**t<double> wipb = new l**t<double>(); l**t<double> wicd = new l**t<double>(); l**t<double> wiphenol = new l**t<double>(); l**t<double> win = new l**t<double>(); l**t<double> wip = new l**t<double>(); wimn = evaluate.weightcalculate(l**tmn, summn); winh3n = evaluate.weightcalculate(l**tnh3n, sumnh3n); wicod = evaluate.weightcalculate(l**tcod, sumcod); wio2 = evaluate.weightcalculate(l**to2, sumo2); wicr = evaluate.weightcalculate(l**tcr, sumcr); wipb = evaluate.weightcalculate(l**tpb, sumpb); wicd = evaluate.weightcalculate(l**tcd, sumcd); wiphenol = evaluate.weightcalculate(l**tphenol, sumphenol); win = evaluate.weightcalculate(l**tn, sumn); wip = evaluate.weightcalculate(l**tp, sump); //记?录?每?种?物?质¨º的ì?熵¬?值¦ì double simn, sinh3n, sicod, sio2, sicr, sipb, sicd, siphenol, sin, sip; simn = evaluate.entropyvalue(wimn); sinh3n = evaluate.entropyvalue(winh3n); sicod = evaluate.entropyvalue(wicod); sio2 = evaluate.entropyvalue(wio2); sicr = evaluate.entropyvalue(wicr); sipb = evaluate.entropyvalue(wipb); sicd = evaluate.entropyvalue(wicd); siphenol = evaluate.entropyvalue(wiphenol); sin = evaluate.entropyvalue(win); sip = evaluate.entropyvalue(wip); //计?算?每?种?物?质¨º的ì?权¨¡§重? double[] sivalue = new double[] { simn, sinh3n, sicod, sio2, sicr, sipb, sicd, siphenol, sin, sip }; double sumsi = 0; double[] weight = new double[10]; for (int i = 0; i < 10; i++) { sumsi = sumsi + sivalue[i]; } for (int j = 0; j < 10; j++) { weight[j] = (1.0 - sivalue[j]) / (10 - sumsi); } //隶¢£¤属º?度¨¨ l**t<double[]> mimn = new l**t<double[]>(); l**t<double[]> minh3n = new l**t<double[]>(); l**t<double[]> micod = new l**t<double[]>(); l**t<double[]> mio2 = new l**t<double[]>(); l**t<double[]> micr = new l**t<double[]>(); l**t<double[]> mipb = new l**t<double[]>(); l**t<double[]> micd = new l**t<double[]>(); l**t<double[]> miphenol = new l**t<double[]>(); l**t<double[]> min = new l**t<double[]>(); l**t<double[]> mip = new l**t<double[]>();mimn = evaluate.membership(fmn, l**tmn1); minh3n = evaluate.membership(fnh3n, l**tnh3n1); micod = evaluate.membership(fcod, l**tcod1); mio2 = evaluate.membership(fo2, l**to21); micr = evaluate.membership(fcr, l**tcr1); mipb = evaluate.membership(fpb, l**tpb1); micd = evaluate.membership(fcd, l**tcd1); miphenol = evaluate.membership(fphenol, l**tphenol1); min = evaluate.membership(fn, l**tn1); mip = evaluate.membership(fp, l**tp1); //定¡§义°?一°?个?泛¤o型¨ª用®?于®¨²存ä?储ä¡é测a站?历¤¨²年¨º来¤¡ä水?质¨º等쨨级?变à?化¡¥情¨¦况? l**t<int> rank = new l**t<int>(); //矩?阵¨®相¨¤乘?计?算?得ì?到ì?水?质¨º等쨨级? for (int m = 0; m < mimn.count; m++) { //声¦¨´明¡â一°?个?10行d5列¢d的ì?二t维?数ºy组ᨦ，ê?用®?于®¨²存ä?放¤?隶¢£¤属º?度¨¨矩?阵¨® double[,] matrixmembership = new double[10, 5]; //for (int i = 0; i < 10; i++) //{ for (int j = 0; j < 5; j++) { matrixmembership[0, j] = mimn[m][j]; matrixmembership[1, j] = minh3n[m][j]; matrixmembership[2, j] = micod[m][j]; matrixmembership[3, j] = mio2[m][j]; matrixmembership[4, j] = micr[m][j]; matrixmembership[5, j] = mipb[m][j]; matrixmembership[6, j] = micd[m][j]; matrixmembership[7, j] = miphenol[m][j]; matrixmembership[8, j] = min[m][j]; matrixmembership[9, j] = mip[m][j]; } //} //声¦¨´明¡â一°?个?1行d5列¢d的ì?一°?维?数ºy组ᨦ用®?于®¨²存ä?放¤?1行d10列¢d的ì?权¨¡§重?数ºy组ᨦ与®?10行d5列¢d的ì?隶¢£¤属º?度¨¨矩?阵¨®的ì?乘?积y double[] result = new double[5]; for (int i = 0; i < 5; i++) { for (int j = 0; j < 10; j++) { result[i] += weight[j] * matrixmembership[j, i]; } } //得ì?到ì?的ì?一°?维?数ºy组ᨦresult[]中d的ì?最á?大䨮值¦ì的ì?下?标਺加¨®1就¨ª该?指?标਺的ì?水?质¨º等쨨级? double temp = result[0]; for (int i = 1; i < 5; i++) { if (result[i] > temp) { temp = result[i]; } } for (int j = 0; j < 5; j++) { if (temp == result[j]) { rank.add(j + 1); break; } } } //用®?柱¨´状á¡ä图ª?来¤¡ä表ਪ示º?某3一°?测a站?历¤¨²年¨º来¤¡ä水?质¨º变à?化¡¥情¨¦况? d**playwindow.v**ibility = v**ibility.v**ible; evaluatechart1.title = evaluatecombobox3.selecteditem.tostring() + "测a站?历¤¨²年¨º来¤¡ä水?质¨º等쨨级?变à?化¡¥"; featurelayer featurelayer = mymap.layers["djkfeature"] as featurelayer; l**t<keyvaluepair<string, int>> lst2 = new l**t<keyvaluepair<string, int>>(); l**t<string> dt = new l**t<string>(); foreach (graphic graphic in featurelayer.graphics) { if (graphic.attributes["name"].tostring() == evaluatecombobox3.selecteditem.tostring()) { datetime shijian = (datetime)graphic.attributes["时º¡à间?"]; dt.add(shijian.year.tostring()); } } for (int i = 0; i < dt.count; i++) { lst2.add(new keyvaluepair<string, int>(dt[i], rank[i])); } columnseries columnseries = new columnseries(); binding keybinding = new binding(); keybinding.path = new propertypath("key"); columnseries.independentvaluebinding = keybinding; columnseries.title = "水?质¨º等쨨级?"; binding valuebinding = new binding(); valuebinding.path = new propertypath("value"); columnseries.dependentvaluebinding = valuebinding; columnseries.horizontalalignment = horizontalalignment.stretch; columnseries.verticalalignment = verticalalignment.stretch; columnseries.itemssource = lst2; th**.evaluatechart1.series.clear(); th**.evaluatechart1.series.add(columnseries); } private void evaluatequit_click(object sender, routedeventargs e) { evaluateborder.v**ibility = v**ibility.collapsed; evaluatewindow.v**ibility = v**ibility.collapsed; d**playwindow.v**ibility = v**ibility.collapsed; } public class binddata { public datetime dt { set; get; } public int grade { set; get; } } //得ì?到ì?每?种?物?质¨º的ì?所¨´有®d记?录? private l**t<double> weightcalculatedata(string element) { featurelayer featurelayer = mymap.layers["djkfeature"] as featurelayer; featurelayer.autos**e = true; l**t<double> elementl**t = new l**t<double>(); foreach (graphic graphic in featurelayer.graphics) { elementl**t.add(convert.todouble(graphic.attributes[element])); } return elementl**t; } 20210311