Re: [討論] 這個mpower該如何解??

看板MATLAB作者zupo (外派越南工作者)時間11年前 (2013/12/22 10:21)推噓0(0推 0噓 0→)

留言0則, 0人參與討論串2/2 (看更多)

※ 引述《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