数字下变频matlab源程序

%main_OFDM.m

%一个相对完整的OFDM通信系统的仿真设计，包括编码，调制，IFFT， %上下变频，高斯信道建模，FFT，PAPR抑制，各种同步，解调和解码等模 %块，并统括系统性能的仿真验证了系统设计的可靠性。

clear all close all clc

%++++++++++++++++++++++++++全局变量++++++++++++++++++++++++++++++ % seq_num 表示当前帧是第几帧 % count_dds_up 上变频处的控制字的累加

% count_dds_down 下变频处的控制字的累加（整整） % count_dds_down_tmp 下变频处的控制字的累加（小数） % dingshi 定时同步的定位

% m_syn 记录定时同步中的自相关平台 global seq_num global count_dds_up global count_dds_down

global count_dds_down_tmp global dingshi global m_syn

%+

%main_OFDM.m

%一个相对完整的OFDM通信系统的仿真设计，包括编码，调制，IFFT， %上下变频，高斯信道建模，FFT，PAPR抑制，各种同步，解调和解码等模 %块，并统括系统性能的仿真验证了系统设计的可靠性。

clear all close all clc

%++++++++++++++++++++++++++全局变量++++++++++++++++++++++++++++++ % seq_num 表示当前帧是第几帧 % count_dds_up 上变频处的控制字的累加

% count_dds_down 下变频处的控制字的累加（整整） % count_dds_down_tmp 下变频处的控制字的累加（小数） % dingshi 定时同步的定位

% m_syn 记录定时同步中的自相关平台 global seq_num global count_dds_up global count_dds_down

global count_dds_down_tmp global dingshi global m_syn

%+

书中Matlab源程序 第1 章 绪论

【例1-1】有一名学生，期末有5门功课要考试，可用的复习时间有18小时。假定这五门课程分别是数学、英语、计算机基础、画法几何和专业概论。如果不复习直接参加考试，这五门功课预期的考试成绩分别为65分、60分、70分、60分和65分。复习以1小时为一单元，每增加1小时复习时间，各门功课考试成绩就有可能提高，每复习1小时各门功课考试成绩提高的分数分别为3分、4分、5分、4分和6分。问如何安排各门功课的复习时间可使平均成绩不低于80分，并且数学和英语成绩分别不低于70分和75分。

解：设分配在数学、英语、计算机基础、画法几何和专业概论这五门功课的复习时间分别为x1,x2,x3,x4,x5，则可列出如下的目标函数和限制条件为：

f(x)?x1?x2?x3?x4?x5 x1?x2?x3?x4?x5?18

(3x1?4x2?5x3?4x4?6x5?320)/5?80

3x1?65?70

本例具体程序如下：

%li_1_1 f=[1 1 1 1 1];

