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

返回列表

小木蟲zb

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時(shí)
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

[求助] matlab程序的問題

小弟用matlab計(jì)算一個(gè)偏微分方程，程序運(yùn)行不出來，求大神相助，程序如下
function DBT_MNA%DBT用1表示,MNA用2表示
clear all;clc
A=importdata('DBT+MNA-RS-03.txt');%數(shù)據(jù)源文件中橫坐標(biāo)為g-MDF/g-AC
t2=A(:,1)/7.03;y2=A(:,2);%橫坐標(biāo)由g-MDF/g-AC化為h-1
q10=zeros(1,20);
c10=zeros(1,20);
q20=zeros(1,20);
c20=zeros(1,20);
y0=[c10 q10 c20 q20];
c0=7.03;
options=odeset('relTol',1e-6);
tspan=[0:0.02:t2(end)];
[t,y]=ode15s(@fangcheng,tspan,y0,options);
plot(t*7.03,y(:,40)/c0,'ro',t*7.03,y(:,80)/c0,'bo','LineWidth',3);hold on
plot(A(:,1),A(:,2),A(:,1),A(:,3),'LineWidth',3)
xlabel('Amount of treated MDF(g-MDF/g-AC)'),ylabel('c/c0'),axis([-0.05 120 -0.05 1.05]),
legend('擬合DBT','擬合MNA','實(shí)驗(yàn)曲線DBT','實(shí)驗(yàn)曲線MNA','Location','best') ,title('DBT+MNA-RS')
grid,hold off
%---------------------------------------------------------
function dydt=fangcheng(t,y)
kf=200;
a=0.363;%空隙率
u=0.361032/a;%床層顆粒間隙中流體流速
c0=7.03;
L=8.31e-2;
N=20;
dz=L/N;
p=436.8/(1-0.363);%顆粒表觀密度對曲線影響較大從434.2該為436.8
q1=y(1:N);c1=y(N+1:2*N);
q2=y(2*N+1:3*N);c2=y(3*N+1:4*N);
[ce1,ce2]=equilibrium(q1,q2);
%-----------定義dq/dt-----------------------------
for j=1:N
dq1dt(j)=kf/p*(c1(j)-ce1(j));
dq2dt(j)=kf/p*(c2(j)-ce2(j));
end
%---------定義dC/dt-----------------------------------
dc1dt(1)=-u*(c1(1)-c0)/dz-(1-a)/a*p*dq1dt(1);
dc2dt(1)=-u*(c2(1)-c0)/dz-(1-a)/a*p*dq2dt(1);
for j=2:N
dc1dt(j)=-u*(c1(j)-c1(j-1))/dz-(1-a)/a*p*dq1dt(j);
dc2dt(j)=-u*(c2(j)-c2(j-1))/dz-(1-a)/a*p*dq2dt(j);
end
dydt=[dq1dt dc1dt dq2dt dc2dt]';
function [ce1,ce2]=equilibrium(q1,q2)
[ce1,ce2]=solve('q1=0.45967ce1^0.3839/(ce1^0.0457+0.3*ce2^0.0003206',...
'q2-0.17356*ce2^(-1.55201)/(ce2^(-2)+0.1*ce1^(-1))');[ Last edited by 小木蟲zb on 2013-4-26 at 09:30 ]

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : DBT+MNA-RS-03.txt

2013-04-26 09:30:48, 990 bytes

» 猜你喜歡

本人考085602 化學(xué)工程專碩已經(jīng)有23人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)

1樓 2013-04-26 09:28:55

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

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

求助MATLAB做三維圖問題已經(jīng)有4人回復(fù)
matlab語句問題已經(jīng)有5人回復(fù)
matlab中遇到的畫圖小問題，程序沒問題。標(biāo)題要長@@@@@@@@@@@@@@@@@@@@@@ 已經(jīng)有8人回復(fù)
求助一個(gè)matlab關(guān)于分類討論的問題已經(jīng)有9人回復(fù)
matlab編程中遇到的問題已經(jīng)有4人回復(fù)
Matlab中save問題已經(jīng)有5人回復(fù)
問下這個(gè)問題怎么解決.....真心啊已經(jīng)有4人回復(fù)
matlab程序問題，總是出錯(cuò)，麻煩大家看看問題在那里？已經(jīng)有6人回復(fù)
請教關(guān)于matlab矩陣轉(zhuǎn)換問題已經(jīng)有6人回復(fù)
關(guān)于Matlab作圖的問題已經(jīng)有14人回復(fù)
matlab的函數(shù)調(diào)用問題已經(jīng)有3人回復(fù)
問個(gè)關(guān)于matlab里find函數(shù)的問題！已經(jīng)有4人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
求助使用MATLAB計(jì)算綠色葉片的覆蓋度的問題已經(jīng)有24人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
MATLAB 畫圖小問題已經(jīng)有11人回復(fù)
同樣的MATLAB程序，為何在別人的電腦上運(yùn)行不了已經(jīng)有21人回復(fù)
關(guān)于matlab 2010b 的一個(gè)小問題已經(jīng)有7人回復(fù)
matlab做函數(shù)圖象的問題已經(jīng)有6人回復(fù)

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

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

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