SSB(滤波法法)调制解调Matlab.m文件源码

%----------------输入参数

%m：信号

%ts:系统时域采样间隔、

%df:所需的采样频率

%fs:系统采样频率

%----------------输出（返回）函数

%M：傅里叶变换后的频谱序列

%m:输入信号参与傅里叶变换后对应序列，需要注意的是，该序列与输入信号m的区别，其长度是不一样的，输入的m长度不一定是2的整数次幂，而傅里叶变换要求输入信号长度为2的整数次幂，

%故傅里叶变换前需要对m信号进行补零操作，其长度有所增加，估输出参数中的m为补零后的输入信号，其长度与输入参数m不一样，但与M,f的长度是一样的，

%并且，其与时间序列t所对应的序列m(1:length(t))与输入参数中的m是一致的。

%df1:返回的频率分辨率

%f:与M相对应的频率序列

[M,m,df1] = fftseq(m,ts,df);

f = [0:df1:df1*(length(m)-1)]-fs/2;

M = M / fs;

End

%低通滤波器

function [H,f] = lp_f(n,f_cutoff,df1,fs,p)

%低通滤波器函数 输入设计的滤波器参数，产生低通滤波器频率特性函数H和频率向量f %------------输入参数

%n 低通滤波器的输入信号长度

%f_cutoff 低通滤波器截止频率

%df1 频率分辨率

%fs 抽样频率

%p 滤波器振幅

%------------输出（返回）函数

%H 低通滤波器频率响应

%f 频率向量

% 滤波器设计

n_cutoff = floor(f_cutoff/df1);

f = [0:df1:df1*(n-1)]-fs/2;

H = zeros(size(f));

H(1:n_cutoff) = p * ones(1,n_cutoff);

H(length(f)-n_cutoff+1:length(f)) = p * ones(1,n_cutoff);

end

%带通滤波器