NotAnyMore_HCWY-CSDN博客

原创编程中的小问题记录

Matlabtopkrows 可以返回最大（小）的k个值。不知道是否用了最小堆排序。

2020-11-10 17:34:22 122

前前后后有两篇英文论文被一区TOP期刊接受了。写的过程中没有给老外润色（好像是武大现在不给报销了，所以老板就默认让我直接投了），审稿人和编辑都没有指出有任何语言问题，因此我想自己算是攒了一些经验了，因此记录下来，希望对有缘人有用，也让自己温习。但我并不打算按逻辑整理，只是想到啥写啥。1、写作之前应该掌握一些最基础的规则，推荐书籍：（1）《MLA科研论文写作规范》（第七版），标点的用法，论文的结构都涉及到，非常专业。（2）An outline of scientific writing for r

2020-11-13 16:55:20 193

原创需要更加努力才配得上过去的自己

很久没有好好写东西了，突然想写点什么，因为昨晚又感觉自己打了鸡血，想要与前半个学期有些堕落的我作一些告别。此刻耳机放着bob dylan 的歌，很和适宜。要写也只能写在这里了，因为没熟悉的人（其实是没人）看。哈哈。　　自己高中是学文科的，每当有feel的时候，总能写个几段，可是现在每天都是论文，代码加矩阵的日子，似乎让自己变得不会抒发感情了，除了哈哈似乎啥子东西都写不出，嗯，再试试往下凭感觉

2017-01-31 01:41:46 451 1

原创 SIFT+RANSAC做图像矫正

做图像矫正时使用了一下sift算法，尽管sift确实很牛，但还是会出现一些误匹配，直接计算两张影像的单应矩阵会出现很大误差，因此可以在计算时使用RANSAC算法在单应矩阵的约束下剔除误匹配，并计算单应矩阵，基本思路为：1、进行sift匹配。2、RANSAC计算单应矩阵。3、使用插值函数进行重采样。基于OpenCV的c++代码：#include <iostream>#include <fstrea

2016-01-08 09:38:24 6431 1

原创高斯牛顿迭代求解非线性回归问题

