【正文】 range=minmax(lamda39。) %計算級比的范圍x1=cumsum(x0)。 %累加運算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。Y=x0(2:n)。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。Dx+a*x=b39。,39。x(0)=x039。)。 %求微分方程的符號解x=subs(x,{39。a39。,39。b39。,39。x039。},{u(1),u(2),x0(1)})。 %代入估計參數(shù)值和初始值yuce1=subs(x,39。t39。,[0:n1])%求已知數(shù)據(jù)的預測值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運算，還原數(shù)據(jù)epsilon=x039。yuce %計算殘差delta=abs(epsilon./x039。) %計算相對誤差rho=1(*u(1))/(1+*u(1))*lamda39。 %計算級比偏差值，u(1)=a程序結果：yuce = +003 * Columns 1 through 5 Columns 6 through 10 epsilon = Columns 1 through 5 0 Columns 6 through 10 delta = Columns 1 through 5 0 Columns 6 through 10 rho = Columns 1 through 5 Columns 6 through 9 程序8煤炭消費量 clc,clearx0=[896396631084913272167791842819952207372100323100]。 %注意這里為列向量n=length(x0)。lamda=x0(1:n1)./x0(2:n) %計算級比range=minmax(lamda39。) %計算級比的范圍x1=cumsum(x0)。 %累加運算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。Y=x0(2:n)。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。Dx+a*x=b39。,39。x(0)=x039。)。 %求微分方程的符號解x=subs(x,{39。a39。,39。b39。,39。x039。},{u(1),u(2),x0(1)})。 %代入估計參數(shù)值和初始值yuce1=subs(x,39。t39。,[0:n1])%求已知數(shù)據(jù)的預測值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運算，還原數(shù)據(jù)epsilon=x039。yuce %計算殘差delta=abs(epsilon./x039。) %計算相對誤差rho=1(*u(1))/(1+*u(1))*lamda39。 %計算級比偏差值，u(1)=a程序結果：yuce = +004 * Columns 1 through 5 Columns 6 through 10 epsilon = +003 * Columns 1 through 5 0 Columns 6 through 10 delta = Columns 1 through 5 0 Columns 6 through 10 rho = Columns 1 through 5 Columns 6 through 9 程序9第一產(chǎn)業(yè)煤炭消費量clc,clearx0=[606363]。 %注意這里為列向量n=length(x0)。lamda=x0(1:n1)./x0(2:n) %計算級比range=minmax(lamda39。) %計算級比的范圍x1=cumsum(x0)。 %累加運算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。Y=x0(2:n)。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。Dx+a*x=b39。,39。x(0)=x039。)。 %求微分方程的符號解x=subs(x,{39。a39。,39。b39。,39。x039。},{u(1),u(2),x0(1)})。 %代入估計參數(shù)值和初始值yuce1=subs(x,39。t39。,[0:n1])%求已知數(shù)據(jù)的預測值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運算，還原數(shù)據(jù)epsilon=x039。yuce %計算殘差delta=abs(epsilon./x039。) %計算相對誤差rho=1(*u(1))/(1+*u(1))*lamda39。 %計算級比偏差值，u(1)=a程序結果：yuce = Columns 1 through 5 Columns 6 through 10 epsilon = Columns 1 through 5 0 Columns 6 through 10 delta = Columns 1 through 5 0 Columns 6 through 10 rho = Columns 1 through 5 Columns 6 through 9 程序10第二產(chǎn)業(yè)煤炭消費量clc,clearx0=[]。 %注意這里為列向量n=length(x0)。lamda=x0(1:n1)./x0(2:n) %計算級比range=minmax(lamda39。) %計算級比的范圍x1=cumsum(x0)。 %累加運算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。Y=x0(2:n)。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。Dx+a*x=b39。,39。x(0)=x039。)。 %求微分方程的符號解x=subs(x,{39。a39。,39。b39。,39。x039。},{u(1),u(2),x0(1)})。 %代入估計參數(shù)值和初始值yuce1=subs(x,39。t39。,[0:n1])%求已知數(shù)據(jù)的預測值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運算，還原數(shù)據(jù)epsilon=x039。yuce %計算殘差delta=abs(epsilon./x039。) %計算相對誤差rho=1(*u(1))/(1+*u(1))*lamda39。 %計算級比偏差值，u(1)=a程序結果：yuce = +004 * Columns 1 through 5 Columns 6 through 10 epsilon = +003 * Columns 1 through 5 0 Columns 6 through 10 delta = Columns 1 through 5 0 Columns 6 through 10 rho = Columns 1 through 5 Columns 6 through 9 程序11第三產(chǎn)業(yè)煤炭消費量clc,clearx0=[]。 %注意這里為列向量n=length(x0)。lamda=x0(1:n1)./x0(2:n) %計算級比range=minmax(lamda39。) %計算級比的范圍x1=cumsum(x0)。 %累加運算B=[*(x1(1:n1)+x1(2:n)),ones(n1,1)]。Y=x0(2:n)。u=B\Y %擬合參數(shù)u(1)=a,u(2)=bx=dsolve(39。Dx+a*x=b39。,39。x(0)=x039。)。 %求微分方程的符號解x=subs(x,{39。a39。,39。b39。,39。x039。},{u(1),u(2),x0(1)})。 %代入估計參數(shù)值和初始值yuce1=subs(x,39。t39。,[0:n1])%求已知數(shù)據(jù)的預測值y=vpa(x,6) %其中的6表示顯示6位數(shù)字yuce=[x0(1),diff(yuce1)] %差分運算，還原數(shù)據(jù)epsilon=x039。yuce %計算殘差delta=abs(epsilon./x039。) %計算相對誤差rho=1(*u(1))/(1+*u(1))*lamda39。 %計算級比偏差值，u(1)=a程序結果：8e1*t)yuce = Columns 1 through 5 Columns 6 through 10 Column 11 epsilon = Columns 1 through 5 0 Columns 6 through 10 Column 11 delta = Columns 1 through 5 0 Columns 6 through 10 Column 11 rho = Columns 1 through 5 Columns 6 through 10 34 / 35