古典显示格式求解抛物型偏微分方程(Matlab实现)

clear
clc
at=[0,1];
ax=[0,1];
t1=0.01/6;
x1=0.1;
th=at(1):t1:at(2);
xh=ax(1):x1:ax(2);
n1=length(xh);
n2=length(th);
u=zeros(n2,n1);
syms x t;
bot=sin(pi*x);
boxr=0;
boxl=0;
%下面求边界值
bt=zeros(1,n1);
for i=1:n1bt(i)=subs(bot,x,xh(i));u(1,i)=bt(i);
end
bxr=zeros(1,n2);
bxl=zeros(1,n2);
for i=1:n2bxr(i)=subs(boxr,t,th(i));bxl(i)=subs(boxl,t,th(i));u(i,1)=bxr(i);u(i,n1)=bxl(i);
end
p=t1/(x1)^2;
for i=2:n2for j=2:n1-1u(i,j)=p*u(i-1,j-1)+(1-2*p)*u(i-1,j)+p*u(i-1,j+1);end
end
uu=zeros(n2,n1);
for i=1:n2for j=1:n1uu(i,j)=exp(-pi^2*th(i))*sin(pi*xh(j));end
end
mesh(xh,th,u);%数值解图
%mesh(xh,th,uu);%真实值图
u2=uu-u;
%mesh(xh,th,u2);%误差图;

古典显示格式求解抛物型偏微分方程(Matlab实现)相关推荐

  1. 解抛物型微分方程c语言程序,偏微分方程的matlab代码_偏微分方程matlab程序,matlab偏微分代码...

    matlab偏微分方程含源代码,大家可以下载运行,里面的都能用 Pie Eais ae h t lat BLep 话k气飞"目基回 古班比比,一生热件身的些懂 1身定 回量1 古典显式格式稳 ...

  2. 差分法数值求解热传导偏微分方程

    代码和推导文档可以在这里下载:差分法求解热传导偏微分方程Matlab代码 1.热传导偏微分方程 热传导方程,是传热数中,经典的数学方程.这里说的热传导方程,不仅仅是指传热学中的方程,而是指与热传导偏微 ...

  3. matlab抛物偏微分方程,抛物型方程差分求解 跪求MATLAB解抛物型偏微分方程的程序...

    为什么抛物线方程与圆方程联立不能使用韦达定理 很容易了解到,抛物线和圆的交点均在X轴上方\"其实这时你应该注意到一点就是,这两个交点的纵坐标是相等的,所以其实对应的是一个y值,也就是你列的一 ...

  4. matlab 离散偏微分,有限差分法求解偏微分方程MATLAB.doc

    有限差分法求解偏微分方程MATLAB.doc 南京理工大学 课程考核论文 课程名称: 高等数值分析 论文题目: 有限差分法求解偏微分方程 姓 名: 罗 晨 学 号: 115104000545 成 绩: ...

  5. matlab pdepe函数边界,科学网-使用MATLAB中pdepe函数求解一维偏微分方程-邓浩鑫的博文...

    由于自己科研水平较低,记录的各种体会更多的是给自己做个小结,错误之处,欢迎大家指正. 使用MATLAB求解偏微分方程或者方程组,大致有三类方法.第一种是使用MATLAB中的PDE Toolbox,PD ...

  6. 傅里叶谱方法-傅里叶谱方法求解基本偏微分方程(一维波动方程、 二维波动方程、一维非线性薛定谔方程)及其Matlab程序实现

    3.2 傅里叶谱方法求解基本偏微分方程 (组) 3.2.1 一维波动方程 对于一根两端固定.没有受到任何外力的弦, 若只研究其中的一段, 在不太长的时间 里, 固定端来不及对这段弦产生影响, 则可以认 ...

  7. 偏微分方程matlab课件,MATLAB偏微分方程

    4.8.2  偏微分方程 在自然科学的很多领域内,都会遇到微分方程初值问题,特别是偏微分方程,它的定解问题是描述自然界及科学现象的最重要的工具.可以说,几乎自然界和各种现象都可以通过微分方程(特别是偏 ...

  8. 用matlab求解jacobi,用jacobi迭代法求解线性方程,求助matlab大师

    问题描述: 用jacobi迭代法求解线性方程,求助matlab大师 用jacobi迭代法求解线性方程 9x1-x2+x3=10 -x1+10x2-2x3=7 -2x1+x2+10x3=6,设迭代初值为 ...

  9. 【单目标优化求解】基于matlab增强型黑猩猩优化器算法求解单目标优化问题【含Matlab源码 2013期】

    ⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[单目标优化求解]基于matlab增强型黑猩猩优化器算法求解单目标优化问题[含Matlab源码 2013期] 点击上面蓝色字体,直接付费下 ...

最新文章

  1. hybbs接口php,HYBBS
  2. 操作系统(三)内存管理
  3. controller,sevices层,java初步了解
  4. Go语言程序的数组初始化
  5. CentOS+NodeJS环境
  6. 计算机设备运输规范,《电子计算机机房设计规范》GB50174-93
  7. 基于51单片机的GPS公交自动报站系统
  8. python爬虫图书信息并存入数据库,以及安装工具库
  9. 机器学习——特征工程——数据的标准化(Z-Score,Maxmin,MaxAbs,RobustScaler,Normalizer)
  10. Android Studio 个性化主题配置。
  11. jQuery中的动画 -- 案例
  12. Python爬取高颜值美女(爬虫+人脸检测+颜值检测)
  13. c语言人民币小写转换成大写,C语言实现人民币小写转大写
  14. 如何使用 CSS flex box 和 Javascript 设计棋盘
  15. 计算机组成原理实验主存储器,计算机组成原理实验课件--存储器.ppt
  16. 7月5号cf服务器维护,cf2015年5月7日维护到几点 cf5月7日更新公告
  17. 硬盘格式化后数据如何恢复
  18. 我的软设、系分考试经验(书)
  19. 【分享】国内著名互联网企业PHP笔试题
  20. C++ 模板与泛型编程简述

热门文章

  1. TIPS001:Mendeley下载全文及其命令逻辑
  2. 计算机视觉常见数据集 转载
  3. html语音输入功能讯飞,讯飞输入法电脑版语音输入功能怎么使用?
  4. linux在终端打开dbi,Linux下安装DBI和DBD
  5. 科学家发现阻碍生发物质 治秃顶药5年内或可问世
  6. Mac上修复Gitee错误Oauth: Access token is expired
  7. TMS320C28x的中断
  8. FAF世链区块链大会|圆桌论坛:NFT虚拟世界的机遇和挑战
  9. Java Swing面板布局之箱式布局BoxLayout
  10. 观察:若干现有的RSS天气预报服务