function o=impulse_response(AR)
%計算加權沖擊響應h(n)；用于通過p(z)，計算碼書域目標信號及矩陣fi
%輸出---h,加權沖擊響應，40個樣點    

changshu=[1 0.75 0.75^2 0.75^3 0.75^4 0.75^5 0.75^6 0.75^7 0.75^8 0.75^9 0.75^10];
[o,t]=impz(1,AR.*changshu,40);
%round(o*2^12)'