Matlab几种卷积的实现与比较conv与filter,conv2与filter2 matlab conv2

最近在做控制算法实现的时候,对于其中参杂的各种差分、卷积很头疼,就在网上搜集了些资料,汇总于此,以做备忘。

在MATLAB中,可以用函数y=filter(p,d,x)实现差分方程的仿真,也可以用函数y=conv(x,h)计算卷积。

(1)即y=filter(p,d,x)用来实现差分方程,d表示差分方程输出y的系数,p表示输入x的系数,而x表示输入序列。输出结果长度数等于x的长度。

实现差分方程,先从简单的说起:
filter([1,2],1,[1,2,3,4,5]),实现y[k]=x[k]+2*x[k-1]
y[1]=x[1]+2*0=1(x[1]之前状态都用0)
y[2]=x[2]+2*x[1]=2+2*1=4

(2)y=conv(x,h)是用来实现卷级的,对x序列和h序列进行卷积,输出的结果个数等于x的长度与h的长度之和减去1。

卷积公式:z(n)=x(n)*y(n)=∫x(m)y(n-m)dm.

程序一:以下两个程序的结果一样

(1)h = [3 2 1 -2 1 0 -4 0 3];%impulse response

x = [1 -2 3 -4 3 2 1];%input sequence

y = conv(h,x);

n = 0:14;

subplot(2,1,1);

stem(n,y);

xlabel('Timeindex n'); ylabel('Amplitude');

title('Output Obtained by Convolution'); grid;

(2)x1 = [x zeros(1,8)];

y1 = filter(h,1,x1);

subplot(2,1,2);

stem(n,y1);

xlabel('Time index n'); ylabel('Amplitude');

title('OutputGenerated by Filtering'); grid;

程序二:filter和conv的不同

x=[1,2,3,4,5];
h=[1,1,1];

y1=conv(h,x)
y2=filter(h,1,x)
y3=filter(x,1,h)

结果:y1=13691295

y2=136912

‍ y3=136

可见:filter函数y(n)是从n=1开始,认为所有n<1都为0;而conv是从卷积公式计算,包括n<1部分。

因此filter 和conv的结果长短不同

程序三:滤波后信号幅度的变化

num=100; %总共1000个数
x=rand(1,num); %生成0~1随机数序列
x(x>0.5)=1;
x(x<=0.5)=-1;
h1=[0.2,0.5,1,0.5,0.2];
h2=[0,0,1,0,0];
y1=filter(h1,1,x);
y2=filter(h2,1,x);
n=0:99;
subplot(2,1,1);
stem(n,y1);
subplot(2,1,2);
stem(n,y2);

MATLAB中提供了卷积运算的函数命令conv2,其语法格式为:
C = conv2(A,B)
C =conv2(A,B)返回矩阵A和B的二维卷积C。若A为ma×na的矩阵,B为mb×nb的矩阵,则C的大小为(ma+mb-1)×(na+nb-1)。

例:
A=magic(5)
A =
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
>> B=[1 2 1 ;0 2 0;3 1 3]
B =
1 2 1
0 2 0
3 1 3
>> C=conv2(A,B)
C =
17 58 66 34 32 38 15
23 85 88 35 67 76 16
55 149 117 163 159 135 67
79 78 160 161 187 129 51
23 82 153 199 205 108 75
30 68 135 168 91 84 9
33 65 126 85 104 15 27
MATLAB图像处理工具箱提供了基于卷积的图象滤波函数filter2,filter2的语法格式为:
Y = filter2(h,X)
其中Y =filter2(h,X)返回图像X经算子h滤波后的结果,默认返回图像Y与输入图像X大小相同。例如:
其实filter2和conv2是等价的。MATLAB在计算filter2时先将卷积核旋转180度,再调用conv2函数进行计算。
Fspecial函数用于创建预定义的滤波算子,其语法格式为:
h = fspecial(type)
h = fspecial(type,parameters)
Matlab几种卷积的实现与比较(conv与filter,conv2与filter2) matlab conv2
参数type制定算子类型,parameters指定相应的参数,具体格式为:
type='average',为均值滤波,参数为n,代表模版尺寸,用向量表示,默认值为[3,3]。
type='gaussian',为高斯低通滤波器,参数有两个,n表示模版尺寸,默认值为[3,3],sigma表示滤波器的标准差,单位为像素,默认值为0.5

  

爱华网本文地址 » http://www.413yy.cn/a/25101016/295903.html

更多阅读

扎头发的方法有几种 精 简单的扎头发方法

今天在网上看到了几种扎头发的方法,真的是简单易学,而且都很漂亮精致,下面小编就为大家一一介绍吧。扎头发的方法有几种 精——第一款扎头发的方法有几种 精 1、先把头顶的头发分成三份扎头发的方法有几种 精 2、把2放在3的下面扎头发

几种花草常见的害虫及防治小妙招 常见害虫

很多人喜欢在室内养些花草,既绿化环境又陶冶情趣。但是花草不仅养起来有些娇气,还很容易受到一些小虫的“青睐”,使人烦恼,又很无奈。给大家介绍几种花草常见的害虫及防治方法。菊花易招蚜虫蚜虫,多是绿色和黑色,

教你几种袜子娃娃的做法 用袜子做娃娃教程

几种袜子娃娃的做法,姐妹们可以来学习一下了。袜子娃娃在一针一线中诞生,缝缝补补中成长,袜子娃娃顾名思义就是用袜子做的娃娃.之所以用袜子,一方面是因为袜子的图案比较立体,另一方面是袜子的弹性比较大,塞进丝绵之后弹性会比一般的布

黑木耳功效和几种吃法 黑木耳的功效与作用

黑木耳蒸吃有营养,但必须隔水蒸小时,然后可以用芝麻酱拌着吃,如果再加黄瓜拌会起到功半事倍的效果。降脂,防血稠,而且还有减肥的作用。黑木耳炒黄花菜木耳含有丰富的蛋白质、铁、磷、维生素B等健脑需要的营养素,其中维生素B2含量较蔬菜高

声明:《Matlab几种卷积的实现与比较conv与filter,conv2与filter2 matlab conv2》为网友劳资很爷们分享!如侵犯到您的合法权益请联系我们删除