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...