关闭

基于多项滤波的数字正交变换MATLAB仿真程序

标签: matlab
789人阅读 评论(1) 收藏 举报
分类:
function filter_emu( Num )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
%num_size = sizeof(Num);
%利用多项滤波器的分支特性获得两个有半个采样点时延的低通滤波器系数</span>
filt1 = zeros(1,8);
filt2 = zeros(1,8);
for m = 1:64
    if(mod(m,8) == 4)
        filt1(i) = Num(m);
    end
    if(mod(m,8) == 0)
        filt2(i) = Num(m);
        i = i + 1;
    end
end
%------生成窄带信号,中频150MHz,带宽不大于20MHz
%------仿真信号x(t)=a(t)*cos[2*pi*f0*t+phi(t)]
f0 = 1.5e8;     %中心频率
fs = 2e8;       %采样频率
N = 1600;       %取的样本点数
n = 0:N-1;      %取的样本序列
t = n/fs;       %获得以1/fs为时间间隔的采样序列
%a = 1+cos(2*pi*1000*t);      
                %获取a(t)的采样点
phi = 2*pi*2e6;
                %获取phi
xt = cos(2*pi*f0*t+phi*t);
                %生成窄带信号并获取其采样点
%------进行2倍抽取并混频
nt = 1:N/20;
xi = xt(2*nt).*((-1).^nt);  %获取同相分量
xq = xt(2*nt-1).*((-1).^nt);%获取正交分量

%------打印同相分量与正交分量
figure(1);
plot(nt,xi,'r');
hold on;
plot(nt,xq,'g');
grid on;
%------滤波
x1 = filter(filt1,1,xi);
x2 = filter(filt2,1,xq);
%------打印恢复后的同相正交分量
figure(2);
plot(nt,x1,'r');
hold on;
plot(nt,x2,'g');
grid on;

end
注:Num为生成的63阶FIR滤波器系数存放的数组,通过fdatool工具得到。
0
0

猜你在找
【直播】机器学习&数据挖掘7周实训--韦玮
【套餐】系统集成项目管理工程师顺利通关--徐朋
【直播】3小时掌握Docker最佳实战-徐西宁
【套餐】机器学习系列套餐(算法+实战)--唐宇迪
【直播】计算机视觉原理及实战--屈教授
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之矩阵--黄博士
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之凸优化--马博士
【套餐】Javascript 设计模式实战--曾亮
查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:4852次
    • 积分:96
    • 等级:
    • 排名:千里之外
    • 原创:5篇
    • 转载:0篇
    • 译文:0篇
    • 评论:1条
    文章存档
    最新评论