
v=ones(10,1);
k=2*v;
f=3*v;
z=4*v;
[A,B,C,D]=zhuiganfa(v,k,f,z)
Q=diag(B)+diag(A(2:end),-1)+diag(C(1:end-1),1);
function [A,B,C,D]=zhuiganfa(v,k,f,z)
f1=cumsum(f);
A=v+[0;f1(1:end-1)]-v(1);
v1=[v(2:end);0];
f2=[0;f1(2:end)];
B=-v.*k-v1-f2+v(1);
C=v.*k;
D=-f.*z;
end
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删