A=[1 1 1 1 1; -3 -4 -5 -4 -6; -3 0 0 0 0; 0 -4 0 0 0; 3 0

因为最近在做宽带数字信号的数字下变频，所以重点看了这一方面的论文。下文重点对IEEE上的一篇论文Real-time wideband DDC based on parallel architecture in synthetic instrument 做了个读书笔记，并对自己在这方面的工作做了一个小结。 主要内容：这篇论文主要解决宽带数字中频信号处理过程中的一些问题，并行处理结构可以有效的降低采样率，多通道并行处理则可以有效降低处理成本。该论文介绍了数字下变频技术，分析了多种数控振荡器导致输出相位的不同，并总结了减少信噪比的误差因素，提供了特殊工程设计的指数分布依据。

图1 DDC总体结构

原理：数字下变频器(DDC)是接收机A／D变换后，首先要完成的处理工作，一般的DDC由本地振荡器(NCO)、混频器、低通滤波器和抽取器组成.主要作用：其一是把中频信号变为零中频信号；其二是降低采样率。从频谱上看，数字下变频将A/D采样后信号从中频变换,到基带。这样的处理由两步完成：首先是将输入信号与正交载波相乘，然后进行数字滤波滤除不需要的频率分量。NCO，混频器，数字滤波器速率要等于采样率，采样率低于600MHz，很难实时的在FPGA中进行处理。

数字

clear all;

clc;

tic;

X=zeros(1000,4);Y=zeros(1000,4);U=zeros(1000,1);

M=3;

a=[28 18 74 74 70 72 60 36 12 18 14 90 78 24 54 62 98 36 38 32]; b=[42 50 34 6 18 98 50 40 4 20 78 36 20 52 6 60 14 58 88 54]; Aeq=zeros(20,80);

for j=1:20

for k=1:20

if j==k

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=21:40

if k==j+20

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=41:60

if k==j+40

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=61:80

if k==j+60

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

end

Aeq;

A=zeros(4,80);

for i=1:4

for k=1:80

if i==1

for k=1:20

A(i,k)=1;

end

end

if i==2

for k=21:40

A(i,k)

%*************************************************************************************

% This function pertains to the addition of AWGN with mean zero and % parameter 'variance' to an input signal. %

% AUTHOR: Wenbin Luo % DATE : 04/12/01 % % SYNOPSIS: y = awgn(x,var) % x ---> input signal % var ---> variance

% y ---> y = x + AWGN

%***********************************************************************************

function y = awgn(x,var) w = randn(1,lengt

电工专业必看

%简单潮流计算的牛顿拉夫逊程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写 %对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点 %编号，将变压器的串联阻抗置于低压侧处理。

%第三列为支路的串列阻抗参数。

%第四列为支路的对地导纳参数。

%第五列为含变压器支路的变压器的变比

%第六列为变压器是否含有变压器的参数，其中“1”为含有变压器，

%“0”为不含有变压器。

%B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。

%X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地 %参数。

n=input('请输入节点数:n=');

n1=input('请输入支路数:n1=');

isb=input('请输入平衡节点号:isb=');

pr=input('请输入误差精度:pr=');

B1=input('请输入支路参数:B1=');

B2=input('请输入节点参数:B2=');

X=input('节点号和对地参数:X=');

Y=zeros(n

%*************************************************************************************

% This function pertains to the addition of AWGN with mean zero and % parameter 'variance' to an input signal. %

% AUTHOR: Wenbin Luo % DATE : 04/12/01 % % SYNOPSIS: y = awgn(x,var) % x ---> input signal % var ---> variance

% y ---> y = x + AWGN

%***********************************************************************************

function y = awgn(x,var) w = randn(1,lengt

龙源期刊网 http://www.qikan.com.cn

数字下变频基于FPGA的软件设计与实现

作者：马纯

来源：《科技视界》2016年第14期

【摘 要】雷达的数字下变频功能主要是将天线接收到的中频回波信号通过A/D变换后，进行数字下变频处理，转为两路I/Q基带数据。本文主要设计了一种基于FPGA的数字下变频方法，通过对数字控制振荡器（NCO）及低通滤波器（FIR）的设计及实现，完成了对于不同频率本振信号的数字下变频处理。结果表明，基于该方案设计的数字下变频功能已在实际系统中得到应用。

【关键词】数字下变频（DDC）；FPGA；数字滤波器

【Abstract】The function of DDC is mainly to convert the signal received after A/D conversion by the antenna to Digital Down Convert（DDC）， finally into two I/Q baseband data. This paper mainly design a method about Digital Down Con

1、Newdon迭代法求解非线性方程

function [x k t]=NewdonToEquation(f,df,x0,eps) %牛顿迭代法解线性方程

%[x k t]=NewdonToEquation(f,df,x0,eps) %x:近似解 %k:迭代次数 %t:运算时间

%f:原函数，定义为内联函数 ?:函数的倒数，定义为内联函数 %x0:初始值 %eps:误差限 %

%应用举例：

%f=inline('x^3+4*x^2-10'); ?=inline('3*x^2+8*x');

%x=NewdonToEquation(f,df,1,0.5e-6) %[x k]=NewdonToEquation(f,df,1,0.5e-6) %[x k t]=NewdonToEquation(f,df,1,0.5e-6)

%函数的最后一个参数也可以不写。默认情况下，eps=0.5e-6 %[x k t]=NewdonToEquation(f,df,1)

if nargin==3 eps=\end tic; k=0; while 1

x=\ k=\

if abs(x-x0) < eps || k >30