博主大毕设关于数字下变频(DDC)的CUDA实现,预计工期比较长,所以留下一些文字记录. 主要分为两部分工作,Matlab仿真部分和CUDA实现. 由于很久没有仿真了,所以先用一个简单的AM调制仿真练手.具体代码如下: 代码都是基础的所以就不解释了(环境matlab2016) clc; fm=; %信号频率 fc=; %载波频率 fs=; %抽样频率 Am=; A=; N=; K=N-; n=:N-; t=(:/fs:K/fs); yt=Am*cos(*pi*fm*t); figure() su…

引用来自:http://blog.csdn.net/sinwel/article/details/8115673 %仿真参数中的含义 % Ts 表示间隔Ts时间采样,这个越小越接近连续信号,而实际上不可能真正的仿真连续信号. % df 表示变化后频率分量的最小间隔,即频谱中相邻的点间隔频率是多大,你也可以叫它 频率分辨力% 为什么df = fs/N-1呢?% 因为采样频率 fs=1/Ts,表示的意义是连续信号用fs去采样,得到了% 离散的信号,因此实际上的连续信号在这里并没有给出,因为MATLA…

要用到FIR滤波器和抽样器下面研究这两个的Matlab实现: Fir滤波器: matlab上fir滤波器的关键字是fir1 在command窗口输入help fir1出现帮助文档: >> help fir1 fir1 FIR filter design using the window method.                                              fir滤波器使用窗函数法设计 B = fir1(N,Wn) designs an N'th order…

在工程实际中遇到的信号经常伴有噪声,为了消除或减弱噪声,提取有用信号,必须进行滤波,能实现滤波功能的系统称为滤波器.严格地讲,滤波器可以定义为对已知的激励提供规定响应的系统,响应的要求可以在时域或频域内给定.滤波器的功能主要由其系统函数决定.滤波技术在通信测量.数据采集.控制系统及电力系统等领域得到了广泛的应用.滤波器的设计在这些领域中必不可少, 甚至是至关重要的环节, 其优劣将直接影响整个系统的性能.滤波器设计的关键是找到一个稳定可实现的逼近函数, 其频率特性满足所需滤波器的所有要求.虽然目前…

clf; fs=32;N=1024;   %采样频率和数据点数 n=0:N-1; t=n/fs;   %时间序列 x=1*sin(t); %信号 y=fft(x,N);    %对信号进行快速Fourier变换 mag=abs(y);     %求得Fourier变换后的振幅 f=n*fs/N;    %频率序列 subplot(2,1,1),plot(t, x);   %绘出随频率变化的振幅 xlabel('频率/Hz'); ylabel('振幅');title('N=128');grid o…

命令和一些特殊的变量 who:查看有哪些变量1. whos:可以查看变量的大小.字节和类型等资料. clear:如果单独使用则是清空所有命令,若后面跟着一个变量名称则为删除此变量. clc:清空命令行窗口. Inf:表示无穷大. eps:2.2004e-016 NaN:not a number pi:3.1415925··· 查看关键字可使用命令:iskeyword Numeric Display ‘Format’ Format:可设置输出格式.例:>>format long https://…

http://blog.sina.com.cn/s/blog_68f3a4510100qvp1.html 注:转载请注明出处——by author. 我们知道Fourier分析是信号处理里很重要的技术,matlab提供了强大的信号处理能力,但是有一些细节部分需要我们注意. 记信号f(t)的起始时间为t_start, 终止时间为t_end, 采样周期为t_s, 可以计算信号的持续时间Duration为 t_end – t_start, 信号离散化造成的采样点数 N = Duration/t_s +…

视频来源:https://www.bilibili.com/video/av51932171?t=628. 博文来源:https://ww2.mathworks.cn/help/matlab/ref/fft.html?searchHighlight=fft&s_tid=doc_srchtitle 视频来源很好的解释了: 1 .傅里叶变换过程,经过傅里叶变化得到了,频率w,振幅a0,相位角φ: 2. 傅里叶变换 主要应用领域: 声音, 图像处理: 博文则很好的解释了: 1.  傅里叶变换在matl…

