Sumujingling-CSDN博客

转载使用 CUBLAS 库给矩阵运算提速

转自http://www.cnblogs.com/scut-fm/p/3756242.html前言　编写 CUDA 程序真心不是个简单的事儿，调试也不方便，很费时。那么有没有一些现成的 CUDA 库来调用呢？　　答案是有的，如 CUBLAS 就是 CUDA 专门用来解决线性代数运算的库。　　本文将大致介绍如何使用 CUBLAS 库，同时演示一个使用 CUBLAS 库进行矩阵乘法的例子。

2016-09-07 21:57:50 1347

原创 Matlab之宽带波束形成

波束形成的学习都从常规的开始。常规宽带波束形成原理上很简单。信号->各频率子带->窄带波束形成->相加。周末的时候场外求助了一种相移的DAS波束形成。果然高手和新手是不一样的。其中体现的矩阵运算思想十分值得学习！关键代码：Phase_shift=exp(2i * pi * Freq. * Tao); Freq 和 Tao 都是 8*data’length； Freq的每行相同，均为信号的

2016-06-13 16:02:35 12613 2

原创 CUDA之窄带常规波束形成

思路现在手上有了cuda的复数矩阵乘法和复数矩阵转置。理论上讲可以做一个简单的波束形成了。按照matlab之并行计算的思想把for循环都变成矩阵来做。复数矩阵定义typedef struct { int width; int height; int stride; double * real; double * imag;} Matrix;复数矩

2016-05-31 22:31:27 2224 1

原创 CUDA之矩阵乘法——复数

做好矩阵乘法和转置之后本来开心得不行的！准备上手做个最基本的波束形成了！突然发现希尔伯特变换完以后需要进行各种复数的运算…所以临时补写了一个复数乘法… 学着学着好像有点感觉了~！还是蛮有意思的。当然前提是能调试成功。用一句傅小姐的名言鼓励一下“只要心甘情愿任何事情都会变得简单！”。

2016-05-30 16:55:25 4125 1

原创 CUDA之学习资料

NYU Prof. Mohamed Zahran 的GPU 教学资料胡文美教授的CUDA讲座

2016-05-30 10:48:48 560

原创 CUDA之矩阵转置

