clear;

sigma=input（‘白噪聲標準差設定為??‘）;
r=input（‘生成數據周期數（1—3）??‘）;

%%%%%%%%%%%%%%生成a=1P=6的M序列%%%%%%%%%%%%%%%%%
M=[0?0?1?0?1?0?0];
Np=2^6-1;
for?k=1:252
?????if?M（7）==0?
????????u（k）=1;end;
?????if?M（7）==1??
????????u（k）=-1;end;
?????for?i=7:-1:2
????????M（i）=M（i-1）;
?????end
?????M（1）=M（6）+M（7）;
?????if?M（1）==2?
????????M（1）=0;end;
end
u=u-mean（u）;
%%%%%%%%%%%%%%%生成高斯白噪聲%%%%%%%%%%%%%%%%%%%%
A=179;??xi=11;??M=2^15;
for?k=1:252????????
????ksai=0;
????for?i=1:12
????????xi=A*xi;
????????xi=mod（xiM）;
????????ksai=ksai+（xi/M）;
????end
v（k）=sigma*（ksai-6）;
end

%%%%%%%%%%%%%%%%%%理論輸出數據%%%%%%%%%%%%%%%%%%%%%%
k=120;?T1=8.3;?T2=6.2;?T0=1;
K=k/（T1*T2）;
E1=exp（-T0/T1）;
E2=exp（-T0/T2）;
x（1）=0;y（1）=0;
for?k=2:252
????x（k）=E1*x（k-1）+T1*K*（1-E1）*u（k-1）+T1*K*（T1*（E1-1）+T0）*（u（k）-u（k-1））/T0;
????y（k）=E2*y（k-1）+T2*（1-E2）*x（k-1）+T2*（T2*（E2-1）+T0）*（x（k）-x（k-1））/T0;
end
z=y+v;

%%%%%%%%%%%%%%%%%%計算互相關函數%%%%%%%%%%%%%%%%%%
for?k=1:Np
????R（k）=0;
for?i=（Np+1）:（r+1）*Np
????R（k）=R（k）+u（i-k）*z（i）;
end
Rmz（k）=R（k）/（r*Np）;
end
%%%%%%%%%%%%%%%%%%脈沖響應估計值%%%%%%%%%%%%%%%%%%
for?k=1:Np
g（k）=（Rmz（k）-Rmz（Np-1））*Np/（Np+1）;
end

%%%%%%%%%%%%%%%%%%脈沖響應理論值%%%%%%%%%%%%%%%%%%
K=120;
for?k=1:Np
g0（k）=K*（exp（-k/T1）-exp（-k/T2））/（T1-T2）;
end
%%%%%%%%%%%%%%%%%%脈沖響應估計誤差%%%%%%%%%%%%%%%%%%%
G=g0-g;
GG=G.*G;
gg=g0.*g0;
E=sqrt（sum（GG）/sum（gg））
figure
plot（1:Npg‘r*‘）
hold?on
plot（1:Npg0‘.‘）
legend（‘脈沖響應估計曲線‘‘脈沖響應理論曲線‘）

%%%%%%%%%%%%%%%%%%%%%%計算噪聲方差%%%%%%%%%%%%%%%%%
a=0;
for?i=1:length（v）
????a=a+v（i）;
end
a=a/252;
varv=0;
for?i=1:252
????varv=varv+（v（i）-a）^2;
end
varv=varv/252;
%%%%%%%%%%%%%%%%%%%計算過程輸出數據方差%%%%%%%%%%%%%%%%
b=0;
for?i=1:length（y）
????b=b+y（i）;
end
b=b/252;
vary=0;
for?i=1:252
????vary=vary+（y（i）-b）^2;
end
vary=vary/252;
%%%%%%%%%%%%%%%%%%%%計算過程信噪比%%%%%%%%%%%%%%%%%%
xinzaobi=sqrt（varv/vary）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????177583??2011-10-17?15:31??實驗報告.docx

?????文件???????1990??2011-04-25?20:37??main.m

-----------?---------??----------?-----??----

???????????????179573????????????????????2