一. 读写图像文件 1. imread imread函数用于读入各种图像文件,如:a=imread('e:/w01.tif') 注:计算机E盘上要有w01相应的.tif文件. 2. imwrite imwrite函数用于写入图像文件,如:imwrite(a,'e:/w02.tif',’tif’) 3. imfinfo imfinfo函数用于读取图像文件的有关信息,如:imfinfo('e:/w01.tif') 二. 图像的显示 1. image image函数是MATLAB提供的最原始的图像显示…

matlab图像处理 转自:http://www.cnblogs.com/lovebay/p/5094146.html 1. 图像和图像数据 缺省情况下,MATLAB将图像中的数据存储为双精度类型(double),64位浮点数,所需存储量很大:MATLAB还支持另一种类型无符号整型(uint8),即图像矩阵中每个数据占用1个字节. 在使用MATLAB工具箱时,一定要注意函数所要求的参数类型.另外,uint8与double两种类型数据的值域不同,编程需注意值域转换. 从uint8到double的转…

转自博客:http://blog.163.com/fei_lai_feng/blog/static/9289962200991713415422/ 一. 读写图像文件 1. imread imread函数用于读入各种图像文件,如:a=imread('e:\w01.tif') 注:计算机E盘上要有w01相应的.tif文件. 2. imwrite imwrite函数用于写入图像文件,如:imwrite(a,'e:\w02.tif',’tif’) 3. imfinfo imfinfo函数用于读取图像文…

前面一篇随笔叙述了关于MATLAB中FFT_HDL_Optimzed模块定点(IEEE754单精度float格式)二进制与十进制转换实现,有需要的往前翻一下即可.这一篇随笔将叙述 FFT HDL Optimized 模块实现 8192 burst Radix 2 结构的定点 FFT (定点格式依照 IEEE-754 single-precision format,小数为 23 位二进制精度),并通过 Xilinx LogiCORE FFT MEX Function interface 进行单精度…

