-A)*x0+bРnorm(x-x1)>eps?Р是Рn>=200?Р例题流程图Р输入系数矩阵AР ↓Р输入初始向量x0Р及常数向量bР[x,n]=richason(A,b,x0)Р ↓Р输出计算结果Р↓Р输出迭代次数Р↓РA*x验证结果Р↓Р解:?Р(1)算法说明Р分析已给方程可知,为拉普拉斯方程,在MATLAB工具箱PDETOOL中可看成椭圆型方程,转化为标准形式如下:Р因此,对应的c=-1,a=0,f=0,然后根据给出的边界约束条件,在微分方程工具箱中选择所需要的条件,РDirichlet条件Р РNeumann条件Р Р其中n是上的单位外法矢量,g,q,h和r是定义在上的函数｡Р(题目中Γ1与Γ2分别代表x+y=2与x-y=2这两条边界线)Р(2)操作流程Р设置坐标限Р选择Options栏中Axes Limits选项,输入坐标范围Р绘制区域图Р点击绘制多边形键画出要求的区域图Р设置边界条件Р选择Boundary中的Boundary Mode,设置为边界模式;双击各条边界线,由方程组中已知边界条件设定Р设置方程参数Р点击,将已知方程对照标准偏微分方程形式,知c=-1,a=0,f=0｡Р Р剖分网格Р按顺序点击两按钮,细分网格｡Р绘制温度分布图Р点击绘制三维示意图:Р(3)简易流程图Р 开始Р绘制要求区域图Р设置边界条件Р设置方程参数Р剖分网格Р绘制温度分布示意图Р结束Р实验1 用GUI方式解下列PDEР解:Р(1)算法说明Р同上题,由已给方程可知,为拉普拉斯方程,在PDETOOL中可看成椭圆型方程,转化为标准形式如下:Р因此,对应的c=1,a=0,f=0,然后根据给出的边界约束条件,在微分方程工具箱中选择所需要的条件,РDirichlet条件Р РNeumann条件Р Р其中n是上的单位外法矢量,g,q,h和r是定义在上的函数｡Р(2)操作流程Р设置坐标限Р绘制区域图
全文预览
