matlab
文章平均质量分 71
enjoy_learn
这个作者很懒,什么都没留下…
展开
-
低通滤波器实现过程解析
首先来看一个低通滤波的例子:给定一组波形如下: 如果想要获取低通0-0.5Hz的波形:需要调用低通滤波器可得到:很简单的一个例子:程序见:点击打开链接低通滤波器实现过程:主函数: y=lowp(x,f1,f3,rp,rs,Fs)X:输入信号;Fs: 采样频率f1:通带截止频率,rp:通带衰减f3:阻带截止频率,rs:原创 2017-01-07 15:34:08 · 11962 阅读 · 0 评论 -
matlab 画图时日期格式坐标轴的处理
引自:http://blog.sina.com.cn/s/blog_6328ac360102wk5n.htmlmatlab作图时经常会遇到X轴为时间的情况,查了下可用用datenum()与datetick(‘x’,6);实现坐标轴的时间刻度绘制;...转载 2019-05-31 09:54:11 · 18688 阅读 · 0 评论 -
修改matlab读取文件夹名字顺序
dirOutput1=dir(fullfile(str_path,’*XX’));%匹配XX;file_path_24=sortObj(dirOutput1);function file=sortObj(file)for i=1:length(file)A{i}=file(i).name;end% [~, ind]=natsortfiles(A);% [~, ind]=natsort...翻译 2019-03-01 10:34:59 · 1412 阅读 · 1 评论 -
相关计算xcorr
参考自:https://blog.csdn.net/u010592995/article/details/79841108N=256;di=20;x1=test(1+di-10:N+di-10);%x1延后x2 10个点x2=test(1+di:N+di)时域相关计算:N =length(x2);[xc,zuobiao]= xcorr(x1,x2);[k,ind] = max(xc...转载 2019-01-03 20:03:43 · 1792 阅读 · 0 评论 -
matlab卷积
参考自:http://sepwww.stanford.edu/data/media/public/sep//prof/index.htmlfree book Earth Sounding Analysis Processing Vers.pdfy0=x0b0;y1=x1b0+x0b1;y2=x2b0+x1b1+x0b2;y3=x3b0+x2b1+x1*b2;…卷积定义:functi...原创 2018-12-18 20:57:21 · 2346 阅读 · 1 评论 -
matlab画一个圆,并均分角度
作图画一个圆,并标注角度;clc;clearclose allR=1;t=0:pi/20:2*pi;x=R*cos(t);y=R*sin(t);plot(x,y),axis equaln=72;a=2*pi/n;% for k=0:n-1% hold on% plot([-cos(pi+k*a),cos(pi+k*a)],[-sin(pi+k*a),sin(pi+k*a)]),...原创 2018-11-27 19:39:46 · 4973 阅读 · 0 评论 -
envelope包络与findpeaks求极值点包络对比
matlab 自带函数:function [upperEnv,lowerEnv] = envelope(x, n, method)函数可使用三种方式计算:如下matlab帮助文档:% [YUPPER,YLOWER] = ENVELOPE(X,N,ENVTYPE) specifies the type of envelope% to return. The default is 'an...原创 2018-12-02 21:24:29 · 7454 阅读 · 0 评论 -
matlab滤波器
懒得写了,贴这function y=bandp(x,f1,f3,fsl,fsh,rp,rs,Fs) %带通滤波 %使用注意事项:通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一半 %即,f1,f3,fs1,fsh,的值小于 Fs/2 %x:需要带通滤波的序列 % f 1:通带左边界100 % f 3:通带右边界 % fs1:衰减截止左边界 % fsh:衰变截止右边界 ...原创 2018-07-20 20:19:08 · 9339 阅读 · 0 评论 -
关于fft相位
参考:这里写链接内容angle与phase的区别:源地址: http://www.ilovematlab.cn/thread-259049-1-1.htmlphase 和 angle 在输入为单个标量数据时,没有差别,二者都是用 atan2 函数来求输入数据的四象限辐角。但是对于向量或矩阵数据输入时,二者差别非常大。phase 只支持标量和一维(行、列)向量输入,不支持二维或高维矩...转载 2018-06-01 12:49:19 · 4881 阅读 · 0 评论 -
滑动平均滤波
滑动平均滤波,窗口长度为6; matlab窗口长度6进行滤波,输入为512点,输出为507点,最后几个点未进行处理; 滑动平均后,波形中的小波动消失; 用于平滑噪声: matlab :codefunction T1 = huadong( d_temp )%UNTITLED Summary of this function goes here% Detailed ex...原创 2018-05-30 14:42:05 · 31942 阅读 · 2 评论 -
时频域能量相等(parseval定理)
帖子-时域和频域能量相等 Parseval 定理 有限上序列x{k}的离散fourier变换是正交变换,满足Parseval能量守恒定理,反映了序列在时域的能量等于其变换域的能量。 关于能量定义:信号幅度平方的积分,如果是数字信号,能量就是各点信号幅度值平方后的求和。 论坛帖子中关于等式关系给出的结论是: 求和 (x(tn)^2)T=RMS^2*Ttotal=求和(P(fn)原创 2017-03-26 18:13:24 · 31873 阅读 · 4 评论 -
关于幅值谱逆变换得原始波形的问题(matlab)
前些天突然碰到一个问题,一直未曾细究过,在剔除频谱中某些频率成分时,时域能量与频域能量本应该是相等的,这就是著名的 Parswel定理嘛,无奈忘了时域有效值和频谱直接算有效值到底是怎么样一个对应关系,后面再细细扒拉扒拉; 尽然忘了换算关系,fft逆变换总可以吧,这时才发现原来这也有问题哈。得到了频谱,但是要从频谱再倒回去得到时域波形,可就是难了,查了下资料看来相要准确的时域波形还原,需要相位信原创 2017-03-09 22:30:07 · 3574 阅读 · 1 评论 -
matlab 批量读取execl(csv)文件
一直没时间整理自己写的垃圾代码,如今代码乱的一团糟,今天把matlab读取excel文件拿出来;需要根据数据格式稍作修改就可以用;读取核心的语句莫过这两句:dir(); csvread();文件下载:%read csvdata%author:enjoy_learn@126.com%str='D:\玛依塔斯X期20#\1V';%确定分析目录function x=batch原创 2017-02-26 18:58:21 · 22283 阅读 · 0 评论 -
小波包分解一例
说到小波,真是痛苦,记得上小波分析这门课的时候,听的是一头雾水,都是说小波是继傅里叶变换后又一重大发现,小波的论文快是被学术界写烂了,但是工作了才发现小波真的是应用的少之又少,还是15年七八月份的事,应用研究了下,现在想想,真是快万忘了。推荐一本书不错:MATLAB 6.5 辅助小波分析与应用 飞思科技.pdf.不是讲理论,而是直接将怎么用的。下载链接:拿一组波形示意如下:对此进行6...原创 2017-02-22 20:41:29 · 12287 阅读 · 7 评论 -
matlab 比较两个文件夹内CSV文件夹
比较两个文件夹的文件名,不同则删除;%比较两个文件夹下文件名,输出文件不同的文件名stra='';strb='';str11=[stra,'\*.csv'];jiegou1=dir(str11);%等价于下面这行geshu1=length(jiegou1);str12=[strb,'\*.csv'];jiegou2=dir(str12);%等价于下面这行geshu2=lengt原创 2017-02-27 22:43:13 · 590 阅读 · 0 评论 -
C语言笔记--代码学习笔记--C语言语法--基本操作运算-basic-logorithm
2017-1-7function1:C语言如何实现一个复数加法,好像C语言里没有matlab的复数变量i,那么怎么办:typedef struct _F_COMPLEX {double real;double img;}F_COMPLEX;然后在需要的地方调用:void add(F_COMPLEX a,F_COMPLEX b,F_COMPLEX *c){原创 2017-01-07 16:48:36 · 861 阅读 · 0 评论 -
matlab结构体数组获取变量名
sheetnames_nf=fieldnames(data.data);name1=sheetnames_nf{1};name2=sheetnames_nf{2};normal_data=getfield(data.data,name1);原创 2019-06-03 13:00:21 · 8641 阅读 · 0 评论