应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 请教一个线性方程组的求解问题
21/1返回列表
查看: 260  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lkbiandou

金虫 (小有名气)

[求助] 请教一个线性方程组的求解问题

请各位帮忙看看这个程序哪里有问题，为什么总是错误呢。
就是当i=1的时候，有一个表达式，但是j，k可以随意变化的。
然后就是当j=1或者n的时候，i，k可以随意变化。
k=1的时候是初始条件，就是T(i,j,1)=483.15固定值。
在上述条件之外，T（i,j,k）有一个表达式。
这个方程组怎么解呢，谢谢了！帮帮忙了，已经为这个问题花费很久时间了
～～～～～～～～～～～～～～～～～～～～～～～～～～～～
clear all ; clc;
T(1:10,1:10,1)=483.15;
m=10; n=10;K=10;
Ramda = 0.34;       % W/(m K)
Kw=25;
rho=1440;Cor=0.9;
Tw=483.15;
R=0.08;
x0=0.3959; dx=0.06041; df=pi/6/10;
dz=0.1;
%反应器入口
P=5;w=5; R0=0.0125;L=10;Dr=1;G=500;Cp=1;M=100;
a5=L*M*Ramda/R0^2/G/Cp; a6=rho*Cor*L*M/G/Cp*(1+1);

for k=1:K;
    for j=1:n;
        for i=1:m;                        
             %k=1，入口位置
       T(i,j,1)=483.15;
             d33=1/(x0+i*dx)^2*(T(i,j+1,k)-2*T(i,j,k)+T(i,j-1,k))/df^2;
             d44=(T(i+1,j,k)-2*T(i,j,k)+T(i-1,j,k))/dx^2+1/(x0+i*dx)*(T(i+1,j,k)-T(i-1,j,k))/(2*dx);
           %除边界以外浓度和温度的表达式
      T(i,j,k+1)=T(i,j,k)+dz*(a5*(d33+d44)+a6);     %当前层计算下一层
      T(1,j,k)=1/(3+2*dx*Kw*R/Ramda)*(4*T(2,j,k)-T(3,j,k)+2*dx*Kw*R*Tw/Ramda);%，i=1
             T(i,1,k)=(4*T(i,2,k)-T(i,3,k))/3;%，j=1            
             T(i,n+1,k)=(4*T(i,n,k)-T(i,n-1,k))/3;%，j=n     
        end            
     end   
end
T
～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～

[ Last edited by lkbiandou on 2013-8-3 at 21:40 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-08-03 21:37:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lkbiandou

金虫 (小有名气)
没有人愿意帮忙吗？
赞一下 回复此楼
2楼2013-08-04 20:17:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lkbiandou 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