自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+

钟工的算法小仓库

多数人在抱怨中成长,少数人在成功之前默默耕耘

  • 博客(153)
  • 收藏
  • 关注

原创 [机器学习]机器学习中的小笔记

岭回归,Lasso回归在拟合回归中回出现过拟合现象,表现为拟合方差的过大,训练系数W的过大,加入正则化L1 L2项训练使得方差减少,偏差增大;Lasso-L1回归特性:某些系数可以为0,变为稀疏特征,但其损失方程不可导,a系数增大收敛比岭回归快;Ridge-L2回归特性:系数变小但不为0,方程可导,收敛没Lasso回归快;kappa系数...

2019-09-04 22:18:26 181

原创 [学术资源]常用论文网址汇集

免费下载论文https://sci-hub.twPython区 ----------------------------------------------------------------Python 配置:如何使用pip安装:https://jingyan.baidu.com/article/466506583fec5af549e5f825.htmlSci...

2019-08-12 16:11:33 272

原创 [信号处理]独立成分分析ICA - 基础

为什么随机变量X和Y不相关却不一定独立?为什么随机变量X和Y不相关却不一定独立? - 孙逍的回答说的是对的,我举个更直观的例子吧(X,Y) 均匀分布在单位元 x^2 + y^2 = 1上X和Y的(线性)相关系数是0。为什么呢?直观来说,因为是个圆,如果你画一条线性回归的线,线的斜率是正的还是负的都不合适,因为是对称的。数学上因为E(X|Y) = E(Y|X) = 0所以E(X) = E(Y) = 0,而且E(XY) = E[E(XY|X)] = E[X E(Y|X) ] = 0

2021-10-04 20:53:59 327

原创 睡眠监测中的一些术语和指标(AASM)

低通气指数apnea hypopnea index (AHI) 呼吸紊乱指数respiratory disturbance index (RDI) 眼动电图electrooculography,EOG睡眠呼吸暂停低通气综合征sleep apnoea hypopnea syndrome中矢状下颌运动midsagittal jaw movem...

2019-11-20 13:49:54 2360

原创 [C&C++]C语言如何动态分配二维数组

#include <stdlib.h>int main(){ // 为二维数组分配3行 int **a = (int **)malloc(3 * sizeof(int *)); // 为每行分配5列 for(int i = 0; i < 3; i++) { a[i] = (int *)mall...

2019-09-25 11:42:10 312

原创 [Matlab]Matlab中改变y轴刻度为特定字符

