xbianpin.m

来自「就是qpsk用matlab实现调制解调以及扩频接扩」· M 代码 · 共 38 行

function k2=xbianpin(R,ft,N,Ts);  
                   
fs=1/Ts;
T=0.001/N;                           %码元时间（宽度）
tdi=0:2*Ts:(N*T)-Ts; 
R1=R.*cos(2*pi*ft*tdi);

[b,a]=ellip(5,0.5,80,4e6/(fs/4),'low'); %得到滤波器的系数矩阵
[H,w]=freqz (b,a,512); %H 为滤波器的系统函数
sf=filter(b,a,R1); %用设计的滤波器过滤原始信号

S1=fft(sf,1024); %求出原始信号的傅立叶变换
w=(0:511)/512*(fs/4);
figure(20)            %画出经过低通滤波器后频谱
plot(w,abs(S1(1:512)'));
title('低通滤波器后的频谱' );
set(gcf,'color' ,'white' );
xlabel(' Frequency(Hz)' );
ylabel(' Mag.of frequency response' ); grid;

figure(21);
plot(w*fs/(2*pi),abs(H));
title('IIR 带通滤波器幅频响应' );
xlabel(' Frequency(Hz)' );
ylabel(' Magnitude of frequencyresponse' );


%==================================================
%设计低通滤波器
%===================================================
figure(4);
subplot(212);
plot(tdi,sf);
axis([0 0.001 -2 2]);
title('去载波后的s(t)图形');
xlabel('时间坐标                             t');
ylabel('amplitude');
end