版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

lkbiandou

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1230.5
帖子: 208
在線: 94.6小時
蟲號: 945069
注冊: 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

[求助] 請教一個matlab求解非線性方程組的問題

問題是：
T(i,j,1)=初值，就是整個第一頁都是一樣的值，請問這個該怎么處理？
現(xiàn)在編的程序根本無法運行起來的。各位幫幫忙，幫我看看怎么辦。
拜謝各位了！

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;
  %反應(yīng)器入口
  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 ]

回復(fù)此樓

» 猜你喜歡

光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有201人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有79人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
matlab如何求解一個非線性微分方程組已經(jīng)有8人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
matlab解非線性方程組解不開已經(jīng)有3人回復(fù)
請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。已經(jīng)有17人回復(fù)
求Matlab高手解決線性方程組的迭代求解問題已經(jīng)有17人回復(fù)
求教：線性系統(tǒng)與線性微分方程的關(guān)系。非線性系統(tǒng)能用傳遞函數(shù)來描述嗎？已經(jīng)有7人回復(fù)
matlab求解非線性方程組已經(jīng)有3人回復(fù)
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復(fù)
請教用matlab求解一個非線性偏微分方程組的數(shù)值解已經(jīng)有8人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
初學(xué)者，請教大家關(guān)于matlab解多元二次方程組已經(jīng)有4人回復(fù)
求助如何用matlab取出直線上的點已經(jīng)有10人回復(fù)
matlab中怎么取方程組的一組解并畫圖已經(jīng)有5人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求助有關(guān)Matlab有約束非線性最優(yōu)化問題已經(jīng)有9人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】請教非齊次常微分方程組的解析解法已經(jīng)有4人回復(fù)
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復(fù)
【求助】問一個關(guān)于matlab中l(wèi)egend的問題已經(jīng)有11人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)
【求助】求一個用牛頓迭代法解非線性方程組的程序已經(jīng)有10人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

去掉子函數(shù)里的clc;clear;

贊一下

回復(fù)此樓

showmethemoney

2樓2013-07-31 14:45:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lkbiandou

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1230.5
帖子: 208
在線: 94.6小時
蟲號: 945069
注冊: 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

2樓: Originally posted by csgt0 at 2013-07-31 14:45:37
去掉子函數(shù)里的clc;clear;

版主，你好。謝謝你的幫助！
可是去掉之后還是出現(xiàn)了錯誤，能幫我再看看嗎？
謝謝您了

贊一下

回復(fù)此樓

3樓2013-07-31 15:18:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

引用回帖:

3樓: Originally posted by lkbiandou at 2013-07-31 15:18:52
版主，你好。謝謝你的幫助！
可是去掉之后還是出現(xiàn)了錯誤，能幫我再看看嗎？
謝謝您了...

function f = TxEquations(T)
的返回的f是哪個？你沒有給返回值

贊一下

回復(fù)此樓

showmethemoney

4樓2013-07-31 15:26:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lkbiandou

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1230.5
帖子: 208
在線: 94.6小時
蟲號: 945069
注冊: 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

4樓: Originally posted by csgt0 at 2013-07-31 15:26:35
function f = TxEquations(T)
的返回的f是哪個？你沒有給返回值...

版主，你好。
我想解決的問題是這樣的，T(i,j,1)=Tin表示的是第一面的溫度，也是初值，然后T(i,j,k)表示第k面上每個點的溫度，他與其他溫度點的關(guān)系在for語句里面了，我現(xiàn)在想求解的就是在初值條件下，每個面上每個點的溫度，也就是T。
您看我應(yīng)該怎么做呢

贊一下

回復(fù)此樓

5樓2013-07-31 15:30:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

引用回帖:

5樓: Originally posted by lkbiandou at 2013-07-31 15:30:04
版主，你好。
我想解決的問題是這樣的，T(i,j,1)=Tin表示的是第一面的溫度，也是初值，然后T(i,j,k)表示第k面上每個點的溫度，他與其他溫度點的關(guān)系在for語句里面了，我現(xiàn)在想求解的就是在初值條件下，每個面上每 ...

就是從第一層T計算剩下的9層T嗎，那你直接修改子程序，不要主程序了，不需要解方程呀。不過你的for循環(huán)里問題很多，跟算法有關(guān)，得仔細調(diào)調(diào)。

贊一下

回復(fù)此樓

showmethemoney

