版塊導(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í)熱門版塊排行榜

返回列表

王冉冉

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 190.5
散金: 5
紅花: 1
帖子: 76
在線: 33.4小時(shí)
蟲號(hào): 1872053
注冊(cè): 2012-06-28
性別: MM
專業(yè): 信息處理方法與技術(shù)

[求助] matlab求解一個(gè)非齊次偏微分方程

方程的具體形式在文件附件中（3）中，用C-N方法解的話我編寫的程序?yàn)?br />
p=2.1;%密度
c=1.17;%比熱容
ke=0.0063;%熱導(dǎo)率
g=0.93;%各向異性系數(shù)
u1=3;%吸收系數(shù)
u2=260;%散射系數(shù)
dx=0.015;
x=0:dx:0.15;
J1=length(x);
dt=1e-4;
N=3001;
t=(0:N-1)*dt;

u3=sqrt(3*u1*(u1+(1-g)*u2));
l=ke*dt/(p*c*(dx^2));
d=dt/(p*c);
I0=180;

S=zeros(N,J1);
I=zeros(N,J1);
for i=2:J1-1
   for m=1:9
         I(m,1)=I0;
         I(m,i)=I0;
         I(m,J1)=I0;
         S(m,1)=u1*I(1,1);
         S(m,i)=u1*I(1,i)*exp(-u3*(i-1)*dx);
         S(m,J1)=u1*I(1,11)*exp(-u3*10*dx);
         S(10:N,

=0;
end
end
u=zeros(N,J1);

ubc_l=S(:,1)/(p*c);
ubc_r=S(:,11)/(p*c);
u_initial=310*ones(1,J1);
u(1,

=u_initial;
a=ones(J1,1)*(l+1);
b=-l/2*ones(J1-1,1);
c=-l/2*ones(J1-1,1);
f=zeros(J1,1);
for n1=2:N
temp1=l/2*ubc_l(n1-1)+(1-l)*u(n1-1,1)+l/2*u(n1-1,2);
f(1)=temp1+l/2*ubc_l(n1)+1000*d*S(n1-1,1);
f(2:J1-1)=l/2*u(n1-1,1:J1-2)+(1-l)*u(n1-1,2:J1-1)+l/2*u(n1-1,3:J1)+1000*d*S(n1-1,2:J1-1);
f(J1)=l/2*u(n1-1,J1-1)+(1-l)*u(n1-1,J1)+l/2*ubc_r(n1-1)+l/2*ubc_r(n1)+1000*d*S(n1-1,J1);
u(n1,

=my_chasing(a,b,c,f);
end
plot(t,u(:,1)-273)
其中S相當(dāng)于Q,u相當(dāng)于T,為什么我畫出的圖形衰減太快，my_chasing是追趕法
另外對(duì)于這類方程我也不清楚它的邊界調(diào)節(jié)如何確定。
請(qǐng)專家們幫幫我萬分感謝�。。。。�

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 新建WPS文字文檔(3).wps

2012-07-02 17:42:22, 12.5 K

» 猜你喜歡

290求調(diào)劑已經(jīng)有7人回復(fù)
一志愿中南化學(xué)（0703）總分337求調(diào)劑已經(jīng)有8人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑一志愿南京航空航天大學(xué)289分已經(jīng)有3人回復(fù)
A區(qū)線材料學(xué)調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩297已過四六級(jí)求調(diào)劑推薦已經(jīng)有11人回復(fù)
316求調(diào)劑已經(jīng)有5人回復(fù)
一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有5人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有8人回復(fù)
329求調(diào)劑已經(jīng)有8人回復(fù)

1樓 2012-07-02 17:50:52

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

這樣一個(gè)偏微分方程有沒有解析解？已經(jīng)有4人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
【求助】matlab求解一個(gè)簡(jiǎn)單的積分（急！�。。� 已經(jīng)有7人回復(fù)
【求助】用pdepe求解偏微分方程已經(jīng)有10人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】如何求解下列偏微分方程已經(jīng)有13人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】用MATLAB求解偏微分方程的數(shù)值解懸賞20 已經(jīng)有8人回復(fù)
【求助】求MAtlab求解程序！�。� 已經(jīng)有20人回復(fù)
【分享】PDECOL-經(jīng)典數(shù)值求解非線性偏微分方程組源程序已經(jīng)有142人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

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