passbandfilter.m

完整的幅相调制信号可k阶非线性功率谱分析

function y=passbandfilter(fs,B,B3dB,fc)
%本函数设计带通滤波器
%函数形式：y=passbandfilter(fs,B,fc)
%参数说明：
%    y    输出滤波器抽头系数
%    fs   采样速率
%    B    信号带宽
%    B3dB    信号3dB带宽
%    fc   载波频率

%函数体
f=[fc-B/2 fc-B3dB/2 fc+B3dB/2 fc+B/2];
% rp = 1;          % Passband ripple
% rs = 80;         % Stopband ripple
a = [0 1 0];       % Desired amplitudes
% Compute deviations
dev = [0.00001 0.01 0.00001]; 
% dev = [(10^(rp/20)-1)/(10^(rp/20)+1)  10^(-rs/20)
% (10^(rp/20)-1)/(10^(rp/20)+1)]; 
[n,fo,ao,w] = firpmord(f,a,dev,fs);
y = firpm(n,fo,ao,w);
% freqz(y,1,1024,fs);
% title('Lowpass Filter Designed to Specifications');