6樓2013-07-31 16:09:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lkbiandou

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1230.5
帖子: 208
在線: 94.6小時
蟲號: 945069
注冊: 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

6樓: Originally posted by csgt0 at 2013-07-31 16:09:39
就是從第一層T計算剩下的9層T嗎，那你直接修改子程序，不要主程序了，不需要解方程呀。不過你的for循環(huán)里問題很多，跟算法有關(guān)，得仔細調(diào)調(diào)。...

是的，我就是想通過第一層來計算其余各層的，那就是去掉解方程那個選項？直接用子程序處理啊？
版主您覺得for循環(huán)主要是那些問題了？

贊一下

回復(fù)此樓

7樓2013-07-31 16:30:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

引用回帖:

7樓: Originally posted by lkbiandou at 2013-07-31 16:30:48
是的，我就是想通過第一層來計算其余各層的，那就是去掉解方程那個選項？直接用子程序處理啊？
版主您覺得for循環(huán)主要是那些問題了？...

不知道是不是這樣的。主要是矩陣維數(shù)超了

CODE:

 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=0.1;

 %反應(yīng)器入口

 a5=1; a6=1; d33=2; d44=3;

 for k=1:K-1; 

    for j=1:n; 

        for i=1:m-1;

             %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);  %變更k層第2行數(shù)據(jù)，第一層不變

             %AB段，j=1的處理

             T(i+1,1,k)=T(i,1,k)+0;           %變更k層第1列數(shù)據(jù)

             %CD段，j=n

             T(i+1,n,k)=T(i,n,k)+0;           %變更k層第n列數(shù)據(jù)

             %除邊界以外濃度和溫度的表達式

             T(i,j,k+1)=T(i,j,k)+dz*(a5*(d33+d44)+a6);     %當(dāng)前層計算下一層     

        end            

     end   

 end

贊一下

回復(fù)此樓

showmethemoney

8樓2013-07-31 17:00:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lkbiandou

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1230.5
帖子: 208
在線: 94.6小時
蟲號: 945069
注冊: 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

8樓: Originally posted by csgt0 at 2013-07-31 17:00:17
不知道是不是這樣的。主要是矩陣維數(shù)超了

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 ...

之前一直有個疑問，就是for循環(huán)的時候，i，j，k的取值范圍該怎么辦呢，。
我自己再嘗試下，今天非常感謝你，總體來說是有溫度數(shù)據(jù)出來了，比起之前那樣全是錯誤的提示，要進步了許多了。

贊一下

回復(fù)此樓

9樓2013-07-31 22:54:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lkbiandou

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1230.5
帖子: 208
在線: 94.6小時
蟲號: 945069
注冊: 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

版主，您好。我現(xiàn)在對方程做了些修改。方程代碼如下：
～～～～～～～～～～～～～～～～～～～～～～～～～～～～
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;
  %反應(yīng)器入口
  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;
～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～
就是當(dāng)i=1的時候，有一個表達式，但是j，k可以隨意變化的。
然后就是當(dāng)j=1或者n的時候，i，k可以隨意變化。
k=1的時候是初始條件，就是T(i,j,1)=483.15固定值。
在上述條件之外，T（i,j,k）有一個表達式。
版主，您看，這個方程組怎么解呢，謝謝了！版主幫幫忙了，已經(jīng)為這個問題花費很久時間了

贊一下

回復(fù)此樓

10樓2013-08-02 14:53:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 lkbiandou 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 085602 289分求調(diào)劑 +3	WWW西西弗斯 2026-03-24	3/150	2026-03-24 14:26 by JourneyLucky
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線 +3	哈哈157349 2026-03-21	3/150	2026-03-24 14:11 by zhyzzh
[考研] 材料292調(diào)劑 +6	橘頌思美人 2026-03-23	6/300	2026-03-24 13:32 by 連洋洋洋
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +12	吃吃吃才有意義 2026-03-19	13/650	2026-03-24 10:06 by NCUT_Yu
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	6/300	2026-03-23 21:45 by 子系博
[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風(fēng)里
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-20	3/150	2026-03-22 16:00 by ColorlessPI
[考研] 354求調(diào)劑 +7	Tyoumou 2026-03-18	10/500	2026-03-22 11:11 by 人來盛
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky

24小時熱門版塊排行榜

[求助] 請教一個matlab求解非線性方程組的問題

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: