计算机模拟第二课.pptVIP

1.3 计算机模拟的步骤 1.明确模拟系统 2.建立数学模型 3.模型转换（数学模型转化为模拟模型） 4.设计模拟程序 5.模型实验 6.模型检验（正确性、稳定性） 7.实验结果的评价与分析 1.4 计算机模拟应用举例 1. 追击敌机问题 设敌机B在80公里高处以VB速度沿直线匀速飞行，我方飞机 F以VF速度进行追击，假设两机始终保持在同一个平面内，我方飞机机头始终对准敌机，并设VB=20km/min， VF=40km/min，设两机相距在10km之内我机即可进行攻击摧毁敌机，若我机在10分钟内没有捕捉到敌机，则认为敌机已逃脱，追击失败。 计算机模拟图示 function sim_pushbutton_Callback(hObject, eventdata, handles) % hObject handle to sim_pushbutton (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) hold on % 图形迭加 axis([0,100,0,100]); % 设置坐标轴 k=0; A=[0,0]; %我机初始位置 B=[0,80]; %敌机初始位置 v=1; dt=1; %离散时间改变量 d=80; %相距距离 tol=str2num(get(handles.tol_edit,String)); while dtol plot(A(1),A(2),‘r.’); %画我机位置 plot(B(1),B(2),‘b+’); %画敌机位置 pause(0.1); k=k+1; B=B+[v*dt,0]; %敌机移动位置 e=B-A; %指向敌机向量 d=norm(e); e0=e/d; %取向量方向（单位化） A=A+2.0*v*dt*e0; %我机追击位置 str=strcat(k=,num2str(k), A(,num2str(A),) B(,num2str(B),) d=, num2str(d)); set(handles.out_edit,String,strvcat(get(handles.out_edit,String),str)); end set(handles.pos_text,String,strcat((,num2str(B(1)),, 80))) 2. 计算中子穿过铅墙的比例问题 原子反应堆由铅墙护卫，设铅墙厚5cm。中子从反应堆垂直进入铅墙，它直线行走每1cm，就会遇到一个铅原子而发生碰撞。碰撞会随机地改变中子的运动方向。几次碰撞以后中子可能因方向改变而返回反应堆。没有返回的中子，如果一直在铅墙里碰撞，则撞过10次以后就没有能量继续运动了。自然有些中子在10次碰撞内，穿过铅墙而越出去，造成危害。 求中子越出的百分比。 第二章 Matlab 语言的基本使用方法  2.1 基本运算 1、数学式的运算就像用计算器进行一样简单方便。 1+2+3+4+5 ans= 15 如果算式是 x=1+2+3+4+5 x=1+2+3+4+5 x= 15 这样我们就可以使用变量来进行更复杂的问题求解。 x=15 y=10 z=7 total=x+y*2+z*5 average=total/3 问题：ans代表什么变量？为何有的有ans而有的没有呢？ 2、一行中可以同时有几个语句，它们只要用逗号或分号隔开就可以了（应注意：它们的功能各不相同）。例如： a=2;b=3; a=2,b=3,c=4 3. 表2-1 matlab基本算术运算符 4.混合运算也涉及到运算次序问题 从左到右、幂次方优先级最高、乘除次之，最后是加减。如果有括号，则括号优先执行。 例2-1 有一个半径r=3的圆，计算其面积。 r=3; area=pi*r^2 Matlab可以将计算结果以