24小时热门版块排行榜

返回列表

lkbiandou

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1230.5
帖子: 208
在线: 94.6小时
虫号: 945069
注册: 2010-01-19
专业: 化学反应工程

[交流] 遇到一个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

回复此楼

» 猜你喜欢

1楼 2013-11-06 11:06:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

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人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 lkbiandou 的主题更新

返回列表