MATLAB下的数字信号处理实现示例

MATLAB  下的数字信号处理实现示例本部分内容是本讲义中数据信号处理实验部分实验项目在  MatLab  下实现代码。之所以提供这些代码，是希望通过研究以下代码，能够更快、更好地掌握用MatLab  进行数据信号处理实验的方法；提高实验质量。希望同学们在阅读代码的时候，注意学习方法，在最短的时间内熟悉MatLab，提高应用能力。示例中有些部分是实验项目中的内容实现，有些是一些典型例题的实现。研究示例代码，倡导个性化编程是我们的目标，希望同学们能在在进行实验项目的过程中提高MatLab  的应用能力；在学习MatLab  编程的同时加强对数字信号处理有关实验项目的理解。以下代码段均在  MatLab5.3  下调试通过，但是由于排版或其他一些原因，可能有部分代码段不能得到正常结果。您可以在“http://202.38.75.33/dsp/matlab/”得到本讲义的修订内容，同时可以在这个网址获取所有代码。附录  C1  信号、系统和系统响应1、理想采样信号序列（1）首先产生信号x(n),0<=n<=50n=0:50;  %定义序列的长度是50A=444.128;  %设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001;  %采样率w0=50*sqrt(2.0)*pi;  %ω符号在MatLab  中不能输入，用w  代替x=A*exp(-a*n*T).*sin(w0*n*T);  %pi  是MATLAB  定义的π，信号乘可采用“.*”close  all  %清除已经绘制的x(n)图形subplot(3,1,1);stem(x);  %绘制x(n)的图形title(‘理想采样信号序列’);  %设置结果图形的标题（2）绘制信号x(n)的幅度谱和相位谱k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X);  %绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’);angX=angle(X);  %绘制x(n)的相位谱subplot(3,1,3);stem(angX)  ;  title  (‘理想采样信号序列的相位谱’)（3）改变参数为：  1,  0.4,  2.0734,  1  0  A  =  α  =  Ω  =  T  =n=0:50;  %定义序列的长度是50A=1;  a=0.4;  w0=2.0734;  T=1;  %设置信号有关的参数和采样率Tx=A*exp(-a*n*T).*sin(w0*n*T);  %pi  是MATLAB  定义的π，信号乘可采用“.*”