Matlab 数字信号处理编程导论 by Source Code.docxVIP

Matlab 数字信号处理编程导论 by Source Code最近在做信号处理和模式识别的相关技术的研究，有感于介绍这方面的入门文章太少，希望记录些有用的知识帮助需要从事相关工作的朋友少走些弯路。Matlab大家应该都不陌生，园子里面应该有相当一部分人是学电信或自动控制出身的，我们从事着不同的行业，我们每天写着不同的软件。如果你在做信号处理或者通信相关的行业，相信如果有一天你需要做些原型设计，本文可以引导你快速的开始。首先让我们来看看matlab的产品定位：The Language of Technical ComputingMATLAB??is a high-level language and interactive environment that enablesyou to perform computationally intensive tasks faster than with traditionalprogramming languages such as C, C++, and Fortran.?过多的也没必要解释了，就是说一个字“快”，使用matlab可以使你想法更快的被实现出来，做设计的效率是其他编程语言说无法比拟的。据本人所知，数字信号处理里面的滤波器设计，频谱分析，通信中的各种调制技术的设计，matlab可以说已经成了标准工具，就像windows上开发C#代码，大部分人都使用visual studio一样。?首先让我们看看如何产生一些简单有用的序列如何编写：---注：源代码来自互联网，仅供学习参考之用单位冲击：function?[x,n]?=?impseq(n0,n1,n2)%?Generates?x(n)?=?delta(n-n0);?n1?=?n,n0?=?n2%?----------------------------------------------%?[x,n]?=?impseq(n0,n1,n2)%if?((n0??n1)?|?(n0??n2)?|?(n1??n2))?error(arguments?must?satisfy?n1?=?n0?=?n2)endn?=?[n1:n2];%x?=?[zeros(1,(n0-n1)),?1,?zeros(1,(n2-n0))];x?=?[(n-n0)?==?0];?单位阶越：????function [x,n] = stepseq(n0,n1,n2)% Generates x(n) = u(n-n0); n1 = n,n0 = n2% ------------------------------------------% [x,n] = stepseq(n0,n1,n2)%if ((n0 n1) | (n0 n2) | (n1 n2))error(arguments must satisfy n1 = n0 = n2)endn = [n1:n2];%x = [zeros(1,(n0-n1)), ones(1,(n2-n0+1))];x = [(n-n0) = 0];?再看看如何实现一些简单算法和计算：序列相加：?function?[y,n]?=?sigadd(x1,n1,x2,n2)%?implements?y(n)?=?x1(n)+x2(n)%?-----------------------------%?[y,n]?=?sigadd(x1,n1,x2,n2)%???y?=?sum?sequence?over?n,?which?includes?n1?and?n2%??x1?=?first?sequence?over?n1%??x2?=?second?sequence?over?n2?(n2?can?be?different?from?n1)%n?=?min(min(n1),min(n2)):max(max(n1),max(n2));?%?duration?of?y(n)y1?=?zeros(1,length(n));?y2?=?y1;%?initializationy1(find((n=min(n1))(n=max(n1))==1))=x1;%?x1?with?duration?of?yy2(find((n=min(n2))(n=max(n2))==1))=x2;%?x2?with?duration?of?yy?=?y1+y2;%?sequence?addition?序列相乘：??function?[y,n]?=?sigmult(x1,n1,x2,n2)%?implements?y(n)?=?x1(n)*x2(n)%?-----------------------------%?[y,n]?=?sigmult(x1,n1,x2,n2