解决的问题：已知函数：，数据X=[1.0,2.0,3.0,4.0,5.0]'; Y=[15.0,21.656,27.928,34,39.944]'; 求参数 b c d。该问题属于非线性回归，不能直接使用最小二乘或梯度下降法求解。但可以使用高斯牛顿法迭代求解。首先将b c d 作为变量泰勒展开得。。。。。。。懒得写了。。。。。代码copy:clcclear%训练数据x=[1.0,2

2016-01-07 23:16:12 3252

原创线性函数的协方差传播

开始着手学习平差啦，第一个要接触的就是协方差传播，而书上最先介绍的便是线性函数的协方差传播，鉴于这个地方十分重要，自己来推一遍。一切都源于统计里的方差计算公式:（1）以及协方差公式：（2） 1、首先是推导单一线性函数传播情况即，其中Z是t维向量，K为t*n的系数矩阵，X为n*t的矩阵，K0为t*1的矩阵。根据公式（1）有

2015-05-27 23:41:18 2054 1

原创通过投影矩阵推导最小二乘法

最小二乘用于参数拟合，说它是统计中最重要的公式都不为过，其线性推导如下。对于一个公式，其中A为已知数据组成的列空间，x为待估计参数，b是x的线性组合。然而，在很多情况下，由于误差或其它原因，使得b不再A的列空间中，这会造成方程无解。于是必须将b投影到A的列空间。其投影的误差为, 其中A为投影矩阵，当的值最小时矩阵e应该位于与A垂直的空间的列空间（这里画几条线段就很好理解），所以有。将式子展

2015-05-26 17:45:41 1588

原创基于特征值的斐波那契数列求解

最近在听Gilbert Strang的线代公开课，再讲矩阵的对角化展开的应用时，提到了利用特征值求斐波那契数列的例子。该方求指定位置的数列值时，其计算量大大缩小。理论基础为，将两端各自相乘会发现，这意味着求矩阵的幂只需计算其特征值。接着便需要把斐波那契数列转换到线代的问题上来。斐波那契数列的简单表达就是，这个数列第0项为0，第一项是1，从第三项开始，每一项都等于前两项之和。代数表达即为，令

2015-05-24 13:25:17 3176

转载 Hash算法

C H A P T E RT W OFUNDAMENTALS Most of the advancements throughout the history of computers have been because of particular problems which required solving. In this chapter, we

2014-06-18 13:05:42 483

转载堆

http://www.benfrederickson.com/2013/10/10/heap-visualization.html

2014-06-18 12:43:55 381

原创排序算法（未完成）

#include#include#include//插入排序int* sort(int* array){ for(int i=1;array[i]!=NULL;i++) { int tmp=array[i]; for(int j=i;j>0&&tmp<array[j-1];j--) array[j]=array[j-1]; array[j]=t

2014-04-01 21:24:52 541

原创平衡树实现（未完成）

基本操作1、清空整棵树。2、查找元素。3、查找最大值和最小值。4、插入。5、删除。插入和删除明天再做#include#includetypedef int elemtype;////////////////////////////////////////////////////////////////////////////定义avl树结构str

2014-03-27 18:45:14 503

原创链表中递归查找元素，非递归查找元素以及基数排序（未完成）josephus问题（未完成）

#ifndef CIRCLECHAIN_H#define CIRCLECHAIN_Hstruct NODE;struct LIST;typedef NODE* node;typedef LIST* list;typedef int elemtype;//单链表struct NODE{ int n; node next;};struct LIST{ node tail;

2014-03-26 22:45:42 725

原创链表实现多项式相加相乘

csdn太差还是自己操作不对代码插入就会卡死#include#include#include"string.h"typedef int elemType;struct Node{elemType coefficient;int highPower;Node *next;}; struct List{Node *heade

2014-03-26 22:41:40 1136

原创简单的单链表实现 c

只实现创建链表头，插入操作，与打印操作，其它的慢慢加#include#include#include"string.h"typedef int elemType;struct Node{ elemType element; Node *next;}; struct List{ Node *header; Node *tail;};typedef Node* n

2014-03-24 19:27:25 541

原创数据结构与算法C语言描述第三章练习

习题3.2 3.4 3.5 3.5的实现是错误的，后面再来改吧，因为部分习题使用的是相同的数据结构，所以就把他们都放在一起，这里的基本数据结构是循环链表头文件包含基本的数据结构实现以及习题里所要求实现的函数

2014-03-20 17:00:13 850

原创 3.6例题实现方法

与书上的实现相比，自己的实现增加了自动判断字符串大小的部分从而动态分配数组#include #include #includeusing namespace std;void reverse(char s[]);/*itoa function,change integer into string,savd in s[];*/char* itoa(int n);

2014-03-20 16:37:11 500

原创 itoa函数：将数字n转换为字符串并保存到s中 3.6例题实现

#include #include #includeusing namespace std;void reverse(char s[]);/*itoa function,change integer into string,savd in s[];*/char* itoa(int n);int main(){ //test and I alway

2014-03-20 16:35:02 844

原创扩展速记符号字符串练习3-3

因为对string类不是很熟悉，不知道如何在string后面直接插入一个字符，所以不得已用指针而为之，缺陷就是一开始必须显示规定扩展字符串的大小。后面有时间了用string来实现以下。不得不得不夸一下K.R这本书，给出的例子或是习题虽然用到的语法规则很简单，但都需要思考一番才能做出来，适合我这种需要提高的菜鸟。 #include #include #includeus

2014-03-20 12:09:23 992

原创 reverse函数倒置字符串各个字符位置 3.5

#include#include #include /* reverse function,reverse every char in the sting s*/void reverse(char s[]){ //remember the strlen() function int j=0; int tmp; printf("\nlll"); //the f

2014-03-20 10:41:15 568

原创 atoi函数：将s转换为整形数

#include#include /* atoi function:change string to integer**//*if there is space in the head,skip it*//*if it is "+" or "-" in the head,read it*/int atoi(char s[]){ int i=0;//counter

2014-03-20 10:14:32 901

原创指针二维数组的传递与函数参数

#include"stdlib.h"float* ReturnMinMax(float *array[],int maxrow,int maxcol);int main(){ float array[3][3]={8,6,2,5,7,3,4,8,3}; float* p[3]; p[0]=array[0]; p[1]=array[1]; p[2]=array

2014-03-12 16:42:09 641

转载动态分配数组

int *a; int N; scanf("%d", &N); a = (int *) malloc(N * sizeof(int)); .... free(a);

2014-02-02 13:58:20 659

原创结合matlab与arcpy的土壤样点聚类与成图

matlab 进行均值聚类与分类for i=1:1clcformat long g%均值聚类opts = statset('Display','final');[Idx,Ctrs,SumD,D] = kmeans(youjizhi,20,'Replicates',2000,'Options',opts, 'Start','uniform');points=[Idx youjizhi

2014-01-17 19:58:05 1736

原创获取arcgis克里格插值后的每个点的插值结果

最近在arcgis克里格插值后需要得到插值后的矩阵。找了很久才找到工具箱中的calculateZvalue，但是一次只能计算一个点，因此需要用脚本去提取。import arcpyarcpy.CheckOutExtension("GeoStats")total=0points=[]for width in range(37578564,37636464,50):

2013-11-27 23:28:27 5709 2

原创 c++search记录

1、c++如何按空格拆分字符串std::vector split(std:: string str,std:: string pattern) { std:: string::size_type pos; std::vector result; str+=pattern; // 扩展字符串以方便操作 int size=str.size

2013-11-27 23:18:49 991

wd171027的专栏