计算机仿真总述应用举例.pptVIP

第六十一页，共一百三十页，2022年，8月28日 例8：方程组求解 第六十二页，共一百三十页，2022年，8月28日 初始值x0=[-5 –5] 第六十三页，共一百三十页，2022年，8月28日 %程序myfun88.m function f=myfun88(x) f=[2*x(1)-x(2)-exp(x(1)); -x(1)+2*x(2)-exp(-x(2))]; %程序fxxfc.m x0=[-5;-5]; [x，F]=fsolve(@myfun88,x0) 第六十四页，共一百三十页，2022年，8月28日 例9： 数值积分 x=0:0.1:10 plot(x , x.^2) 第六十五页，共一百三十页，2022年，8月28日 程序：只需一条语句area=quad(x.^2,4,7) 结果： area = 93 验算： 第六十六页，共一百三十页，2022年，8月28日 例10：微分方程的数值解法 例：解经典范得波（Van der Pol)微分方程 令：ω=2 第六十七页，共一百三十页，2022年，8月28日 第六十八页，共一百三十页，2022年，8月28日 [t,y]=ode23(myfun5,[0 30],[1;0],[],3); y0=[1;0] y1=y(:,1); y2=y(:,2); plot(t,y1,:b,t,y2,-r) legend(‘速度’,‘位移’) function yp=myfun5(t,y); yp=[y(2) 2*(1-y(1)^2)*y(2)-y(1)]; 第六十九页，共一百三十页，2022年，8月28日 偏微分方程数值解法工具箱可应用于弹性力学、静电学、电磁学、热传导和发散等许多领域。 第七十页，共一百三十页，2022年，8月28日 例11：优化设计 第七十一页，共一百三十页，2022年，8月28日 问题1 某车间有甲乙两台机床，可用于加工三种工件，这两台机床的可用台时数分别为700和800，三种工件的数量分别为300、500和400，已知机床加工情况如下表： 第七十二页，共一百三十页，2022年，8月28日 问如何分配机床的加工任务，才能既满足加工工件的要求，又使总加工费用最低？ 解：设在甲机床上加工工件1、2和3的数量分别是x1、 x2和x3，在乙机床上加工工件1、2和3的数量分别是x4、 x5和x6，根据条件可建立数学模型： 第七十三页，共一百三十页，2022年，8月28日 第七十四页，共一百三十页，2022年，8月28日 程序： f=[13;9;10;11;12;8]; A=[0.4 1.1 1 0 0 0 0 0 0 0.5 1.2 1.3]; b=[700;800]; aeq=[1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1]; beq=[300 500 400]; lb=zeros(6,1); [x,fval,exitflag,output,lambda]=linprog(f,A,b,aeq,beq,lb) 第七十五页，共一百三十页，2022年，8月28日 运行结果： x = 0.0000 500.0000 0.0000 300.0000 0.0000 400.0000 fval = 1.1000e+004 exitflag = 1 第七十六页，共一百三十页，2022年，8月28日 即：在甲机床上加工500个工件2，在乙机床上加工300个工件1，加工400个工件3，可在满足条件情况下使总加工费用最小，为11000元。且收敛正常。 第七十七页，共一百三十页，2022年，8月28日 问题2 求： 在区间[2 5]间的最小值 在区间[0 3]间的最大值 第七十八页，共一百三十页，2022年，8月28日 程序： %y5.m x=0:0.1:8; y=2*exp(-x).*sin(x); plot(x,y) fn=2*exp(-x)*sin(x); [xmin,fmin]=fminbnd(fn,2,5) fm=-2*exp(-x)*sin(x); [xmax,fmax]=fminbnd(fm,0,3) 第七十九页，共一百三十页，2022年，8月28日 结果： xmin = 3.9270 fmin = -0.0279 xmax = 0.7854 fmax = -0.6448 第八十页，共一百三十页，2022年，8月28日 例12：二维图形处理 Demo/Graphics/2-D plots 第八十一页，共一百三十页，2022年，8月28日 例13： 三维图形