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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

百里道

木蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1797
散金: 5
帖子: 46
在線: 11.1小時
蟲號: 1070800
注冊: 2010-08-07
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[求助] 有限差分的實現(xiàn)

n=200;
L=1;
t=0.1;
Kn=1;
%輸入值
dx=L/n;
w=dx/Kn;
dt=w;
m=t/dt;
%規(guī)定數(shù)組
x=linspace(0,1,n);
e1=[n,m];
e2=[n,m];
e=[n,m];
%計算過程
%初始值
for j=2:n;
e1(j,1)=0;
e2(j,1)=0;
end
%邊界值
for k=1:m;
e1(1,k)=0.5;
e2(1,k)=0.5;
e1(n,k)=0;
e2(n,k)=0;
end
%內(nèi)部過程
for j=2:n-1;
for k=2:m;
      e1(j+1,k+1)=(1-w/2)*e1(j,k)+w/2*e2(j,k);
      e2(j-1,k+1)=(1-w/2)*e2(j,k)+w/2*e1(j,k);
end
end
for j=1:n;
for k=1:m;
      e(j,k)=e1(j,k)+e2(j,k);
end
end
plot(x,e,'k-');
grid on
請教一下規(guī)定數(shù)組這么寫合適嗎？

回復(fù)此樓

» 猜你喜歡

288求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有19人回復(fù)
材料考研調(diào)劑已經(jīng)有3人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有6人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有9人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有12人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有5人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有6人回復(fù)
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有7人回復(fù)

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

Matlab求解二階偏微分方程組，希望能給出相應(yīng)的求解代碼已經(jīng)有10人回復(fù)
請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。已經(jīng)有17人回復(fù)
有限差分離散二階偏微分方程后如何編譯matlab程序解方程組？已經(jīng)有3人回復(fù)
求助，含有常數(shù)項的差分方程怎么轉(zhuǎn)換成狀態(tài)空間方程？已經(jīng)有5人回復(fù)
matlab有限差分法求解復(fù)雜的偏微分方程組已經(jīng)有19人回復(fù)
matlab 積分求助已經(jīng)有8人回復(fù)
積分方程和微分方程已經(jīng)有4人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】有人用有限差分方法計算過G點聲子譜嗎？已經(jīng)有10人回復(fù)
【求助】求微弱差分電流轉(zhuǎn)差分電壓芯片型號已經(jīng)有4人回復(fù)
100金幣下載論文求助已經(jīng)有3人回復(fù)
【求助】請教蟲友關(guān)于matlab差分方程的求解和畫圖程序已經(jīng)有4人回復(fù)

1樓 2012-04-08 19:24:54

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

ltcc2000

木蟲 (著名寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 4074.5
散金: 365
紅花: 3
帖子: 2330
在線: 154.5小時
蟲號: 1647979
注冊: 2012-02-27
性別: GG
專業(yè): 控制理論與方法

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

★
感謝參與，應(yīng)助指數(shù) +1
臭水溝: 金幣+1, 謝謝交流~~ 2012-04-09 09:41:09

你規(guī)定的這個數(shù)組只有兩個數(shù)，需要修改一下

贊一下

回復(fù)此樓

海闊憑魚妖

2樓2012-04-09 08:37:43

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

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