24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lkbiandou

金虫 (小有名气)

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

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

问题是：
T(i,j,1)=初值，就是整个第一页都是一样的值，请问这个该怎么处理？
现在编的程序根本无法运行起来的。各位帮帮忙，帮我看看怎么办。
拜谢各位了！

function myEquations_test
clear all ; clc;
Tin(1:10,1:10,1)=483.15;
x0=Tin(1:10,1:10,1)
T=fsolve(@TxEquations,x0);
function f = TxEquations(T)
global n F G rc dz M  i j k
clc;clear;
m=10; n=10;K=10;
Ramda = 0.45;    % W/(m K)
Kw=1.2;
Tw=483.15;
R=0.08;
x0=0.3959; dx=0.06041;
dz=0.1;
Tin=483.15;
  %反应器入口
  a5=1; a6=1; d33=2; d44=3;
%  T(:,:,

;
for k=1:K;
for j=1:n;
      for i=2:m;
         %k=1，入口位置
         T(i,j,1)=483.15;
         %初始条件AD段的处理,i=1
         T(2,j,k)=T(1,j,k)-dx*Kw*R/Ramda*(T(1,j,k)-Tw);
         %AB段，j=1的处理
         T(i+1,1,k)=T(i,1,k)+0;
         %CD段，j=n
         T(i+1,n,k)=T(i,n,k)+0;
         %除边界以外浓度和温度的表达式
         T(i,j,k+1)=T(i,j,k)+dz*(a5*(d33+d44)+a6);
      end
end
end
T
%--------------------------------------------------------------------------

[ Last edited by lkbiandou on 2013-7-30 at 23:04 ]

回复此楼

» 猜你喜欢

🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队已经有0人回复
第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有271人回复
求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有2人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
matlab如何求解一个非线性微分方程组已经有8人回复
MATLAB求解非线性方程组已经有5人回复
matlab解非线性方程组解不开已经有3人回复
请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。已经有17人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
求教：线性系统与线性微分方程的关系。非线性系统能用传递函数来描述吗？已经有7人回复
matlab求解非线性方程组已经有3人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
请教用matlab求解一个非线性偏微分方程组的数值解已经有8人回复
求助matlab---fsolve解非线性方程组已经有6人回复
初学者，请教大家关于matlab解多元二次方程组已经有4人回复
求助如何用matlab取出直线上的点已经有10人回复
matlab中怎么取方程组的一组解并画图已经有5人回复
matlab求解非线性方程组已经有16人回复
求助有关Matlab有约束非线性最优化问题已经有9人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】请教非齐次常微分方程组的解析解法已经有4人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】问一个关于matlab中legend的问题已经有11人回复
【求助】求教matlab解非线性方程组已经有9人回复
【求助】求一个用牛顿迭代法解非线性方程组的程序已经有10人回复

1楼2013-07-30 16:10:51

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

lkbiandou

金虫 (小有名气)

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

引用回帖:

8楼: Originally posted by csgt0 at 2013-07-31 17:00:17
不知道是不是这样的。主要是矩阵维数超了

clear all ; clc;
T(1:10,1:10,1)=483.15;
m=10; n=10;K=10;
Ramda = 0.45; % W/(m K)
Kw=1.2;
Tw=483.15;
R=0.08;
x0=0.3959; dx=0.06041;
dz ...

版主，您好。我现在对方程做了些修改。方程代码如下：
～～～～～～～～～～～～～～～～～～～～～～～～～～～～
m=10; n=10;K=10;
Ramda = 0.45;    % W/(m K)
Kw=1.2;
Tw=483.15;
R=0.08;
x0=0.3959; dx=0.06041; df=pi/6/10;
dz=0.1;
Tin=483.15;
  %反应器入口
  a5=1; a6=1;
  %d33=2; d44=3;
%  T(:,:,

;
for k=2:K;
for j=1:n;
      for i=1:m;
         while i<=1;
            %T(i,j,1)=483.15;
            T(1,j,k+1)=T(1,j,k)+dz*(a5*(1/(x0+i*dx)^2*(T(1,j+1,k)-2*T(1,j,k)+T(1,j-1,k))/df^2+...
            (T(2,j,k)-2*T(1,j,k)+T(1,j,k)+dx*Kw*R/Ramda*(T(1,j,k)-Tw))/dx^2+...
            1/(x0+i*dx)*(T(2,j,k)-T(1,j,k)-dx*Kw*R/Ramda*(T(1,j,k)-Tw)/(2*dx)))+a6); %i=1的处理
         end
         while j<=1;
            T(i,1,k+1)=T(i,1,k)+dz*(a58(1/(x0+i*dx)^2*(T(i,2,k)-2*T(i,1,k)+T(i,2,k))/df^2+...
            (T(i+1,1,k)-2*T(i,1,k)+T(i-1,1,k))/dx^2+1/(x0+i*dx)*(T(i+1,1,k)-T(i-1,1,k))/(2*dx))+a6);%j=1的处理
         end
         while j>=n;
            T(i,n,k+1)=T(i,n,k)+dz*(a5*(1/(x0+i*dx)^2*(T(i,n-1,k)-2*T(i,n,k)+T(i,n-1,k))/df^2+...
            (T(i+1,n,k)-2*T(i,n,k)+T(i-1,n,k))/dx^2+1/(x0+i*dx)*(T(i+1,n,k)-T(i-1,n,k))/(2*dx))+a6);%j=n的处理
         end
         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);
      end
end
end
f=T;
～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～
就是当i=1的时候，有一个表达式，但是j，k可以随意变化的。
然后就是当j=1或者n的时候，i，k可以随意变化。
k=1的时候是初始条件，就是T(i,j,1)=483.15固定值。
在上述条件之外，T（i,j,k）有一个表达式。
版主，您看，这个方程组怎么解呢，谢谢了！版主帮帮忙了，已经为这个问题花费很久时间了