发明名称 一种免耕开沟器结构优化与作业效果预测方法
摘要 本发明涉及一种免耕开沟器结构优化与作业效果预测方法,包括以下步骤:1)决策变量与指标的确定;2)数据采集:收集成熟地区的决策变量与指标的数据;3)构建开沟器结构优化及开沟效果预测的多目标优化模型;4)收集待免耕播种地区土壤条件和作业条件等数据以重构多目标优化模型,得到开沟器的最优结构参数;5)将4)中所述的待免耕播种地区土壤数据、作业条件数据和开沟器结构参数,代入到各指标目标函数中对作业效果进行预测。本发明为待免耕播种地区的开沟器结构参数进行优化和作业效果进行预测,整个过程均通过计算机模拟完成,无需进行田间试验且省时、省力,提高了开沟器结构优化效率和作业效果,促进保护性耕作技术的推广。
申请公布号 CN105225001A 申请公布日期 2016.01.06
申请号 CN201510604248.3 申请日期 2015.09.21
申请人 中国农业大学 发明人 王庆杰;张祥彩;李洪文;何进;李问盈;郑智旗;王宪良;尤炳晓;郑侃
分类号 G06Q10/04(2012.01)I;G06Q50/02(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 陈波
主权项 一种免耕开沟器结构优化与作业效果预测方法,其特征在于:该方法主要包括以下步骤:1)决策变量和指标的选择:根据成熟地区和待免耕播种地区开沟器的结构特点和开沟效果的要求,确定多目标优化模型的决策变量和指标;所述决策变量包括开沟器结构参数、土壤数据、作业条件数据,所述指标包括作业效果、作业能耗、工作阻力;2)数据采集:按照步骤1)中确定的决策变量和指标,收集成熟地区的决策变量和指标的数据;3)构建开沟器结构优化及开沟效果预测的多目标优化模型:将步骤2)中采集的各决策变量数据进行标准化处理,根据标准化的决策变量数据,确定各指标对应的目标函数及权重系数,构建开沟器结构优化及开沟效果预测的多目标优化模型,按以下步骤进行:a)数据标准化:首先根据步骤2)中采集的决策变量数据构成原始数据矩阵X<sup>*</sup>,将原始数据矩阵X<sup>*</sup>转化为标准化矩阵X,其中x<sub>i</sub>为标准化矩阵X中的元素,表示标准化后第i个决策变量值,i=1,2,…,n;b)构造各指标的目标函数:通过标准化处理的决策变量数据,由Matlab软件对数据进行拟合,分别建立以决策变量为变量、各指标为因变量的拟合方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>F</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>&lsqb;</mo><msub><mi>a</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>c</mi><mrow><mi>k</mi><mi>i</mi><mi>j</mi></mrow></msub><msubsup><mi>x</mi><mi>i</mi><mi>j</mi></msubsup><mo>&rsqb;</mo><mo>,</mo><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>m</mi><mo>;</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000807692430000011.GIF" wi="1419" he="198" /></maths>根据各决策变量的要求,确定各决策变量的约束条件:g(x<sub>i</sub>)≤0(i=1,2,…,n);得到各指标的目标函数为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mo>{</mo><mtable><mtr><mtd><mrow><mi>min</mi><mi> </mi><msub><mi>F</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>&lsqb;</mo><msub><mi>a</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>c</mi><mrow><mi>k</mi><mi>i</mi><mi>j</mi></mrow></msub><msubsup><mi>x</mi><mi>i</mi><mi>j</mi></msubsup><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mi>g</mi><mi>h</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&le;</mo><mn>0</mn></mrow></mtd><mtd><mrow><mo>(</mo><mi>h</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mtd></mtr></mtable></mtd></mtr></mtable><mo>,</mo><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>m</mi><mo>;</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000807692430000021.GIF" wi="1611" he="373" /></maths>同时,各指标的向量目标函数为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>min</mi><mi>U</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>min</mi><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>F</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>F</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>F</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mi>g</mi><mi>h</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&le;</mo><mn>0</mn></mrow></mtd><mtd><mrow><mo>(</mo><mi>h</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mtd></mtr></mtable></mtd></mtr></mtable></mfenced><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000807692430000022.GIF" wi="1601" he="276" /></maths>其中F<sub>k</sub>(x)被同等地极小化;式中,F<sub>k</sub>(x)表示第k个指标的目标函数中决策变量所拟合的函数部分,a<sub>ki</sub>表示第k个指标中第i个决策变量的系数,<img file="FDA0000807692430000023.GIF" wi="70" he="90" />表示第k个指标中第i个决策变量的j次方数,c<sub>kij</sub>表示<img file="FDA0000807692430000024.GIF" wi="73" he="97" />的系数,j表示第k个指标中第i个决策变量的方数,p为最高次数;s.t.g<sub>h</sub>(x<sub>i</sub>)≤0表示决策变量应满足所有的约束条件;c)确定各指标的权重系数:依照F<sub>1</sub>(x),F<sub>2</sub>(x),…,F<sub>m</sub>(x)在整个多目标优化模型中的重要程度,相应地给出一组加权系数w<sub>1</sub>,w<sub>2</sub>,…,w<sub>m</sub>;采用层次分析法,按各指标的顺序构造两两比较的m阶判断矩阵B:<img file="FDA0000807692430000025.GIF" wi="1619" he="346" />其中,b<sub>αβ</sub>﹥0;b<sub>αβ</sub>=1/b<sub>βα</sub>;α=β时,b<sub>αβ</sub>=1;b<sub>αβ</sub>是第α个指标对第β个指标的相对重要性,其值是根据资料数据、经验以及实际要求来确定;由特征向量W=[w<sub>1</sub>,w<sub>2</sub>,…,w<sub>m</sub>]<sup>T</sup>得出各指标的权重系数为w<sub>1</sub>,w<sub>2</sub>,…,w<sub>m</sub>;取其界限为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>w</mi><mi>k</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo><msub><mi>w</mi><mi>k</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000807692430000031.GIF" wi="905" he="179" /></maths>d)构建开沟器结构优化及开沟效果预测的多目标优化模型:运用向量目标函数最小化对目标函数进行优化,要使各指标的目标函数都达到最优的效果,采用多目标优化方法来解决,用F<sub>k</sub>(x)与w<sub>k</sub>(k=1,2,…,m)的线性组合构成一个多目标优化模型,其表达式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>{</mo><mtable><mtr><mtd><mrow><mi>min</mi><mi> </mi><msup><mi>U</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>min</mi><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>w</mi><mi>k</mi></msub><msub><mi>F</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mi>g</mi><mi>h</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&le;</mo><mn>0</mn></mrow></mtd><mtd><mrow><mo>(</mo><mi>h</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mtd></mtr></mtable></mtd></mtr></mtable><mo>,</mo><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>m</mi><mo>;</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000807692430000032.GIF" wi="1600" he="357" /></maths>即将多目标优化问题转化为单目标优化问题;4)开沟器结构优化:根据待免耕播种地区作业效果及作物生长的要求,查阅待免耕播种地区开沟器结构参数和作业条件数据,测量待免耕播种地区的土壤数据,将收集的数据代入到步骤3)中的多目标优化模型,重新构建待免耕播种地区的开沟器结构优化及开沟效果预测的多目标优化模型,利用Matlab软件求解待免耕播种地区所用开沟器的最优结构参数;5)开沟作业效果预测:将步骤4)中所述的待免耕播种地区的土壤数据、作业条件数据及开沟器最优结构参数,代入到步骤3)中各指标的目标函数中,预测开沟作业效果。
地址 100094 北京市海淀区圆明园西路2号