古典显示格式求解抛物型偏微分方程(Matlab实现)
古典显示格式求解抛物型偏微分方程(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实现)相关推荐
- 解抛物型微分方程c语言程序,偏微分方程的matlab代码_偏微分方程matlab程序,matlab偏微分代码...
matlab偏微分方程含源代码,大家可以下载运行,里面的都能用 Pie Eais ae h t lat BLep 话k气飞"目基回 古班比比,一生热件身的些懂 1身定 回量1 古典显式格式稳 ...
- 差分法数值求解热传导偏微分方程
代码和推导文档可以在这里下载:差分法求解热传导偏微分方程Matlab代码 1.热传导偏微分方程 热传导方程,是传热数中,经典的数学方程.这里说的热传导方程,不仅仅是指传热学中的方程,而是指与热传导偏微 ...
- matlab抛物偏微分方程,抛物型方程差分求解 跪求MATLAB解抛物型偏微分方程的程序...
为什么抛物线方程与圆方程联立不能使用韦达定理 很容易了解到,抛物线和圆的交点均在X轴上方\"其实这时你应该注意到一点就是,这两个交点的纵坐标是相等的,所以其实对应的是一个y值,也就是你列的一 ...
- matlab 离散偏微分,有限差分法求解偏微分方程MATLAB.doc
有限差分法求解偏微分方程MATLAB.doc 南京理工大学 课程考核论文 课程名称: 高等数值分析 论文题目: 有限差分法求解偏微分方程 姓 名: 罗 晨 学 号: 115104000545 成 绩: ...
- matlab pdepe函数边界,科学网-使用MATLAB中pdepe函数求解一维偏微分方程-邓浩鑫的博文...
由于自己科研水平较低,记录的各种体会更多的是给自己做个小结,错误之处,欢迎大家指正. 使用MATLAB求解偏微分方程或者方程组,大致有三类方法.第一种是使用MATLAB中的PDE Toolbox,PD ...
- 傅里叶谱方法-傅里叶谱方法求解基本偏微分方程(一维波动方程、 二维波动方程、一维非线性薛定谔方程)及其Matlab程序实现
3.2 傅里叶谱方法求解基本偏微分方程 (组) 3.2.1 一维波动方程 对于一根两端固定.没有受到任何外力的弦, 若只研究其中的一段, 在不太长的时间 里, 固定端来不及对这段弦产生影响, 则可以认 ...
- 偏微分方程matlab课件,MATLAB偏微分方程
4.8.2 偏微分方程 在自然科学的很多领域内,都会遇到微分方程初值问题,特别是偏微分方程,它的定解问题是描述自然界及科学现象的最重要的工具.可以说,几乎自然界和各种现象都可以通过微分方程(特别是偏 ...
- 用matlab求解jacobi,用jacobi迭代法求解线性方程,求助matlab大师
问题描述: 用jacobi迭代法求解线性方程,求助matlab大师 用jacobi迭代法求解线性方程 9x1-x2+x3=10 -x1+10x2-2x3=7 -2x1+x2+10x3=6,设迭代初值为 ...
- 【单目标优化求解】基于matlab增强型黑猩猩优化器算法求解单目标优化问题【含Matlab源码 2013期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[单目标优化求解]基于matlab增强型黑猩猩优化器算法求解单目标优化问题[含Matlab源码 2013期] 点击上面蓝色字体,直接付费下 ...
最新文章
- hybbs接口php,HYBBS
- 操作系统(三)内存管理
- controller,sevices层,java初步了解
- Go语言程序的数组初始化
- CentOS+NodeJS环境
- 计算机设备运输规范,《电子计算机机房设计规范》GB50174-93
- 基于51单片机的GPS公交自动报站系统
- python爬虫图书信息并存入数据库,以及安装工具库
- 机器学习——特征工程——数据的标准化(Z-Score,Maxmin,MaxAbs,RobustScaler,Normalizer)
- Android Studio 个性化主题配置。
- jQuery中的动画 -- 案例
- Python爬取高颜值美女(爬虫+人脸检测+颜值检测)
- c语言人民币小写转换成大写,C语言实现人民币小写转大写
- 如何使用 CSS flex box 和 Javascript 设计棋盘
- 计算机组成原理实验主存储器,计算机组成原理实验课件--存储器.ppt
- 7月5号cf服务器维护,cf2015年5月7日维护到几点 cf5月7日更新公告
- 硬盘格式化后数据如何恢复
- 我的软设、系分考试经验(书)
- 【分享】国内著名互联网企业PHP笔试题
- C++ 模板与泛型编程简述