[转帖]laplas到z转换

suffer

1. 
   
2. for i=1:N+1,
   
3. for j=1:N+1,
   
4. V(i,j)=0;
   
5. for k=0:j-1,
   
6. V(i,j)=V(i,j)+(-1)^(i-1-k)*C(j-1-k,N-(i-1))*C(k,i-1);
   
7. end
   
8. end
   
9. end
   
10. V
    
11. 
    
12. Az0=0;
    
13. for i=1:N+1,
    
14. Az0=Az0+At(N+1-(i-1))*(2/T)^(i-1)*V(i,N+1);
    
15. end
    
16. K0=1/Az0
    
17. 
    
18. for j=1:N+1,
    
19. Bz(N+2-j)=0;
    
20. for i=1:M+1,
    
21. Bz(N+2-j)=Bz(N+2-j)+K0*Bt(M+2-i)*(2/T)^(i-1)*V(i,j);
    
22. end
    
23. end
    
24. Bz
    
25. for j=1:N,
    
26. Az(N+2-j)=0;
    
27. for i=1:M+1,
    
28. Az(N+2-j)=Az(N+2-j)+K0*At(N+2-i)*(2/T)^(i-1)*V(i,j);
    
29. end
    
30. end
    
31. Az(1)=1;
    
32. Az

复制代码

[ 本帖最后由 suffer 于 2006-10-9 20:43 编辑 ]