| 查看: 274 | 回复: 1 | ||
[求助]
请教一个线性方程组的求解问题
|
|
请各位帮忙看看这个程序哪里有问题，为什么总是错误呢。 就是当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 ] |
回复此楼» 猜你喜欢
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有278人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
-
德国亥姆霍兹Hereon中心招收两位医用镁合金腐蚀与LPSO相变方向2026公派博士生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一个matlab求解非线性方程组的问题
已经有9人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- MATLAB求解非线性方程组
已经有5人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 拟牛顿法求解复杂（带根号的）非线性方程组
已经有17人回复
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 第5章_解线性方程组的直接方法.ppt
已经有13人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 非线性方程组的解法
已经有5人回复
- 请教用matlab求解一个非线性偏微分方程组的数值解
已经有8人回复
- 解线性代数方程组
已经有10人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求教 多元指数方程如何解 有软件否？
已经有6人回复
- 【求助】求解非线性方程
已经有8人回复
- 【讨论】关于线性试验中线性方程的标准
已经有8人回复
2楼2013-08-04 20:17:13