鏈�鏂扮簿鍝佹枃妗ｏ紝鐭ヨ瘑鍏变韩锛?/p>
绗簲娆′綔涓氾紙鍓嶄笁棰樺啓鍦ㄤ綔涓氱焊涓婏級
涓�銆?/p>
鐢ㄦ湁闄愬樊鍒嗘柟娉曟眰瑙ｄ竴缁撮潪瀹氬父鐑紶瀵兼柟绋嬶紝
鍒濆鏉′欢鍜岃竟鐣屾潯浠惰璇存槑
鏂囦欢锛?/p>
鐑墿鏁ｇ郴鏁拔?/p>
=const
锛?/p>
2
2
T
T
t
x
飦?/p>
飩?/p>
飩?/p>
飥?/p>
飩?/p>
飩?/p>
1.
鐢?/p>
Tylaor
灞曞紑娉曟帹瀵煎嚭
FTCS
鏍煎紡鐨勫樊鍒嗘柟绋?/p>
2.
璁ㄨ璇ユ柟绋嬬殑鐩稿鎬у拰绋冲畾鎬э紝骞惰鏄庣ǔ瀹氭�ц姹傚姹傝В宸垎鏂圭▼鐨勫奖鍝嶃�?/p>
3.
璇存槑璇ユ柟绋嬬殑绫诲瀷鍜屽畾瑙ｆ潯浠讹紝濡備綍鍦ㄧ▼搴忎腑瀹炵幇杩欎簺瀹氳В鏉′欢銆?/p>
4.
缂栧啓
M
鏂囦欢姹傝В涓婅堪鏂圭▼锛屽苟鐢ㄩ�傚綋鐨勬枃瀛楀绋嬪簭鍋氬嚭璇存槑銆?/p>
锛堥儴鍒嗙敱缃戠粶鎼滅储寰?
鍒帮紝娣诲姞锛屼慨鏀瑰悗寰楀埌銆?/p>
锛?/p>
function rechuandaopde
%
浠ヤ笅鎵�鐢ㄦ暟鎹紝闄や簡
t
鐨勮寖鍥存垜鏍规嵁棰樼洰瑕佹眰鍙栧埌浜?/p>
20000
锛屽叾浣欏潎浠?/p>
涓緱鏉?/p>
a=0.00001;%a
鐨勫彇鍊?/p>
xspan=[0 1];%x
鐨勫彇鍊艰寖鍥?/p>
tspan=[0 20000];%t
鐨勫彇鍊艰寖鍥?/p>
ngrid=[100 10];%
鍒嗗壊鐨勪唤鏁帮紝鍓嶉潰鐨勬槸
t
杞寸殑锛屽悗闈㈢殑鏄?/p>
x
杞寸殑
f=@(x)0;%
鍒濆�?/p>
g1=@(t)100;%
杈圭晫鏉′欢涓�
g2=@(t)100;%
杈圭晫鏉′欢浜?/p>
[T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%
璁＄畻鎵�璋冪敤鐨勫嚱鏁?/p>
[x,t]=meshgrid(x,t);
mesh(x,t,T);%
鐢诲浘锛屽苟涓旀妸鍧愭爣杞村悕绉版敼涓?/p>
x
锛?/p>
t
锛?/p>
T
xlabel('x')
ylabel('t')
zlabel('T')
T%
杈撳嚭娓╁害鐭╅樀
dt=tspan(2)/ngrid(1);%t
姝ラ暱
h3000=3000/dt;