(转于它处,仅供参考) 1.. Ctrl+C 中断正在执行的操作 如果程序不小心进入死循环,或者计算时间太长,可以在命令窗口中使用Ctrl+c来中断.MATLAB这时可能正疲于应付,响应会有些滞后. 2. figure命令新建一个绘图窗口 figure 可以打开一个空的绘图窗口,接下的绘图命令可以将图画在它里面,而不会覆盖以前的绘图窗口.当有多个figure窗口时,在命令窗口中执行如Plot等命令将覆盖当前figure窗口中的对象.所谓的当前figure窗口,也就是最后一次查看的窗口(可以用命令…

原文地址:matlab并行计算,大家共同学习吧,涉及到大规模数据量处理的时候还是效果很好的 今天搞了一下matlab的并行计算,效果好的出乎我的意料. 本来CPU就是双核,不过以前一直注重算法,没注意并行计算的问题.今天为了在8核的dell服务器上跑程序才专门看了一下.本身写的程序就很容易实现并行化,因为beamline之间并没有考虑相互作用.等于可以拆成n个线程并行,要是有550核的话,估计1ms就算完了... 先转下网上找到的资料. 一.Matlab并行计算原理梗概 Matlab的并行计算实…

1.FFT算法概要: FFT(Fast Fourier Transformation)是离散傅氏变换(DFT)的快速算法.即为快速傅氏变换.它是根据离散傅氏变换的奇.偶.虚.实等特性,对离散傅立叶变换的算法进行改进获得的. 2.FFT算法原理: 离散傅里叶变换DFT公式: FFT算法(Butterfly算法) 设x(n)为N项的复数序列,由DFT变换,任一X(m)的计算都需要N次复数乘法和N-1次复数加法,而一次复数乘法等于四次实数乘法和两次实数加法,一次复数加法等于两次实数加法,即使把一次复数…

matlab矩阵内存预分配就意味着,划定一个固定的内存块,各数据可直接按"行.列指数"存放到对应的元素中.若矩阵中不预配置内存.则随着"行.列指数"的变大.MATLAB就必须不断地为矩阵找到新的"空的内存",从而导致"建造矩阵"的速度大大下降. Hilbert矩阵的数学描写叙述 Hilbert矩阵是著名的"坏条件"矩阵,当中其第(i,j)元素的表达式是a(i,j)=1/(i+j-1). 代码比較: % 20…

一.Matlab并行计算原理梗概 Matlab的并行计算实质还是从主从结构的分布式计算.当你初始化Matlab并行计算环境时,你最初的Matlab进行自动成为主节点,同时初始化多个(具体个数手动设定,详见下文)Matlab计算子节点.Parfor的作用就是让这些子节点同时运行Parfor语句段中的代码.Parfor运行之初,主节点会将Parfor循环程序之外变量传递给计算子节点.字节点运算过程时互补干扰,运算完毕,则应该有相应代码将各自子节点得到的结果组合到同一个数组变量中,并返回到Matlab…

Roberts边缘检测算子:根据一对互相垂直方向上的差分可用来计算梯度的原理,采用对角线方向相邻两像素之差. 小波变换的方法比较适用于展现夹带在正常信号中的瞬间反常现象,具有方向敏感性.所以可以边缘检测. Roberts边缘检测算子: clear; I=imread('D:\文件及下载相关\图片\gray2.png'); I=rgb2gray(I); grayPic=mat2gray(I); [m,n]=size(grayPic); newGrayPic=grayPic; robertsNum=…

本来CPU就是双核,不过以前一直注重算法,没注意并行计算的问题.今天为了在8核的dell服务器上跑程序才专门看了一下.本身写的程序就很容易实现并行化,因为beamline之间并没有考虑相互作用.等于可以拆成n个线程并行,要是有550核的话,估计1ms就算完了... 先转下网上找到的资料. 一.Matlab并行计算原理梗概 Matlab的并行计算实质还是主从结构的分布式计算.当你初始化Matlab并行计算环境时,你最初的Matlab进程自动成为主节点,同时初始化多个(具体个数手动设定,详见下文)M…

首先,在写这篇博客之前,我还没有完全学会FFT. 先把会的部分打好,加深一下记忆(也可以说是做笔记吧). 初三了,还不会FFT,要退役喽-- 多项式乘法 点开这篇博客之前,你就应该知道,FFT是用来求多项式乘法的. 什么是多项式,什么是多项式乘法? 不讲.初一内容. 如果要求多项式乘法,有一个非常显然的做法,就是暴力. 时间复杂度是O(N2)O(N^2)O(N2)的,很朴素. 然而FFT这个东西可以将其复杂度优化到O(Nlg⁡N)O(N\lg N)O(NlgN). 点值表示法 对s于一个多项式A…

转自:http://silencethinking.blog.163.com/blog/static/911490562008928105813169/ AX=B或XA=B在MATLAB中,求解线性方程组时,主要采用前面章节介绍的除法运算符"/"和"\".如: X=A\B表示求矩阵方程AX＝B的解: X＝B/A表示矩阵方程XA=B的解. 对方程组X＝A\B,要求A和B用相同的行数,X和B有相同的列数,它的行数等于矩阵A的列数,方程X＝B/A同理. 如果矩阵A不是方阵…

R的优点:免费,开源,体积小.缺点:对大文本处理差,另外一个也在于开源,package如果出错,烦死你.当你跑比较大的simulation,对效率有要求的时候,有时还是不得不用C,这可能是10小时和10分钟的差别,毫不夸张.SAS流行于公司,R流行于研究机构和大学数据分析不是单纯的靠软件来做的,需要很好的数学基础. 统计学工具各有千秋.https://englianhu.wordpress.com/statistics/学了R,可以免去学spss,matalab,ucinet等等众多的软件,可以…

01 概述 Greedy Randomized Adaptive Search,贪婪随机自适应搜索(GRAS),是组合优化问题中的多起点元启发式算法,在算法的每次迭代中,主要由两个阶段组成:构造(construction)和局部搜索( local search). 构造(construction)阶段主要用于生成一个可行解,而后该初始可行解会被放进局部搜索进行邻域搜索,直到找到一个局部最优解为止. 02 整体框架 如上面所说,其实整一个算法的框架相对于其他算法来说还算比较简单明了,大家可以先看以…

题目 题目大意 也懒得解释题目大意了-- 正解 正解居然是\(FFT\)? 不要看题目的那个式子这么长,也不要在那个式子上下手. 其实我们会发现,不同的\((x_i-x_j,y_i-y_j,z_i-z_j)\)并不多. 如果我们求出每个三元组的出现次数,后面的就好做了. 那怎么求呢? 祭出我们的大杀器--\(FFT\). 考虑只有一个维怎么做.设两个多项式分别为\(A\)和\(B\). 对于\(x_i\),就在\(A\)的\(x_i\)这一位上的系数加一: 对于\(x_j\),就在\(B\)的\…

教学目标:1.使学生了解JAVA课程的性质.定位.作用:为什么要学习JAVA?让学生知道如何学好JAVA: 教学内容: 一.        问几个问题 1.             你们到这里来干什么的? 来学习JAVA程序设计 为什么要来学习JAVA呢? 找个好工作,拿到高薪水 2.             怎么样才能达到你们的目标呢? 有的同学说学好JAVA就可以了,但是如何才能学好呢? 学好JAVA的几个关键 1.             认认真真上课 2.             按时按…

FFT和功率谱估计 用Fourier变换求取信号的功率谱---周期图法 clf; Fs=1000; N=256;Nfft=256;%数据的长度和FFT所用的数据长度 n=0:N-1;t=n/Fs;%采用的时间序列 xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N); Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);%Fourier振幅谱平方的平均值,并转化为dB f=(0:length(Pxx)-1)*Fs/length(Pxx…

实验的过程中,经常需要对所采集的数据进行频谱分析,软件的选择对计算速度影响挺大的.我在实验过程中,通常使用Origin7.5来进行快速傅里叶变换,因为方便快捷,计算之后,绘出来的图也容易编辑.但是当数据容量太大,达到100M大小,这时候使用Origin7.5进行快速傅里叶变换,运算速度非常慢,甚至运算不出来. 对大容量的数据进行快速傅里叶变换,我使用Matlab,运算速度比Origin7.5快很多.但是使用Matlab进行FFT时,需要进行一些小的处理,才能使运算结果与使用Origin7.5进行…

参考文章:http://www.cnblogs.com/adgk07/p/9314892.html 首先根据他这个代码和我之前手上已经拥有的那个代码,编写了一个适合自己的代码. 首先模仿他的代码,测试成功. 思路: 短时傅里叶变换,其实还是傅里叶变换,只不过把一段长信号按信号长度(nsc).重叠点数(nov)重新采样. % 结合之前两个版本的stft,实现自己的周期图法,力求通俗易懂,代码分明.% 该代码写的时候是按照输入信号为实数的思路写的,在每个片段fft时进行前一半行的转置存储.后续代码思…

普通调幅波AM的频谱,大信号包络检波频谱分析 u(t)=Ucm(1+macos t)cos ct ma称为调幅系数 它的频谱由载波,上下边频组成 , 包络检波中二极管截去负半周再用电容低通滤波,可以得到基带信号,那么,截去负半周后的AM信号必定包含基带信号的频谱.我们可以通过matlab来验证. %已知基带信号为1hz,载波为64hz,调制系数ma=0.3,采样频率1024hz,FFT变换区间N为2048 clear; fs=; f=; %1hz基带信号 fc=; %64hz载波 t=:/f…

一:FFT原理 1.1 DFT计算 在一个周期内的离散傅里叶级数(DFS)变换定义为离散傅里叶变换(DFT). \[\begin{cases} X(k) = \sum_{n=0}^{N-1}x(n)W_N^{kn}, & 0 \le k \le {N-1} \\ x(n) = \frac{1}{N} \sum_{k=0}^{N-1}X(k)W_N^{-kn}, & 0 \le n \le {N-1} \\ \end{cases} \] 其中,\(W_N = e^{-j\frac{2\pi}…