紫卓执守-CSDN博客

原创 matlab实现小波图像去噪

本文主要内容1.简单置零操作实现小波域去噪2.基于简单机器学习思想的小波去噪下面是本次内容的全部代码代码后面会有解释如果你需要直接使用下面的代码，别忘了带走几个脚本中要调用到的函数clcclearnoiseAmp=2.1;%噪声的幅度sourceImage=imread('colorWoman.jpg');%读取图片noiseFreqCut=0.42...

2019-09-05 19:57:28 12315 5

原创 VGA编程接口讲解

1VGA接口和管脚定义的介绍VGA是一种常见的显示器接口，接口大概是下面这样想必大家之前都见过。上面有15个针脚，通过给这15个针脚发送信号脉冲就可以控制显示屏显示我们需要的图像。下面仅仅提一下简单的几个本文需要用到的针脚定义1 VGA-R 表示红色的亮度2 VGA-G 表示绿色的亮度3 VGA-B 表示蓝色的亮度RGB三个亮度电平组合起来就是我们最终看到的颜色，每...

2018-12-14 19:44:43 8976 2

原创用verilog HDL实现数字基带信号的2FSK调制

2FSK的介绍可以参考：https://blog.csdn.net/qq_39148922/article/details/84337730下面介绍verilog HDL的源代码module FSK( input clk,//时钟信号 output reg [7:0]sigOut,//输出已调信号 output reg [7:0]carryWave1,//1对应的载波幅度 ...

2018-11-25 09:56:32 10050 20

感觉eigen使用<<实现对矩阵的赋值用起来特别的方便，所以也想在自己做的类里面实现这个功能。自己根据功能实现出来了一个版本，目前还没有添加异常检查，仅仅实现了基本的功能，也不知道eigen底层是不是这么实现的，但总之是聊胜于无了。其实比较关键的一个操作就是在类里面设置了一个私有类，这个私有类充当一个缓冲区的功能，并且重载了逗号操作符，从而实现了功能。#include<iostream>#include<vector>class MyVector{priv

2021-09-20 15:36:12 823

原创 StandardScaler中fit_transform和transform的区别

在使用StandardScaler的时候需要先新建一个它的对象from sklearn.preprocessing import StandardScalerstdScaler=StandardScaler()fit_transform是对数据做归一化，归一化后的数据的每1列的均值为0,方差为2.例如下面这样的一个数据[[1,2,0],[-1,-2,4]]对它进行fit_transform操作后会得到[[ 1. 1. -1.] [-1. -1. 1.]]这个时候StandardS

2021-08-06 13:52:51 1690

原创 cuda核函数传入device函数指针的使用例子

cuda的global函数里面可以调用__device__函数，在有特殊需要的时候，还可以把__device__函数作为参数传入到一个__global__函数中在cuda里面不能像c++那样简单地传入函数的指针，需要在传入前对函数的指针做一些包装。例如typedef double(*funcFormat)(int,char);这里面double表示函数的返回值，int,char是函数的参数列表，所有满足这种格式的函数都可以用这种函数类型指代。上面的funcFormat是一种函数类型，在这里可以把

2021-07-08 22:41:41 3271

原创编译cuda代码时出现unsupported GNU version! gcc versions later than 8 are not supported!

ubuntu20.04系统中，正常使用nvcc编译cuda代码用得还是很方便快乐的。但有一天想写一个像模像样的工程的时候，尝试用c++调用cuda代码的时候，却出现了下面这样的报错。In file included from /usr/include/cuda_runtime.h:83, from <command-line>:/usr/include/crt/host_config.h:138:2: error: #error -- unsupporte

2021-06-14 18:42:36 3883

原创用cuda实现对较大内存数据的并行处理

刚刚学cuda的时候以为调用核函数的时候可以开辟任意大小的线程，所以当自己尝试用cuda实现卷积运算的时候，调用核函数的时候用了这样的写法cudaConv<<<1,memSize>>>(gpuMem);实际上cuda每个block的线程数是有上限的。下面这个知乎答案说明了怎样获取cuda可开辟的最大线程数。https://www.zhihu.com/question/462801801/answer/1925184931当需要开辟多于maxThreadsPerBl

