matlab模型求解.ppt

科学计算与MATLAB实验讲义;第三讲 ＭＡＴＬＡＢ模型求解 ——（常微分、偏微分）方程（组），概率统计;内容提要;1、 实验目的;2、实验内容;; x = fzero(FUN,x0) %x0可以是数，或区间 x = fzero(FUN,x0,options) [x,fval]= fzero(FUN,x0,options) [x,fval,exitflag] = fzero(FUN,x0,options);;;单因素方差分析 函数 anova1 格式 p = anova1(X) %X的各列为彼此独立的样本观察值，其元素个数相同，p为各列均值相等的概率值，若p值接近于0，则原假设受到怀疑，说明至少有一列均值与其余列均值有明显不同。 p = anova1(X,group) %X和group为向量且group要与X对应 p = anova1(X,group,displayopt) % displayopt=on/off表示显示与隐藏方差分析表图和盒图 [p,table] = anova1(…) % table为方差分析表 [p,table,stats] = anova1(…) % stats为分析结果的构造 说明 anova1函数产生两个图：标准的方差分析表图和盒图。;方差分析表中有6列： 第1列(source)显示：X中数据可变性的来源； 第2列(SS)显示：用于每一列的平方和； 第3列(df)显示：与每一种可变性来源有关的自由度； 第4列(MS)显示：是SS/df的比值； 第5列(F)显示：F统计量数值，它是MS的比率； 第6列显示：从F累积分布中得到的概率，当F增加时，p值减少。;2、 双因素方差分析 函数 anova2 格式 p = anova2(X,reps) p = anova2(X,reps,displayopt) [p,table] = anova2(…) [p,table,stats] = anova2(…) 说明 执行平衡的双因素试验的方差分析来比较X中两个或多个列（行）的均值，不同列的数据表示因素A的差异，不同行的数据表示另一因素B的差异。如果行列对有多于一个的观察点，则变量reps指出每一单元观察点的数目，每一单元包含reps行，如： reps=2 其余参数与单因素方差分析参数相似。;3、概率统计图;（3）绘制威布尔(Weibull)概率图形 函数 weibplot 格式 weibplot(X) %若X为向量，则显示威布尔(Weibull)概率图形，若X为矩阵，则显示每一列的威布尔概率图形。 h = weibplot(X) %返回绘图直线的柄 说明 绘制威布尔(Weibull)概率图形的目的是用图解法估计来自威布尔分布的数据X，如果X是威布尔分布数据，其图形是直线的，否则图形中可能产生弯曲。;（5）附加有正态密度曲线的直方图 函数 histfit 格式 histfit(data) %data为向量，返回直方图 和正态曲线。 histfit(data,nbins) % nbins指定bar的个数， 缺省时为data中数据个数的平方根。;3、 上机实践;对于等间隔节点;在tn节点上，微分方程可以写为;具体求解过程为：;简单欧拉方法程序 function [outx,outy]=MyEuler(fun,x0,xt,y0,PointNum) %MyEuler 用前向差分的欧拉方法解微分方程 %fun 表示f（x，y） %x0,xt表示自变量的初值和终值 %y0表示函数在x0处的值,其可以为向量形式 %PointNum表示自变量在[x0,xt]上取的点数 if nargin5 | PointNum=0 %如果函数仅输入4个参数值，则PointNum默认值为100 PointNum=100; end if nargin4 %y0默认值为0 y0=0; end h=(xt-x0)/PointNum;%计算步长h x=x0+[0:PointNum]*h;%自变量数组 y(1,:) = y0(:);%将输入存为行向量，输入为列向量形式 for k = 1:PointNum f=feval(fun,x(k),y(k,:));%计算f(x,y)在每个迭代点的值 f=f(:); y(k + 1,:) =y(k,:) +h*f; %对于所取的点x迭代计算y值 end outy=y; outx=x; %plot(x,y)%画出方程解的函数图;2、不动点迭代法求解非线性方程 迭代法是一种逐次逼近的方法，用某个固定公式反复校正根的近似值，使之逐步精确化，最后得到满足精度要求的结果。;;将连续函数方程f