按照通信原理书上的公式，写进去就行了
频域的　用fft就行了

FFT函数具体怎么用啊？我查了点资料，写了程序但是运行出来不知道波形对不对，请各位帮忙看下：
t=[0.0001:0.0001:0.30]
Fs=500;%采样频率
a=ones(1,500)
s=[1*a,-2*a,0*a,1*a,-2*a,0*a]
f=100*sqrt(34)
v=2*cos(2*pi*f*t)
figure(1)
plot(t,v)
v1=s.*v
figure(2)
plot(t,v1)
l_s=length(t);
n=ceil(l_s/2); 
Fs=500;%采样频率
ps=fft(v1,l_s); 
A=2*abs(ps)/l_s;   
A(1)=A(1)/2;  
tt=Fs*(0:n-1)/l_s; 
figure(3)
plot(tt,A(1:n));