- 博客(114)
- 收藏
- 关注
RSS订阅原创 自然数平方的倒数和为多少
简单数学，做起来不简单，需要用泰勒公式、求根变换和伟达定理！故1＋1/2²＋1/3²＋ …[在这里插入图片描述]由韦达定理,常数项为1时,根的倒数和＝一次项系数的相反数。这个首先是由欧拉推出来的,要用到泰勒公式,属于大学范围。故方程sin√y/√y＝0的根为π²,(2π)²,…而方程sinx＝0的根为0,±π,±2π,…＝0的根为π²,(2π)²,…即1/π²＋1/(2π)²＋…自然数平方的倒数和为多少？于是sinx/x＝1－x。
2023-07-07 19:57:35
393
原创 python的reshape和matlab的reshape，我们不一样
python的reshape和matlab的reshape，形似而神不同，千万当心
2022-07-27 18:11:40
626
1
原创 软件移植问题：Linux上好好程序到Windows就跑不通？
问题描述1 Linux代码跑到windows平台上一跑就崩溃，细看下发现是所有路径都寻找错误。于是定位到路径操作filename = path.split("/")[-1]原来是linux上路径分隔用斜杠/，而windows用反斜杠,无奈全改过来跑通。不再写死斜杠，用Os.sep替代最好。filename = path.split(os.sep)[-1]教训跨平台的路径分隔符，路径操作一定注意。os.sep替代写死的分隔符，斜杠/。...
2022-05-11 23:33:14
417
原创 【矩阵运算】矩阵的迹以及迹对矩阵求导总结
矩阵求导最终都是华为标量求导，迹就是最简单的衡量标量 一定要掌握迹求导原因总结技巧迹求导原因矩阵求导最终都是华为标量求导，迹就是最简单的衡量标量 一定要掌握总结Tr(AB)=Tr(BA)Tr(AB) =Tr(BA)Tr(AB)=Tr(BA)Tr(ABC)=Tr(BCA)=Tr(CAB)Tr(ABC) = Tr(BCA) = Tr(CAB)Tr(ABC)=Tr(BCA)=Tr(CAB)Tr(A)=Tr(A′)Tr(A) = Tr(A')Tr(A)=Tr(A′)d(Tr(XB)）=d(Tr(BX)
2022-03-24 17:40:06
2896
原创 musi算法使用简记
musicMUSIC是超分辨算法中最有工程使用价值的一种。原理很简洁明了，难点在工程应用。文章《Multiple Emitter Location and Signal Parameter Estimation》非常经典。公式就是那个经典的PMUSIC=1/(aH(θ)UNUNHa(θ))P_{MUSIC} = 1 / (a^H(\theta)U_N U_N^H a(\theta))PMUSIC=1/(aH(θ)UNUNHa(θ))其中a(θ)a(\theta)a(θ)是导向矢量， UNU_NU
2022-03-03 10:48:37
1023
原创 matlab强制转换向量或矩阵为行向量，列向量
有时候需要将向量或矩阵转为行向量，方便输出或操作，可以不用判断，利用(😃，所有元素“：”来自动获取拉平的值（注意matlab是按列优先取值）。免于判断烦恼。a = [1,2,3,4]a_row = a(:)'a2 = a'a2_row = a(:)'A1 = [1,2,3;4,5,6]A1_row = A1(:)'A1b = A1'A1_row = A1b(:)'resultsa = 1 2 3 4a_row = 1 2 3 4
2022-03-02 10:24:38
9656
原创 os.path.join最大好处是兼容windows和linux系统，频繁需要的操作
#os.path.join最大的意义是兼容windows和linux系统，不用管是"/"还是“\"连接各个路径；os.getcwd() pt = os.listdir("../")absp = os.path.abspath("../")for it in pt: pa = os.path.join(absp,it) if(os.path.isdir(pa)): print(os.path.join(absp,it))...
2022-02-17 11:09:48
952
原创 python的nan生成、表示、判决和忽略nan的统计运算
"""nan 三种生成方法，判断方法，来源"""#Gen nanxf = float('nan')import numpy as npxn = np.nanfrom decimal import Decimalxd = Decimal('nan')# souce of nan, 'inf' -'inf', inf/infx1 = float('inf') + 2*float('inf')x2 = float('inf') - float('inf')x3 = float('i
2022-02-11 17:33:32
1878
原创 发现一本C++面向对象编程的好书
C++好书目标目标适合 大型程序 的C++语言，要想获得兼顾灵活性和高效率，难免的是面向对象和模板模式，这里想发现一本实操的好书不容易。最近发现一本杜茂康《面向对象程序设计》讲类、对象、友元、继承、多态性、虚函数、重载、I/O流类库、文件、模板与STL、C++ Windows程序的结构、消息驱动、MFC应用程序框架、GDI、菜单、对话框、工具栏、文档与视图等内容比较好。这里可以下载：下载链接：面向对象程序设计-杜茂康https://cn.mianfei-dianzi-shu.site/leixin
2022-01-28 09:14:27
1001
原创 拉格朗日 优美典范
%syms： good lagrange 朗格朗日好样的； %两种方式1：x=x0, eval(f1); 2:使用{x,y},{x,y} subs替换函数，属性-值的花括号（元组对）即可；syms x1 x2 x3 x4 x x5 x6exp = (x-x2)*(x-x3)*(x-x4)/((x1-x2)*(x1-x3)*(x1-x4))f1 = subs(exp, {x,x1,x2,x3,x4},{0,-2,-1,1,2})f2 = subs(exp, {x,x1,x2,x3,x4},{0,-1,
2022-01-20 21:59:14
207
原创 MU-MIMO有意思的用户配对/调度
简单compare相关性来逐个选取R = V1*V2’；ui = smallest®;code%MU_MIMO_Schedule.m MU配对逻辑% core: smallest relative users schedule together; R = V1V2', min norm_f(R)% 1 find strong user0; 2 find smallest relative R(Sta0,Sta1); 3 find% samllest relative R(Sta(n-1), S
2022-01-07 10:36:52
1105
1
原创 Linux libreOffice的找不到Language spell check拼写检查问题
Linux上用Office系统安装方法去掉spell check拼写安装方法整体上掠过，下载libreOffice，解压后，打开Program，找到soffice后打开，我用的时xlx，所以都以soffice做说明。word之类也有对应的执行文件。去掉spell check拼写保密系统里面不能拷贝出来数据，所以希望直接截图统计数据，但是有个麻烦，系统自动spell check拼写检测，有很多下滑波浪线非常碍眼。看了好多，找language》options》去掉勾languae spell chec
2021-12-24 09:43:29
953
1
原创 PAR 相控阵方法，阵元间隔原理
这里写自定义目录标题Motivation搞不懂相控阵波束如何变化，为何Massive MIMO的阵元都设定为0.5λ，为何波束有旁瓣？于是查到一个图，搞了2个公式，1个程序初步探索了这个玩意儿！基于ULA线性阵列的相控阵变波束情况。引子我們觀察到這是一個極座標下線陣的波束掃描圖，其中外圈標註的是度數0~360°，0°表示垂直於線陣的方向。其中，一圈一圈的環代表的是幅度，標出的100,80,60…是dB值。波束增益最大的是指向0°的方向，最大值100dB，旁瓣比主瓣低接近40dB，說明進行了降低旁瓣
2021-11-25 15:32:13
1276
2
原创 cpp虚函数，基类命名指向子类实例必备佳选
简单虚函数//============================================================================// Name : p1022VirtualDeconstructor.cpp// Author : perfey// Version :// Copyright : Your copyright notice// Description : Hello World in C++, Ansi-st
2021-10-22 22:02:35
144
原创 Python sympy做代数运算、解Cholesky求逆的L和L逆矩阵示例
提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档文章结构sympy包问题一、代码二、运行结果总结sympy包问题在1.6.x的包中发现用“from sympy import *"遇上频繁的Warning问题？解决方法：shell power >> pip install -U sympy 更新到了1.9系列，问题解决！提示：sympy不仅仅有symbol系列的代数运算功能，其Matrix系列的矩阵数值运算也比较方便，比如Matrix（[[1,2],[3,4]]）这
2021-10-22 12:11:48
907
原创 余弦相似性和相关
person相关首先的概念是相关性是衡量线性关系，而非“非线性”关系。如上图，二维平面点集合第一行可以发现是围绕直线的(x,y)点，比如(1,1.1) 和（1.1,1）几乎为1，反之最右侧的是（-1,1.1） (-1.1,1)，在直观上，他们两队坐标是各自相关的，却是-1和1差距，在两个极端。但实际还是一回事儿，反应了两个变量x,y的"线性相关"程度。—同向增长，反向增长。 介于中间的0.8/0.4这些点则是弱线性相关，增长反向不总一致或反向。0处则是完全无规律。第二行，发现相关方向和斜率无关（线
2021-09-30 21:44:23
1648
原创 小浮点大名堂，1.0元累加2000万次等于2000万吗？ float的下溢
//============================================================================// Name : p930Add2000wTimes.cpp// Author : perfey// Version :// Copyright : Your copyright notice// Description : 将浮点数1，累加2000万次，得到的还是2000万吗？, Ansi-style/
2021-09-30 18:35:21
135
原创 matlab中builtin的regexp用法
提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档regexp @ matlab前言一、怎么用？二、Example1.run2. result3.简化返回总结前言regexp是Regular expression 各取前三个字母组合成的命令。含义一目了然，做regular expression match（正则表达匹配，类似的python中是re包）一、怎么用？– [S, E, TE, M, T, NM, SP] = regexp (STR, PAT)– […] = regex
2021-09-24 10:51:15
882
原创 【硬件资源节约】技巧：利用复数乘模块实现共轭复数乘-不经过取反加1
复数乘简单实现共轭乘问题解决问题已经有一个复数乘模块 Cmultiply(A,B)=C，实现了复数乘法功能：A*B = (a+bj)(c+dj) = ac-bd + (bc+ad)j=C则基于此模块实现复数共轭乘法 A′∗B′HA'*B'^HA′∗B′H?(a+bj)∗(c−dj)=(ac+bd)+(bc−ad)j(a+bj)*(c-dj) = (ac+bd) + (bc-ad)j(a+bj)∗(c−dj)=(ac+bd)+(bc−ad)j要求,不用对B取共轭（因为取反加1，涉及2次硬件计算）是
2021-09-17 09:37:01
389
原创 cpp虚函数、纯需函数样例
//============================================================================// Name : prac0908virtualFunction.cpp// Author : perfey// Version :// Copyright : Your copyright notice// Description : practice virtual function and pure
2021-09-08 15:38:00
86
原创 python自动脚本执行——之子进程管理包 subprocess
自动化脚本魅力subprocess包很管用subprocess包很管用最常用的方法有：subprocess.run .call （注意 shell = True）老的版本有：subprocess.getstatusoutpu(“cmd”)
2021-08-30 15:44:04
144
原创 学生信息转换系统，doc读入，一键录入xlxs
'''逐个读取学生信息，转换为格式表'''import osimport docximport xlrd, xlwtfrom xlutils.copy import copyimport re#fileName = "一二班信息"for d in os.listdir(): if d.endswith(".doc") or d.endswith(".docx"): fileName = '.'.join(d.split(".")[:-1]) pri
2021-08-29 12:03:58
146
原创 HTML开坑
Html不仅仅做网页，做属性输入和配置文件也不错配置个visual studiohelloworld配置个visual studio因为也不是新手，跳过用txt文本来写 之类的东西了！配置了一个VS的编译器，好用啊！html各种提示，只需要简单下一个open browser插件就行兼容cpp等语言，非常便利发现所谓标记语言和编程语言是有差异的，html天生不是为了实现功能，实现需要加载js（javascript）之类的插件。helloworld 基本结构 标题、段落、对齐、背景 超
2021-07-25 12:54:15
80
1
原创 WLAN的一个典型实现
规格时钟频率： 160MHz带宽： 10、20 或 40MHz10MHz 和 20MHz 带宽符合 802.1140MHz 实现为双时钟 20MHz，不符合 802.11 HT40 波形OFDM 格式： 64 个子载波，16 个样本循环前缀NONHT 波形 (11ag)：48 个数据子载波，4 个导频音HTMF 波形 (11n)：52 个数据子载波，4 个导频音帧格式：根据 802.11-2012 的第 18 和 20 节的规定：前导码（10 次重复 16 样本短训练符号，2.5 次重复
2021-07-03 12:03:28
314
原创 ofstream输出文件流简易操作
/* * fstreamSteps.cpp * * Created on: 2021年6月26日 * Author: Perfey */#include<iostream>#include<fstream>using namespace std;int main(){ std::ofstream confg_; if (!confg_ .good()) { //std::cout <<strerror(errno)<
2021-06-26 09:28:56
496
原创 实部虚部合并和提取，巧用位运算符＜＜,＞＞, |,宏告警
#include <iostream>using namespace std;#define WARN(eq, str)\{\ if (!(eq))\ {\ cout<<"Warning, eq is not satisfied!" << str <<endl;\ }\}int main(){ cout << "Hello Complex Digital World!"<<endl; // c
2021-06-19 20:40:34
335
原创 cpp namespace使用，全局变量，嵌套namespace
//============================================================================// Name : prac0616NameSpace.cpp// Author : perfey// Version :// Copyright : Your copyright notice// Description : namespace//=============================
2021-06-16 15:04:57
977
原创 html的页面内超链接，为何跳不动
基本格式“点我跳到末尾”想要跳转的正文部分（内部链接锚点位置）：<a id="chap1">正文中：内部链接锚点位置</a>索引目录部分（点击这里可以跳转到预设锚点）：<a href="#chap1">目录中可以跳转的位置</a>坑点编程时候，一定要加上id = “” xxx” 上下索引标号"“点我调到末尾”...
2021-06-15 22:54:19
2782
原创 电力线频率-wifi的荧光效应
荧光效应解释室内常开的荧光灯其交流电是有个频率（比如50hz，60hz），人眼不可见，但是在环境反射散射会造成一种遍布室内的100hz或120hz左右调幅信号（倍频干涉的较强值，其它50hz，200hz之类较为弱一点）。这里高能预警，光的相干、干涉和混频。这些信号会叠加在wifi射频信号上，造成一种电磁波快速小幅度波动现象。—这里涉及到环境曼散射小信号直接调制了Radio信号（现象基带信号x在环境中，神奇地加载、变频到Radio上的场景—光是神奇的，也有干涉、调制、变频等通信中电磁波的一切现象。毕竟本质都
2021-06-11 10:21:00
687
1
原创 IQ复数信号和实数信号处理差别和殊途同归
实数iq调制和复数的对比，晕了好多年，就是不明白。知道现在记住信道只传实数信号，世界上没有复数信号，只是一种数学手段，处理和展示方法，实数处理完全可替代 只是没有相位信息，不如复数直观。这才豁然开朗。...
2021-06-10 21:26:41
2603
原创 matlab的cell索引差异、mkdir参数形式差异
matlab简洁语言和多变形式，容易有些啼笑皆非错误。序号名称举例说明1cell的索引(),{}有差异a = {‘first’,‘second’}创建一个cell类型--a[1] vs a{1}前面依然是cell类型，后面是cell（小蜂窝）内部的类型；即1个返回cell，1个返回字符串类型2mkdir两种形参var = “result”建立一个字符串变量--mkdir var这时候建立的不是预期的result文件夹，而是名称为var的文件
2021-06-07 18:08:27
218
原创 近似计算-简单举例 复数乘法和strassen算法
原因计算是第一生产力，节约计算等于提升了计算效率，自然就有巨大价值。尤其是AI、通信、并发这些耗计算大户下。方法以前从来没接触过，所以没法长篇大论，只举例我自己知道的2个。一维的复数计算（或者类似（a+b）*（c+d）这种多项式乘法），二维的矩阵计算一维复数计算，每4次乘法节约1次，提效25%原本计算（a+jb）(c+jd) = (ac -bd) + j(bc + ad), 4次乘法一目了然。节约计算：先算x= (a+b)(c-d) = (ac-bd)+bc-ad; 左侧只有1次乘法（加法耗计
2021-06-05 18:48:31
547
原创 ofdm小问题汇总
1 非理想因素 1）频偏CFO（时钟+其它）、2）采样定时偏差、3）相噪（高频更不可忽视）、4）PA非线性、5）DC偏置2 加扰目的 信元随机化；避免由长0长1引发高脉冲（时域的长1频域是尖峰）3 星座图的各点错误概率相同吗？ 1）4QAM以上高阶都不同（外围星座点概率低，包含无线空间）；2）应用重传上；3）研发非QAM的高阶调制方式4 FFT变换两种形式 1）由0N-1和N/2N/2-1 (时序起始不同，前者matlab，后者协议）； 2）两者可以等效（一种等效法是圆周移位一半N/2，并有F矩阵式.
2021-05-29 12:00:35
472
原创 预处理三剑客（宏定义，文件包含和条件编译）之第三种——条件编译。 #ifndef #ifdef
预处理三剑客（宏定义，文件包含和条件编译）之第三种——条件编译。 #ifndef #ifdef用法 #ifdef 标识符 #ifndef 标识符 程序段1 程序段1 #else #else 程序段2
2021-05-19 10:51:02
72
原创 位数统计问题：获取一系列数据，统计bit数，并画出直方图
"""查看所有数据有多少bit"""import collectionsimport numpy as npimport matplotlib.pyplot as pltimport matplotlibwith open("f.txt") as f: x = f.readlines()digs = []for i in x: if len(i.split(" ")) > 1 : #print(i.split(" ")[1].split("\n")
2021-05-18 21:30:40
110
原创 python随着脚本传入系统参数
easy way# 简便方法:sys.argv输入一个列表，第0个为函数名本身，依次空格隔开各个字符串输入import syssysInput = sys.argvprint("input type", type(sysInput))for i in sysInput: print("Input para : ", i) #即第一个参数是 sysInput[1]。。。run(base) D:\Code\Python\CodeVerify\BasicSyntax>pyt
2021-05-17 16:03:30
141
原创 Wifi开坑
缘起：读书期间，大量接触通信基础，无线快速发展，目光关注点自然是高高的基站塔，以及小小的手机（这个信息社会浓缩操作系统、通信、开源和智能化（智能掌机）的结晶当然也不简单）。 至于短距离通信里面的形形色色，眼花缭乱的Wifi、BLE（蓝牙）、Zigbee、体域网则是感知中盲区，异类。一直到博士后期到职业才隐隐感知到，天下大事必出于细的道理。云动：就如“小小”体域网（Body area network）中除了通用的天线设计、收发、同步、多通道处理外还有特别关注的鲁棒性设计，如何高可靠无中断地始终获取无线信
2021-05-17 15:07:05
28566
6
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人