2021-06-11 19:01:54 430

原创 cuda共享内存的作用和功能验证

在cuda中定义共享内存需要使用关键字__shared__在一个核函数里面声明一个共享内存，其它的线程也可以访问它。共享内存是不能被初始化的，只能先定义再赋值。下面是验证共享内存功能的一个例子。#include<iostream>#define MEM_SIZE 10__global__ void shareVisit(int* gpuMem){ //定义共享变量 __shared__ int shareValue; //当前位置的线程编号 int

2021-05-31 20:33:12 578

原创 Cmakelists.txt添加.h和.cpp文件

文件目录结构是这样代码都是最基本的代码,仅仅是个小的democmakelists.txt里面的内容#cmake的版本要求cmake_minimum_required(VERSION 3.6)#工程的名称project(algorTest)#添加可执行文件add_executable(main main.cpp)#第2个选项必须要有,可以选PUBLIC|PRIVATE|INTERFACEtarget_sources(main PUBLIC "solution.cpp")如果需

2020-12-11 21:30:09 8877 2

原创 vistual studio2015 cuda编程初体验 gpu速度的简单验证

做了一个小的cuda编程的例程代码，验证了一个gpu的运算速度。vs新建工程。新建工程后会有一个自带的例程代码，这个代码是直接调用cuda做运算的，这里对代码做了简单的修改，让cpu和gpu做相同的运算，然后比较它们两个的运算时间，最后时间比较结果如下：第2个是gpu的速度代码如下：...

2020-12-06 11:59:56 165

原创 cmakelists.txt 遍历文件夹下所有的文件名

cmakelists.txt内容如下,仅仅是一个参考.代码中的file指定目录下所有的符合命名规则的文件名写入到变量LIB_NAMES中,foreach遍历变量中的每个名字,message则是输出每个文件名代码中的命名规则表示的是所有以.lib结尾的文件....

2020-11-05 12:07:26 5264

原创 linux环境下matlab调用plot报错无法加载libmwosgserver.so(can‘t load libmwosgserver.so)

matlab R2016b的linux版本使用hold on或者plot等涉及到图形的命令的时候，会出现以下报错。Can’t reload ‘/usr/local/matlab-res.2015b/bin/glnxa64/libmwosgserver.so’在usr/lib下面搜索找到libstdc++.so.6复制到matlab下的bin/glnxa64目录下，并重启matlab从而解决问题。复制的命令（路径不一定完全相同）cp /usr/lib64/libstdc++.so.6 /usr/loc

2020-11-02 11:31:06 1156

原创 hp惠普暗影精灵6在ubuntu 18.04下无法打开wifi的问题

暗影精灵6使用的wifi模块的驱动对linux内核版本要求较高，ubuntu18.04默认的内核版本无法支持。所以才会在软件更新里面找不到对应的驱动。手动升级内核也是可以的，不过比较麻烦。直接重装u系统buntu 20.04从而解决问题，毕竟新电脑里面还没装什么东西，损失不大~...

2020-10-28 22:10:23 1000

原创递归实现八皇后问题

研究八皇后问题只是为了学习回溯，可惜还没有开始看回溯就想到了解法，也许这个解法远不如最优解，却磨灭了我学习回溯的动力。先贴上代码，以后可能会补充。#include<iostream>#include<vector>#define CHESS_LOAD 1//棋子落子的时候对势力范围上的操作#define CHESS_LEAVE -1//棋子离开的时候对势力范围...

2020-02-11 17:20:14 257

原创并不优雅的背包算法

问题描述：有若干物品，物品的质量和体积已知。有一个容积固定的背包，怎样在背包中装入质量尽可能大的物品代码如下，注释应该足够详细，只是代码风格不太优雅（捂脸）#include<iostream>#include<vector>using namespace std;//用来定义物品，需要提供质量和体积class BagObj{public: boo...

