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

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

24小時(shí)熱門版塊排行榜

返回列表

lkbiandou

金蟲 (小有名氣)

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

[求助] 請(qǐng)教一個(gè)matlab求解非線性方程組的問題

問題是：
T(i,j,1)=初值，就是整個(gè)第一頁都是一樣的值，請(qǐng)問這個(gè)該怎么處理？
現(xiàn)在編的程序根本無法運(yù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;
         %除邊界以外濃度和溫度的表達(dá)式
         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論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有112人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有79人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)

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

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

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

【答案】應(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小時(shí)
蟲號(hào): 945069
注冊(cè): 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

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

引用回帖:

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

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

贊一下

回復(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小時(shí)
蟲號(hào): 945069
注冊(cè): 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

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

引用回帖:

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

就是從第一層T計(jì)算剩下的9層T嗎，那你直接修改子程序，不要主程序了，不需要解方程呀。不過你的for循環(huán)里問題很多，跟算法有關(guān)，得仔細(xì)調(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小時(shí)
蟲號(hào): 945069
注冊(cè): 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

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

引用回帖:

7樓: Originally posted by lkbiandou at 2013-07-31 16:30:48
是的，我就是想通過第一層來計(jì)算其余各層的，那就是去掉解方程那個(gè)選項(xiàng)？直接用子程序處理�。�
版主您覺得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ù)

             %除邊界以外濃度和溫度的表達(dá)式

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

        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小時(shí)
蟲號(hào): 945069
注冊(cè): 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 ...

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

贊一下

回復(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小時(shí)
蟲號(hào): 945069
注冊(cè): 2010-01-19
專業(yè): 化學(xué)反應(yīng)工程

引用回帖:

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

贊一下

回復(fù)此樓

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

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

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

返回列表

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

[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike
[考研] 080500求調(diào)劑 +3	zzzzfan 2026-03-24	3/150	2026-03-24 16:38 by barlinike
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 336求調(diào)劑 +4	收到VS 2026-03-20	4/200	2026-03-23 19:02 by macy2011
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 297求調(diào)劑 +11	戲精丹丹丹 2026-03-17	12/600	2026-03-21 17:47 by ColorlessPI
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[基金申請(qǐng)] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky

24小時(shí)熱門版塊排行榜

[求助] 請(qǐng)教一個(gè)matlab求解非線性方程組的問題

» 猜你喜歡

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

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

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

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

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

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