| 查看: 517 | 回复: 0 | |||
[交流]
遇到一个matlab的问题
|
最近碰到一个问题，在求解哦偏微分方程的时候，采用了有限差分法。但是结果总是不如人意。能不能麻烦各位帮我按照你们的思路写一下呢？不甚感激，当然如果有人是有偿代写，可以跟我联系，我们交流一下！ 然后我自己写了点， function crack_2D clear all; clc Ramda = 0.45; % W/(m K) %d33=2;d44=3; %T(1:11,1:11,1:11)=0; T(1:10,1:10)=483; x0=T(1:10,1:10); %X=fsolve(@temperature,x0); %function f=temperature(T) m=10; n=10;K=10; Kw=1.2;Tw=513;R=0.08; Ramda = 0.45; % W/(m K) x0=0.3959; dx=0.06041;df=pi/6/10;dz=0.1; P=5;w=34.572; R0=0.0125;L=10;Dr=1;G=500;Cp=1;M=8.2636;rho=1440;Cor=0.9; a5=L*M*Ramda/R0^2/G/Cp; a6=rho*Cor*L*M/G/Cp*0.32; %(-r1(i,j,k)*HCO(i,j,k)-r2(i,j,k)*HCO2(i,j,k)); %*************************************************************************% for h=1:100 for i=1:m; T(i,1)=483; %这个是入口条件 %T(i,n)=(4*T(i,n-1)-T(i,n-2))/3; T(i,n+1)=(4*T(i,n)-T(i,n-1))/3; %这些表达式都是参考别的文献得来 end for j=2:n; T(1,j)=1/(1+Ramda/(Kw*R*dx))*(Tw+Ramda/(Kw*R*dx)*(T(2,j))); %这个是紧挨加热管的温度表达式 T(m+1,j)=(4*T(m,j)-T(m-1,j))/3; %边界条件 end for j=1:n; for i=2:m; d44=(T(i+1,j)-2*T(i,j)+T(i-1,j)+T(i+1,j+1)-2*T(i,j+1)+T(i-1,j+1))/(2*dx^2)+1/(x0+i*dx)*(T(i+1,j)-T(i-1,j)+T(i+1,j+1)-T(i-1,j+1))/(4*dx); T(i,j+1)=T(i,j)+dz*a5*d44+dz*a6; %这个是除边界外，也就是微分方程表达式 %T(i,n+1)=T(i,n-1); end end end f=T 希望高手能按照你们的思路帮我看看，不要直接看我的程序，我的很垃圾的。 谢谢各位了。 未命名.jpg |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
matlab的plot画图（如何画出与取值数目不一样的标记符）
已经有4人回复
- matlab的m文件关联问题
已经有17人回复
- Matlab数值积分结果是虚数
已经有6人回复
- 一个matlab 多项式问题~~
已经有3人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- 请教大家一个反卷积的问题~~帮忙推导或者matlab编程计算~~
已经有5人回复
- matlab编程中遇到的问题
已经有4人回复
- 求助Matlab中如何把几个比特合为一个symbol
已经有4人回复
- 最近用matlab采集数据遇到一个问题，求助！
已经有6人回复
- matlab 积分求助
已经有8人回复
- 求matlab 解方程的问题
已经有5人回复
- 急求一个matlab编码，各位如果知道一定要帮我一把啊
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- Matlab中RBF建模，仿真是遇到的问题
已经有6人回复
- 关于matlab 2010b 的一个小问题
已经有7人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- hysys模拟一个例题
已经有6人回复
- matlab做函数图象的问题
已经有6人回复
- 【求助】matlab如何只读取图像的一部分？
已经有11人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30