[转帖]laplas到z转换
for i=1:N+1,
for j=1:N+1,
V(i,j)=0;
for k=0:j-1,
V(i,j)=V(i,j)+(-1)^(i-1-k)*C(j-1-k,N-(i-1))*C(k,i-1);
end
end
end
V
Az0=0;
for i=1:N+1,
Az0=Az0+At(N+1-(i-1))*(2/T)^(i-1)*V(i,N+1);
end
K0=1/Az0
for j=1:N+1,
Bz(N+2-j)=0;
for i=1:M+1,
Bz(N+2-j)=Bz(N+2-j)+K0*Bt(M+2-i)*(2/T)^(i-1)*V(i,j);
end
end
Bz
for j=1:N,
Az(N+2-j)=0;
for i=1:M+1,
Az(N+2-j)=Az(N+2-j)+K0*At(N+2-i)*(2/T)^(i-1)*V(i,j);
end
end
Az(1)=1;
Az
[ 本帖最后由 suffer 于 2006-10-9 20:43 编辑 ]
页:
[1]