matlab下的数字信号处理示例new

《matlab下的数字信号处理示例new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、MATLAB下的数字信号处理实现示例附录一信号、系统和系统响应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;x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”closeall%清除已经绘制的x(n)图形subplot(3,1,1);stem(x);%绘制x(n)的图形title(‘理想采样信号序列

2、’);（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）改变参数为：A=1,α=0.4,Ω=2.0734,T=10n=0:50;%定义序列的长度是50A=1;%设

3、置信号有关的参数a=0.4;T=1;%采样率w0=2.0734;x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”closeall%清除已经绘制的x(n)图形subplot(3,1,1);stem(x);%绘制x(n)的图形title(‘理想采样信号序列’);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(‘理想采样信号

4、序列的幅度谱’);angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(‘理想采样信号序列的相位谱’)2、单位脉冲序列在MatLab中，这一函数可以用zeros函数实现：n=1:50;%定义序列的长度是50x=zeros(1,50);%注意：MATLAB中数组下标从1开始x(1)=1;closeall;subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n’*k);magX=abs(X);

5、%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(‘单位冲击信号的相位谱’)3、矩形序列n=1:50x=sign(sign(10-n)+1);closeall;subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);k=-25:25;X=x*(exp(-j*pi/25)).^(n’*k);magX=abs(X);%绘制x(n)的幅度谱subp

6、lot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);angX=angle(X);%绘制x(n)的相位谱subplot(3,1,3);stem(angX);title(‘单位冲击信号的相位谱’)4、特定冲击串：x(n)=δ(n)+2.5δ(n−1)+2.5δ(n−2)+δ(n−3)n=1:50;%定义序列的长度是50x=zeros(1,50);%注意：MATLAB中数组下标从1开始x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;closeall;subplot(3,1,1);stem(x);tit

7、le(‘单位冲击信号序列’);k=-25:25;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(‘单位冲击信号的相位谱’)+∞5、卷积计算：y(n)=x(n)∗h(n)=∑x(m)h(n−m)m=−∞在MATLAB中。提供了卷积函数conv，即y=conv(x,h)，调用十分方便。例

8、如：系统：h(n)=δ(n)+2.5δ(n−1)+2.5δ(n−2)+δ(n−3)b−αnT信号：x(t)=Aesin(ΩnT),0≤n<50a0n=1:50;%定