%filter?tap:k
%data?length:N

g=100;??%統計仿真次數為g
N=1024;??%輸入信號抽樣點數N
k=128;???%時域抽頭LMS算法濾波器階數
pp=zeros（gN-k）;??%將每次獨立循環的誤差結果存于矩陣pp中，以便后面對其平均
u=0.0001;

for?q=1:g
????t=1:N;
????a=1;
????s=a*sin（0.05*pi*t）;?%input?signal
????figure（1）;
????subplot（311）
????plot（treal（s））;??%信號s的時域波形
????title（‘信號s時域波形‘）;
????xlabel（‘n‘）;
????ylabel（‘s‘）;
????axis（[0N-a-1a+1]）;
????xn=awgn（s5）?????%加入零均值高斯白噪聲，信噪比為3dB
????
????%設置初值
????y=zeros（1N）;???%output?signal
????y（1:k）=xn（1:k）;??%將輸入信號xn的前k個值作為輸出y的前k個值