Re: [問題]在計算中有關多次數列疊代的問題!!!!

看板MATLAB作者 (小胖)時間13年前 (2012/08/24 10:10), 編輯推噓1(100)
留言1則, 1人參與, 最新討論串2/2 (看更多)
把D1 改成D1(i) 或者是直接改成下面程式碼 clear clc Ea=0.45; D0=7*10^(-6); k=8.617*10^(-5); i=0:1:5400 ; y1=873-i*(873/5400) D=D0*exp(-Ea./(k*y1)); % Diffusion coefficient D1=sqrt(2*D.*i); % Diffusion length plot(i,D1,'color',[0 0 1]); ※ 引述《johnny4003 ()》之銘言: : 目前在計算有關擴散的問題,由於時間係數i為0~5400, : 代入y1(i)的值再代入D(y1),而所得的D(y1)需要再乘上 : i,才能得到結果的D1。 : 例如i=10 : 先計算出y1(10),將y1(10)代入D方程式得D(y1(10)),然 : 後再將其D(y1(10))*10開根號求結果。 : 請問下列方程式哪邊有問題?謝謝~~ : clear : clc : Ea=0.45; : D0=7*10^(-6); : k=8.617*10^(-5); : for i=0:1:5400 ; //時間 : y1=873-i*(873/5400) : D=D0*exp(-Ea/(k*y1)); // Diffusion coefficient : D1=sqrt(2*D*i); // Diffusion length : end : plot(i,D1,'color',[001]); -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 1.162.16.86

08/25 01:01, , 1F
感謝感謝~~
08/25 01:01, 1F
文章代碼(AID): #1GDkAZ2z (MATLAB)
文章代碼(AID): #1GDkAZ2z (MATLAB)