这是非稳态一维热传导的方法,也叫古典显格式。

如果是做数学建模,就别用了,这种方法计算量比较大,算的很慢,而且收敛不好。

但是如果实在没办法也能凑合用。

该改的地方我都用???代替了。

给个详细解释https://wenku.baidu.com/view/78a359d43b3567ec112d8a77.html?qq-pf-to=pcqq.group

  1. function rechuandao() %
  2.  
  3. Llist = ???
  4. N = ; % 空间点数
  5. M = ;
  6. alfa = ??? % 导热 / (密度*比热)
  7. lambda = 0.4; % 稳定条件
  8. %************
  9. zone=Llist()/N;%空间步长
  10. z=:zone:Llist();%空间点坐标数组
  11. z=z';%矩阵转置
  12.  
  13. n = lambda * (Llist()/N)^ / alfa();
  14.  
  15. TM = M * n; % 总时间
  16.  
  17. t=:n:TM; %时间点坐标数组
  18. t=t'; %矩阵转置
  19.  
  20. %计算初值和边值
  21. T=zeros(N+,M+); %赋温度矩阵初值
  22. Ti=init_fun(z); %各空间点处的初始条件
  23. To=border_funo(t); %内边界条件
  24. Te=border_fune(t); %外边界条件
  25. T(:,)=Ti; %(初始条件 温度都为37)
  26. T(,:)=To; %(边界条件x=0处温度为)
  27. T(N+,:)=Te; %(边界条件x=L 处温度为)
  28.  
  29. %用差分法求出温度T与杆长L、时间t的关系
  30. for k=:M %时间点数
  31. n=;
  32. while n<=N %空间点数
  33. T(n,k+)=lambda*(T(n+,k)+T(n-,k))+(-*lambda+)*T(n,k);
  34. n=n+;
  35. end
  36. end
  37. %设置立体网格
  38. for i=:M+
  39. X(:,i)=z;
  40. end
  41. for j=:N+
  42. Y(j,:)=t;
  43. end
  44. mesh(X,Y,T);
  45. view([ - ]);
  46. xlabel('长度');
  47. ylabel('时间');
  48. zlabel('温度');
  49. csvwrite('csv.csv',T)
  50.  
  51. function y=init_fun(z)%初值条件
  52. y=???
  53. return
  54.  
  55. function y=border_funo(t)%z=0的边界条件
  56. y=???
  57. return
  58.  
  59. function y=border_fune(t)%z=L的边界条件
  60. y= ???
  61. return
  62. %

求解热传导方程matlab的更多相关文章

  1. 龙哥库塔法or欧拉法求解微分方程matlab实现

    举例:分别用欧拉法和龙哥库塔法求解下面的微分方程 我们知道的欧拉法(Euler)"思想是用先前的差商近似代替倒数",直白一些的编程说法即:f(i+1)=f(i)+h*f(x,y)其 ...

  2. matlab中求解线性方程组的rref函数

    摘自:http://www.maybe520.net/blog/987/ matlab中怎么求解线性方程组呢? matlab中求解线性方程组可应用克拉默法则(Cramer's Rule)即通过det( ...

  3. matlab 解方程组

    1.解方程 最近有多人问如何用matlab解方程组的问题,其实在matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法:(1)x=in ...

  4. 【数学建模】day01-线性规划问题

      线性规划问题是在一组线性约束条件下,求线性目标函数最大/最小值的问题.这些约束条件有不等式约束.等式约束以及边界约束,这和中学讲的线性规划无异. 此类问题的MATLAB标准形式为: 其中,max问 ...

  5. 数学建模三剑客MSN

    前言 不管是不是巴萨的球迷,只要你喜欢足球,就一定听说过梅西(Messi).苏亚雷斯(Suarez)和内马尔(Neymar)这个MSN组合.在众多的数学建模辅助工具中,也有一个犀利无比的MSN组合,他 ...

  6. 模拟退火算法Python编程(2)约束条件的处理

    1.最优化与线性规划 最优化问题的三要素是决策变量.目标函数和约束条件. 线性规划(Linear programming),是研究线性约束条件下线性目标函数的极值问题的优化方法,常用于解决利用现有的资 ...

  7. 人工智能与智能系统2-> 机器人学2 | 时间与运动

    <Robotics, Vision and Control>学习到第三章,我才发现这本书是有配套视频的,第二章看的好辛苦,很多地方生硬理解了一下,现在打算把视频再好好看一看,作为补充,也会 ...

  8. MATLAB求解代数方程、微分方程的一些常用指令

    MATLAB版本:R2015b 1.求解符号矩阵的行列式.逆.特征值.特征向量 A = sym('[a11, a12; a21, a22]');deltaA = det(A)invA = inv(A) ...

  9. 蚁群算法求解旅行商问题(附c和matlab源代码)

    前几天写了个模拟退火算法的程序,然后又陆陆续续看了很多群智能算法,发现很多旅行商问题都采用蚁群算法来求解,于是开始写蚁群算法的模板.网上关于蚁群算法的理论很多就不再这里赘述了,下面直接上代码和进行简单 ...

随机推荐

  1. Vue项目开启步骤

    ## Build Setup # install cnpm //安装淘宝CNPM镜像npm install -g cnpm --registry=http://registry.npm.taobao. ...

  2. CodeForces - 510B Fox And Two Dots (bfs或dfs)

    B. Fox And Two Dots time limit per test 2 seconds memory limit per test 256 megabytes input standard ...

  3. Oracle data guard学习

    Oracle data guard学习:三思笔记 Data guard 1data guard结构: data guard是一个集合,由一个primary数据库(生产数据库)和一个或多个standby ...

  4. java中try{}catch{}和finally{}的执行顺序问题

     今天我给大家讲解一下java的的错误和异常处理机制以及相关异常的执行顺序问题.如有不足的地方,欢迎批评指正~ 1.首相简单介绍一下java中的错误(Error)和异常(Exception) 错误和异 ...

  5. c# 设置水印,消除水印

    方案1: 图像处理 opencv etc 方案2: 开源框架,直接使用,已经优化 https://github.com/itext/itextsharp https://github.com/itex ...

  6. 【转】link_to中delete无效的问题解决

    问题一 1.问题描述:点击[删除]链接之后,跳到了show页面,完全没有弹出框提示 <%= link_to "删除", product_path(product.id), : ...

  7. Drools学习笔记4—Consequence/RHS

    Right Hand Side,当LHS所有条件满足才会执行 可以使用LHS部分定义的绑定变量.全局变量.或者直接编写JAVA代码. 提供宏函数操作working memory fact对象,如ins ...

  8. scrapy-redis源码抛析

    #scrapy-redis--->queue.py-->class FifoQueue 队列 LifoQueue(lastinfirstout栈) #self.server父类Base中链 ...

  9. IRedisClient

    事实上,IRedisClient里面的很多方法,其实就是Redis的命令名.只要对Redis的命令熟悉一点就能够非常快速地理解和掌握这些方法,趁着现在对Redis不是特别了解,我也对着命令来了解一下这 ...

  10. Java之——java.lang.NoSuchMethodException: [org.springframework.web.multipart.MultipartFile;.()

    转自:https://blog.csdn.net/l1028386804/article/details/65449355 ava.lang.NoSuchMethodException: [org.s ...