计算方法大作业讲解.docx
《计算方法大作业讲解.docx》由会员分享，可在线阅读，更多相关《计算方法大作业讲解.docx（14页珍藏版）》请在冰豆网上搜索。
计算方法大作业讲解
计算方法大作业
一、计算方法学后感
计算方法学习主要内容为函数逼近论，数值微分，数值积分，误差分析等。
常用方法有迭代法、差分法、插值法、有限元素法等，印象深刻的是三次样条
差值函数，大概是做了作业的缘故。
计算也称数值分析，数学不好的人也没有学会太多东西，不过还是有些小
体会的：
1>主要矛盾和次要矛盾之间的关系。
现实问题中有很多约束条件，需
要我们有侧重的保留摒弃，辨析主要矛盾和次要矛盾，从而提出合理假设；
2>尺有所短寸有所长。
没有完美无缺的算法，虽然我们看到有不断地改进
优化算法，但这些往往都是以牺牲某些优点为代价的。
比如提高精度，往往会
导致格式复杂，产生较大运算量；
3>原则不能变。
算法也是要讲原则的，比如要谈算法的优劣性前提是要保
证算法的可靠性（相容、收敛、稳定等）。
门外汉的感受也就这水准了。
至于计算方法思想在实际生活中的应用，我
想来想去就这些了。
（1）天气预报：
天气会受各种因素的影响，稍微一些因素发生改变就会产
生很大的变化，所以天气预报其实是一件比较困难的工作。
古代人们用占卜或
者经验总结等方式来预计天气状况，这倒更像是统计学。
而有了计算机，我们
就可以通过数值模拟来预报天气。
具体过程就是：
首先根据大气运动列出数学
物理方程（偏微分），其次对空间进行网格划分，然后通过观测数据给出初值
条件，最后通过数值方法求解这些偏微分方程得到网格点处的数值解。
这也是
为什么主持人总是说大概在...地区，大致在...时段，可能有...量级的降水...因
为时空是连续的，而网格划分不可能无限密，所得的数值解也存在误差。
（2）等离子体：
对等离子体现有的理论描述中，磁流体力学、符拉索夫方
程、福克-普朗克方程等都是微分方程，包含很多参量，如果要求出解析解，物
理模型往往需要过分简化以至于无法精确和全面的包罗各种效应，所以需要数
值计算，这也是等离子体物理学研究中很重要的一个方面。
比如最简单的单粒
子模型，它的牛顿洛伦兹方程是这样：
。
这就
是一个二阶微分方程，可以用数值分析里的隐式欧拉，显式欧拉，中点格式，
龙格库塔等方法来求解。
前不久有个应用可以计算出家族人物的关系，感觉挺有意思的。
我就考虑
了一下它的算法步骤，编程是硬伤，所以只能先看看。
1.先计算目标人物和我在家族中的辈分差
'我，哥哥，姐姐，弟弟，妹妹，丈夫，妻子‘辈分=0
'爸爸，妈妈‘辈分=1
'儿子，女儿‘辈分=-1
比如：
'妈妈的爸爸的哥哥‘然后计算辈分差+1+1+0=2，说明
目标人物比我大了两辈。
2.确定目标人物的性别
比如：
”哥哥的儿子的女儿“然后，目标人物只跟xxxxx的**,性别只跟
相关，于是得到目标人物性别。
**
确定目标人物在家族中处在什么分支中。
3.
大概只有，和爸爸一边，和妈妈一边，和丈夫一边，和妻子一边，和女儿
一边，和儿子一边。
”哥哥的爸爸的爸爸“=比如：
”哥哥的哥哥的哥哥的爸爸的爸爸“
最后确定目标人物和爸爸一边。
确定的三个参数好像就能确定了吧，根据，12，3
水槽、21水槽的原理分
#探#探挡
泵#2泵#1水槽#2#1水槽2出口1出口
图耦合水槽
系统由两个相同的水槽组成。
两个水槽通过阀门连接在一起，当阀门Q3
关闭时，两个水槽都能够相互独立的工作，互不影响而当阀门打开时，Q3Q3;
两个水槽的液位之间相互耦合，水槽之间的流量（）受两个水槽液位的影响。
Q12流入水槽的流量受泵的控制，对泵的控制是通过调节电动机的电枢电压来DC
实现的。
水槽的液位通过压电传感器测量，所测得的输出电压值与液位成比例
关系。
互耦水槽系统虚拟实验室软件界面图
软件界面分为两部分：
上半部分为互耦水槽系统实验模拟图，如图所示，
下半部分为示波器显示和设置界面。
互耦水槽系统实验模拟图
在图中，设定值是通过信号发生器（图中和Signal1SignalGenerator
）提供的，信号发生器可以提供五种类型的信号源，分别为阶跃信Generator2
号、正弦波信号、方波信号、脉冲信号和嗓声信号。
阀门的开关通过双击阀门
上的红色开关来实现，红色开关与管子垂直时，表示阀门关，反之，表示阀门
开。
双击控制器可以设置控制器的三个参数值（、、）。
另外，kkPIDkPIDdpi整个实验可以在开环和闭环两种状态下进行，如果要进行闭环控制，，需将反
馈线路上的开关指向非的那一端。
图中用处不同颜色菱形框标注的地方是60V
示波器所测的参数，下面分别介绍一下：
①为水槽的液位值;②为被控参数设1
定值;③为被控参数设定值与实际值的偏差;④为水槽的液位值;⑤为控制PID2
器输入;⑥为控制器输出。
PID
下半部分的左面为示波器的显示界面，整个界面被横轴和纵轴又分为四个
，纵轴为电压轴，单位为伏特，示波器可以显部分，横轴为时间轴，单位为S
所示为①－⑥，每个部分的运行个部分的运行状态，如图示互耦水槽系统36
状态在示波器上通过不同颜色的曲线显示。
在下半部分的右测还可以对示波器
的参数进行设置，如设置时间轴和电压轴的每一间隔所表示的时间和电压值，
这样可以将曲线进行相应的放大和缩小，另外在示波器设置界面上还有一些功
能按钮，如等等，使示波器处于运行状态，采集RUNRUN,STOP,RESET,DISPLAY
数据使示波器停止运行，这时前面采集到的数据会停在显示界面上，然;STOP
后可以通过按钮，将图形导出成文件，以便进行数据分析。
EXPORT
对耦合水槽建模：
2.dd
d（3-39d号号水槽的截面积其中分别是两个水c36.5
分别是两个水泵向水槽中注水的体积流速的液面高度c
分别是对应的比例常数，通过升高或者降低挡板可以控
两个水槽之间的流速，关号水泵，通号水泵号水槽中注水以控号2121
槽液面高度，因此
（3-40
初始化
5.618
1（3-41设置最大控制信号，ma7（3-42
采样时间秒1
分别使用上一节中的三个模糊控制系统对耦合水槽进行仿真控制，号水槽2
目标液面，一型模控制系统手动调节参数如下PIref=15c
80.070.75
模控制系统选取参数如下PIIT2
0.181.20.070.75
依据参考输入，建模程序如下
h1=0;
h2=0;
h=2;
A=36.52;
a1=5.6186;
a2=10;
q=[37.25677672
69.66088955
75
75
75
75
75
75
75
75
75
75
75
75
74.97903823
73.64093728
70.67536298
66.29923609
60.82554485
54.64709663
48.21781909
42.02959175
36.58216821
32.34352313
29.69841269
28.88575668
29.93409086
32.62119253
36.49739777
40.98847746
45.53065207
49.66539648
53.07066222
55.55247629
57.0241398
57.48593192
57.00837713
55.71847404
53.78749167
51.41903223
48.8363526
46.26828262
43.93350228
42.02352006
40.6855027
40.00711842．
40.00649753
40.63064664
41.76438138
43.24890457
44.90581916
46.56075347
48.06191265
49.29163498
50.17148078
50.66251846
50.76249748
50.501139
49.93428356
49.13729457
48.19794443
47.20897529
46.26058646
45.43321907
44.79114445
44.37746537
44.21114619
44.28654244
44.57558428
45.03233373
45.5992105
46.21391345
46.81604602
47.35266365
47.78229055
48.0772753
48.22458908
48.22529684
48.09297482
47.85134442
47.53137371
47.16808156
46.79727248
46.45242905
46.16198588
45.94719374
45.82074512
45.78626877
45.83871298
45.96553696．
46.14853649
46.3660615
46.59535563
46.81476172
47.0055865
47.1534859
47.24930399
47.28936198
47.27524227
47.21314639
]
fori=1:
1:
100
y（i）=h2;
q1=q（i）;
k1=h*f（q1,h1,h2）;
m1=h*g（h1,h2）;
k2=h*f（q1,h1+0.5*k1,h2+0.5*m1）;
m2=h*g（h1+0.5*k1,h2+0.5*m1）;
k3=h*f（q1,h1+0.5*k2,h2+0.5*m2）;
m3=h*g（h1+0.5*k2,h2+0.5*m2）;
k4=h*f（q1,h1+0.5*k3,h2+0.5*m3）;
m4=h*g（h1+0.5*k3,h2+0.5*m3）;
h1=h1+1/6*（k1+2*k2+2*k3+k4）;
h2=h2+1/6*（m1+2*m2+2*m3+m4）;
end
x=1:
1:
100;
plot（x,y）
holdon
plot（x,15）
模型验证：
3.
内的仿真结果如图所示，其中：
3.8200s
）表示一型模糊控制系统，PIT1PIFS1
）表示使用迭代降型的区间二型模糊控制系统，PIIT2PIFS12KM
）表示使用近似法降型的区间二型模糊控制系统，PIWu-Mendel3IT2PIFS2
三个控制系统对应的性能指标如表中所示，使用（时间乘绝对误ITAE3.1
差的积分）作为衡量控制系统的性能标准。
从图及表的性能指标中，可3.13.8
以看出，使用迭代降型法的模糊控制系统性能最优，然后为使用PIKMIT2
近似降型法的模糊控制系统，一型模糊控制系统的性能最PIWu-MendelIT2PI
差。
与一型模糊控制系统相比，使用迭代降型的模糊控制系统在PIIT2KMPI
所需时间）不增加的情况下，减小了上升至保持上升时间（从稳态值90%10%
控制系统在上升模糊超调量和调节时间，而使用近似降型的PIWu-MendelIT2
时间不明显增加时，超调量和调节时间只是略微下降。
模型建立好以后，应自行验证。
验证指标为使1tdt（t）y（t）?
y2fin?
?
Itm越小越好。
p为起始时间，t?
ttiniiniinifin
t为模型响应。
为实际响应，为终止时间，）y（y（t）tfinm即使模型响应曲线越接近于实际响应曲线越好。
20
15实实实实
实1实实实T1PIFS实实实实实实IT2PIFS1实0实实实实IT2PIFS2实5200100150050/s实0实
耦合水槽控制仿真控制系统仿真对比图3.1
不同控制系统的性能对比3.1
模糊系IT2
模糊系系统类Wu-T1KMMendel
上升时）313030
超调4.2686%2.4333%4.6289%
调节时）
81
9373
=2
6384.4
5927.96747.6
性能指（ITAE）
仿真结果如下：
感受：
虽然很多东西都不是自己想出来的，经过很多的资料搜
索阅读，我感受到了这门课的重要性，整个仿真界都用这门课知识，
整个没有解析解的地方都用它，极度重要，后悔没有好好学数学，
MATLAB的程序是借鉴同学的，学霸挨句讲了一遍，半懂吧！
学霸意
味深长地说，MATLAB即使不参加数模也很重要，默默决定抽时间去
学一下。
升级会员 