求解热传导方程matlab
这是非稳态一维热传导的方法,也叫古典显格式。
如果是做数学建模,就别用了,这种方法计算量比较大,算的很慢,而且收敛不好。
但是如果实在没办法也能凑合用。
该改的地方我都用???代替了。
给个详细解释https://wenku.baidu.com/view/78a359d43b3567ec112d8a77.html?qq-pf-to=pcqq.group
function rechuandao() % Llist = ???
N = ; % 空间点数
M = ;
alfa = ??? % 导热 / (密度*比热)
lambda = 0.4; % 稳定条件
%************
zone=Llist()/N;%空间步长
z=:zone:Llist();%空间点坐标数组
z=z';%矩阵转置 n = lambda * (Llist()/N)^ / alfa(); TM = M * n; % 总时间 t=:n:TM; %时间点坐标数组
t=t'; %矩阵转置 %计算初值和边值
T=zeros(N+,M+); %赋温度矩阵初值
Ti=init_fun(z); %各空间点处的初始条件
To=border_funo(t); %内边界条件
Te=border_fune(t); %外边界条件
T(:,)=Ti; %(初始条件 温度都为37)
T(,:)=To; %(边界条件x=0处温度为)
T(N+,:)=Te; %(边界条件x=L 处温度为) %用差分法求出温度T与杆长L、时间t的关系
for k=:M %时间点数
n=;
while n<=N %空间点数
T(n,k+)=lambda*(T(n+,k)+T(n-,k))+(-*lambda+)*T(n,k);
n=n+;
end
end
%设置立体网格
for i=:M+
X(:,i)=z;
end
for j=:N+
Y(j,:)=t;
end
mesh(X,Y,T);
view([ - ]);
xlabel('长度');
ylabel('时间');
zlabel('温度');
csvwrite('csv.csv',T) function y=init_fun(z)%初值条件
y=???
return function y=border_funo(t)%z=0的边界条件
y=???
return function y=border_fune(t)%z=L的边界条件
y= ???
return
%
求解热传导方程matlab的更多相关文章
- 龙哥库塔法or欧拉法求解微分方程matlab实现
举例:分别用欧拉法和龙哥库塔法求解下面的微分方程 我们知道的欧拉法(Euler)"思想是用先前的差商近似代替倒数",直白一些的编程说法即:f(i+1)=f(i)+h*f(x,y)其 ...
- matlab中求解线性方程组的rref函数
摘自:http://www.maybe520.net/blog/987/ matlab中怎么求解线性方程组呢? matlab中求解线性方程组可应用克拉默法则(Cramer's Rule)即通过det( ...
- matlab 解方程组
1.解方程 最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法:(1)x=in ...
- 【数学建模】day01-线性规划问题
线性规划问题是在一组线性约束条件下,求线性目标函数最大/最小值的问题.这些约束条件有不等式约束.等式约束以及边界约束,这和中学讲的线性规划无异. 此类问题的MATLAB标准形式为: 其中,max问 ...
- 数学建模三剑客MSN
前言 不管是不是巴萨的球迷,只要你喜欢足球,就一定听说过梅西(Messi).苏亚雷斯(Suarez)和内马尔(Neymar)这个MSN组合.在众多的数学建模辅助工具中,也有一个犀利无比的MSN组合,他 ...
- 模拟退火算法Python编程(2)约束条件的处理
1.最优化与线性规划 最优化问题的三要素是决策变量.目标函数和约束条件. 线性规划(Linear programming),是研究线性约束条件下线性目标函数的极值问题的优化方法,常用于解决利用现有的资 ...
- 人工智能与智能系统2-> 机器人学2 | 时间与运动
<Robotics, Vision and Control>学习到第三章,我才发现这本书是有配套视频的,第二章看的好辛苦,很多地方生硬理解了一下,现在打算把视频再好好看一看,作为补充,也会 ...
- MATLAB求解代数方程、微分方程的一些常用指令
MATLAB版本:R2015b 1.求解符号矩阵的行列式.逆.特征值.特征向量 A = sym('[a11, a12; a21, a22]');deltaA = det(A)invA = inv(A) ...
- 蚁群算法求解旅行商问题(附c和matlab源代码)
前几天写了个模拟退火算法的程序,然后又陆陆续续看了很多群智能算法,发现很多旅行商问题都采用蚁群算法来求解,于是开始写蚁群算法的模板.网上关于蚁群算法的理论很多就不再这里赘述了,下面直接上代码和进行简单 ...
随机推荐
- jquery-post 异常
报错:Client closed connection before receiving entire response 前端请求: $.ajax({ type : 'POST' , url : '/ ...
- hihoCoder#1079(线段树+坐标离散化)
时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi和小Ho在回国之后,重新过起了朝7晚5的学生生活,当然了,他们还是在一直学习着各种算法~ 这天小Hi和小Ho所在的学 ...
- 机器学习:评价分类结果(ROC 曲线)
一.基础理解 1)定义 ROC(Receiver Operation Characteristic Curve) 定义:描述 TPR 和 FPR 之间的关系: 功能:应用于比较两个模型的优劣: 模型不 ...
- HashSet中是如何判断元素是否重复的
HashSet不能添加重复的元素,当调用add(Object)方法时候, 首先会调用Object的hashCode方法判hashCode是否已经存在,如不存在则直接插入元素: 如果已存在则调用Obje ...
- Oracle 归档开启切换和归档日志删除(单实例和RAC)
Oracle默认安装后,是没有开启归档模式的,需要手动开启. 开启归档--单实例如果archive log模式下不能正常startup,则先恢复成noarchive log,startup成功后,再s ...
- Celery-4.1 用户指南: Concurrency (并发)
简介 Eventlet 的主页对它进行了描述:它是一个python的并发网络库,可以让你更改如何运行你的代码而不是怎么编写代码. 对高可扩展非阻塞IO操作,它使用 epoll或者libevent. C ...
- [MySQL]修改mysql数据库的root密码的方法
方法1: 用SET PASSWORD命令 mysql -u root mysql> SET PASSWORD FOR 'root'@'localhost' = PASSWORD('newpass ...
- sqlplus 设置显示格式
使用sqlplus查询显示结果,显示很乱,下面有种方法可以让她显示的更好看些.1.设置显示的宽度:设置前可以先查看当前宽度: SQL> show linesize;linesize 100SQL ...
- javascript 中的JSON.stringify - 将对象和数组转换为json格式(来源于网络)
JSON.stringify 函数 (JavaScript) 将 JavaScript 值转换为 JavaScript 对象表示法 (Json) 字符串. JSON.stringi ...
- ios配合iTuns提取应用Documents下的文件到本地
出处:http://blog.csdn.net/jianandjan/article/details/50442988 有一些App需要通过使用iTunes让用户上传和下载文档.要让iOS程序支持iT ...