2020-01-18 13:36:06 153

原创基于opencv识别二值图像特定连通图形的算法介绍

本文主要介绍一些涉及到的编程思想。至于源代码，这个程序功能本身并没有什么现实意义，而且由于程序几乎没有调用外部接口，全是独立实现，程序复杂冗长，也与当今编程方法相背，所以不再展示。这是程序运行的效果：将两个正方形的描边颜色去掉就是原图。提取图像特征在进行图像识别工作的时候，最重要的就是找到要识别的目标区别于其它图形的特征。而且该特征对于识别目标具有稳定性....

2019-09-06 08:02:17 429

原创用matlab制作数字图像高频噪声

高频噪声的制作思想可参考https://blog.csdn.net/qq_39148922/article/details/86609629本文介绍数字图像高频噪声的制作，原始噪声用正态分布随机数生成源代码如下%数字图像高频噪声的制作%运行参数mLength=800;%噪声图像的宽度nLength=600;%噪声图像的高度ampi=0.46;%噪声的幅度freqCut=0...

2019-01-24 13:48:58 2822

转载 matlab dftfilt源代码

源代码如下function g=dftfilt(f,H)%DFTFILT Performs frequency domain filtering.% G=DFTFILT(F,H) filters F in the frequency domain using the% filter transfer function H. The output, G, is the filte...

2019-01-24 13:48:18 2760 2

转载 dftuv源代码

源代码如下function [U, V] = dftuv(M, N)%DFTUV Computes meshgrid frequency matrices.% [U, V] = DFTUV(M, N) computes meshgrid frequency matrices U and% V. U and V are useful for computing frequenc...

2019-01-24 13:39:25 1733

原创用matlab制作高频噪声的制作思想

要制作高频噪声，就要先搞清楚高频噪声的特点。在matlab中，可以中正态随机数比较方便地生成白噪声。白噪声的频谱如下图所示关于在matlab中研究数字图像频域问题的若干注意事项在如下链接文章中讲解↓https://blog.csdn.net/qq_39148922/article/details/86609686要得到高频噪声就要抑制频谱中的低频成分所以对低频的成分乘上...

2019-01-24 13:32:48 4327

原创在matlab中研究数字图像频域问题的注意事项

在进行数字图像去噪、压缩等过程中经常需要分析数字图像的频谱对于二维图像，分析数字图像的频谱主要工具是二维离散傅里叶变换在matlab中对应的主要函数是fft2（本文仅讨论"fft2(image)"这种调用形式) 对于一幅无噪声的，以低频为主的数字图像取傅里叶变换，再取绝对值，画出它的三维图像如下这样的图像即表现为频率的分布但我们在实际处理图像频谱的时...

2019-01-24 13:16:52 1228

原创数字图像中的频率概念

常见的数字图像是一个数值矩阵，所谓数字图像中频率的概念也就是二维矩阵中频率的概念。下图是在数字图像中呈现出来的一个由中心向外扩展的正弦波，哪个是高频哪个是低频一望而知图像中亮度突变的位置会表现出高频可以从下面的例子中感受一下图像中的高频部分由哪些因素产生一幅含有噪声的数字图像将图像通过数字低通滤波器后的结果大部分的噪声被滤除，带有信息的图像...

2019-01-22 06:58:37 1977

原创 matlab对非平稳一维信号的小波变换

matlab实验代码如下clccleart=linspace(0,2*pi,1024);%信号采样时间点signal=sin(10*t).*stepfun(t,pi);%实验信号%进行一维离散小波变换，获得近似分量和细节分量[sigWaveApp,sigWaveDet]=dwt(signal,'db4');%提取小波基的滤波器...

2019-01-17 14:42:14 8971 1

原创一维离散小波变换过程

