上计51_ 2013年最新Matlab语言及其在工程领域中的应用课件.pptVIP

MATLAB上机练习5;定积分数值求解;基本思想：将被积函数在积分区间上离散节点(足够多)处的函数值与其所在小区间长度的乘积之和作为定积分的近似值或者构造一个函数P(x)逼近原函数，使得 数值积分方法主要有Newton-Cotes数值积分法，Gauss积分法和Romberg积分法。;n=1，梯形公式 n=2，复化simpson公式 n=3，Cotes公式;n＝1 梯形公式(trapezoid rule);Matlab中quad用递归自适应的外插Simpon公式 q=quad(@fun, a, b) q=quad(@fun,a , b, tol) q=quad(@fun, a, b,t ol, trace, p1, p2,…) q=quad(@fun, a, b, [ ], [ ], p1, p2,…) ％ p1,p2：直接传递给函数fun的已知参数 ;syms x Isym=vpa(int(exp(-x^2),x,0,1)) Isym = 0.74682413281242702539946743613185 format long d=0.001;x=0:d:1; Itrapz=d*trapz(exp(-x.*x)) Itrapz = 0.74682407149919 fx=exp(-x.^2); Ic=quad(fx,0,1,1e-8) Ic = 0.74682413285445 ;quad8:：自适应的Cotes公式 q=quad8(@fun,a,b) q=quad8(@fun,a,b,tol);;常微分方程（组）数值解;;;求解方法：Runge-Kutta法 ;辛普森公式展开;二阶、三阶Runge-Kutta函数 ode23 [t,y]=ode23(@fun, tspan,y0) [t,y]=ode23(@fun, tspan,y0,options,p1,p2);例1：求解一阶常微分方程：;例2：求解一阶常微分方程组：;例3：求解含积分的一阶常微分方程组：;function dd3 clear all;clc [t,y]=ode45(@fun1,[0,2],[0,0]) plot(t,y(:,1),t,y(:,2)) function dydt=fun1(t,y) a=1; if t1 dydt(1)=1.*sin(y(1))-3.*t.^2-10; elseif t=0.5 dydt(1)=quad(@fun2,0,1,[],[]).*t; else dydt(1)=y(1).^1-t-y(2); end dydt(2)=y(2).^2-2.*t; dydt=dydt(:); function q=fun2(x) q=sqrt(cos(x))-exp(x);;Exercise 1: 练习并熟练掌握用m文件编辑器(editor)编写脚本程序(script file);;例1：求解二阶常微分方程：;例2：求解一阶常微分方程组：;Exercise 2: 练习并熟练掌握用m文件编辑器(editor)编写无输入输出的函数文件(function file);Exercise 3: 练习并熟练掌握用m文件编辑器(editor)编写有输入输出的函数文件(function file); 将下列方程组编写为一个有输入输出(待调用)的函数文件; 将下列方程组编写为一个有输入输出(待调用)的函数文件;;IF;IF