| 查看: 2010 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
请教一个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 ] |
回复此楼» 猜你喜欢
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有136人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有120人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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人回复
|
版主，您好。我现在对方程做了些修改。方程代码如下： ～～～～～～～～～～～～～～～～～～～～～～～～～～～～ 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）有一个表达式。 版主，您看，这个方程组怎么解呢，谢谢了！版主帮帮忙了，已经为这个问题花费很久时间了 |
10楼2013-08-02 14:53:00
2楼2013-07-31 14:45:37
3楼2013-07-31 15:18:52
4楼2013-07-31 15:26:35