在波束形成中除了矩阵乘法，还需要一个矩阵转置。其实用C语言来实现矩阵转置很简单，但是速度可能比较慢。参考博文：CUDA矩阵转置const int TILE_SIZE = 4; __global__ void TransposeKernel(const

2016-05-30 09:26:59 2432

原创 CUDA之矩阵乘法——非方阵计算

说明A矩阵为M * N，B矩阵为N * M，求C = A * B;

2016-05-26 18:12:27 2237

原创 MATLAB之并行算法学习

MATLAB的加速1.在占用时间较长的代码上进行算法的优化 2.打开MATLAB并行池使用parfor等运算 3.将循环转化为矩阵0526更新本来在做CUDA的矩阵乘法，在设计参数的时候突然想到matlab这边可以把一个循环转化成矩阵。也就是3的思想。波束形成的主要部分：data = hilbert(data);R1 = data'*data;for i = 1:length(theta)

2016-05-26 10:39:38 1200

原创 CUDA之矩阵乘法——TILE&sharedmemory

CUDA 矩阵乘法将输入数据分成很多个TILE使用shared memory进行并行计算

2016-05-25 10:12:02 3223

原创 CUDA之矩阵乘法——globalmemory

根据CUDA讲座调试矩阵乘法报错错误 17 error : no instance of overloaded function “cudaMalloc” matches the argument list E:\Niki\MVDR_BTR\MVDR_BTR\MatrixMulOnDevice.cu 16 修正：把CUDAcudaMalloc(&Nd, size);改成cudaMa

2016-05-12 08:17:45 1093

原创坐标系的转换

今天的任务是需要把一个以距离r为横坐标，以角度theta为纵坐标的彩色图像，转换为一个以theta和r为元素的扇形图。由于没有找到可以直接画出极坐标彩图的函数，所以最终需要给直角坐标系下的像素点赋值。两种思路，一种是用x,y来表示theta和r，以读取原图中的彩色像素值。即， r = sqrt(x^2+y^2); theta = atan(y/x

2016-05-05 08:51:47 880

原创水声成像之基本概念

（水）声学一直整得一知半解，遇到两个概念有点不太明白了。有选择性散射叫做散射，无选择性散射叫做漫射。百度百科：前向散射：散射时，如果相对于电子入射方向，散射角小于90°，则为前向散射，否则为背向散射。后向散射：在合成孔径雷达及散射计等遥感器中，所观测的散射波的方向是入射方向，这个方向上的散射就称作后向散射。我的理解：根据定义，那么，双/多基地声呐是前

2016-04-29 09:45:08 2102

原创 MFC炼狱之旅（二）画图像

teechart莫名其妙搞不出来（我造是我道行不够）所以我想转投opencv的怀抱噜！看别人的代码中，使用CvvImage类，由于是OpenCV的问题，新版本没有了CvvImage，查阅网上资料，只要自己把CvvImage.h和CvvImage.cpp加到工程中就行了。而我却找不到这两个文件，不过这里有。博客：http://blog.chinaunix.net/uid-2602

2016-03-31 11:14:04 508

原创 MFC炼狱之旅（一）文件读写及曲线图

1.文件读写读FILE *fd1;errno_t err1=fopen_s(&fd1,"E:\\DATA_20150831\\DATA1\\DATA1_1\\Record9.bin","rb");fread(beamdoa,sizeof(double),4*180,fd1);fclose(fd1);notice：在控制台工程里面用fopen打开就好好的，在MFC里面他就

2016-03-29 15:50:31 859

原创 Arrayfire学习笔记（三）常规窄带波束形成

经过两天的熟悉以及小伙伴的指导，蹒跚着开始了程序的实（测）现（试）有了一些常用的矩阵运算，基本上就可以实现一个煎蛋的波束形成代码了。在VS上写矩阵比matlab痛苦很多，matlab：R = signal_d'*signal_d; for i = 1:length(theta) a_s = exp(-jay*2*pi*f0*[0:N-1]

2016-03-17 13:15:14 3375

原创 Arrayfire学习笔记（二）开始滑翔吧！

（虽然写的是原创，但是好像我的功能就是提取我需要的信息然后摘录翻译下来？好方。）1.支持的数据类型：If no data type is specified, ArrayFire defaults to f32.b8 8-bit boolean values (bool)f32 real single-precision (float)c32 complex single-pr

2016-03-15 08:30:13 1929

原创 Arrayfire学习笔记（一）安装

1.安装 arrayfire的官网上有教程。按照步骤基本就可以完成安装。http://arrayfire.org/docs/installing.htm安装的过程就是很简单了，用户需要做的就是设置一个环境变量。即，In edit mode, append AF_PATH%/lib. NOTE: Ensure that there is a semi-colon separating AF_

2016-03-15 08:09:22 4894

原创 MATLAB之极坐标绘图

目标是要绘制一个二维的极坐标彩色图。输入参数有三个，一个是角度，一个是半径，一个是颜色。说到极坐标绘图，第一个想到的就是polar啦~那就先试试吧！1.polar绘图polar函数用来绘制极坐标图，调用格式为：polar（theta，rho，选项）其中，theta为极坐标极角，rho为极径，选项的内容和plot函数相似。例：绘制的极坐标图clccl

2016-03-14 10:17:24 26540 2

原创 MATLAB之波束形成学习

MathWorks官网上有一个Phased Array System Toolbox： http://cn.mathworks.com/help/phased/index.html其中有一项叫做Beamforming and Direction of Arrival Estimation： http://cn.mathworks.com/help/phased/beamforming-and

2016-03-03 09:30:03 21520 4

摸爬滚打琦琦子的博客