作者zupo (外派越南工作者)
標題Re: [討論] 這個mpower該如何解??
時間Sun Dec 22 10:20:55 2013
※ 引述《NChewinggum (鈉鎂壬)》之銘言：
: 我輸入程式碼後出現這個錯誤提示：
: Undefined function 'mpower' for input arguments of type 'cell'
: 初學Ｍａｔｌａｂ幾乎甚麼都不懂，請教神人大大們這是甚麼意思呢？？
: 如果是說我的斜率ｍ沒有定義，後來定義了ｍ＝－１,１還是一樣
: 改了好久都沒辦法解決
: 在此請求各位幫小弟解惑，謝謝～
: 附上程式碼
: A=2.53*10^8, B=5.42*10^3, Cp=1004, Lv=2.5*10^6, g=9.8, Rd=287,dZ=1
: for k=1:10000, Z(1)=0, T(1)=305, P(1)=1000, m(1)=-1.1;
: es(T)= A*2.7183^{-[B/T(k)]}
: qvs=0.622*es(T)/P(k)
: Cp(k)=g*{1+[Lv*qvs/Rd*T(k)]}/{Cp+[Lv*Lv*qvs/(Rd*T(k)*T(k))]}
: m(k)=-g/Cp(k)
: P(k+1)=P(k)*exp(Z(k)/k)
: Z(k)=(Rd*T(k)/g)*ln(P(k)/P(k+1))
: T(k+1)=T(k)+m(k)*Z(k)
: end
程式有點在亂寫,講真的,初學者沒錯. 錯的是應該要先學到基本再來
這樣寫是不會有結果的
我幫你debug了一些. 但很明顯數學式跟邏輯應該是錯的
你在自己修一下
--------------------------------------------------------
function pttex143
%
es = zeros(length(1:10000),1);
qvs = zeros(length(1:10000),1);
Cp = zeros(length(1:10000),1);
m = zeros(length(1:10000),1);
P = zeros(length(1:10000),1);
Z = zeros(length(1:10000),1);
T = zeros(length(1:10000),1);
%
A=2.53*10^8;
B=5.42*10^3;
Cp(1,1)=1004;
Lv=2.5*10^6;
g=9.8;
Rd=287;
dZ=1;
Z(1,1)=0;
T(1,1)=305;
P(1,1)=1000;
m(1,1)=-1.1;
for k=1:10000
es(k,1)= A.*2.7183.^(-B./T(k,1));
qvs(k,1)=0.622.*es(k,1)./P(k,1);
Cp(k+1,1)=g.*(1+(Lv.*qvs(k,1)./Rd.*T(k,1))./(Cp(k,1)+(Lv.*Lv.*qvs(k,1)...
./(Rd.*T(k,1).*T(k,1)))));
m(k,1)=-g./Cp(k,1);
P(k+1,1)=P(k,1).*exp(Z(k,1)./k);
Z(k+1,1)=(Rd.*T(k,1)./g).*log(P(k,1)./P(k+1,1));
T(k+1,1)=T(k,1)+m(k,1).*Z(k,1);
end
k = 1:10001;
i = 1:10000;
figure(1)
plot(i,es,i,qvs)
legend('es','qvs')
------------------------------------------------------------------
所有數值都是一樣那表示不是迴圈有錯就是算式有錯
自己再看看哪裡有問題吧
--
1.MATLAB programming 2.ASPEN process simulation package
3.FORTRN programming 4.Advance Engineering Mathematics
5.Process Control Theory
6.Chemical Engineering Basic Theory(Kinetic.thermodynamics.transport)
7.Numerical Method and Analysis
8.MATLAB Toolbox.Simulink system basic design
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 113.162.89.91