一类严格线性不等式组的解法.pdfVIP

独创性说明 作者郑重声明：本硕士学位论文是我个人在导师指导下进行的研究工 作及取得研究成果。尽我所知，除了文中特别加以标注和致谢的地方外， 论文中不包含其他人已经发表或撰写的研究成果，也不包含为获得大连理 工大学或者其他单位的学位或证书所使用过的材料。与我一同工作的同志 对本研究所做的贡献均已在论文中做了明确的说明并表示了谢意。 、， 作者签名： Z远 日期： 鲤￡!渔 大连理工大学硕士研究生学位论文 大连理工大学学位论文版权使用授权书 本学位论文作者及指导教师完全了解“大连理工大学硕士、博士学位 论文版权使用规定”，同意大连理工大学保留并向国家有关部门或机构送 交学位论文的复印件和电子版，允许论文被查阅和借阅。本人授权大连理 工大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，也 可采用影印、缩印或扫描等复制手段保存和汇编学位论文。 作者签名： 主堑坠 导师签名．一丕亟二 且年』月盟日 大连理工大学硕士学位论文 1 绪论 线性不等式组问题在数学、物理、计算机应用以及管理等多种领域中都有广泛的应 用。很多实际问题都可以应用到或者转化为线性不等式组问题，如[1]中的微分包含关 于非光滑区域生存性的判别问题，[2]中的移动通信系统中的功率最优控制问题，[3]中 的钻井布局模型的构造问题，[4]中的医药工作中调配辅料食谱、选方配药、选择医疗 方案等问题，以及[5]中的基于效用最大化的投资组合模型等等。 经过一系列处理，弹簧振动系统的反二次特征值问题 MXA2+CXA+KX=0 (1．1) 其中，A、X为给定的特征值与特征向量，要求的M、C、K分别为质量、阻尼、弹性系 数矩阵(c≠0)，可转化为如下线性不等式组(MoodyChu，2006年8月在大连理工大学 的讲座)： (≯枷。 m：， 其中Y，z为给定的上三角矩阵。不同于一般的线性不等式组问题，(1．2)要求所有的不等 式均为严格不等式。 解决线性不等式组问题有以下两种思路：～是考虑线性不等式组的相容性，即先判 断线性不等式组是否有解，在此基础上再给出适当的算法求解：一是不考虑不等式组的 相容性，直接给出算法求解，由计算结果来判断原不等式组是否有解。基于第一种思路， 关于不等式组相容性的判别已经有了很多结果，如[6]中的ABs～SG算法，[7]中的关于s 矩阵的～些结论，等等；但这些方法只是对某些特殊结构的线性不等式组给出了相容性 的充分或必要条件，对结构一般的线性不等式组则没有统一的结果。而基于第二种思路， 比较常用的是将线性不等式组问题看作线性规划问题，用单纯形方法[8]或者内点法[9] 求解。 将(1．2)看作线性规划问题来求解容易得到边界点。鉴于这种情形，我们希望找到 别的方法来避免产生上述情况。我们不针对(1．2)中的特殊矩阵，而是考虑更一般的情况， 将(1．2)抽象为一般形式： 问题(P)：给定A∈Rm，求x∈R”， S．t．J出刈．．{ (1．3)(．) l x0 由于(P)中的不等式均为严格不等式，故若将其看作线性规划问题用单纯形方法来解 f将目标函数设为常数，式(1．3)作为约束条件)，根据其最优解必在可行域顶点处得到可 知，容易得到边界点，即，使得(P)中某些不等式变为等式，故单纯形方法不可行，具体 例子将在第四章给出。因此，我们考虑将其转化为Minimax问题来求解。 本文的安排如下：下一章我们介绍相关理论，第三章给出解严格线性不等式组问题 口)的几个解法，最后，第四章我们将给出数值实验。