- 博客(25)
- 收藏
- 关注
RSS订阅原创 多项式平滑：Savitzky-Golay滤波
多项式平滑%五点三次平滑法n = length(x);a = x;for k = 1:50 b(1) = (69*a(1)+4*(a(2)+a(4))-6*a(3)-a(5))/70; b(2) = (2*(a(1)+a(5))+27*a(2)+12*a(3)-8*a(4))/35; for j = 3:n-2 b(j) = (-3*(a(j-2...
2020-04-17 17:08:11
4393
原创 word 中 ctrl+v不能使用问题
1.问题描述最近打开word，发现ctrl+v不管用了，每次都要选择粘贴类型，使用默认之后还不能解决问题2.解决方法文件->选项->加载项->管理下拉框选择模板->转到->取消勾选MathType->确定尝试了一下，成功解决...
2020-04-15 13:59:27
10209
5
原创 matlab 出现问题：未定义函数或变量
1）键入函数或变量名时的拼写错误（例如，当您要键入小写字母 l 时键入了数字1）。同时，MATLAB 的新版本提供“是不是想输入:”来解决这些拼写错误。例如：foo = 42;fo0复制代码未定义函数或变量 ’fo0’是不是想输入:foo2）函数或变量名的错误大小写。同时，MATLAB 的新版本提供“是不是想输入:”来解决这些拼写错误。3）改变了路径，使得被调用的函数不在 ...
2019-12-05 18:19:54
22350
原创 matlab画图，均值±标准差
Average=[12,11,7,7,6,5];Variance=[0.5,0.4,0.3,1,0.3,0.5];Time=0:1:5;errorbar(Time,Average,Variance,'-or')
2019-10-30 09:00:47
9881
原创 visio查找希腊字母及下标
base visio查找希腊字母插入——符号——TIMES——Times New Roman——子集希腊语上下标是这样的：选中要成为上标的文字,ctrl+shift+"="选中要成为下标的文字,ctrl+"="...
2019-10-19 11:14:32
2760
原创 关于matlab 画水平线
1）画水平线 x=get(gca,'xlim'); y=0.001; hold on plot(x,[y y]) 2） 垂直线：如 x = 5 line([5 5],[0 100]) 或者plot([5 5],[0 100]) 水平线：如 y = 5 line([0 100]，[5 5])...
2019-10-18 22:14:39
10259
原创 c语言的绝对值
fabs与abs的区别：fabs的参数为double型，返回值也是double型。abs的参数为int型，返回值也是int型。abs是求一个整数的绝对值，而fabs是求一个实数的绝对值。最后记得加上#include<math.h>...
2019-07-05 09:43:58
4738
转载 C语言实现，MATLAB生成巴特沃斯滤波器
环境：matlab2014a工具：Filter Design & Analysis Tool参数设置如下图：这里用的是250HZ的采样频率，低通截止频率为20HZ，6阶巴特沃斯滤波器。设置好参数后，点最下方的“Design Filter”，生成参数。然后在菜单File下Export导出参数：点Export后即可在工作区看到导出的数据Hd.最后使用命令行将参数转换为C语言可...
2019-07-04 15:28:27
4148
1
原创 卡尔曼求roll
#卡尔曼滤波#求取rollimport mathQ_angle = 0.001 #陀螺仪噪声的协方差Q_gyro = 0.005 #陀螺仪漂移噪声的协方差C_0 = 1R_angle = 0.5 #加速度计的协方差q_bias = 0.1 #陀螺仪漂移dt = 0.005P = [0, 0, 0, 0]PP = ...
2019-06-21 15:34:55
192
原创 卡尔曼求pitch
#卡尔曼滤波#求取pitchimport mathQ_angle = 0.001 #陀螺仪噪声的协方差Q_gyro = 0.005 #陀螺仪漂移噪声的协方差C_0 = 1R_angle = 0.5 #加速度计的协方差q_bias = 0.1 #陀螺仪漂移dt = 0.005P = [0, 0, 0, 0]PP =...
2019-06-21 15:34:24
225
原创 四元数法求姿态角，九轴
coding:utf-8四元数法9轴计算#第一版import mathIMU算法更新Kp = 2.0 # 比例增益控制加速度计/磁强计的收敛速度Ki = 0.005 # 积分增益控制陀螺偏差的收敛速度halfT = 0.01 # 采样周期的一半def Update_IMU(ax, ay, az, gx, gy, gz, mx, my, mz):global q0, q1, ...
2019-06-21 15:32:57
3038
3
原创 四元数法求姿态角，六轴
coding:utf-8四元数法6轴计算#修改版import mathIMU算法更新#Kp = 100 # 比例增益控制加速度计/磁强计的收敛速度#Ki = 0.002 # 积分增益控制陀螺偏差的收敛速度#halfT = 0.001 # 采样周期的一半Kp = 2.0 # 比例增益控制加速度计/磁强计的收敛速度Ki = 0.005 # 积分增益控制陀螺偏差的收敛速度...
2019-06-21 15:32:05
2947
3
转载 四元数AHRS姿态解算和IMU姿态解算分析
AHRS是自动航向基准系统(Automatic Heading Reference System)的简称。目前，使用四元数来进行AHRS姿态解算的算法被广泛采用于四轴飞行器上。该算法源自英国Bristol大学的Ph.D Sebastian Madgwick，他在2009年开发并发布了该算法。下面我们来对该算法的代码进行详细分析。AHRS是自动航向基准系统(Automatic Heading Ref...
2019-06-17 10:27:01
1228
3
原创 deconv 函数介绍
deconvDeconvolution and polynomial divisionSyntax[q,r] = deconv(u,v)Description[q,r] = deconv(u,v) 使用长除法将向量v从向量u中解卷积，并返回商q和余数r，使得u = conv（v，q）+ r。 如果u和v是多项式系数的向量，则对它们进行去卷积相当于将由u表示的多项式除以由v表示的多项式。...
2019-02-21 11:04:49
5989
转载 功率谱密度
（如果不先对数据都减去一个均值，那么接近0Hz部分的功率会非常非常大~这是有害的，一般做数据处理的时候如果只关注交流量的话，就要先使得整个序列的均值为0，方法就是给每个数减去该序列的一个均值2012-01-12）from http://www.ymlib.net/article/sort010/info-279.htmlmatlab求功率谱密度代码实例clear;Fs=1000; %采样频...
2019-01-11 10:08:10
3066
转载 [转载]matlab中[Pxx,f] = periodogram(x,window,nfft,fs)
原文地址：matlab中[Pxx,f] = periodogram(x,window,nfft,fs)作者：海上冲浪手http://zhidao.baidu.com/question/403225309.htmlmatlab中[Pxx,f] = periodogram(x,window,nfft,fs)periodogram是用来计算功率谱密度的，参数中，X：所求功率谱密度的信号；win...
2019-01-07 17:43:11
11025
3
原创 频谱相谱功率谱
第一：频谱一.调用方法X=FFT(x)；X=FFT(x，N)；x=IFFT(X);x=IFFT(X,N)用MATLAB进行谱分析时注意：（1）函数FFT返回值的数据结构具有对称性。例：N=8;n=0:N-1;xn=[4 3 2 6 7 8 9 0];Xk=fft(xn)→Xk =39.0000 -10.7782 + 6.2929i ...
2019-01-07 15:46:21
1189
原创 fft用法
matlab 中fft的用法一.调用方法X=FFT(x)；X=FFT(x，N)；x=IFFT(X);x=IFFT(X,N)用MATLAB进行谱分析时注意：（1）函数FFT返回值的数据结构具有对称性。例：N=8;n=0:N-1;xn=[4 3 2 6 7 8 9 0];Xk=fft(xn)→Xk =39.0000 -10.7782 + 6.2929i...
2019-01-02 14:44:20
1209
原创 EEGLAB 使用记录
EEGLAB 使用记录EEGLAB 使用记录1 EEGLAB 安装2 导入数据2.1 导入MATLAB矩阵数据（*.mat）3 预处理3.1 滤波记录一下个人使用的需求和过程。推荐两本比较全面详细的EEGLAB 使用手册！1 EEGLAB 安装下载安装包，在MATLAB中点击设置路径，将EEGLAB文件和子文件夹设置完成即可，然后在命令行输入：eeglab12 导入数...
2018-12-28 17:31:21
886
原创 gradle
Gradle是一个基于Apache Ant和Apache Maven概念的项目自动化建构工具。它使用一种基于Groovy的特定领域语言来声明项目设置，而不是传统的XML。当前其支持的语言限于Java、Groovy和Scala，计划未来将支持更多的语言。上面是维基上对Gradle的解释,相信一个没有接触过构建的人是不大能看明白的,当初我也是.下面是我对Gradle通俗的理解:软件开发讲究代码复用,通...
2018-12-14 11:21:01
84
原创 时域和频域
1.最简单的解释频域就是频率域，平常我们用的是时域，是和时间有关的，这里只和频率有关，是时间域的倒数。时域中，X轴是时间，频域中是频率。频域分析就是分析它的频率特性！图像处理中：空间域，频域，变换域，压缩域等概念！只是说要将图像变换到另一种域中，然后有利于进行处理和计算比如说：图像经过一定的变换（Fourier变换，离散yuxua DCT 变换），图像的频谱函数统计特性：图像的...
2018-11-26 15:15:45
5682
转载 快速傅里叶变换
一、快速傅里叶介绍傅立叶原理表明：任何连续测量的时序或信号，都可以表示为不同频率的余弦（或正弦）波信号的无限叠加。FFT是离散傅立叶变换的快速算法，可以将一个信号变换到频域。那其在实际应用中，有哪些用途呢？1.有些信号在时域上是很难看出什么特征的，但是如果变换到频域之后，就很容易看出特征（频率，幅值，初相位）；2.FFT可以将一个信号的频谱提取出来，进行频谱分析，为后续滤波准备；3.通过对...
2018-11-22 17:00:55
430
1
转载 LSB概念
最近做了一块板子，当然考虑到元器件的选型了，由于指标中要求精度比较高，所以对于AD的选型很慎重。我们所讲的精度通常是指它的精确度，其实这是错误的。精度又叫做精密度，是跟准确度相对应的一个概念。就像打靶一样，打的准，那就说它的准确度比较高；而每两个靶之间能打出的偏移越小，那它的精密度就越高。精密度与准确度合起来称为精确度。但是鉴于大家都将精度指代了精确度，那以下所说的精度如无特别指出，都是指精确度...
2018-11-15 19:48:49
9225
原创 matlab快捷操作
MATLAB中怎样在一段程序的每一行前面一次性加上百分号？将要加百分号的程序段选上，然后 ctrl+R 就可以实现全部加上百分号； ctrl+T 可以取消所选语句前面的百分号另外，还有一个比较好用的快捷键 ctrl+i 可以将一段程序整理为按规范的层次结构...
2018-11-15 09:31:28
515
1
原创 dps定义
dps：全量程，L3G4200D可以设置3种dps，分别最大可测量：250/500/2000度/秒 的角速度。mdps/digit：灵敏度，指输出的每数字1代表的毫角速度 。当dps设置为 250dps 时，灵敏度是8.75mdps/digit，则：输出是2000时，角速度 = 2000*8.75/1000=17.5度/秒...
2018-11-14 15:50:54
2407
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人