




付费下载
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
微分方程数值解大作业双曲型方程求解如下双曲型方程其中关于时间和空间的二阶偏导数均采用二阶中心差商逼近.1.,且初边条件如下:存在精确解为:2.,且初边条件如下:存在精确解为:第一题数值解的源程序为:functionu=hyperbolic(minx,maxx,mint,maxt,M,N,a)h=(maxx-minx)/M;%横轴步长t=(maxt-mint)/N;%纵轴步长r=t/h;u=zeros(M+1,N+1);%先定义u为M+1*N+1维的矩阵u(1,:)=0;%左边界为0u(M,:)=0;%右边界为0u(:,1)=bottoma(h,M)';%下边界由已知函数bottoma计算fori=1:M-1%对下边界上的导数g1的差分格式进行定义g1(i)=g(i*h);endforj=2:M%由导数的差分格式求第2行的u值u(j,2)=0.5*(a^2*r^2*(u(j+1,1)+u(j-1,1))+2*(1-a^2*r^2)*u(j,1)+t*g1(j-1));endfork=2:N%从第三行开始,一行一行地求u值forj=2:Mu(j,k+1)=a^2*r^2*(u(j+1,k)+u(j-1,k))+2*(1-a^2*r^2)*u(j,k)-u(j,k-1);endend下边界的导数的定义functiong=g(x)g=0下边界函数的定义functiona=bottoma(h,m)pi=3.1415926;a(1)=0;a(m+1)=0;fori=2:ma(i)=sin(pi*(i-1)*h);end精确解的源程序:functionu=hy(minx,maxx,mint,maxt,M,N)pi=3.1415926;x=minx+(maxx-minx)*(0:M)/M;t=mint+(maxt-mint)*(0:N)/N;forj=1:N+1fori=1:M+1u(i,j)=cos(pi*t(j))*sin(pi*x(i));endend稳定时的误差sum(sum(hyperbolic(0,1,0,2,4,8,1)-hy(0,1,0,2,4,8)))ans=2.587544328334521e-007误差很小第二题数值解代码functionu=hyperbolic(minx,maxx,mint,maxt,M,N,a)h=(maxx-minx)/M;%横轴步长t=(maxt-mint)/N;%纵轴步长r=t/h;u=zeros(M+1,N+1);%先定义u为M+1*N+1维的矩阵u(1,:)=0;%左边界为0u(M,:)=0;%右边界为0u(:,1)=bottoma(h,M)';%下边界由已知函数bottoma计算fori=1:M-1g1(i)=g(i*h);%对下边界上的导数g1的差分格式进行定义endforj=2:M%由导数的差分格式求第2行的u值u(j,2)=0.5*(a^2*r^2*(u(j+1,1)+u(j-1,1))+2*(1-a^2*r^2)*u(j,1)+t*g1(j-1));endfork=2:N%从第三行开始,一行一行地求u值forj=2:Mu(j,k+1)=a^2*r^2*(u(j+1,k)+u(j-1,k))+2*(1-a^2*r^2)*u(j,k)-u(j,k-1);endend下边界函数的定义functiona=bottoma(h,m)a=0;下边界的导数的定义functiong=g(x)pi=3.141592653589793;g=sin(4*pi*x);稳定的精确解functionu=hy(minx,maxx,mint,maxt,M,N)pi=3.141592653589793;x=minx+(maxx-minx)*(0:M)/M;t=mint+(maxt-mint)*(0:N)/N;forj=1:N+1fori=1:M+1u(i,j)=sin(t(j))*sin(4*pi*x(i));endend稳定性条件和Courant条件分析。我们采用分离变量的方法对格式进行稳定性分析,可得传播矩阵:G而这个矩阵的特征方程为:λ他的根λ1,2按模不大于1的充要条件为2-ar≤1.亦即,当上式成立时,vonNeuma
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 二零二五年度崇明区危化品运输车辆租赁合同范本
- 2025年度房产抵押贷款买卖合同书
- 2025春国家开放大学教育教学改革研究项目合同
- 二零二五年度铜管材购销合同模板
- 2025年防雷接地施工及维护一体化服务合同
- 2025版风机租赁与销售一体化服务合同模板
- 2025版个人二手房买卖协议含房屋质量保证期限及维修责任
- 2025版葡萄酒线上线下联合促销销售合同
- 2025版柴油市场调查与分析合同模板
- 2025年度海上货物运输合同-集装箱运输管理及安全协议
- 股权代持协议范本:股权代持与股权质押
- 青春期生理讲课件
- 2025-2030年中国财税服务行业市场深度调研及发展前景与投资研究报告
- 贵州省2024年高考真题政治试卷(含答案)
- 2024年咸阳市社区工作者计划招聘真题
- 2025年四川酒业茶业投资集团有限公司及下属子公司招聘笔试参考题库含答案解析
- 新能源汽车热管理系统的能量优化与梯级利用策略探讨
- 工程控制面试题及答案
- (高清版)DB13 5808-2023 餐饮业大气污染物排放标准
- 【泡泡玛特营销策略研究的文献综述】3100字
- 2025农业银行笔试题目及答案
评论
0/150
提交评论