美赛论文写作
美赛:全国大学生数学建模竞赛论文格式规范---论文写作规范

全国大学生数学建模竞赛论文格式规范●本科组参赛队从A、B题中任选一题,专科组参赛队从C、D题中任选一题。
(全国评奖时,每个组别一、二等奖的总名额按每道题参赛队数的比例分配;但全国一等奖名额的一半将平均分配给本组别的每道题,另一半按每题论文数的比例分配。
)●论文用白色A4纸打印;上下左右各留出至少2.5厘米的页边距;从左侧装订。
●论文第一页为承诺书,具体内容和格式见本规范第二页。
●论文第二页为编号专用页,用于赛区和全国评阅前后对论文进行编号,具体内容和格式见本规范第三页。
●论文题目、摘要和关键词写在论文第三页上(无需译成英文),并从此页开始编写页码;页码必须位于每页页脚中部,用阿拉伯数字从“1”开始连续编号。
注意:摘要应该是一份简明扼要的详细摘要,请认真书写(但篇幅不能超过一页)。
●从第四页开始是论文正文(不要目录)。
论文不能有页眉或任何可能显示答题人身份和所在学校等的信息。
●论文应该思路清晰,表达简洁(正文尽量控制在20页以内,附录页数不限)。
●引用别人的成果或其他公开的资料(包括网上查到的资料) 必须按照规定的参考文献的表述方式在正文引用处和参考文献中均明确列出。
正文引用处用方括号标示参考文献的编号,如[1][3]等;引用书籍还必须指出页码。
参考文献按正文中的引用次序列出,其中书籍的表述方式为:[编号] 作者,书名,出版地:出版社,出版年。
参考文献中期刊杂志论文的表述方式为:[编号] 作者,论文名,杂志名,卷期号:起止页码,出版年。
参考文献中网上资源的表述方式为:[编号] 作者,资源标题,网址,访问时间(年月日)。
●在论文纸质版附录中,应给出参赛者实际使用的软件名称、命令和编写的全部计算机源程序(若有的话)。
同时,所有源程序文件必须放入论文电子版中备查。
论文及源程序电子版压缩在一个文件中,一般不要超过20MB,且应与纸质版同时提交。
(如果发现程序不能运行,或者运行结果与论文中报告的不一致,该论文可能会被认定为弄虚作假而被取消评奖资格。
美赛论文(最终版)

For office use onlyT1________________ T2________________ T3________________ T4________________Team Control Number 46639Problem ChosenCFor office use onlyF1________________F2________________F3________________F4________________2016 MCM/ICM Summary SheetAn Optimal Investment Strategy ModelSummaryWe develop an optimal investment strategy model that appears to hold promise for providing insight into not only how to sort the schools according to investment priority, but also identify optimal investment amount of a specific school. This model considers a large number of parameters thought to be important to investment in the given College Scorecard Data Set.In order to develop the required model, two sub-models are constructed as follows: 1.For Analytic Hierarchy Process (AHP) Model, we identify the prioritizedcandidate list of schools by synthesizing the elements which have an influence on investment. First we define the specific value of any two elements’ effect on investment. And then the weight of each element’s influence on investment can be identified. Ultimately, we take the relevant parameters into the calculated weight, and then we get any school’s recommended value of investment.2.For Return On Investment M odel, it’s constructed on the basis of AHP Model.Let us suppose that all the investment is used to help the students to pay tuition fee.Then we can see optimal investment as that we help more students to the universities of higher return rate. However, because of dropout rate, there will be an optimization investment amount in each university. Therefore, we can change the problem into a nonlinear programming problem. We identify the optimal investment amount by maximizing return-on-investment.Specific attention is given to the stability and error analysis of our model. The influence of the model is discussed when several fundamental parameters vary. We attempt to use our model to prioritize the schools and identify investment amount of the candidate schools, and then an optimal investment strategy is generated. Ultimately, to demonstrate how our model works, we apply it to the given College Scorecard Data Set. For various situations, we propose an optimal solution. And we also analyze the strengths and weaknesses of our model. We believe that we can make our model more precise if more information are provided.Contents1.Introduction 21.1Restatement of the Problem (2)1.2Our Approach (2)2.Assumptions 23.Notations 34.The Optimal Investment Model 44.1Analytic Hierarchy Process Model (4)4.1.1Constructing the Hierarchy (4)4.1.2Constructing the Judgement Matrix (5)4.1.3Hierarchical Ranking (7)4.2Return On Investment Model (8)4.2.1Overview of the investment strategy (8)4.2.2Analysis of net income and investment cost (9)4.2.3Calculate Return On Investment (11)4.2.4Maximize the Total Net Income (11)5.Test the Model125.1Error Analysis (12)5.2Stability Analysis (13)6.Results136.1Results of Analytic Hierarchy Process (13)6.2Results of Return On Investment Model (14)7.Strengths and Weaknesses157.1Strengths (15)7.2Weaknesses (16)References16 Appendix A Letter to the Chief Financial Officer, Mr. Alpha Chiang.171.Introduction1.1Restatement of the ProblemIn order to help improve educational performance of undergraduates attending colleges and universities in the US, the Goodgrant Foundation intends to donate a total of $100,000,000 to an appropriate group of schools per year, for five years, starting July 2016. We are to develop a model to determine an optimal investment strategy that identifies the school, the investment amount per school, the return on that investment, and the time duration that the organization’s money should be provided to have the highest likelihood of producing a strong positive effect on student performance. Considering that they don’t want to duplicate the investments and focus of other large grant organizations, we interpret optimal investment as a strategy that maximizes the ROI on the premise that we help more students attend better colleges. So the problems to be solved are as follows:1.How to prioritize the schools by optimization level.2.How to measure ROI of a school.3.How to measure investment amount of a specific school.1.2Our ApproachWe offer a model of optimal investment which takes a great many factors in the College Scorecard Data Set into account. To begin with, we make a 1 to N optimized and prioritized candidate list of school we are recommending for investment by the AHP model. For the sake that we invest more students to better school, several factors are considered in the AHP model, such as SAT score, ACT score, etc. And then, we set investment amount of each university in the order of the list according to the standard of maximized ROI. The implement details of the model will be described in section 4.2.AssumptionsWe make the following basic assumptions in order to simplify the problem. And each of our assumptions is justified.1.Investment amount is mainly used for tuition and fees. Considering that theincome of an undergraduate is usually much higher than a high school students, we believe that it’s necessary to help more poor students have a chance to go to college.2.Bank rates will not change during the investment period. The variation ofthe bank rates have a little influence on the income we consider. So we make this assumption just to simplify the model.3.The employment rates and dropout rates will not change, and they aredifferent for different schools4.For return on investment, we only consider monetary income, regardlessof the intangible income.3.NotationsWe use a list of symbols for simplification of expression.4.The Optimal Investment ModelIn this section, we first prioritize schools by the AHP model (Section 4.1), and then calculate ROI value of the schools (Section 4.2). Ultimately, we identify investment amount of every candidate schools according to ROI (Section 4.3).4.1Analytic Hierarchy Process ModelIn order to prioritize schools, we must consider each necessary factor in the College Scorecard Data Set. For each factor, we calculate its weight value. And then, we can identify the investment necessity of each school. So, the model can be developed in 3 steps as follows:4.1.1Constructing the HierarchyWe consider 19 elements to measure priority of candidate schools, which can be seen in Fig 1. The hierarchy could be diagrammed as follows:Fig.1AHP for the investment decisionThe goal is red, the criteria are green and the alternatives are blue. All the alternatives are shown below the lowest level of each criterion. Later in the process, each alternatives will be rated with respect to the criterion directly above it.As they build their hierarchy, we should investigate the values or measurements of the different elements that make it up. If there are published fiscal policy, for example, or school policy, they should be gathered as part of the process. This information will be needed later, when the criteria and alternatives are evaluated.Note that the structure of the investment hierarchy might be different for other foundations. It would definitely be different for a foundation who doesn't care how much his score is, knows he will never dropout, and is intensely interested in math, history, and the numerous aspects of study[1].4.1.2Constructing the Judgement MatrixHierarchy reflects the relationship among elements to consider, but elements in the Criteria Layer don’t always weigh equal during aim measure. In deciders’ mind, each element accounts for a particular proportion.To incorporate their judgments about the various elements in the hierarchy, decision makers compare the elements “two by two”. The fundamental scale for pairwise comparison are shown in Fig 2.Fig 2Right now, let's see which items are compared. Our example will begin with the six criteria in the second row of the hierarchy in Fig 1, though we could begin elsewhere if we want. The criteria will be compared as to how important they are to the decisionmakers, with respect to the goal. Each pair of items in this row will be compared.Fig 3 Investment Judgement MatrixIn the next row, there is a group of 19 alternatives under the criterion. In the subgroup, each pair of alternatives will be compared regarding their importance with respect to the criterion. (As always, their importance is judged by the decision makers.) In the subgroup, there is only one pair of alternatives. They are compared as to how important they are with respect to the criterion.Things change a bit when we get to the alternatives row. Here, the factor in each group of alternatives are compared pair-by-pair with respect to the covering criterion of the group, which is the node directly above them in the hierarchy. What we are doing here is evaluating the models under consideration with respect to score, then with respect to Income, then expenditure, dropout rate, debt and graduation rate.The foundation can evaluate alternatives against their covering criteria in any order they choose. In this case, they choose the order of decreasing priority of the covering criteria.Fig 4 Score Judgement MatrixFig 5 Expenditure Judgement MatrixFig 6 Income Judgement MatrixFig 7 Dropout Judgement MatrixFig 8 Debt Judgement MatrixFig 9 Graduation Matrix4.1.3 Hierarchical RankingWhen the pairwise comparisons are as numerous as those in our example, specialized AHP software can help in making them quickly and efficiently. We will assume that the foundation has access to such software, and that it allows the opinions of various foundations to be combined into an overall opinion for the group.The AHP software uses mathematical calculations to convert these judgments to priorities for each of the six criteria. The details of the calculations are beyond the scope of this article, but are readily available elsewhere[2][3][4][5]. The software also calculates a consistency ratio that expresses the internal consistency of the judgments that have been entered. In this case the judgments showed acceptable consistency, and the software used the foundation’s inputs to assign these new priorities to the criteria:Fig 10.AHP hierarchy for the foundation investing decision.In the end, the AHP software arranges and totals the global priorities for each of the alternatives. Their grand total is 1.000, which is identical to the priority of the goal. Each alternative has a global priority corresponding to its "fit" to all the foundation's judgments about all those aspects of factor. Here is a summary of the global priorities of the alternatives:Fig 114.2 ROI Model4.2.1 Overview of the investment strategyConsider a foundation making investment on a set of N geographically dispersed colleges and university in the United States, D = {1, 2, 3……N }. Then we can select top N schools from the candidate list which has been sorted through analytic hierarchy process. The total investment amount is M per year which is donated by the Goodgrant Foundation. The investment amount is j m for each school j D ∈, satisfying the following balance constraint:j j D mM ∈=∑ (1)W e can’t invest too much or too little money to one school because we want to help more students go to college, and the student should have more choices. Then the investment amount for each school must have a lower limit lu and upper limit bu as follows:j lu m bu ≤≤ (2)The tuition and fees is j p , and the time duration is {1,2,3,4}j t ∈. To simplify ourmodel, we assume that our investment amount is only used for freshmen every year. Because a freshmen oriented investment can get more benefits compared with others. For each school j D ∈, the number of the undergraduate students who will be invested is j n , which can be calculated by the following formula :,jj j j m n j D p t =∈⨯ (3)Figure12The foundation can use the ROI model to identify j m and j t so that it canmaximize the total net income. Figure1 has shown the overview of our investment model. We will then illustrate the principle and solution of this model by a kind of nonlinear programming method.4.2.2 Analysis of net income and investment costIn our return on investment model, we first focus on analysis of net income and investment cost. Obviously, the future earnings of undergraduate students are not only due to the investment itself. There are many meaning factors such as the effort, the money from their parents, the training from their companies. In order to simplify the model, we assume that the investment cost is the most important element and we don’t consider other possible influence factors. Then we can conclude that the total cost of the investment is j m for each school j D ∈.Figure 13For a single student, the meaning of the investment benefits is the expected earnings in the future. Assuming that the student is not going to college or university after graduating from high school and is directly going to work. Then his wage base is 0b as a high school graduate. If he works as a college graduate, then his wage base is 0a . Then we can give the future proceeds of life which is represented symbolically by T and we use r to represent the bank rates which will change over time. We assume that the bank rates will not change during the investment period. Here, we use bank rates in 2016 to represent the r . The future proceeds of life of a single undergraduate student will be different due to individual differences such as age, physical condition environment, etc. If we consider these differences, the calculation process will be complicated. For simplicity’s sake, we uniform the future proceeds of life T for 20 years. Then we will give two economics formulas to calculate the total expected income in the next T years for graduates and high school graduates:40(1)Tk k a u r +==+∑(4) 40(1)T kk b h r +==+∑(5) The total expected income of a graduate is u , and the total expected income of a highschool graduate is h .Then, we continue to analyze the net income. The net income can be calculated by the following formula:os NetIncome TotalIncome C t =- (6) For each school j D ∈, the net income is j P , the total income is j Q , and the cost is j m . Then we will get the following equation through formula (6):j j j P Q m =- (7)Therefore, the key of the problem is how to calculate j Q . In order to calculate j Q, weneed to estimate the number of future employment j ne . The total number of the invested is j n , which has been calculated above. Considering the dropout rates j α and the employment rates j β for each school j , we can calculate the number of future employment j ne through the following formula:(4)(1)jt j j j j n e n βα-=⨯⨯- (8)That way, we can calculate j Q by the following formula:()j j Q ne u h =⨯- (9)Finally, we take Eq. (2) (3) (4) (7) (8) into Eq. (6), and we will obtain Eq. (9) as follows:4(4)00400(1)()(1)(1)j TT t j j j j j k kk k j jm a b P m p t r r βα+-+===⨯⨯-⨯--⨯++∑∑ (10) We next reformulate the above equation of j P for concise presentation:(4)(1)j t j jj j j jc m P m t λα-⨯⨯=⨯-- (11)where jj j p βλ= and 400400(1)(1)TT k kk k a b c r r ++===-++∑∑ .4.2.3 Calculate Return On InvestmentROI is short of return on investment which can be determined by net income andinvestment cost [7]. It conveys the meaning of the financial assessment. For each schoolj D ∈ , the net income is j P , and the investment cost equals to j m . Then the j ROIcan be calculated by the following formula:100%j j jP ROI m =⨯ (12)We substitute Eq. (10) into Eq. (11), and we will get a new formula as follows:(4)((1)1)100%j t j j j jc ROI t λα-⨯=⨯--⨯ (13)4.2.4 Maximize the Total Net IncomeGiven the net income of each school, we formulate the portfolio problem that maximize the total net income, S=Max(4)((1))j t j jj j j j Dj Djc m P m t λα-∈∈⨯⨯=⨯--∑∑ (14)S. T.jj DmM ∈=∑,{1,2,3,4}t = ,j lu m bu ≤≤ ,Considering the constraint jj DmM ∈=∑, we can further simplify the model,S is equivalent to S’=Max(4)((1))j t j jj j j Dj Djc m P t λα-∈∈⨯⨯=⨯-∑∑ (15)S. T.jj DmM ∈=∑,{1,2,3,4t = ,j l u m b u ≤≤. By solving the nonlinear programming problem S’, we can get the sameanswer as problem S.5. Testing the Model 5.1 Error AnalysisSince the advent of analytic hierarchy process, people pay more attention to it due to the specific applicability, convenience, practicability and systematization of the method. Analytic hierarchy process has not reached the ideal situation whether in theory or application level because the results depend largely on the preference and subjective judgment. In this part, we will analyze the human error problem in analytic hierarchy process.Human error is mainly caused by human factors. The human error mainly reflects on the structure of the judgment matrix. The causes of the error are the following points:1. The number of times that human judge the factors’ importance is excessive.2. The calibration method is not perfect.Then we will give some methods to reduce errors:1. Reduce times of human judgment. One person repeatedly gave the samejudgment between two factors. Or many persons gave the same judgment between two factors one time. Finally, we take the average as result.2. Break the original calibration method. If we have defined the ranking vector111121(,...)n a a a a =between the factor 1A with others. Then we can get all theother ranking vector. For example : 12122111(,1...)na a a a a =.5.2 Stability AnalysisIt is necessary to analyze the stability of ranking result [6], because the strong subjectivefactors. If the ranking result changed a little while the judgment changed a lot, we can conclude that the method is effective and the result is acceptable, and vice versa. We assume that the weight of other factors will change if the weight of one factor changed from i ξ to i η:[8](1)(,1,2...,)(1)i j j i i j n i j ηξηξ-⨯==≠- (16)And it is simple to verify the equation:11nii η==∑ (17)And the new ranking vector ω will be:A ωη=⨯ (18)By this method, the Relative importance between other factors remain the same while one of the factor has changed.6. Results6.1 Results of Analytic Hierarchy ProcessWe can ranking colleges through the analytic hierarchy process, and we can get the top N = 20 schools as follows6.2 Results of Return On Investment ModelBased on the results above, we next use ROI model to distribute investment amountj m and time duration j t for each school j D ∈ by solving the following problem:Max (4)((1))j t j jj j j Dj Djc m P t λα-∈∈⨯⨯=⨯-∑∑S. T.jj DmM ∈=∑,{1,2,3,4t = , j l u m b u≤≤ . In order to solve the problem above, we collected the data from different sources. Inthe end, we solve the model with Lingo software. The program code is as follows:model: sets:roi/1..20/:a,b,p,m,t;endsets data:a = 0.9642 0.9250 0.9484 0.9422 0.9402 0.9498 0.90490.9263 0.9769 0.9553 0.9351 0.9123 0.9410 0.98610.9790 0.9640 0.8644 0.9598 0.9659 0.9720;b = 0.8024 0.7339 0.8737 0.8308 0.8681 0.7998 0.74920.6050 0.8342 0.8217 0.8940 0.8873 0.8495 0.87520.8333 0.8604 0.8176 0.8916 0.7527 0.8659;p = 3.3484 3.7971 3.3070 3.3386 3.3371 3.4956 3.22204.0306 2.8544 3.1503 3.2986 3.3087 3.3419 2.78452.9597 2.92713.3742 2.7801 2.5667 2.8058;c = 49.5528;enddatamax=@sum(roi(I):m(I)/t(I)/p(I)*((1-b(I))^4)*c*(1-a(I)+0.05)^(4-t(I)));@for(roi:@gin(t));@for(roi(I):@bnd(1,t(I),4));@for(roi(I):@bnd(0,m(I),100));@sum(roi(I):m(I))=1000;ENDFinally, we can get the investment amount and time duration distribution as follows:7.Strengths and Weaknesses7.1Strengths1.Fixing the bank rates during the investment period may run out, but it will haveonly marginal influences.2.For return on investment, we only consider monetary income, regardless of the3.intangible income. But the quantization of these intangible income is very importantand difficult. It needs to do too much complicated technical analysis and to quantify 4.too many variables. Considering that the investment persists for a short time, thiskind of random error is acceptable.5.Due to our investment which is freshmen oriented, other students may feel unfair.It is likely to produce adverse reaction to our investment strategy.6.The cost estimation is not impeccable. We only consider the investment amount andignore other non-monetary investment.5. AHP needs higher requirements for personnel quality.7.2Weaknesses1.Our investment strategy is distinct and clear, and it is convenient to implement.2.Our model not only identifies the investment amount for each school, but alsoidentifies the time duration that the organization’s money should be provide d.3.Data processing is convenient, because the most data we use is constant, average ormedian.4.Data sources are reliable. Our investment strategy is based on some meaningful anddefendable subset of two data sets.5.AHP is more simple, effective and universal.References[1] Saaty, Thomas L. (2008). Decision Making for Leaders: The Analytic Hierarchy Process for Decisions in a Complex World. Pittsburgh, Pennsylvania: RWS Publications. ISBN 0-9620317-8-X.[2] Bhushan, Navneet, Kanwal Rai (January 2004). Strategic Decision Making: Applying the Analytic Hierarchy Process. London: Springer-Verlag. ISBN 1-8523375-6-7.[3] Saaty, Thomas L. (2001). Fundamentals of Decision Making and Priority Theory. Pittsburgh, Pennsylvania: RWS Publications. ISBN 0-9620317-6-3.[4] Trick, Michael A. (1996-11-23). "Analytic Hierarchy Process". Class Notes. Carnegie Mellon University Tepper School of Business. Retrieved 2008-03-02.[5] Meixner, Oliver; Reiner Haas (2002). Computergestützte Entscheidungs-findung: Expert Choice und AHP – innovative Werkzeuge zur Lösung komplexer Probleme (in German). Frankfurt/Wien: Redline Wirtschaft bei Ueberreuter. ISBN 3-8323-0909-8.[6] Hazelkorn, E. The Impact of League Tables and Ranking System on Higher Education Decision Making [J]. Higher Education Management and Policy, 2007, 19(2), 87-110.[7] Leslie: Trainer Assessment: A Guide to Measuring the Performance of Trainers and Facilitors, Second Edition, Gower Publishing Limited, 2002.[8] Aguaron J, Moreno-Jimenea J M. Local stability intervals in the analytic hierarchy process. European Journal of Operational Research. 2000Letter to the Chief Financial Officer, Mr. Alpha Chiang. February 1th, 2016.I am writing this letter to introduce our optimal investment strategy. Before I describe our model, I want to discuss our proposed concept of a return-on-investment (ROI). And then I will describe the optimal investment model by construct two sub-model, namely AHP model and ROI model. Finally, the major results of the model simulation will be showed up to you.Considering that the Goodgrant Foundation aims to help improve educational performance of undergraduates attending colleges and universities in the US, we interpret return-on-investment as the increased income of undergraduates. Because the income of an undergraduate is generally much higher than a high school graduate, we suggest all the investment be used to pay for the tuition and fees. In that case, if we take both the income of undergraduates’ income and dropout rate into account, we can get the return-in-investment value.Our model begins with the production of an optimized and prioritized candidate list of schools you are recommending for investment. This sorted list of school is constructed through the use of specification that you would be fully qualified to provided, such as the score of school, the income of graduate student, the dropout rate, etc. With this information, a precise investment list of schools will be produced for donation select.Furthermore, we developed the second sub-model, ROI model, which identifies the investment amount of each school per year. If we invest more money in a school, more students will have a chance to go to college. However, there is an optimal investment amount of specific school because of the existence of dropout. So, we can identify every candidate school’s in vestment amount by solve a nonlinear programming problem. Ultimately, the result of the model simulation show that Washington University, New York University and Boston College are three schools that worth investing most. And detailed simulation can be seen in our MCM Contest article.We hope that this model is sufficient in meeting your needs in any further donation and future philanthropic educational investments within the United States.。
美赛论文模版

摘要:第一段:写论文解决什么问题1.问题的重述a. 介绍重点词开头:例1:“Hand move” irrigation, a cheap but labor-intensive system used on small farms, consists of a movable pipe with sprinkler on top that can be attached to a stationary main.例2:……is a real-life common phenomenon with many complexi t ies.例3:An (effective plan) is crucial to………b. 直接指出问题:例1:We find the optimal number of tollbooths in a highway toll-plaza for a given number of highway lanes: the number of tollbooths that minimizes average delay experienced by cars.例2:A brand-new university needs to balance the cost of information technology security measures wi t h the potential cost of attacks on its systems.例3:We determine the number of sprinklers to use by analyzing the energy and motion of water in the pipe and examining the engineering parameters of sprinklers available in the market.例4: After mathematically analyzing the …… problem, our modeling group would like to present our conclusions, strategies, (and recommendations )to the …….例5:Our goal is... that (mini mizes the time )……….2.解决这个问题的伟大意义反面说明。
美赛写作材料模板及参赛经验分享

Tittle of paperSummary/Abstract Key words:I.Introduction(引言)Organ transplantation is a preferable treatment for the most serious forms of end-stage diseases. In recent years, advances in medical science and technology have made solid organ transplantation an increasingly successful and common medical procedure, a literal ''second chance at life". Not only does it offer the best hope for complete rehabilitation, but it has also proved to be the most cost-effective of all treatment options, including dialysis. Consequently, more and more people are benefiting from organ transplants and their survival rates are steadily improving. The surgical techniques involved have been mastered for half a century and are now considered as routine. The two main sources of kidneys for transplantation are deceased-donor kidneys and live-donations from family and friends. However, unfortunately, there is a considerable shortage of donor organs, compared to demands. As a matter of fact, efficient matching and allocation of organs donated has become an exigent problem.The United Network for Organ Sharing (UNOS), as the operator of the Organ Procurement and Transplantation Network (OPTN), is responsible for transplant organ distribution in the United States. UNOS oversees the allocation of many different types of transplants, including liver, kidney, pancreas, heart, lung, and cornea.Focusing on kidney transplantation, based on UNOS Kidney Allocation Model, we develop a mathematical model for US transplant networks. First, incomingorgans are matched with waiting candidates by medical institutions considering the factors as ABO blood compatibility, the degree of recipient major HLA mismatch in order to obtain a matching degree applied on the allocation part. After that, from the patients’perspective, on the basis of linear regression, priority weight is defined by pondering age, disease severity, time on waiting, PRA level, and region. Applying this mechanism of ranking, we realize MWBM (Maximum Weight Bipartite-graph Matching) and SMGS (Stable Matching based on Gale-Shapley algorithm). MWBM focuses on the optimal assignment of donors following the algorithm of bipartite-graph maximum weight matching; SMGS emphasizes the process of exchanges in order to obtain the stable exchanges between donors and candidates on the waiting list.II.The Description of Problem(问题重述)III.Basic Assumptions●The level of mismatch is only relative to the number of antigens.●The data and information are accurately registered according to the medicalmeasures●The data and information are refreshed in time according to the status of thepatients●No differences in the quality of the donor kidneys●The quality of the donor kidney is constantIV.Definitions and Notations●Kidney transplantation: A kidney transplant is a surgical procedure to implant ahealthy kidney into a patient with kidney failure.●Prioritization●MD: Matching Degree●PW: Prioritization weight●MWSM: Maximum Weight Bipartite Matching●SMGS: Stable Matching based on Gale-Shapley algorithm或V.ModelsThrough the investigation of US transplantation network, we draw a general picture of the mechanism. With reference to some resources available on the website of UNOS, a flow chart (Figure 1) is developed showing the procedure of the network.Currently, the initial waiting list is composed of patients who are waiting for a kidney or combined kidney-pancreas transplant. For the first time, the patients arerequested to show the correct and scientific information to the US kidney transplant network which is needed for donor-recipient matching, the ranking of patients on the waiting list, and determining the outcome of those transplanted. The patients’waiting lists are composed of initial patients, historical patients and unsuccessful recipient after transplantation. Historical patients refer to registered patients whose status have changed and have an influence on the procedure. A patient is taken off the waiting list when a graft is offered and accepted by that patient or the patient is dead while waiting for a transplant. Unsuccessful recipients refer to the patients who have a bad result of transplantation calling for transplantation again, as it is so-called relistFigure 1. A schematic depicting the steps occurring in the transplantation networks......Table 1.Survival rate involving HLA mismatchVI.Conclusions.Our model for the optimal allocation of the donor organs is established by three modules, procurement of MD and PW, optimal assignment by MWBM model and Stable Matching of Gale-Shapley algorithm. The model has offered a convincing procedure of the allocation with the ……VII.Strengths and weaknesses(模型优缺点)Strengths●……Weaknesses●VIII.References注意文献的积累,不要等到文章写完再去重新寻找文献。
美赛数模论文

MCM 2015 Summary Sheet for Team 35565For office use onlyT1________________ T2________________ T3________________ T4________________Team Control Number35565Problem ChosenBFor office use onlyF1________________F2________________F3________________F4________________ SummaryThe lost MH370 urges us to build a universal search plan to assist searchers to locate the lost plane effi-ciently and optimize the arrangement of search plans.For the location of the search area, we divided it into two stages, respectively, to locate the splash point and the wreckage‟s sunk point. In the first stage, we consider the types of crashed aircraft, its motion and different position out of contact. We also consider the Earth‟s rotation, and other factors. Taking all these into account, we establish a model to locate the splash point. Then we apply this model to MH370. we can get the splash point in the open water is 6.813°N 103.49°E and the falling time is 52.4s. In the second stage, considering resistances of the wreckage in different shapes and its distribution affected by ocean currents, we establish a wreckage sunk point model to calculate the horizontal displacement and the angle deviation affected by the ocean currents. The result is 1517m and 0.11°respectively. Next, we extract a satellite map of submarine topography and use MATLAB to depict seabed topography map, determining the settlement of the wreckage by using dichotomy algorithm under different terrains. Finally, we build a Bayesian model and calculate the weight of corresponding area, sending aircrafts to obtain new evidence and refresh suspected wreckage area.For the assignment of the search planes, we divide it into two stages, respectively, to determine the num-ber of the aircraft and the assignment scheme of the search aircraft. In the first stage, we consider the search ability of each plane and other factors. And then we establish global optimization model. Next we use Dinkelbach algorithm to select the best n search aircrafts from all search aircrafts. In the second stage, we divide the assignment into two cases whether there are search aircrafts in the target area. If there is no search aircraft, we take the search area as an arbitrary polygon and establish the subdivision model. Considering the searching ability of each plane, we divide n small polygons into 2n sub-polygons by using NonconvexDivide algorithm, which assigns specific anchor points to these 2n sub-polygons re-spectively. If there exist search aircrafts, we divide the search area into several polygons with the search aircrafts being at the boundary of the small polygons. To improve search efficiency, we introduce” ma x-imize the minimum angle strategy” to maximize right-angle subdivision so that we can reduce the turning times of search aircraft. When we changed the speed of the crashed plane about 36m/s, the latitude of the splash point changes about 1°.When a wreck landing at 5.888m out from the initial zone, it will divorce from suspected searching area, which means our models are fairly robust to the changes in parameters. Our model is able to efficiently deal with existing data and modify some parameters basing the practical situation. The model has better versatility and stability. The weakness of our model is neglect of human factors, the search time and other uncontrollable factors that could lead to deviation compared to practical data. Therefore, we make some in-depth discussions about the model, modifying assumptions establish-Searching For a Lost PlaneControl#35565February 10, 2014Team # 35565 Page 3 of 47 Contents1 Introduction (5)1.1 Restatement of the Problem (5)1.2 Literature Review (6)2 Assumptions and Justifications (7)3 Notations (7)4 Model Overview (10)5 Modeling For Locating the Lost Plane (10)5.1 Modeling For Locating the Splash Poin t (11)5.1.1 Types of Planes (11)5.1.2 Preparation of the Model—Earth Rotation (12)5.1.3 Modeling (13)5.1.4 Solution of The Model (14)5.2 Modeling For Locating Wreckage (15)5.2.1 Assumptions of the Model (16)5.2.2 Preparation of the Model (16)5.2.3 Modeling (21)5.2.4 Solution of the Model (25)5.3 Verification of the Model (26)5.3.1 Verification of the Splash Point (26)5.3.2 Verification of the binary search algorithm (27)6 Modeling For Optimization of Search Plan (29)6.1 The Global Optimization Model (29)6.1.1 Preparation of the Model (29)6.1.2 Modeling (31)6.1.3 Solution of the Model (31)6.2 The Area Partition Algorithm (33)6.2.1 Preparation of the Model (33)6.2.2 Modeling (34)6.2.3 Solution of the Model (35)6.2.4 Improvement of the Model (36)7 Sensitivity Analysis (38)8 Further Discussions (39)9 Strengths and Weaknesses (41)9.1 Strengths (41)9.2 Weaknesses (42)10 Non-technical Paper (42)1 IntroductionAn airplane (informally plane) is a powered, fixed-wing aircraft that is propelled for-ward by thrust from a jet engine or propeller. Its main feature is fast and safe. Typi-cally, air travel is approximately 10 times safer than travel by car, rail or bus. Howev-er, when using the deaths per journey statistic, air travel is significantly more danger-ous than car, rail, or bus travel. In an aircraft crash, almost no one could survive [1]. Furthermore, the wreckage of the lost plane is difficult to find due to the crash site may be in the open ocean or other rough terrain.Thus, it will be exhilarating if we can design a model that can find the lost plane quickly. In this paper, we establish several models to find the lost plane in seawater and develop an op-timal scheme to assign search planes to model to locate the wreckage of the lost plane.1.1 Restatement of the ProblemWe are required to build a mathematical model to find the lost plane crashed in open water. We decompose the problem into three sub-problems:●Work out the position and distributions of the plane‟s wreckage●Arrange a mathematical scheme to schedule searching planesIn the first step, we seek to build a model with the inputs of altitude and other factors to locate the splash point on the sea-level. Most importantly, the model should reflect the process of the given plane. Then we can change the inputs to do some simulations. Also we can change the mechanism to apply other plane crash to our model. Finally, we can obtain the outputs of our model.In the second step, we seek to extend our model to simulate distribution of the plane wreckage and position the final point of the lost plane in the sea. We will consider more realistic factors such as ocean currents, characteristics of plane.We will design some rules to dispatch search planes to confirm the wreckage and de-cide which rule is the best.Then we attempt to adjust our model and apply it to lost planes like MH370. We also consider some further discussion of our model.1.2 Literature ReviewA model for searching the lost plane is inevitable to study the crashed point of the plane and develop a best scheme to assign search planes.According to Newton's second law, the simple types of projectile motion model can work out the splash point on the seafloor. We will analyze the motion state ofthe plane when it arrives at the seafloor considering the effect of the earth's rotation,After the types of projectile motion model was established, several scientists were devoted to finding a method to simulate the movement of wreckage. The main diffi-culty was to combine natural factors with the movement. Juan Santos-Echeandía introduced a differential equation model to simplify the difficulty [2]. Moreover,A. Boultif and D. Louër introduced a dichotomy iteration algorithm to circular compu-ting which can be borrowed to combine the motion of wreckage with underwater ter-rain [3]. Several conditions have to be fulfilled before simulating the movement: (1) Seawater density keeps unchanged despite the seawater depth. (2) The velocity of the wreck stay the same compared with velocity of the plane before it crashes into pieces.(3) Marine life will not affect our simulation. (4) Acting forceof seawater is a function of the speed of ocean currents.However the conclusion above cannot describe the wreckage zone accurately. This inaccuracy results from simplified conditions and ignoring the probability distribution of wreckage. In 1989, Stone et.al introduced a Bayesian search approach for searching problems and found the efficient search plans that maximize the probability of finding the target given a fixed time limit by maintaining an accurate target location probabil-ity density function, and by explicitly modeling the target‟s process model [4].To come up with a concrete dispatch plan. Xing Shenwei first simulated the model with different kinds of algorithm. [5] In his model, different searching planes are as-sessed by several key factors. Then based on the model established before, he use the global optimization model and an area partition algorithm to propose the number of aircrafts. He also arranged quantitative searching recourses according to the maxi-mum speed and other factors. The result shows that search operations can be ensured and effective.Further studies are carried out based on the comparison between model andreality.Some article illustrate the random error caused by assumptions.2 Assumptions and JustificationsTo simplify the problem, we make the following basic assumptions, each ofwhich is properly justified.●Utilized data is accuracy. A common modeling assumption.●We ignore the change of the gravitational acceleration. The altitude of anaircraft is less than 30 km [6]. The average radius of the earth is 6731.004km, which is much more than the altitude of an aircraft. The gravitational accele-ration changes weakly.●We assume that aeroengine do not work when a plane is out of contact.Most air crash resulted from engine failure caused by aircraft fault, bad weather, etc.●In our model, the angle of attack do not change in an air crash and thefuselage don’t wag from side to side. We neglect the impact of natural and human factors●We treat plane as a material point the moment it hit the sea-level. Thecrashing plane moves fast with a short time-frame to get into the water. The shape and volume will be negligible.●We assume that coefficient of air friction is a constant. This impact is neg-ligible compared with that of the gravity.●Planes will crash into wreckage instantly when falling to sea surface.Typically planes travel at highly speed and may happen explosion accident with water. So we ignore the short time.3 NotationsAll the variables and constants used in this paper are listed in Table 1 and Table 2.Table 1 Symbol Table–ConstantsSymbol DefinitionωRotational angular velocity of the earthg Gravitational accelerationr The average radius of the earthC D Coefficient of resistance decided by the angle of attack ρAtmospheric densityφLatitude of the lost contact pointμCoefficient of viscosityS0Area of the initial wrecking zoneS Area of the wrecking zoneS T Area of the searching zoneK Correction factorTable 2 Symbol Table-VariablesSymbol DefinitionF r Air frictionF g Inertial centrifugal forceF k Coriolis forceW Angular velocity of the crash planev r Relative velocity of the crash planev x Initial velocity of the surface layer of ocean currentsk Coefficient of fluid frictionF f Buoyancy of the wreckagef i Churning resistance of the wreckage from ocean currents f Fluid resistance opposite to the direction of motionG Gravity of the wreckageV Volume of the wreckageh Decent height of the wreckageH Marine depthS x Displacement of the wreckageS y Horizontal distance of S xα Deviation angle of factually final position of the wreckage s Horizontal distance between final point and splash point p Probability of a wreck in a given pointN The number of the searching planeTS ' The area of sea to be searched a i V ˆ The maximum speed of each planeai D The initial distance from sea to search planeai A The search ability of each plane is),(h T L i The maximum battery life of each plane isi L The mobilized times of each plane in the whole search )1(N Q Q a a ≤≤ The maximum number of search plane in the searching zone T(h) The time the whole action takes4 Model OverviewMost research for searching the lost plane can be classified as academic and practical. As practical methods are difficult to apply to our problem, we approach theproblem with academic techniques. Our study into the searching of the lost plane takes several approaches.Our basic model allows us to obtain the splash point of the lost plane. We focus on the force analysis of the plane. Then we We turn to simple types of projectile motion model. This model gives us critical data about the movement and serves as a stepping stone to our later study.The extended model views the problem based on the conclusion above. We run diffe-rential equation method and Bayesian search model to simulate the movement of wreckage. The essence of the model is the way to combine the effect of natural factors with distribution of the wreckage. Moreover, using distributing conditions, we treat size of the lost plane as “initial wreckage zone” so as to approximately describe the distribution. Thus, after considering the natural factors, we name the distribution of wreckage a “wreck zone” to minimize searching zone. While we name all the space needed to search “searching zone”.Our conclusive model containing several kinds of algorithm attempts to tackle a more realistic and more challenging problem. We add the global optimization model and an area partition algorithm to improve the efficiency of search aircrafts according to the area of search zone. An assessment of search planes consisting of search capabili-ties and other factors are also added. The Dinkelbach and NonConvexDivide algo-rithm for the solutions of the results are also added.We use the extended and conclusive model as a standard model to analyze the problem and all results have this two model at their cores.5 Modeling For Locating the Lost PlaneWe will start with the idea of the basic model. Then we present the Bayesian search model to get the position of the sinking point.5.1 Modeling For Locating the Splash PointThe basic model is a academic approach. A typical types of projectile behavior con-sists of horizontal and vertical motion. We also add another dimension consider-ing the effect of the earth's rotation. Among these actions, the force analysis is the most crucial part during descent from the point out of contact to the sea-level. Types of plane might impact trajectory of the crashing plane.5.1.1 Types of PlanesWe classify the planes into six groups [7]:●Helicopters: A helicopter is one of the most timesaving ways to transfer be-tween the city and airport, alternatively an easy way to reach remote destina-tions.●Twins Pistons: An economical aircraft range suitable for short distance flights.Aircraft seating capacity ranging from 3 to 8 passengers.●Turboprops: A wide range of aircraft suitable for short and medium distanceflights with a duration of up to 2-4 hours. Aircraft seating capacity ranging from 4 to 70 passengers.●Executive Jets:An Executive Jet is suitable for medium or long distanceflights. Aircraft seating capacity ranging from 4 to 16 passengers●Airliners:Large jet aircraft suitable for all kinds of flights. Aircraft seatingcapacity ranging from 50 to 400 passengers.●Cargo Aircrafts:Any type of cargo. Ranging from short notice flights carry-ing vital spare parts up to large cargo aircraft that can transport any volumin-ous goods.The lost plane may be one of these group. Then we extract the characteristics of planes into three essential factors: mass, maximum flying speed, volume. We use these three factors to abstract a variety of planes:●Mass: Planes of different product models have their own mass.●Maximum flying speed: Different planes are provided with kinds of me-chanical configuration, which will decide their properties such as flying speed.●Volume: Planes of distinct product models have different sizes and configura-tion, so the volume is definitive .5.1.2 Preparation of the Model —Earth RotationWhen considering the earth rotation, we should know that earth is a non-inertial run-ning system. Thus, mobile on the earth suffers two other non-inertial forces except air friction F r . They are inertial centrifugal force F g and Coriolis force F k . According to Newton ‟s second law of motion, the law of object relative motion to the earth is:Rotational angular velocity of the earth is very small, about .For a big mobile v r , it suffers far less inertial centrifugal force than Coriolis force, so we can ignore it. Thus, the equation can be approximated as follows:Now we establish a coordinate system: x axis z axis pointing to the east and south re-spectively, y axis vertical upward, then v r , ω and F r in the projection coordinate system are as follows:⎪⎪⎩⎪⎪⎨⎧++=⋅⋅-⋅⋅=++=kdt dz j dt dy i dt dx m v k j w kF j F i F F r rz ry rx r φωφωcos sinφis the latitude of the lost contact point of the lost plane. Put equation 1-3 and equa-tion 1-2 together, then the component of projectile movement in differential equation is:ma FF F k g r=++srad ⋅⨯=-5103.7ωmamv F r r =+ω2⎪⎪⎪⎩⎪⎪⎪⎨⎧+⋅=+⋅=+⎪⎭⎫ ⎝⎛+⋅-=m F dt dx w dt z d m F dt dx w dt y d m F dt dz dt dy w dtx d rz ry rx φφφφsin 2cos 2sin cos 22222225.1.3 ModelingConsidering the effect caused by earth rotation and air draught to plane when crashing to sea level, we analyze the force on the X axis by using Newton ‟s second law, the differential equation on x y and axis, we can conclude:In conclusion, we establish the earth rotation and types of projectile second order dif-ferential model:()⎪⎩⎪⎨⎧+-⋅'⋅⋅=''-⋅'+⋅'⋅⋅-=''-⋅'⋅⋅=''m gf y w m z m f z x w m y m f y w m x m obj 321cos 2cos sin 2sin 2.φφφφAccording to Coriolis theorem, we analyze the force of the plane on different direc-tions. By using the Newton ‟s laws of motion, we can work out the resultant accelera-tion on all directions:⎪⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎪⎨⎧'+'+'⋅'⋅⋅⨯+='+'+'⋅'⋅⋅⨯+='+'+'⋅'⋅⋅⨯+=⋅⋅-⋅⋅=⋅⨯=⋅'''⋅===-2222222225)()()(21)()()(21)()()(21cos sin 103.704.022z y x z c F f z y x y c F f z y x x c F f k j w s rad S y x F c D rz D ryD rx D ρρρφωφωωμφC D is the angle of attack of a plane flew in the best state, w is the angular speed of a moving object, vector j and k are the unit vector on y and z direction respectively,μisrx F y w m x m -⋅'⋅⋅⨯=''φsin 2()ry F z x w m y m -'+⋅'⋅⨯-=''φφcos sin 2mg F y w m z m rz +-⋅'⋅⋅⋅=''φcos 2the coefficient of viscosity of the object.5.1.4 Solution of the ModelWhen air flows through an object, only the air close to layer on the surface of the ob-ject in the laminar airflow is larger, whose air viscosity performance is more noticea-ble while the outer region has negligible viscous force [8]. Typically, to simplify cal-culation, we ignore the viscous force produced by plane surface caused by air resis-tance.Step 1: the examination of dimension in modelTo verify the validity of the model based on Newton ‟s second theorem, first, we standardize them respectively, turn them into the standardization of dimensionless data to diminish the influence of dimensional data. The standard equation is:Step 2: the confirmation of initial conditionsIn a space coordinate origin based on plane, we assume the earth's rotation direc-tion for the x axis, the plane's flight heading as y axis, the vertical downward di-rection for z axis. Space coordinate system are as follows:Figure 1 Space coordinate systemStep 3: the simplification and solutionAfter twice integrations of the model, ignoring some of the dimensionless in thesxx y i -=integral process, we can simplify the model and get the following:⎪⎪⎪⎩⎪⎪⎪⎨⎧+'⋅⋅⋅-⋅'⋅⨯='''-⋅⋅⋅-⋅'⋅⨯-=''⋅'⋅⨯=''g z m s c y w z y v m s c z w y y w x D D 220)(2cos 2)(2cos 2sin 2ρφρφφWe can calculate the corresponding xyz by putting in specific data to get the in-formation about the point of losing contact.Step 4: the solution of the coordinateThe distance of every latitude on the same longitude is 111km and the distance ofevery longitude on the same latitude is 111*cos (the latitude of this point) (km). Moreover, the latitude distance of two points on the same longitude is r ×cos(a ×pi/180) and the longitude distance of two points on the same latitude is: r ×sin(a ×pi/180)[9].We assume a as the clockwise angle starting with the due north direction and r as the distance between two points; X 、Y are the latitude and longitude coordinates of the known point P respectively; Lon , Lat are the latitude and longitude coordi-nates of the unknown point B respectively.Therefore, the longitude and latitude coordinates of the unknown point Q is:⎪⎪⎩⎪⎪⎨⎧⨯⨯+=⨯⨯⨯⨯+=111)180/cos()180/cos(111)180/sin(pi a r Y Lat pi Y pi a r X LonThus, we can get coordinates of the point of splash by putting in specific data.5.2 Modeling For Locating WreckageIn order to understand how the wreckage distributes in the sea, we have to understand the whole process beginning from the plane crashing into water to reaching the seaf-loor. One intuition for modeling the problem is to think of the ocean currents as astochastic process decided by water velocity. Therefore, we use a differential equation method to simulate the impact on wreckage from ocean currents.A Bayesian Searching model is a continuous model that computing a probability dis-tribution on the location of the wreckage (search object) in the presence of uncertain-ties and conflicting information that require the use of subjective probabilities. The model requires an initial searching zone and a set of the posterior distribution given failure of the search to plan the next increment of search. As the search proceeds, the subjective estimates of the detection will be more reliable.5.2.1 Assumptions of the ModelThe following general assumptions are made based on common sense and weuse them throughout our model.●Seawater density keeps unchanged despite the seawater depth.Seawater density is determined by water temperature, pressure, salinity etc.These factors are decided by or affected by the seawater density. Considering the falling height, the density changes slightly. To simplify the calculation, we consider it as a constant.●The velocity of the wreck stay the same compared with velocity of theplane before it crashes into pieces. The whole process will end quickly witha little loss of energy. Thus, we simplify the calculation.●Marine life will not affect our simulation.Most open coast habitats arefound in the deep ocean beyond the edge of the continental shelf, while the falling height of the plane cannot hit.●Acting force of seawater is a function of the speed and direction of oceancurrents. Ocean currents is a complicated element affected by temperature, wide direction, weather pattern etc. we focus on a short term of open sea.Acting force of seawater will not take this factors into consideration.5.2.2 Preparation of the Model●The resistance of objects of different shapes is different. Due to the continuityof the movement of the water, when faced with the surface of different shapes, the water will be diverted, resulting in the loss of partial energy. Thus the pressure of the surface of objects is changed. Based on this, we first consider the general object, and then revise the corresponding coefficients.●Ocean currents and influencing factorsOcean currents, also called sea currents, are large-scale seawater movements which have relatively stable speed and direction. Only in the land along the coast, due to tides, terrain, the injection of river water, and other factors, the speed and direction of ocean currents changes.Figure 2Distribution of world ocean currentsIt can be known from Figure 2 that warm and cold currents exist in the area where aircraft incidences happened. Considering the fact that the speed of ocean currents slows down as the increase of the depth of ocean, the velocity with depth sea surface currents gradually slowed down, v x is set as the initial speed of ocean currents in subsequent calculations.●Turbulent layerTurbulent flow is one kind of state of the fluid. When the flow rate is very low, the fluid is separated into different layers, called laminar flow, which do not mix with each other. As the flow speed increases, the flow line of the fluid begins to appear wavy swing. And the swing frequency and amplitude in-creases as the flow rate increases. This kind of stream regimen is called tran-sition flow. When the flow rate becomes great, the flow line is no longer clear and many small whirlpools, called turbulence, appeared in the flow field.Under the influence of ocean currents, the flow speed of the fluid changes as the water depth changes gradually, the speed and direction of the fluid is un-certain, and the density of the fluid density changes, resulting in uneven flow distribution. This indirectly causes the change of drag coefficient, and the re-sistance of the fluid is calculated as follows:2fkvGLCM texture of submarine topographyIn order to describe the impact of submarine topography, we choose a rectan-gular region from 33°33…W, 5°01…N to 31°42‟W , 3°37‟N. As texture is formed by repetitive distribution of gray in the spatial position, there is a cer-tain gray relation between two pixels which are separated by a certain dis-tance, which is space correlation character of gray in images. GLCM is a common way to describe the texture by studying the space correlation cha-racter of gray. We use correlation function of GLCM texture in MATLAB:I=imread ('map.jpg'); imshow(I);We arbitrarily select a seabed images and import seabed images to get the coordinate of highlights as follows:Table 1Coordinate of highlightsNO. x/km y/km NO. x/km y/km NO. x/km y/km1 154.59 1.365 13 91.2 22.71 25 331.42 16.632 151.25 8.19 14 40.04 18.12 26 235.77 13.93 174.6 14.02 15 117.89 14.89 27 240.22 17.754 172.38 19.23 16 74.51 12.29 28 331.42 24.455 165.71 24.82 17 45.6 8.56 29 102.32 19.486 215.75 26.31 18 103.43 5.58 30 229.1 18.247 262.46 22.96 19 48.934 3.51 31 176.83 9.188 331.42 22.34 20 212.42 2.85 32 123.45 3.239 320.29 27.55 21 272.47 2.48 33 32.252 11.7910 272.47 27.55 22 325.85 6.45 34 31.14 27.811 107.88 28.79 23 230.21 7.32 35 226.88 16.0112 25.579 27.05 24 280.26 9.93 36 291.38 5.46Then we use HDVM algorithm to get the 3D image of submarine topography, which can be simulated by MATLAB.Figure 3 3D image of submarine topographyObjects force analysis under the condition of currentsf is the resistance, f i is the disturbance resistance, F f is the buoyancy, G isgravity of object.Figure 4Force analysis of object under the conditions of currentsConsidering the impact of currents on the sinking process of objects, wheninterfered with currents, objects will sheer because of uneven force. There-。
美赛范文

24914
Problem Chosen
B
2014 M athem atical Contest in M odeling (M CM ) Sum m of evaluating the best college coach ever lies largely in evaluating the performance of a college coach reasonably, especially when we consider both genders and the time factor among various sports. Existing models, however, do not take into account the above factors thoroughly. This problem can be attributed to an evaluation problem. We present a novel solution called “Weighted Topological Sorting Model” to setting up a scientific assessment system of indexes to choose the top 5 coaches in the history of three professional sports. We firstly obtain the data of all college coaches. In particular, we select five representative indicators tightly related the state of coaches and their own backgrounds by means of clustering and analysis. Then we formulate the Analytical Hierarchy Process (AHP) to get the ranking of coaches. In order to reduce the subjectivity of AHP, we combined AHP with Triangular Fuzzy Function. In particular, considering the time line horizon, we make the weight of different periods of coaching based on Life Cycle Theory. We conclude a more rational solution. We show that this strategy is not optimal for its instability but can be improved by establishing a Fuzzy Neutral Network Model using the above fuzzy rules we have got. It also can reduce evaluation errors. Moreover, we develop a new evaluation without traditional and external individual evaluation. This is called Weighted Topological Sorting Model. By distinguishing teams of different genders, we compare each two coaches to obtain the directed and weighted network of contacts. We compare the data with The ESPN Rankings to verify our models. Finally, based on the analysis about the strengths and weaknesses of the proposed models, we have also considered the refinement of our models, taking unobservable factors into account.
美赛 文献综述 开头

美赛文献综述开头
文献综述是科研工作中的重要环节,通过对已有文献的整理、分析和总结,可以帮助研究者更好地把握前人工作的现有成果和不足之处,为自己的研究提供指导和启示。
本文将结合最新的研究进展,就美国数学建模大赛(MCM/ICM)的相关文献进行综述,希望能为广大参赛者提供一些有益的启发。
首先,美赛作为一个知名的数学建模竞赛,每年吸引了来自全球各地的数千名参赛者。
参赛者需要在规定的时间内,针对所提出的问题进行深入的研究和分析,最终给出合理的建模方案和解决方案。
因此,对于参赛者而言,对前人的研究成果进行充分的了解是至关重要的。
其次,通过对美赛历年获奖作品的文献进行分析,可以发现一些共性和规律。
比如,优秀的团队往往具有良好的团队协作能力和高效的分工合作方式;在问题建模和解决过程中,他们通常能够充分挖掘数据的潜在规律,结合数学工具和计算机算法,提出创新性的解决方案。
这些都是参赛者可以借鉴和学习的经验和教训。
最后,希望参赛者在阅读文献综述的过程中,不仅能够了解前人的研究成果和方法,还要结合自身的实际情况和研究对象,找到适合自己的解题思路和技巧。
同时,要注重团队协作能力的培养和提升,不断学习和改进自己的数学建模能力,为未来的参赛之路奠定坚实的基础。
综上所述,文献综述在美赛准备过程中具有重要的指导意义,通过对前人的研究成果进行深入的分析和总结,可以帮助参赛者更好地把握竞赛的要点和难点,为取得优异的成绩提供有力的支持。
希望广大参赛者能够认真对待文献综述,积极学习和借鉴前人的经验,不断提升自身的数学建模能力,为美赛的成功参赛打下坚实的基础。
《美赛论文写作》课件

掌握PPT制作技巧,制作高 效且展示专业化的幻灯片。
通过学习和实践,在美赛中 展示专业、高效的文学写作 和幻灯片制作技巧。
问题解答
1 答疑解惑
2 课程反馈
回答学员在学习和实践中遇到的问题,保障学员 实现学习效果。
收集学员对该课程的反馈,反馈结果将用于后续 课程的改善和优化。
5
论文基本要素
介绍论文的基本组成部分和写作要求,引导 选手开始论文写作。
语言表达与修辞
着重讲解论文细节部分,从选题到修辞手法, 给出一些可供参考的写作技巧。
参考文献引用与写作规范
介绍如何引用文献、如何写参考文献和常见 的写作规范。
PPT制作技巧
PPT总体设计理念
从PPT设计的整体思路出发,介绍PPT主题、设计风格 和模板选择等方面的技巧。
给出竞赛的具体时间节点,包括摸底赛和正式比赛。
竞赛准备
列一些准备工作和注意事项,包括选题、下载数据 等。
美赛论文写作策略
1
内容组织与结构分析
2
讲解论文结构的设计、关键信息的呈现,并
且给出一些有效性的写作技巧。
3
展示数据与图表设计
4
有效、清晰的表达数据是论文的重要组成部
分,此节介绍如何呈现数据及其图表设计。
图形设计与视频/音频插入
介绍图形设计、视频、音频、图片等多媒体素材的插 入和处理。
制作专业的动画效果
讲解如何利用PPT自带的动画效果和过渡效果,设计制
如何在团队中协作
介绍如何利用PPT自带的分享功能,在团队中协作
3 展示专业化水平
掌握论文写作技巧,将论文 写得更加精简、清晰、合理。
《美赛论文写作》PPT课 件
这份PPT课件为你准备了美赛论文写作的攻略和PPT制作技巧。从文章结构和 细节到PPT设计和协作分享,掌握技巧,制作专业的论文和幻灯片。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
0 假设条件与解释:
Eg: 2010年的MCM“centroids, clusters, and crime: anchoring the geographic profiles of serial criminals” 中,参赛小组在“假设条件”这一节中逐条列出了建模中的 假设条件,并给出了假设条件的合理性论证。 “criminal’s movement is unconstrained :because of the difficulty of finding real-world distance data, we invoke the ‘Manhattan assumption ’: there are enough streets and sidewalks in a sufficiently grid-like pattern that movements along real-world movement routes is the same as ‘straight-line’ movement in space discretized into city blocks…” (罪犯的活动不受限制:由于实际数据不易获得,本文将采 用“曼哈顿假设”来描述罪犯的活动方式,也就是说,假设 市区内有足够多的街道及人行道纵横交错地构成一个与离散 方格相似的网状布局,罪犯在此网状布局中的活动路径(即 沿着方格的边作运动)与现实世界中的“直线运动”相同…)
论文结构-小节划分
0 重述并澄清赛题 0 列出建模所用的所有前提条件及假设,并给出清晰的
解释。 0 分析赛题,给出建模动机或论证建模的合理性。 0 模型设计 0 讨论如何验证模型,包括误差分析和稳定性测试(如 对条件,敏感度等因素进行分析和测试) 0 讨论模型的优缺点。 0 书写论文的摘要。(摘要必须按要求写在特定的摘要 页上,长度不超过一页,提交时作为参赛论文的首页)
0 即便是已经表述的很精确的概念,仍然可以给出有利于
解题的解读方式。 0 Eg: 对精确到2°以内的解释: “we interpret the error of 2°as a normal distribution,…with standard deviation of 1°” (我们将精确到±2°这个概念理解成一个正态分布,。。。 其标准偏差为1°)
0 一个问题可能存在多种解读的方式,所以书写时必须明确表
述他们是如何解读赛题的。 0 Eg:2010年的MCM中,有一组题要求参赛小组算出棒球棍上 是最佳击球点。可是“最佳击球点”这个概念有多种解释。 在题为“the sweet spot: a wave model of baseball bats”的引 言部分,参赛小组是这样解释最佳击球点的: 0 “there are at least two notions of where the sweet spot should be --- an impact location on the bat that either minimizes the discomfort to the hands, or maximizes the outgoing velocity of the ball We focus exclusively on the second definition. ” (最佳击球点至少存在两种不同的解释: 一种是让击球者的手感觉最舒适的击球位置。 一种是使棒球飞出速度最快的击球位置。 本文只考虑第二种解释)
0 无论解答什么样的赛题,参赛小组都应该在论文中明确
列出所有用到的假设条件,并解释其合理性。如果对某 个假设无法给出满意的解释,则应重新考虑这个假设是 否合理,并经行修改。 模型设计: 在所有能够解决问题的模型中,最简单的模型也许就是 最好的模型。下面是介绍两种建立模型的思路。 系列模型 设计模型时,可以尝试从简单模型开始,逐步加工,修 改及完善,一次比一次更接近现实,最终得到能够圆满 解决问题的模型。 Eg:2008年的MCM竞赛要求建立一个数学模型,预测在全 球温度升高而导致北极冰盖融化的假设下,美国的沿海 地区将会受到怎样的影响。“the impending effects of north polar ice cap melt”,就是采取建立一系列模型的 方法最终给出满意的解答。具体做法如下:
① 首先建立模型1:常温模型。这个模型假设全球温度不变,
冰盖的融化速度不变,以及海洋水量不变。 ② 接着建立模型2:变温模型。这个模型假设全球温度在不 断变化。 ③ 最后建立模型3:气候变暖下的海洋水量模型。这个模型 将前面建立的模型中所忽视的问题考虑进来,包括南北 两半球的相对海洋水面积。 建模时,参赛小组应集中精力设计一个模型,或者最终能导 出一个较好的模型的一系列子模型。不要分散精力,拼凑出 几个很一般的模型,企图用这些一般的模型碰运气。这种做 法会分散评委的注意力,不被看好。
0 摘要应达到吸引读者进一步阅读论文的目的,不要用
“first we…then we…” 这样干巴巴的句子重复论文中 的内容,也不要在开头使用类似“this paper will solve..”这样的句子。摘要的第一句话尤其重要,应该 以吸引人的语言激发读者的兴趣。Eg: 2008年MCM关 于数独游戏的赛题“taking the mystery out of sudoku difficulty: an oracular model”中摘要的第一句话就达到 这样的效果:“ in the last few years, the 9-by-9 puzzle grid known as sudoku has gone from being a popular Japanese puzzle to a global craze.” 0 (在过去的几年中,一个在9乘9方阵上玩的游戏,即 数独游戏,在日本流行后风靡全球)
论文结构 – 如何写结论 conclusions
0 论文的结论部分应描述模型的测试方法以及所用到的
结果。建模是解决问题的手段,即使赛题没有标准答 案,参赛小组也应该在所建立模型的基础上进行分析 并且得出明确的结论。结论应该在结论部分的开始就 明确的写出来。(写的好的论文不需要读者自己去找 结论)。此外,这一节还应对结论经行适当的讨论, 并指明那些是作者的观点,那些是得到的结果。、 0 参赛小组应该对所建立的数学模型经行敏感度分析和 稳定性测试,使模型更具说服力。模型通常会用到一 些参数,在结论中应该讨论这些参数值的轻微变化对 模型及结论产生的影响。
论文结构--写好引言(introduction)
0 引言应该包括以下内容:对赛题的解读,对现有研究结果的综
述与评论以及对解题思路和主要方法的简要介绍。 0 时间安排:应该在解题时就开始思考如果写引言,并在解题过 程中反复修改,直到满意为止。 0 引言的第一句话是全文最重要的句子,用于激发读者阅读论文 的兴趣。同理,引言的第一段话也是最重要的一个段落,应该 写的浅显易懂,尽量少用或者不用数学表达式,使得即使读者 对所要解决的数学问题不了解,也能看懂第一段话并产生继续 阅读的兴趣。 0 赛题选定后,首先要用自己的语言重述赛题,明确解题目标, 并澄清原题叙述上可能出现的模糊概念。
论文结构 –写好摘要
0 摘要应该包含以下内容: ① 赛题重述与阐明:用自己的语言描述将要解决的问 ② ③ ④ ⑤
题。 解释假设条件及其合理性:强调建模所用的假设, 而且清楚地列出建模所需的所有变量。 模型设计及合理性论证:指出所用模型的类型或构 造新的模型。 描述模型的测试情况及灵敏度分析:包括误差分析 等检测项目。 有缺点讨论:包括模型及解题方法的优点与不足。
论文结构 – 论文主体
0 论文主体是用于描述模型的设计,包括列出和论证建
立数学模型所做的假设,以及对赛题的分析。(论文 主体通常分为若干小节。) 0 在描述模型之前,应该将模型设计所用的假设条件一 一列出并解释清楚。不要有未经说明的假设,此外, 还应该对建模的初衷和动机适当地加以讨论。 0 下面是对如果解释假设条件,并讨论模型设计的两种 思路:一种是通过建立一系列的子模型一步一步地构 造最终模型。另一种是以赛题为特例构构造普通模型。
0 普遍模型:
模型设计的另一个要注意的问题就是,不要只针对赛题 给出的参数值设计模型。高水平的论文通常会把赛题看 成一个普遍问题的特例,首先探讨普遍问题的求解,然 后再对赛题这一特例给出具体的解答。 Eg: 1988年的MCM要求参赛小组设计一个算法,根据火 车货运车厢的长度,宽度,高度以及载重量,计算如何 在两节车厢装载大小不一的货物。“loading two railroad cars” 就是这样处理的,论文中作者明确指出: “we have produced a general algorithm to solve this type of problem, but for our problem a relationship exists that greatly simplified the algorithm.” (本文给出了一个解决此类问题的通用算法。利用赛题 给定的参数数值,本算法将大为简化)
美国大学生数学建模竞赛等 级评审标准
0 评审要点 0 是否对赛题给出了满意的解读方式,并对赛题中可能出现的概 0 0 0 0 0 0 0
念给予了必要的澄清; 是否明确列出了建模用到的所有条件及假设,并对其合理性给 出了解释或论证; 是否通过对赛题的分析给出了建模的动机或论证了建模的合理 性; 是否设计出了能有效地解答赛题的模型; 是否对模型给出了稳定性测试; 是否讨论了模型的优缺点,并给出了清晰的结论; 是否给出了圆满的摘要。 没有全部完成解答的论文不但是可接受的,而且如果在某些方 面有创新并有独到之处,仍然可能获得较好的评审结果。
0 在描述结论的时候,应设法使读者认同论文给出的解答,
尽管该解答不一定是最好的。有些问题的解答可能是一个 计算机程序,在这种情况下,应该对程序进行多次测试, 并且每次的测试数据应稍有不同。测试时应尽可能的使用 真实的数据,避免因为人造数据而引起读者对结论的怀疑。 0 在描述结果时一定要给出足够的信息,如果某些结果是由 计算机程序设计计算出来的,则必须指明所使用的计算机 程序。 0 在写结论的时候应该明确的指出模型中的缺陷,表明参赛 小组不但知道模型所含的缺陷,而且也思考过如何改进和 修补。论文缺陷可能在前面描述模型的时候已经提到过, 在这里应该再次指出来。 0 由于时间和计算资源的限制,参赛小组可能无法完成有能 力解决的一些问题。所以,写结论的时候应该明确的指出, 假如有充足的时间和计算资源,参赛小组能够解决的问题。