利用MATLAB编写卷积函数myconv

《利用MATLAB编写卷积函数myconv》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、一、实验目的1.了解MATLAB的基本操作。2.利用MATLAB实现正弦信号采样；3.利用MATLAB编写卷积函数myconV。二、实验条件PC机，MATLAB7.0三、实验内容(-)函数文件与脚本文件的定义，正弦信号的釆样函数文件代码：functionseqs=mysampling(FszL,Fx)t=l/Fs:l/Fs:L/Fs;t1=0:0.001:L/Fs;x=sin(2*pi*Fx*t);xl=sin(2*pi*Fx*tl);plot(t1,xl),holdon;stem(t,x,1r1);holdoff;

2、脚本文件代码：Fs=40;Fx=20;L=15;forFx=20:10:80,mysampling(Fs,L,Fx);pause,End运行结果：1080604020-02-0406-08-100.050.1015020.250.30.350.4实验屮遇到的问题及改正：functionseqs=mysampling(Fs,L,Fx)t=1/Fs:1/Fs:L/Fs;x=sin(2*pi*Fx*t);plot(tzx)zholdon;stem(tzx,1r1);holdoff;一开始并没有添加tl=0:0.001:L/

3、Fs;xl=sin(2*pi*Fx*tl);两行代码，想利用t=l/Fs:l/Fs:L/Fs;x=sin(2*pi*Fx*t);生成的函数直接作为需要抽样的原函数，但是t的采样间隔太大，没法生成平滑的正弦曲线，所以导致最后仿真岀来的图像中岀现不光滑的正弦曲线。经过思考找出问题后，另外编写了一个正弦函数，设定的采样间隔为0.001后，仿真出来的图像中就出现了光滑的正弦曲线。(二)编写卷积函数myconvConv函数代码：functiony=myconvl(hfx)n=length(x);k=length(h);%定义序

4、列长度form=l:1:(n-1)A(m)=0;endform=n:1:n+k-1A(m)=h(m-n+l);%给h序列添(n-l)个0endB=hankel(A);fori=l:1:nC(:,i)=B(:,i);%抽取B矩阵屮从笫1列到第n列的矩阵endD=fliplr(C);%C矩阵左右翻转E=D';%转置y(n+k-1)=0;fora=l:1:n+k-1forb=l:1:nnumble(b)=x(b)*E(b,a);y(a)=numble(b)+y(a);endnumble=[000];貂青零endystem(

5、y);脚木代码：x=[l-231];h=[-32-14];%定义x,h序列myconvl(hAx);运行结果：»UntitledGy=■38-U9-911arts=-38-149-911A»I四、实验结论和讨论通过本次实验运用MATLAB实现采样以及自行编写卷积函数，基本能熟悉运用MATLAB来进行实验，并且对于采样以及卷积的认识更加的深入。在进行本次卷积实验之前，我一直只知道利用MATLAB本身自有的conv函数来进行卷积的运算，并没有想过conv函数的卷积运算功能是如何实现的。通过这次自己编写卷积函数的代码，不仅

6、熟知了hankel,fliplr以及转置等MATLAB中及其有用的用法，而且对于卷积函数有了更深一步的了解。并J4通过自己所写的myconvl函数与系统自带函数conv的比较，发现自己所编写的函数的运行时间要多于corw函数运行的时间，意识到MATLAB的深奥之处，也发现自己的能力还远远不足，需要在进一步的努力。