function [h]=Rayleigh_Clarke(fm,N,fs);
%*************************************
%參考《無線通信原理與應用》P135
% fm:多譜勒頻移
% N :采樣點數
% fs:輸出線譜的采樣頻率,為擴頻之前的符號速率
%*************************************
%頻率間隔
deltaf=2*fm/(2*N-1);
%譜線數
M_Comp=fs/deltaf;
M=2^(ceil(log2(M_Comp)));%ceil數是正無取整

%產生兩路高斯復噪聲
randn('state',sum(100*clock));
h1=randn(1,N).*exp(j*rand(1,N)*2*pi);
%使左右共軛對稱
h1=[conj(fliplr(h1)),h1];

randn('state',sum(100*clock)+500);
h2=randn(1,N).*exp(j*rand(1,N)*2*pi);
h2=[conj(fliplr(h2)),h2];

%生成濾波器幅度特性
n=0:N-2;
Sezf(n+1)=1.5./(pi*fm.*sqrt(1-((n+1/2).*deltaf./fm).^2));
Sezf(N)=2*Sezf(N-1)-Sezf(N-2);
Sezf=[fliplr(Sezf),Sezf];
%高斯噪聲通過多譜勒濾波器
h1=sqrt(Sezf).*h1;
h2=sqrt(Sezf).*h2;
%進行IFFT變換
h1=ifft(h1,M);
h2=ifft(h2,M);

%模平方,相加,開方
h=sqrt((abs(h1)).^2+(abs(h2)).^2);
%幅度歸一化
h=M*h./sum(h);