GM(1,N)動態模型的建立,用灰色系統理論的matlab怎麼樣編寫程式 5

時間 2025-03-12 20:10:13

gm(1,n)動態模型的建立,用灰色系統理論的matlab怎麼樣編寫程式

1樓:網友

clc,clear

x0=[n=length(x0);

lamda=x0(1:n-1)./x0(2:n)range=minmax(lamda)

x1=cumsum(x0)

for i=2:n

z(i)=;

endb=[-z(2:n)',ones(n-1,1)];

y=x0(2:n)';

u=b\yx=dsolve('dx+a*x=b','x(0)=x0');

x=subs(x,,)

yuce1=subs(x,'t',[0:n-1]);

digits(6),y=vpa(x) %為提高**精度,先計算**值,再顯示微分方程的解。

yuce=[x0(1),diff(yuce1)]epsilon=x0-yuce %計算殘差。

delta=abs(epsilon./x0) %計算相對誤差。

rho=1-( 計算級比偏差值。

改變x0的輸入值就行了。

2樓:塗顏帥

描述的不是很明白。

灰色系統**gm(1,1)模型的matlab源** 詳細用法啊

3樓:網友

clc,clear

x0=[n=length(x0);

lamda=x0(1:n-1)./x0(2:n)range=minmax(lamda)

x1=cumsum(x0)

for i=2:n

z(i)=;

endb=[-z(2:n)',ones(n-1,1)];

y=x0(2:n)';

u=b\yx=dsolve('dx+a*x=b','x(0)=x0');

x=subs(x,,)

yuce1=subs(x,'t',[0:n-1]);

digits(6),y=vpa(x) %為提高**精度,先計算**值,再顯示微分方程的解。

yuce=[x0(1),diff(yuce1)]epsilon=x0-yuce %計算殘差delta=abs(epsilon./x0) %計算相對誤差rho=1-( 計算級比偏差值。

你自己改變x0的輸入值就行了。

灰色系統**gm(1,1)模型的matlab詳細源**

4樓:網友

%灰色**模型。

function y=gm11(x,n)

x為行向量資料。

做一次累加。

x1=zeros(size(x));

for i=1:size(x1,2)

x1(i)=sum(x(1:i));

end%x1的均值數列。

z1=zeros(size(x));

for i=1:size(x1,2)-1

z1(i+1)=;

endyn=x(2:size(x,2))';

b=-z1(2:size(z1,2))';

b(:,2)=1;

u=inv((b'*b))*b'*yn;

a=u(1);

b=u(2);

**x2=zeros(1,n);

x2(1)=x(1);

for i=1:n-1

x2(1+i)=(x(1)-b/a)*exp(-a*i)+b/a;

endx2=[0 x2];

y=diff(x2);

用法:假設數列1 2 3 4 6 為已知資料,你要**後面3項,那麼儲存上面的**後輸入命令:

gm11([1 2 3 4 6 ,10)

5樓:第三滴血

clear;

clc;x0=[127627 128453 129227 129988 130756];

n=length(x0);

lamda=x0(1:n-1)./x0(2:n)range=minmax(lamda)

x1=cumsum(x0)

for i=2:n

z(i)=;

endb=[-z(2:n)',ones(n-1,1)];

y=x0(2:n)';

u=b\yx=dsolve('dx+a*x=b','x(0)=x0');

x=subs(x,,)

yuce1=subs(x,'t',[0:n-1]);

digits(6),y=vpa(x) %為提高**精度,先計算**值,再顯示微分方程的解。

yuce=[x0(1),diff(yuce1)]epsilon=(x0-yuce)./x0 %計算殘差delta=abs(epsilon./x0) %計算相對誤差rho=1-( 計算級比偏差值。

急求灰色**模型gm(1,n)matlab源**!

6樓:匿名使用者

急求灰色**模型gm(1,n)matlab源**!

躲營家北釐叔舅第京兼。

灰色系統**gm(1,1)模型的matlab源**

7樓:祿美俞訪天

把資料改一下應該就行了。

clcclear

x0=[1827

道路交通事故財產損失原始資料。

fori=2:8

x1(1)=x0(1);

x1(i)=x1(i-1)+x0(i);

endx1對原始資料進行累加。

fori=1:7

yn(i)=x0(i+1);

endynfori=1:7

b(i)=(;

endbfori=1:7

b(i,1)=b(i);

b(i,2)=1;

endb求得資料矩陣b

c=inv(b'*b)*b'*yn'

a=c(1,1)

u=c(2,1)

通過最小二乘法求出變數a,u

fort=1:8

x11(t)=u/a+exp(-a*t)*(u+1827*a)/exp(-a)/a;

endx11

建立時間響應函式。

q0=x1-x11

fori=2:8

q1(1)=q0(1);

q1(i)=q1(i-1)+q0(i);

endq1殘差第一次累加。

fori=2:8

q2(1)=q1(1);

q2(i)=q2(i-1)+q1(i);

endq2殘差第二次累加。

fori=1:7

b1(i)=(;

endb1fori=1:7

b1(i,1)=b1(i);

b1(i,2)=1;

endb1c1=inv(b1'*b1)*b1'*yn'

a1=c1(1,1)

u1=c1(2,1)

symstq2t=u1/a1+exp(-a1*t)*(u1+q2(1)*a1)/exp(-a1)/a1;

d=diff(q2t)

對二次殘差求灰導。

x=u/a+exp(-a*t)*(u+1827*a)/exp(-a)/a+d;

c=1:5;

y=subs(x,t,c)

z2012=y(3)-y

2)z2013=y(4)-y

3)z2014=y(5)-y

4)%求得2012-2014年道路交通事故財產損失**值。

急求gm(1,n)的matlab程式!!!!

8樓:落地生根發芽

難度太大了 我不會。

n 1an n 1bn都發散則,級數 n 1, ,an 和 n 1, ,bn 都發散 則級數 n 1, ,an bn 發散,為什麼

不一定吧,如果第一個級數裡邊,an n,第二個級數裡邊bn n,這樣級數當然都是發散的,但是每一項是an bn 0這樣的級數顯然不發散。例子不太好。一般的講,應該是考慮an和bn的絕對值,這樣有絕對發散性。級數 cn求和 如果每一項都比已知發散的級數絕對值大,那cn也必然發散。這個可能是叫柯西比較法...

1 求lim n1 2 n 3 n 1 n 20sin x 26求解,需過程

lim n 1 2 n 3 n 1 n e lim n ln 1 2 n 3 n n e lim n 2 n ln2 3 n ln3 1 2 n 3 n e lim n 2 3 n ln2 ln3 1 3 n 2 3 n 1 e 0 ln2 ln3 0 0 1 e ln3 3 0 sin x 2 d...

n 1是什麼意思, n 1 n 表示什麼呢?

懂梗菌 n 1中的n指的是每工作1年支付1個月工資的經濟補償。所謂的 1,是指用人單位解除勞動合同但是未提前30天書面通知勞動者時,額外支付的1個月工資。 n 當n為偶數時等於n n 2 n 4 2當n為奇數時等於n n 2 n 4 3回答者 necing 中級經理 八級 2009 10 23 14...