Matlab二维热传导方程求解
摘要
本文利用有限差分法来求二维热传导方程的数值解,通过Matlab编程求解并作图,进而与解析 解做出的图进行比较,画出误差图。
引言
对于一维热传导方程,通过差分法将方程离散为方程组,利用追赶法求解系数矩阵为三对角矩 阵的方程组来得到数值解。本文采用和一维热传导方程求解类似的方法,通过对x方向和y方向采用 交替隐格式的方法来求解,一个时间层涉及到解两次三对角矩阵,本方法简单易行,而且具有较好 的精度。
求解二维热传导方程的基本步骤 2.1 将定解区域离散 用有限差分方法求解偏微分方程问题必须把连续问题进行离散化。为此首先要对求解区域给出 网格划分,由于求解的问题不同,因此求解区域也不尽相同。下面用例子来说明不同区域的剖分离 散。如果求解区域是 :
把x区域分割成Mx个小区间,把y区域分割成My个小区间,把时间t分割成T个时间层。 即
Matlab二维热传导方程求解相关推荐
- 隐式差分格式matlab,二维热传导方程隐式差分格式系数矩阵特征值的求取
二维热传导方程隐式差分格式系数矩阵特征值的求取 大庆石油学院学报 JOURNALOFDAQINGPETROLEUMINSTITUTE第29卷Vol.29第3期No.32005年6月Jun.2005 二 ...
- 机器视觉学习之--贝叶斯学习 MATLAB二维正态分布二维图
1.贝叶斯介绍 我个人一直很喜欢算法一类的东西,在我看来算法是人类智慧的精华,其中蕴含着无与伦比的美感.而每次将学过的算法应用到实际中,并解决了实际问题后,那种快感更是我在其它地方体会不到的. 一直想 ...
- 二维概率密度求解边缘密度
二维概率密度求解边缘密度 @(概率论) 已知f(x,y)f(x,y),求解fX(x),fY(y)f_X(x),f_Y(y)时,用的是下面的公式: fX(x)=∫+∞−∞f(x,y)dyfY(y)=∫+ ...
- MATLAB二维绘图(二)向图中添加标题,坐标轴,图标和文字信息
MATLAB二维绘图(二)向图中添加标题,坐标轴,图标和文字信息 1.添加标题.图例.x轴信息和y轴信息,示例: %% 添加标题 clear; clc; close all; x = 0:0.1:2* ...
- Matlab 二维网格图pcolor和imagesc区别
Matlab 二维网格图pcolor和imagesc区别 在速度或者衰减层析成像反演之后会得到地下介质每个网格点的速度/衰减值,在进行呈现的时候可以使用maltab的pcolor和imagesc进行二 ...
- matlab 二维矩阵变成一维矩阵
matlab 二维矩阵变成一维矩阵 1.一维变二维: https://blog.csdn.net/qq_40584593/article/details/90691276 reshape 2.a(:) ...
- Matlab二维图形绘制与图形处理
Matlab二维图形绘制与图形处理 一.二维图形绘制 1.极坐标图 2.散点图 3. 平面等值线图 二.图形处理 1.添加格栅,图例 和标注 2.定制坐标 3.在之前基础上继续作图 4.新建图形置于当 ...
- matlab二维绘图部分
matlab二维绘图部分 X,Y是向量,分别表示点集的横坐标和纵坐标 PLOT(X,Y,S) 符号函数(显函数.隐函数和参数方程)画图 (1) ezplot ezplot('f(x)',[a,b]) ...
- matlab二维三维图形绘制和坐标轴范围设置
matlab二维绘图 一. 二维图形(Two dimensional plotting) 1. 基本绘图函数(Basic plotting function):Plot, semilogx, se ...
最新文章
- R语定义函数对宽分布(wide distribution )、有偏分布(skew distribution)的数据进行对数变换(符号对数变换函数、signed log transformation)
- python流程控制语句-4.python流程控制语句介绍
- ios 中的小技巧 - 总有你想要的 一
- 【人脸识别】人脸识别必读论文
- Hyperledger Fabric 1.0 实战开发系列 第一课 系统环境搭建
- 支持批任务的Coscheduling/Gang scheduling
- 有机晶体数据库_Cambridge Structural Database 2017 晶体结构分析软件分享
- leetcode —— 589. N叉树的前序遍历 (使用到遍历多叉树的方法)
- java 日志 异步_log4j 详解异步日志的配置和测试
- 学习Linux让我进入了知名企业 原
- VeraCrypt文件硬盘加密使用教程
- Stuck Stack成过去时,OpenStack升级还可以这样玩?!
- TCP之报文格式解析
- python3笔记5--循环语句
- iEx.ec演示DApp的初步了解,快来“尝鲜”iEx.ec分布式云的强大算力
- about unit test in android ppt,选修5 unit4Learning about language.ppt
- js+jq实现坦克大战
- RGB与LAB优化互转
- Android音乐播放器开发(6)—ListView组件创建歌曲播放列表(内含原理分析)
- 慕课网微信点餐系统之商品信息加载不出来解决方案
热门文章
- 融云 IM 和 RTC 服务,「助攻」智能物流等客户打通链路、完善生态
- 十进制的转换与任意进制的互转
- 计算机教师的名言是,教师的名言
- java班组长竞聘演讲稿_班组长竞聘的演讲稿
- 学生成绩管理系统数据库
- 复旦情商课魅力女教师陈果上课语录(复制)
- rootfs 制作ubuntu_如何构建一个arm64 AArch64的Ubuntu rootfs
- python答案大全_大学慕课2020用Python玩转数据答案大全
- 用友开发者社区 - http://udn.yonyou.com
- 蓝桥杯刷题冲刺 | 倒计时20天