多元微积分及应用.pptVIP

二、实验目的 1.熟习Mathematica编程命令． 2.学会用计算机模拟真实系统模型，并通过计算了解一个系统随向量变化的行为或特性，一便对系统进行评价或预测． 三、预备知识 1.学习有关Mathematica软件流程控制命令．如： If[condition,t]；If[condition,t,t] 功能：若条件condition成立，则执行t，否则返回NULL． Which[test1,value,test2,value,…] 功能：若test1为真,执行value1；若test2为真,则执行value2；… Do[expr,{i,imax}]；Do[expr,{I,imin,imax,di}] 功能：从i到imax，步长为1，计算含有i的表达式expr ． While[test,body]；For[start,test,incr,body] 功能：求test的值，若为True，则计算循环体body，直到得到test的值不是True． 2.学习有关Mathematica软件输入与输出命令．如： Print[expr1,expr2,…] 功能：在同一行中打印表达式expri的值，中间无空格． 四、实验的简单操作程序 由于水的注入与流出，使池中的水的体积，含盐量与含盐率均随时间边变化．初始时刻的含盐率确定以后，每分钟流入含盐率为0.5kg/立方米的水6立方米，并流出混合均匀后的盐水4立方米，直至水的含盐率达到0.2kg/立方米．将此过程转化为流程图，如图1所示： 四、实验的简单操作程序 图1: 四、实验的简单操作程序 根据所作框图可以用Mathematica命令编程．Volume用来表示水的体积，total表示含盐量，rate表示含盐率，而用tt表示时间，其原程序为： time = 0; tt = 0; volume = 2000; total = 2; rate = N[total/volume, 7]; Print[time, volume, salt, saltrate]; While[rate 0.2, time = time + 1; volume = volume + 2; total = total + 3 - 4rate; rate = total/volume; If[time = 10, Print[tt + time, , volume, , total, , rate]; tt = tt + time; time = 0]]; Print[tt + time, , volume, , total, , rate]; 四、实验的简单操作程序 time volume salt saltrate 10 2020 31.6937 0.0156899 20 2040 60.8106 0.0298098 30 2060 89.3761 0.0433843 40 2080 117.4 0.0564421 50 2100 144.914 0.0690067 60 2120 171.935 0.0811012 70 2140 198.479 0.0927474 80 2160 224.566 0.103966 90 2180 250.212 0.114776 100 2200 275.432 0.125196 110 2220 300.242 0.135244 120 2240 324.657 0.144936 130 2260 348.69 0.154287 140 2280 372.355 0.163313 150 2300 395.664 0.172028 160 2320 418.631 0.180444 170 2340 441.266 0.188575 180 2360 463.58 0.196432 185 2370 474.621 0.200262 程序运行结果： 五、思考与提示 1.该计算机模拟是连续模拟还是离散模拟？ 2.用计算机模拟实际问题其实质是什么？ 六、练习内容 如图所示，正方形ABCD的四个顶点各有一人．在某一时刻（设为t0=0），四人同时出发以匀速v走向顺时针方向的下一个人．如果他们始终保持对准目标，则最终将安螺旋状曲线汇合于中