% 画图 h = figure(1); clf(h); ax1 = subplot(2,1,1);bar(BCG_data);title("BCG实时睡眠分期算法结果"); set(gca,'linewidth',0.5,'fontsize',15,'fontname','宋体'); set(gca,'ytick',[0 1 3]); set(g...

2019-09-18 10:22:21 1435

原创 [数字信号处理]信号幅度概率特征算法设计 原创

随机信号处理的一种特殊的方法是通过记录信号幅度的出现概率,来生成特征。如下图所示为一种基本设计:概率分布的范围设定为该段信号的最大最小值,然后分成若干个区间进行分别求和并进行归一化,这种方法对于平稳信号比较适用,但对于幅度变化较大的随机信号,得出的概率模型可能波动很大,并不实际;如下图为一种改进的幅度特征设计方法:该方法首先统计出信号中数据点的中位值,这样做的好处是获得了该信号...

2019-08-17 08:29:26 895

原创 [Python]时间字符串转标准时间戳

# 字符类型的时间tss1 = '2013-10-10 23:40:00'# 转为时间数组timeArray = time.strptime(tss1, "%Y-%m-%d %H:%M:%S")print timeArray # timeArray可以调用tm_year等print timeArray.tm_year # 2013# 转为时间戳timeStamp ...

2019-08-13 20:30:00 466

原创 [原创]数字信号处理中的卷积网络深度学习猜想(3) - 信号动态周期包络算法

基本原理请参考这篇文章:https://blog.csdn.net/xeonmm1/article/details/95735261这次我们希望通过全尺度卷积算法实现一个随机信号的包络计算,包络计算的难题是:如果使用幅度逼近,则可能缺失时域,如果使用时域逼近,则可能缺失幅度,之所以出现这种情况是因为随机信号的信号成分不是固定的,里面包含了从低到高的频率,和不同的幅度信号,为了尽可能缓...

2019-07-23 17:50:16 544

原创 [原创]数字信号处理中的卷积网络深度学习猜想(2) - 信号主周期估计算法 - 一种快速实现

转载请注明出处,Writer:Zhonlihao钟硕 @ 20190713https://blog.csdn.net/xeonmm1全尺度卷积算法的时间复杂度全尺度卷积虽然有非常强的细节提取能力,但其是以大量增加运算量为代价的,下面的伪代码给出了其基本的运算过程:for i = 尺度遍历 n for j = 时域遍历 k for m = 特征遍...

2019-07-14 09:25:33 332

原创 [原创]数字信号处理中的卷积网络深度学习猜想(1) - 信号主周期估计算法

转载请注明出处,Writer:Zhonlihao钟硕 @ 20190713***注意!本文只是猜想,请勿过于认真要作者给出与深度学习的必要联系***https://blog.csdn.net/xeonmm1引言 这是一篇计算机学与信号处理跨界借鉴的文章,涉及电信算法领域,以下内容尽可能不贴图讲解,相关人士老司机可以继续阅读,适逢深度学习的大热,我花了点基本的时间了解深度...

2019-07-13 13:52:59 1367

原创 [个人生活]我生命中超低概率会发生的事件(但发生了)集合

1.在南山公司认识一个同事,某周末在深圳宝安壹方城星巴克看见他买东西喝,但我认为相遇的概率太低,没有过去打招呼,后来周一回到公司,的确是他本人,他也看到我在看他,但大家没有认真打招呼!-_-#2.在深圳打球认识一本地校友华农妹子,她居然认识我在广州的高中同学!3.20190716 四个广州高中同学在深圳相聚吃饭,其中3人在深圳工作,包括我。...

2019-06-14 14:17:31 242

原创 [深度学习]卷积神经网络CNN - 一张图理解滤波器/层数/深度

1.单个卷积滤波器的尺寸一般为3x3 或 5x5;如上图所示的kernel2.滤波器的层数与图像输入的色彩层有关 如上图所示层数为33.滤波器的深度是指布设多少个卷积特征提取器,越深的深度有利于多角度提取相关识别的特征,如上图所示深度为2最终经过合成后特征的图层的个数与上一级的卷积深度一致,从输入的图层个数 转化为特征层的个数;...

2019-05-04 12:24:46 4780

转载 [深度学习]Tensorflow教学 - 手写数字识别例程MNIST代码实现

import tensorflow as tfimport osimport urllibimport tensorflow.examples.tutorials.mnist.input_data as input_data# from tensorflow.examples.tutorials.mnist import input_data# 手动从Lecun乐村大牛的网站下载图片...

2019-05-04 11:48:15 527

原创 [最优化方法]动态规划 - 矩阵链乘法乘积问题

详细请看算法导论教程// MatrixChainOrder.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include <stdlib.h> #include "stdio.h"void MatrixChainOrder(int*p, int n);void printMatrixChainOrder(int** s_mat,...

2019-05-02 21:52:00 539

原创 [最优化方法]动态规划 - 钢条切割问题

# len 0 1 2 3 4 5 6 7 8 9 10Length_Price_Table = [0,1,5,8,9,10,17,17,20,24,30]# 子问题最优解的表Solve_Table = [0,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1]Cut_Table = [[0,-1,-1,-1,-1,-1,-1,-1,-...

2019-05-02 21:45:00 1592

原创 [数值分析]开方根算法 - 牛顿迭代法/二分法/割线法

// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include "math.h"#include "stdlib.h"#define err_limit 0.0001 //注意float的精度,不能要求太高,否则float无法区分#define iter_limit 1000000float ...

2019-05-02 21:33:31 1736

原创 [C&C++]Visual Studio生成DLL动态链接库的注意事项

因为不是所有平台都有DLL的依赖文件所以编译的时候这样选

2019-04-30 09:54:27 392

原创 [C&C++]本地内存管理代码

呵呵 不会告诉你怎么用~skr~#include <stdio.h>#include <stdlib.h>#include <stdint.h>#include <string.h>// config#define mem_debug printf//#define MD_MEM_STATS//#define MD_D...

2019-04-26 14:28:54 174

原创 [深度学习]Tensorflow实现Label一定范围内损失为0的代码实现

介绍:在很多应用中我们不一定要回归到某个Label [a],而是a的某个范围[a-b,a+b]内,均可认为训练成功,下面给出实现代码import tensorflow as tfimport numpy as np# zhonglihao 2019.04# 正负1范围没有lossoutput = [1.0,2.0,3.0,4.0,5.0,6.0]label = [3.5,3.5...

2019-04-23 10:01:29 421

原创 [深度学习]使用Tensorflow实现一个简单异或逻辑的学习

import tensorflow as tfimport numpy as np# zhonglihao 2019.04# 训练集x_train = [[0.0,0.0],[0.0,1.0],[1.0,0.0],[1.0,1.0]]y_train = [[1],[0],[0],[1]]# 网络参数 输入向量为2 隐藏层为4个cell 输出1/0单个节点# zhongliha...

2019-04-22 16:01:00 476

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(8)

第二章:异步并行编程使用parfeval,你可以在后台跑一个特定的功能函数并取出结果而不需要等待其完整跑完。在很多的情况下,提早从循环中释放代码继续运行下去是有利的,例如,对于一个优化的程序来说,就可以提前结束循环当确认了结果是足够好的,你可以在一个或者多个并行池上完成这个功能,请使用parfeval或者parfevalOnall。同样地,这对实时绘制图形也是有利的。请注意,这跟parfor是...

2019-04-10 19:48:53 4619 4

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(7)

Arranged By Zhonglihao @ 2018**请确认Matlab安装时点选了并行计算工具箱并行池(parallel pool)的设置点击右下角中并行图标中的“Parallel preference”选项进入设置界面:在设置中可以选择开启多少个核心做并行计算,如下图所示开启的核心为4个,若填写的数量大于处理器所拥有的最大核心数,则以开启最大核心数为准。下面两个...

2019-04-10 19:44:47 3365

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(6)

Arranged By Zhonglihao @ 2018**请确认Matlab安装时点选了并行计算工具箱改善parfor的性能Matlab中parfor的性能受到很多因素的影响,包括:·数组的大小·创建数组的时间·核心处理数组所有部分的时间·循环次数的大小当数组很大的时候,我们可以创建parallel.pool.Constant字段来加...

2019-04-10 19:37:07 1693 2

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(5)

Arranged By Zhonglihao @ 2018**请确认Matlab安装时点选了并行计算工具箱变量问题观察下面三个parfor索引语法:parfor i = 0:0.2:1 parfor j = 1:2:11 parfor k = 12:-1:1 三句语法在parfor使用中均错误,第一句不能用小数作为递...

2019-04-10 19:32:30 2076

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(4)

嵌套的parfor循环和for循环如下面的代码所示:Code #4-11.parfori=1:102.MyFun(i)3.end4.5.functionMyFun(i)6.parforj=1:57. ...8.end9.end代码#4-1为在循环里面调用一个函数,同时该函数也使用parfor做并行计算,然而,在语句p...

2019-04-06 15:20:54 3338

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(3)

Arranged By Zhonglihao @ 2018**请确认Matlab安装时点选了并行计算工具箱确定parfor中的循环代码片计算是独立的Code #3-11.clearA2.fori=1:83.A(i)=i;4.end5.A6.7.%改写成parfor并行8.clearA9.parfori=1:81...

2019-04-06 14:30:18 3791

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(2)

在改为parfor的过程中代码容易出现的问题Code #2-11.forx=0:0.1:12.fork=2:103.x(k)=x(k-1)+k;4.end5.x6.end7.8.parforx=0:0.1:19.parfork=2:1010.x(k)=x(k-1)+k;11.end12.x...

2019-04-06 14:24:47 1233

原创 [并行计算]Matlab并行计算工具箱(Parallel Computing Toolbox)官方文档教程中文版(1)

Arranged By Zhonglihao @ 2018**请确认Matlab安装时点选了并行计算工具箱第一章:parfor循环并行计算parfor循环介绍parfor循环是Matlab并行计算工具箱用于并行计算单循环的工具,并行度与CPU所具备的核心数有关,一般来说,经过处理的for循环都能够改为parfor循环以提高并行处理速度。Matlab中的CPU并行处理功能如下图所示进...

2019-04-05 21:27:47 22212

原创 [神经网络]Matlab神经网络工具箱提取网络参数进行实现的注意事项

1.将输入特征归一化到浮点[0,1]区间;2.譬如有4个特征[x1 x2 x3 x4]必须对角点进行封边以防止Matlab归一化时将幅度改变,所以应加入训练集两个点[0 0 0 0] =正例,[1 1 1 1] =反例这样;3.注意训练器图示的激活函数,有变换的为sigmoid或tanh网络,一条直线为线性输出以0为简单的分界点;4.从net.IW /net.LW /ne...

2019-03-05 09:52:15 1647 1

原创 [数值分析]插值 - 待定系数法

clc;clear all;close all;n = 1:1:10000;sig = sin(n/100);k = 0;for i = 1:4:10000 k = k + 1; x_arr = n(i:i+3); y_arr = sig(i:i+3); arr_new_y(k) = MUCInterpolation(x_arr,y_arr,4,mea...

2019-03-04 07:40:40 1877

原创 [数值分析]插值 - 牛顿插值法

clc;clear all;close all;n = 1:1:10000;sig = sin(n/100);k = 0;for i = 1:4:10000 k = k + 1; x_arr = n(i:i+3); y_arr = sig(i:i+3); arr_new_y(k) = NewtonInterpolation(x_arr,y_arr,4,...

2019-03-04 07:36:22 896

原创 [数值分析]插值 - 拉格朗日插值法

clc;clear all;close all;n = 1:1:10000;sig = sin(n/100);k = 0;for i = 1:4:10000 k = k + 1; x_arr = n(i:i+3); y_arr = sig(i:i+3); arr_new_y(k) = LagrangeInterpolation(x_arr,y_arr,...

2019-03-04 07:34:04 497

原创 [数值分析]数值积分 - 梯形/辛普森Simpson/柯特斯Newton-cotes方法

clear;close all;clc;format long% matlab 中精度较高的积分f = @(x) sin(x)/(x);x = 0.0000001:0.0000001:2*pi;y = sin(x)./(x);figure(1);plot(x,y);grid on;% 采用高采样率进行积分 [0,1]sum = 0;for i ...

2019-03-04 07:26:50 3597 2

原创 [数值分析]线性方程组求解:列主元高斯消元法

// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include &lt;stdlib.h&gt; #include "stdio.h"void MatrixPrint(double* arr, const int row, const int col);double* MatrixSolve(dou...

2019-03-03 08:55:50 833

原创 [数值分析]线性方程组求解:Jacobi迭代法

// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include &lt;stdlib.h&gt; #include "stdio.h"#include "math.h"void MatrixPrint(double* arr, const int row, const int col);doubl...

2019-03-03 08:52:29 989

原创 [数值分析]线性方程组求解:Gauss-Seidel迭代法

// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include &lt;stdlib.h&gt; #include "stdio.h"#include "math.h"void MatrixPrint(double* arr, const int row, const int col);doubl...

2019-03-03 08:47:46 1614

原创 [线性代数]矩阵求逆

// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。//#include "stdafx.h"#include &lt;stdlib.h&gt; #include "stdio.h"void MatrixPrint(double* arr, const int row, const int col);double* MatrixInverse(d...

2019-03-03 08:39:54 1045

原创 [机器学习]一维聚类:最大类计算

% 一维数组聚类,找出规模最大的类clc;clear all;close all;% 待聚类的一维数据 data_1 = 20*ones(15,1) + randn(15,1);data_2 = 5*ones(5,1) + randn(5,1);data_3 = 15*ones(5,1) + randn(5,1);data_4 = 50*ones(2,1) + randn(2,1)...

2019-03-03 08:36:08 1776

原创 [C&C++]多线程基础编程

1.基本的阻塞型实现// multi_thread.cpp : 此文件包含 "main" 函数。程序执行将在此处开始并结束。//#include "pch.h"#include &lt;iostream&gt;#include &lt;thread&gt;#include &lt;Windows.h&gt;using namespace std;void thread01...

2019-03-02 22:55:44 157

空空如也

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除