精华内容
下载资源
问答
  • 本帖最后由 芒果樱桃 于 2012-4-8 09:05 编辑希望实现功能是通过加自卷积窗函数,获得正确的相位、频率和幅值根据丁康的相位差校正,试着编了一个程序,将信号分成两段,然后分别加窗,通过相位差校正法求出信号信息...

    本帖最后由 芒果樱桃 于 2012-4-8 09:05 编辑

    希望实现功能是通过加自卷积窗函数,获得正确的相位、频率和幅值

    根据丁康的相位差校正,试着编了一个程序,将信号分成两段,然后分别加窗,通过相位差校正法求出信号信息,

    不知程序这样写行不行,新手求解答

    %相位差校正

    clc

    clear all

    fs=1024;

    N=1024;

    M=N/2;

    L=512;

    t1=0:(N-1);  %序列一段长度

    t2=0:(N+L-1);  %序列二段长度(本来是从L到L+N)

    df=fs/N;%频谱分辨率

    %信号函数

    ph1=0.1;%相位1

    ph3=0;%相位3

    f1=50.2;

    f3=99.8;

    x=0.2+6*sin(2*pi*f1*t1/N+ph1*pi/180)+sin(2*pi*f3*t1/N+ph3*pi/180);

    x2=0.2+6*sin(2*pi*f1*t2/N+0.1*pi/180)+sin(2*pi*f3*t2/N+ph3*pi/180);

    %plot(linspace(0,pi,N),angle(fft(x,N)));

    %原始信号函数的输出

    X=fft(x);  %做FFT变换

    Ayy=abs(X)/(N/2);  %取模,换算成实际的幅度

    Ayy(1)=Ayy(1)/2;  %直流分量是(1),上句话多乘了2,所以这句要除以2

    f=(1:N-1)*fs/N;  %换算成实际的频率值

    for k=1:N/2

    ph(k)=phase(X(k))*180/pi+90; %计算相位,换算为角度

    end;

    %加blackmanharris自卷积窗

    window=blackmanharris(M)';

    win=conv(window,window);  %blackmanharris窗函数的时域相乘,即为频域卷积

    win1=ifft(fft(win,N));

    win2=[zeros(1,length(x2)-N),win];%对窗函数补零,序列2段本来应该也是长度为N

    win2=ifft(fft(win2,N+L));

    X1=fft(x.*win1);  %信号1段加窗

    A1=abs(X1);

    Ph1=angle(X1)*180/pi+90;  %求相位

    Ph1=[Ph1,zeros(1,L)];

    X2=fft(x2.*win2);

    A2=abs(X2)/(N/2);

    Ph2=angle(X2)*180/pi+90;

    dph=Ph2-Ph1;

    delta=dph-2*pi*L/N;

    delta=mod(delta,2*pi);

    for k=1:round(f/df);%小数四舍五入,这里返回值是整数,k为谱线号

    if  delta(k)>pi %delta为两序列谱峰线间的相位差,(0~2*pi)

    delta(k)=delta(k)-2*pi;

    elseif delta(k)

    delta(k)=delta(k)+2*pi;

    end

    end

    dfy=(dph-2*pi*L*k/N)/(2*pi*L/N);

    fconr=(k+dfy)*df;

    sita=atan(imag(X)/real(X))+dfy;

    展开全文
  • 采用快速傅里叶变换(FFT)对电力系统进行谐波分析时,由于非同步采样和数据截断...通过MATLAB仿真分析,验证了基于三次样条函数的加Rife-vincent自卷积窗插值FFT算法能够有效抑制频谱泄露和栅栏效应,提高谐波检测精度。
  • matlab函数与滑

    2017-05-11 15:33:00
    在处理图像时,为了提取特征,经常用各种核函数和图像进行卷积,其实就是通过一个矩阵以滑的形式与原图像进行点乘求和,可以看作对一个像素和附近像素进行了加权平均。 比如经常用3x3的近似高斯卷积核 0 1 0 ...

    在处理图像时,为了提取特征,经常用各种核函数和图像进行卷积,其实就是通过一个矩阵以滑窗的形式与原图像进行点乘求和,可以看作对一个像素和附近像素进行了加权平均。

    比如经常用3x3的近似高斯卷积核

    0 1 0

    1 2 1

    0 1 0

    其实就是这幅图的俯视图中z值(高度)的近似,中间最高,边缘低。可见原来的一个像素现在经过卷积由周围的像素加权构成,所以以高斯为例,高斯有平滑效果。

    但是刚才的举例是近似的于真正的高斯分布还有差距,我们怎么用matlab产生真正的高斯核呢?

    fspecial函数:

    可以在matlab的命令行输入help fspecial来查看用法

    使用 h = fspecial('gaussian', hsize,sigma); 就可以产生一个高斯核,其中hsize是卷积核的大小,默认3x3,sigma是高斯核的方差

    在用 filteredI = imfilter(I, h); 对原图像用刚才生成的和进行卷积就可以得到卷积后的图像。

    fspecial还可以生成其他核,都很常用,各有用途。

    'average'  'disk' 'laplacian' 'log' 'motion' 'prewitt'  'sobel'  'unsharp' 。

     

    但有的时候需要自定义一个函数来通过滑窗处理图像,matlab也有相关的函数 -------- blockproc

    B = blockproc(A,[M N],fun)

    A是要处理的矩阵,[M N]是滑窗大小,fun是处理时执行的函数

    转载于:https://www.cnblogs.com/maskmtj/p/6841159.html

    展开全文
  • matlab区分卷积和相关

    千次阅读 2017-11-12 16:21:24
    区分卷积和相关 图像处理中常常需要用一个滤波器做...非线性滤波常见的有中值滤波、最大值滤波等,相当于自定义一个函数,在数学上由于不满足线性变换因此叫做非线性滤波。这里不细研究它。 线性滤波则通常是:将模

    区分卷积和相关

    图像处理中常常需要用一个滤波器空间滤波操作。空间滤波操作有时候也被叫做卷积滤波,或者干脆叫卷积(离散的卷积,不是微积分里连续的卷积);滤波器也有很多名字:卷积模版、卷积核、掩模、窗口等。

    空间滤波可以分为线性滤波和非线性滤波。非线性滤波常见的有中值滤波、最大值滤波等,相当于自定义一个函数,在数学上由于不满足线性变换因此叫做非线性滤波。这里不细研究它。

    线性滤波则通常是:将模版覆盖区域内的元素,以模版中对应位置元素为权值,进行累加。看起来挺简单的,但是要区分相关(cross-correlation)卷积(convolution)两种模式。为什么呢?因为在MATLAB里是有所区分的,而且不少中文书里面把它们混淆了。

    我们最容易的理解是:将模版中元素从左到右、从上到下,作为使用顺序,那么卷积操作的结果,就是模版第一个元素乘以它覆盖的元素,加上模版第二个元素乘以它覆盖的元素,再加上模版第三个元素乘以它覆盖的元素,...,一直加到模版最后一个元素乘以它覆盖的元素。好吧,其实就是:模版覆盖区域内,元素逐一相乘然后累加,此时的对应位置就是上下投影后被覆盖的位置。

    上面这个理解确实是看起来最容易理解的,因而很多中文书把它叫做"卷积"。然而这个概念其实叫做相关,而卷积则相当于:将同样的模版旋转180°后,再做"相关"操作。当然,如果模版是180°对称的那么卷积和相关是相同的。但是并不是所有的模版都对称。因此,我建议,在滑窗操作、计算图像梯度等场合,不要使用“卷积”,而要使用“滤波”或者“相关”。因为,我们通常讲的卷积,其实是相关,那就不要用卷积这个词以免引起混淆。


    imfilter(f, w, filtering_mode, boundary_options, size_options)

    f:图像

    w:滤波模版

    filtering_mode:滤波模式 'corr':相关滤波。[默认值] 'conv':卷积滤波。

    boundary_options:边界选项 P:(没有引号)边界外围补充0。[默认值] 'replicate':边界外围复制边界值 'symmetric':边界外围使用边界镜像 'circular':图像的大小通过讲图像处理为二维周期函数的一个周期来扩展

    size_options:大小选项 'same':输出大小与输入图像f大小相同。[默认值] 'full':输出与扩展(填充)后的图像大小相同。

    展开全文
  • matlab常用内置函数

    千次阅读 2018-09-30 15:00:29
    2.conv() --- 卷积操作 3.subplot() --- 将figure中的图片划分为几块 income = [3.2 4.1 5.0 5.6]; outgo = [2.5 4.0 3.35 4.9]; subplot(2,1,1); plot(income) subplot(2,1,2); plot(outg...

    1.figure() ---  创建一个用来显示图形输出的一个窗口对象

    2.conv() --- 卷积操作

    3.subplot() --- 将figure中的图片划分为几块

    income = [3.2 4.1 5.0 5.6];
    outgo = [2.5 4.0 3.35 4.9];
    subplot(2,1,1); plot(income)
    subplot(2,1,2); plot(outgo)
    先用subplot划分Figure为2*1的区域,上面一个显示income,下面一个显示outgo。
    subplot的格式为(行数,列数,第几个区域的位置),不加逗号也是可以的。

    4.ones() --- 创建元素为0的矩阵  zeros() --- 创建元素为1的矩阵

    5.load(filename) 从 filename 加载数据。

     

    展开全文
  • matlab函数大全-matlab函数大全.doc 比较白痴的东西,不过对初学者也许有用,特分享一下。不要喷我哦!函数具体用法可以用help查一下。附件内容与下面一样的。 Aabs 绝对值、模、字符的ASCII码值 acos 反余弦 ...
  • 或者去掉窗函数但引入时延tao,波形都会畸变。 %h=sin(wc*(n2))./(pi*(n2)); %line14,若注释掉上两行,只采用理想滤波的序列,即此行,效果还行。 len_h=length(h); ...
  • 实验要求,matlab实现卷积运算(不调用函数和方法)1.1 卷积原理讲解1.1.1 卷积计算步骤(图解法)1.1matlab画出两个窗函数1.2 使用c语言画出该图像1.3 matlab卷积实现(非调用函数) 1.实验要求,matlab实现卷积...
  • MATLAB函数参考

    2012-03-19 20:43:56
    matlab函数网上搜集并精心整理,添加了目录,便于查阅。共23页。 文档内容举例: 函数名 功能描述 addpath 增加一条搜索路径 目 录 1 常用命令 1.1 管理用命令 1.2管理变量与工作空间用命令 1.3文件与操作系统...
  • matlab常用函数大全

    2010-06-20 10:28:18
    1.5 指令控制 1.6 操作系统指令 2、运算符和特殊运算符 2.1 算术运算符 2.2 关系运算符 2.3 逻辑操作 2.4 特殊运算符 3、编程语言结构 3.1 控制语句 3.2 计算运行 3.3 脚本文件、函数及变量 ...
  • 传递函数卷积.rar

    2019-06-21 11:11:02
    压缩文件中给出了已知系统传递函数及输入信号的情况下,求解系统输出响应的四种方法,即差分方程法,卷积法,DFT法以及滤波器设计方法;此外还包含创建及其频率响应的matlab设计。
  • matlab 数字信号处理函数

    热门讨论 2011-06-07 20:42:42
    % blackman - 布莱克曼窗函数(第230页) % freqz_m - 改进型的freqz 子程序(第233页) % Hr_Type1 - 计算1型FIR低通滤波器(第215页) % Hr_Type2 - 计算2型FIR低通滤波器(第216页) % Hr_Type3 - 计算3型FIR低通滤波器...
  • 卷积定理可知:时域相乘等于频域卷积,例如矩形窗R5(n)(形式比较简单一种窗函数)的频域幅度谱如图1左子图所示,将它卷积信号的频谱图,就造成“ 拖尾现象 ” ,称之为频谱泄漏。 图1 若序列 的长度为...
  • 将高斯核分为水平方向和垂直方向,先后对图像进行卷积处理,可提高处理速度。 function output=gaussFilter(I,sigma) output=I; ksize=double(uint8(3*sigma)*2+1);%窗口大小一半为3*sigma window = fspecial('...
  • 半年多的坑了,当时遇到的问题都块忘差不多了。原始代码copy的matlab实现鼠标画图,但是在我电脑上跑不了,稍作改动。...%要在函数里load训练好的数字识别卷积神经网络。如果在函数外load,函数执行时,他的工作区...
  • 图像处理中的卷积---2.高斯卷积

    千次阅读 2015-11-26 16:36:27
    那么图像处理中的卷积可以理解为(2D/3D)的函数卷积,也就是在某个移动窗口内的加权(权值由kernel给定)求和。 下面我们给出gaussian kernel做卷积的例子,另附MATLAB实现方法。 sigma = 3;% ...
  • Matlab画冲激信号,stem方法,同一个窗口画多个坐标系,0.8n∗[u...MATLAB提供了conv函数可用来计算下面的卷积和 这里假设 和 都是有限长序列。如果 仅在 区间内为非零, 仅在 区间内为非零,那么 就仅在区间内为零...
  • MATLAB图形图像处理

    热门讨论 2011-01-03 12:20:11
    11.2.6 窗函数法 第十二章 MATLAB 图像运算 12.1 图像点的运算 12.1.1 线性点运算 12.1.2 非线性点运算 12.2 图像的算术运算 12.2.1 加法运算 12.2.2 减 法运算 12.2.3 乘法 运算 12.2.4 除法 运算 12.2.5...
  • VC++ matlab图像处理

    2013-09-15 15:17:54
    11.2.6 窗函数法 第十二章 MATLAB 图像运算 12.1 图像点的运算 12.1.1 线性点运算 12.1.2 非线性点运算 12.2 图像的算术运算 12.2.1 加法运算 12.2.2 减 法运算 12.2.3 乘法 运算 12.2.4 除法 运算 12.2.5...
  • Matlab设计

    2011-12-20 09:44:01
    (提示:Matlab函数length(x)可以计算出x(n)的长度) (2)验证傅立叶变换的时移和频移性质。 3.运用FFT分析信号频谱 (1)对矩形序列x1(n)=R8(n)进行频谱分析,截取长度N=64,记录幅频特性曲线和相频特性曲线;...
  • Matlab fir滤波(conv)

    千次阅读 2020-02-05 17:10:07
    在混合波形中提取目标波形,会经常使用到滤波器,这里通过matlab中的conv函数(向量卷积运算)实现滤波输出。 一:设置滤波器参数 在MATLAB中输入fdatool(Filer Dsign and analysis tool)即可调出滤波器设置界面...
  • 例4.12 由离散传递函数求脉冲过度函数 ex041300 由差分方程求系统函数及脉冲响应 ex050200 宽度L周期N的周期性方波的离散付利叶曲线 ex050500 不同的离散付利叶采样密度对应的时域曲线 ex050600 离散付利叶变换计算...
  • MATLAB下的数字信号处理实现示例 一、 信号、系统和系统响应 ...三、 窗函数法设计 FIR 滤波器 四、 IIR 滤波器的实现 程序下载链接,保证真实 https://mianbaoduo.com/o/bread/mbd-ZZabmpg= ...
  •  7.5基于MATLAB窗函数的FIR滤波器设计  7.6Kaiser窗函数  7.7截尾型傅里叶变换设计方法  7.8频率采样设计法 第8章最小均方设计方法  8.1有限冲激响应综述  8.2最小二乘法  8.3最小二乘FIR滤波器设计  8.4...
  • 数字信号处理课程设计基于matlab 文件内容 /conv(卷积部分) dupconv 重叠保留法 freqcirconv 频域圆周卷积 freqconv 频域卷积 timecirconv 时域圆周卷积 timeconv 时域卷积 /fft(频域...
  • 11.2.6 窗函数法 第十二章 MATLAB 图像运算 12.1 图像点的运算 12.1.1 线性点运算 12.1.2 非线性点运算 12.2 图像的算术运算 12.2.1 加法运算 12.2.2 减 法运算 12.2.3 乘法 运算 12.2.4 除法 运算 12.2.5...
  • 11.2.6 窗函数法 第十二章 MATLAB 图像运算 12.1 图像点的运算 12.1.1 线性点运算 12.1.2 非线性点运算 12.2 图像的算术运算 12.2.1 加法运算 12.2.2 减 法运算 12.2.3 乘法 运算 12.2.4 除法 运算 12.2.5...
  • 24.2.1窗函数 24.2.2数字滤波器频率响应函数 24.3特殊滤波器的MATLAB函数实现 第25章Subsystem模块创建子系统 25.1子系统 25.1.1创建子系统 25.1.2浏览下层子系统 25.1.3条件执行子系统 25.1.4控制流系统 ...
  • 11.2.6 窗函数法 第十二章 MATLAB 图像运算 12.1 图像点的运算 12.1.1 线性点运算 12.1.2 非线性点运算 12.2 图像的算术运算 12.2.1 加法运算 12.2.2 减 法运算 12.2.3 乘法 运算 12.2.4 除法 ...

空空如也

空空如也

1 2 3
收藏数 59
精华内容 23
关键字:

matlab卷积窗函数

matlab 订阅