%產生白噪聲（normal）
nt=randn（11000）;
figureplot（nt）title（‘白噪聲‘）;
xlabel（‘n‘）;
ylabel（‘x（n）‘）;

%s（t）
w0=1000000;
a=1;
k=10000000000;
t=linspace（-0.0001?0.0001?1000）;
st=a*cos（w0*t+0.5*k*t.*t）;
figureplot（tst）title（‘信號‘）;
xlabel（‘t（s）‘）;
ylabel（‘s（t）‘）;

%x（t）
xt=nt+st;
figure?plot（t?xt）?title（‘帶噪聲信號‘）;
xlabel（‘t（s）‘）;
ylabel（‘x（t）‘）;

%y（t）
stt1=fft（st）;
stt2=conj（stt1）;%共軛
st2=ifft（stt2）;
N=size（st22）;
for?n=1:N;
????ht（n+1000）=st2（n）;
end;
ht111=ht（11001:2000）;
yt=conv（ht111xt）;%卷積
t1=10^（-6）*（-999:999）/1000;
figure?plot（t1yt）?title（‘輸出信號‘）;
xlabel（‘t（s）‘）;
ylabel（‘y（t）‘）;

%H（W）
HW=ifft（ht111）;
f=10^7*（1:1000）/1000;
figure?plot（fHW）?title（‘傳遞函數‘）;
xlabel（‘f（Hz）‘）;
ylabel（‘H（W）‘）;

%Y（W）
YW=ifft（yt）;
f=10^7*（-999:999）/1000;
Pyy=YW.*?conj（YW）?/?1000;
figure?plot（f?Pyy）?title（‘輸出頻譜‘）;
xlabel（‘f（Hz）‘）;
ylabel（‘Pyy‘）;

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

?????文件????????928??2007-05-23?02:39??pipeilvbo.m

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

?????????????????1146????????????????????2