计算方法大作业可研报告.docVIP

计算方法大作业 一、计算方法学后感 计算方法学习主要内容为函数逼近论，数值微分，数值积分，误差分析等。常用方法有迭代法、差分法、插值法、有限元素法等，印象深刻的是三次样条差值函数，大概是做了作业的缘故。 计算也称数值分析，数学不好的人也没有学会太多东西，不过还是有些小体会的：1主要矛盾和次要矛盾之间的关系。现实问题中有很多约束条件，需要我们有侧重的保留摒弃，辨析主要矛盾和次要矛盾，从而提出合理假设； 2尺有所短寸有所长。没有完美无缺的算法，虽然我们看到有不断地改进优化算法，但这些往往都是以牺牲某些优点为代价的。比如提高精度，往往会导致格式复杂，产生较大运算量； 3原则不能变。算法也是要讲原则的，比如要谈算法的优劣性前提是要保证算法的可靠性（相容、收敛、稳定等）。 门外汉的感受也就这水准了。至于计算方法思想在实际生活中的应用，我想来想去就这些了。 天气预报：天气会受各种因素的影响，稍微一些因素发生改变就会产生很大的变化，所以天气预报其实是一件比较困难的工作古代人们用占卜或者经验总结等方式来预计天气状况，这倒更像是统计学。而有了计算机，我们就可以通过数值模拟来预报天气。具体过程就是：首先根据大气运动列出数学物理方程（偏微分），其次对空间进行网格划分，然后通过观测数据给出初值条件，最后通过数值方法求解这些偏微分方程得到网格点处的数值解。这也是为什么主持人总是说大概在...地区，大致在...时段，可能有...量级的降水...因为时空是连续的，而网格划分不可能无限密，所得的数值解也存在误差等离子体：对等离子体现有的理论描述中，磁流体力学、符拉索夫方程、福克-普朗克方程等都是微分方程，包含很多参量，如果要求出解析解，物理模型往往需要过分简化以至于无法精确和全面的包罗各种效应，所以需要数值计算，这也是等离子体物理学研究中很重要的一个方面。比如最简单的单粒子模型，它的牛顿洛伦兹方程是这样：就是一个二阶微分方程可以用数值分析里的隐式欧拉，显式欧拉，中点格式，龙格库塔等方法来求解。 水槽的原理分析 图 耦合水槽 2.对耦合水槽建模： (3-39) 其中，为1号和2号水槽的截面积，，分别是两个水槽的液面高度，,分别是两个水泵向水槽中注水的体积流速（），,,分别是与,,对应的比例常数，通过升高或者降低挡板可以控制两个水槽之间的流速，关闭2号水泵，通过1号水泵向1号水槽中注水以控制2号水槽液面高度，因此： ， (3-40) 初始化： , , (3-41) 设置最大控制信号，即 (3-42) 采样时间为1秒。 分别使用上一节中的三个模糊控制系统对耦合水槽进行仿真控制，设2号水槽中目标液面为ref=15cm，一型模糊PI控制系统手动调节参数如下： ， IT2模糊PI控制系统选取参数如下： ，， 模型验证： 200s内的仿真结果如图3.8所示，其中： 1）T1 PI FS表示一型模糊PI控制系统， 2）IT2 PI FS1表示使用KM迭代降型的区间二型模糊PI控制系统， 3）IT2 PI FS2表示使用Wu-Mendel近似法降型的区间二型模糊PI控制系统， 三个控制系统对应的性能指标如表3.1中所示，使用ITAE（时间乘绝对误差的积分）作为衡量控制系统的性能标准。从图3.8及表3.1的性能指标中，可以看出，使用KM迭代降型法的IT2模糊PI控制系统性能最优，然后为使用Wu-Mendel近似降型法的IT2模糊PI控制系统，一型模糊PI控制系统的性能最差。与一型模糊PI控制系统相比，使用KM迭代降型的IT2模糊PI控制系统在保持上升时间（从稳态值10%上升至90%所需时间）不增加的情况下，减小了超调量和调节时间，而使用Wu-Mendel近似降型的IT2模糊PI控制系统在上升时间不明显增加时，超调量和调节时间只是略微下降越小越好。为起始时间，为终止时间，为实际响应，为模型响应。 即使模型响应曲线越接近于实际响应曲线越好。 图3.8 耦合水槽控制仿真控制系统仿真对比 表3.1不同控制系统的性能对比 系统类型 T1模糊系统 IT2模糊系统 KM Wu-Mendel 上升时间(秒) 30 30 31 超调量 4.6289% 2.4333% 4.2686% 调节时间(秒) 93 73 81 性能指标(ITAE) 6747.6 5927.9 6384.4 感受：虽然很多东西都不是自己想出来的，经过很多