matlab上机实验第4次课.pptVIP

考试安排 考试时间：40分钟。每人2题。 在E盘上建立文件夹，以“班级+学号+姓名”命名。 开卷考试，只能带书本，不可拷贝. * * * * 常微分方程的求解 p158 统计作图 p180 多项式拟合 p123 随机试验 p169 设微分方程初值问题： 常微分方程的求解 分为：符号解法和数值解法 符号解法： 命令形式1： dsolve（’eqution’,’var’） 命令形式2： dsolve（’eqution’ , ’cond1,cond2,…’ , ’var’） P158 初始条件 符号解法 dsolve(Dy=y^2,x) 自变量 大写 ?dsolve(Dy=y^2,x) ans = -1/(x-C1) 求解 ?dsolve(x*D2y-3*Dy=x^2,y(1)=0,y(5)=0,x) 大写 对应求导阶数 condition ans =-1/3*x^3+125/468+31/468*x^4 求解常微分方程 [t,y]=ode23(‘fun’,tspan,y0) tspan=[t0 tf]表示积分的起始值和终止值； fun是定义函数的文件名； [t,y]=ode45(‘fun’,tspan,y0) y0是初始状态列向量。 常 微 分 方 程 的 数 值 解 法 高阶微分方程需转化一阶微分方程组 注：只适用于一阶微分方程 例7-45：用数值积分的方法求解微分方程： 设初始时间 ；终止时间 初始条件 分析：求解 令： （化为一阶微分方程）即原微分方程化为： 写成矩阵形式为 （化为一阶微分方程） xdot= 放入函数exf.m中 u function xdot=exf(t,x) u=1-(t.^2)/(2*pi); xdot=[0,1;-1,0]*x+[0 1]*u; 编写函数文件 exf.m clf; t0=0; tf=3*pi; x0t=[0;0]; 定义另外一个函数为主函数 初始和终止时间 初始条件 [t,x]=ode23(exf,[t0,tf],x0t) exf为已定义的子函数 该函数用来记录一阶微分方程 xdot= y=x(:,1); %[t,x]中求出的x是按列排列， 故用ode23求出x后 只要第 一列即为y y2=-1/2*(-2*pi-2+t.^2)/pi-(pi+1)/pi*cos(t); clf, plot(t,y,-, t,y2,o) 解析解为： dsolve(D2y+y=1-t^2 /(2*pi),y(0)=0,Dy(0)=0,t) ans = -1/2*(-2*pi-2+t^2)/pi-(pi+1)/pi*cos(t) P=polyfit（x，y，n） 多项式数据拟合 功能： 将给定向量 x，y对应的（x[i]，y[i]）作为数据点，拟合成n次多项式； 向量x，y具有相同的维数； n为正整数，n值越大则拟和的精度越好； p为多项式的系数向量。 poly2str(p, ‘x’) 将多项式表示成习惯的 形式 p是多项式系数，字符‘x’为自变量 polyval（p，X） 按数组规则计算X处多 项式的值 功能： 计算多项式p的变量在点阵X处的值 ； x可以为向量或矩阵，计算结果是与x同维的向量或矩阵。 多项式的拟合 P123 例?6-17：现有一组实验数据：x的取值是从1到2之间的数，间隔为0.1，y的取值为2.1,3.2, 2.1,2.5,3.2,3.5,3.4,4.1,4.7,5.0,4.8。要求分别用二次、七次拟合曲线来拟合这组数据，观察这2组拟合曲线那个效果更好？ 解：建立M文件 x=1:.1:2; y=[2.1,3.2,2.1,2.5,3.2,3.5,3.4,4.1,4.7,5.0,4.8]; p2=polyfit(x,y,2), %多项式拟合,阶数是2 ,p2为拟合多项式的系数 p7=polyfit(x,y,7); disp(二阶拟合函数),f2=poly2str(p2,x) disp(七阶拟合函数),f7=poly2str(p7,x) x1=1:.01:2; y2=polyval(p2,x1); %多项式p2在x1处的值 y7=polyval(p7,x1); plot(x,y,rp,x1,y2,--, x1,y7); legend(拟合点,二次拟合,七次拟