小波变换的本质不过是一种数学变换在这里仅仅讨论小波的小波变换过程中对输入的信号进行了怎样的操作，尽量不涉及内部细节和数学原理 1.一维小波变换的输入变量是一个【1×n】的矩阵，你也可以把它理解为信号、函数等等2.进行离散小波变换需要预先指定两个滤波器，一个是高通滤波器、另一个是低通滤波器3.将输入的一维向量和滤波器的系统函数卷积得到两个卷积的结果4.得到的两个结果分别进行系...

2019-01-17 12:05:34 9900 5

原创下采样

下采样，通俗地讲就是隔几个数采样一个点下采样过程令y(n)=x(Dn) 例如一个离散序列本来长这样如果系数D=2，那么下采样之后得到的结果就是： ...

2019-01-15 12:06:02 1229

原创 verilog HDL编程驱动VGA接口显示图像代码全注释

关于VGA的编程接口可参考https://blog.csdn.net/qq_39148922/article/details/85005271 本代码驱动屏幕显示红绿蓝三道条纹，具体的代码解释可参考注释module vgaTest( clk,horSync,verSync,rForOut,gForOut,bForOut,reset); input clk;//输入的系统时钟...

2018-12-14 22:38:20 3148

原创 matlab对符号函数赋值——subs inline @

matlab里面的符号函数通常指的是用syms创建出来一个变量，然后对这些变量进行的运算。比如：syms xy=x.^2matlab可以比较方便地对一个符号函数求导或者积分，傅里叶变换拉普拉斯变换等。但是经过运算完成后有时候是需要把运算得到的结果转化回矩阵进行后续的操作的，这个时候就需要通过一个矩阵对符号函数赋值从而达到转化出来一个新的矩阵。subs这个函数常用的情况下有三个参...

2018-12-05 17:27:46 25866

原创 matlab数字图像处理初体验——图片的读取和操作

我们电脑上的以.jpg为后缀名的图片就是最常见的数字图像，这样的图像是可以被matlab读取的。i=imread('woman1.jpg');woman1.jpg是处在matlab当前工作环境的一张图片。运行之后可以发现它得到的是一个三维的矩阵。3层的矩阵每一层代表矩阵的一个颜色，在图像中叫做通道。iDouble=im2double(i);通过上面的代码可以把图片生...

2018-12-02 17:09:18 1976

原创数字带通传输系统——键控法调制

数字带通传输系统，是通信系统中用来传输数字信号的基本模型。而键控法是对数字信号进行调制的一种基本方法。数字基带信号数字基带信号是通信过程中包含信息的原始信号，通常用二进制数字来表示，例如01000101010101010101010101这样……然而如果直接用高电平表示1，低电平表示0，产生出来的信号几乎不可能通过天线发射出去。具体的物理原因是因为信号频率太低。因此需要某种手段让信...

2018-11-22 20:14:23 3740

原创基于matlab产生数字信号的时域波形——将数字信号当作特殊的模拟信号

在使用matlab的时候，可以很方便的表示二进制的序列或者其他形式的数字信号。用一个行向量表示即可。但是有时候需要对数字信号的时域波形进行某些操作，例如对时域信号进行调制。将正弦波和一个向量相乘是不现实的，因此本文讲说明如何将matlab中的向量恢复成细小的采样点组成的函数采样序列写贴上整个的代码%用来产生二进制码的波形%待转换的二进制序列x=round(rand(1,20))...

2018-11-22 16:43:13 6584

原创噪声系数与灵敏度

噪声系数是通信系统中的一个概念。一个信号经过了一个通信系统（通常是接收机），信噪比的相对变化程度就可以用噪声系数来衡量。具体来说，噪声系数指的是输入信噪比和输出信噪比的比值。且通常用分贝来表示有时候对接收机输出信号的信噪比有特定的要求。要达到一个特定的输出信噪比的最小的输入信号功率，就是这种情境下的灵敏度。灵敏度的表达式中有一个参数是没有确定的——输入信号的噪声功率。...

2018-11-15 19:41:10 8044 1

qq_39148922的博客